Incompressible lagrangian fluid flow with thermal coupling

Abstract

A method is presented for the solution of an incompressible viscous fluid flow with heat transfer and solidification using a fully Lagrangian description of the motion. The originality of this method consists in assembling various concepts and techniques which appear naturally due to the Lagrangian formulation. First of all, the Navier-Stokes equations of motion coupled with the Boussinesq approximation must be reformulated in the Lagrangian framework, whereas they have been mostly derived in an Eulerian context. Secondly, the Lagrangian formulation implies to follow the material particles during their motion, which means to convect the mesh in the case of the Finite Element Method (FEM), the spatial discretisation method chosen in this work. This provokes various difficulties for the mesh generation, mainly in three dimensions, whereas it eliminates the classical numerical difficulty to deal with the convective term, as much in the Navier-Stokes equations as in the energy equation. Even without the discretization of the convective term, an efficient iterative solver, which constitutes the only viable alternative for three dimensional problems, must be designed for the class of Generalized Stokes Problems (GSP), which could be able to behave well independently of the mesh Reynolds number, as it can vary greatly for coupled fluid-thermal analysis.

Moreover, it offers a natural framework to treat free-surface problems like wave breaking and rough fluid-structure contact. On one hand, the convection of the mesh during one time step after the resolution of the non-linear system provides explicitly the locus of the domain to be considered. On the other hand, fluid-to-fluid and fluid-to-wall contact, as well as the update of the domain due to the remeshing, must be accurately and efficiently performed. Finally, the solidification of the fluid coupled with its motion through a variable viscosity is considered An efficient overall algorithm must be designed to bring the method effective, particularly in a three dimensional context, which is the ambition of this monograph. Various numerical examples are included to validate and highlight the potential of the method.