Skip to Main content Skip to Navigation
Journal articles

A Newton method with adaptive finite elements for solving phase-change problems with natural convection

Frédéric Hecht 1 Ionut Danaila 2 Moglan Raluca 2 Stéphane Le Masson 3
1 ALPINES - Algorithms and parallel tools for integrated numerical simulations
LJLL - Laboratoire Jacques-Louis Lions, INSMI - Institut National des Sciences Mathématiques et de leurs Interactions, Inria de Paris
Abstract : We present a new numerical system using finite elements with mesh adaptivity for the simulation of solid–liquid phase change systems. In the liquid phase, the natural convection flow is simulated by solving the incompressible Navier–Stokes equations with Boussinesq approximation. A variable viscosity model allows the velocity to progressively vanish in the solid phase, through an intermediate mushy region. The phase change is modeled by introducing an implicit enthalpy source term in the heat equation. The final system of equations describing the liquid–solid system by a single domain approach is solved using a Newton iterative algorithm. The space discretization is based on a P2–P1 Taylor–Hood finite elements and mesh adaptivity by metric control is used to accurately track the solid–liquid interface or the density inversion interface for water flows. The numerical method is validated against classical benchmarks that progressively add strong non-linearities in the system of equations: natural convection of air, natural convection of water, melting of a phase-change material and water freezing. Very good agreement with experimental data is obtained for each test case, proving the capability of the method to deal with both melting and solidification problems with convection. The presented numerical method is easy to implement using FreeFem++ software using a syntax close to the mathematical formulation.
Document type :
Journal articles
Complete list of metadatas
Contributor : Frédéric Hecht <>
Submitted on : Wednesday, August 31, 2016 - 5:35:51 PM
Last modification on : Friday, March 27, 2020 - 4:05:31 AM



Frédéric Hecht, Ionut Danaila, Moglan Raluca, Stéphane Le Masson. A Newton method with adaptive finite elements for solving phase-change problems with natural convection. Journal of Computational Physics, Elsevier, 2014, 274 (1), pp.826-840. ⟨10.1016/⟩. ⟨hal-01358437⟩



Record views