An inverse model for a free-boundary problem with a contact line: Steady case

• Published in: Journal of computational physics Vol. 228; no. 13; pp. 4893 - 4910
• Main Authors: Volkov, Oleg, Protas, Bartosz
• Format: Journal Article
• Language: English
• Published: Kidlington Elsevier Inc 20.07.2009; Elsevier
• Subjects: BOUNDARY CONDITIONS; Computational techniques; Contact line; Exact sciences and technology; Free-boundary problem; ITERATIVE METHODS; MATHEMATICAL METHODS AND COMPUTING; Mathematical methods in physics; MATHEMATICAL SOLUTIONS; PARTIAL DIFFERENTIAL EQUATIONS; Physics; Shape calculus; Sobolev gradients; SOLIDIFICATION; Stefan conditions; Sobolev gradients; 64.70.D; 02.30.Zz; Stefan conditions; Shape calculus; 44.05.+e; Free-boundary problem; Contact line; Partial differential equations; Stefan problem; Free boundary problem; Calculation methods; Closure model; Inverse problems; Functionals; Positions; Calculation; Iterative methods; Performance
• ISSN: 0021-9991; 1090-2716
• DOI: 10.1016/j.jcp.2009.03.042

This paper reformulates the two-phase solidification problem (i.e., the Stefan problem) as an inverse problem in which a cost functional is minimized with respect to the position of the interface and subject to PDE constraints. An advantage of this formulation is that it allows for a thermodynamically consistent treatment of the interface conditions in the presence of a contact point involving a third phase. It is argued that such an approach in fact represents a closure model for the original system and some of its key properties are investigated. We describe an efficient iterative solution method for the Stefan problem formulated in this way which uses shape differentiation and adjoint equations to determine the gradient of the cost functional. Performance of the proposed approach is illustrated with sample computations concerning 2D steady solidification phenomena.