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

JOURNAL OF COMPUTATIONAL PHYSICS

卷:228

Article

Volkov, Oleg¹  Protas, Bartosz¹ 

[1] McMaster Univ, Dept Math & Stat, Hamilton, ON L8S 4K1, Canada

关键词: Free-boundary problem;    Stefan conditions;    Contact line;    Shape calculus;    Sobolev gradients;   

DOI  :  10.1016/j.jcp.2009.03.042

来源: Elsevier

【 摘 要 】

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. (c) 2009 Elsevier Inc. All rights reserved.

【 授权许可】

Free