A CutFEM method for two-phase flow problems

In this article, we present a cut finite element method for two-phase Navier-Stokes flows. The main feature of the method is the formulation of a unified continuous interior penalty stabilisation approach for, on the one hand, stabilising advection and the pressure-velocity coupling and, on the other hand, stabilising the cut region. The accuracy of the algorithm is enhanced by the development of extended fictitious domains to guarantee a well defined velocity from previous time steps in the current geometry. Finally, the robustness of the moving-interface algorithm is further improved by the introduction of a curvature smoothing technique that reduces spurious velocities. The algorithm is shown to perform remarkably well for low capillary number flows, and is a first step towards flexible and robust CutFEM algorithms for the simulation of microfluidic devices

Multiphase modeling of melting : solidification with high density variations using XFEM

La modélisation de la cryolite, utilisée dans la fabrication de l’aluminium, implique plusieurs défis, notament la présence de discontinuités dans la solution et l’inclusion de la difference de densité entre les phases solide et liquide. Pour surmonter ces défis, plusieurs éléments novateurs ont été développés dans cette thèse. En premier lieu, le problème du changement de phase, communément appelé problème de Stefan, a été résolu en deux dimensions en utilisant la méthode des éléments finis étendue. Une formulation utilisant un multiplicateur de Lagrange stable spécialement développée et une interpolation enrichie a été utilisée pour imposer la température de fusion à l’interface. La vitesse de l’interface est déterminée par le saut dans le flux de chaleur à travers l’interface et a été calculée en utilisant la solution du multiplicateur de Lagrange. En second lieu, les effets convectifs ont été inclus par la résolution des équations de Stokes dans la phase liquide en utilisant la méthode des éléments finis étendue aussi. Troisièmement, le changement de densité entre les phases solide et liquide, généralement négligé dans la littérature, a été pris en compte par l’ajout d’une condition aux limites de vitesse non nulle à l’interface solide-liquide pour respecter la conservation de la masse dans le système. Des problèmes analytiques et numériques ont été résolus pour valider les divers composants du modèle et le système d’équations couplés. Les solutions aux problèmes numériques ont été comparées aux solutions obtenues avec l’algorithme de déplacement de maillage de Comsol. Ces comparaisons démontrent que le modèle par éléments finis étendue reproduit correctement le problème de changement phase avec densités variables.The modelling of the cryolite bath, used in the smelting of aluminum, offers multiple challenges, particularly the presence of discontinuities in the solution and a difference in density between the solid and liquid phases. To over come these challenges, several novel elements were developed in this thesis. First of all, the phase change problem, commonly named the Stefan problem, was solved in two dimensions using the extended finite element method. A specially designed Lagrange multiplier formulation, using an enriched Lagrange multiplier solution, was implemented to impose the melting temperature on the interface. The interface velocity is determined by the jump in the heat flux across the interface and was calculated using the Lagrange multiplier values. Secondly, convection was included by solving the Stokes equations in the liquid phase using the extended finite element method as well. Thirdly, the density change between solid and liquid phases, usually neglected in the literature, was taken into account by the addition of a non-zero velocity boundary condition at the solid-liquid interface to maintain mass conservation in the system. Benchmark analytical and numerical problems were solved to validated the various components of the model and the coupled system of equations. The solutions to the numerical problems were compared to the solutions obtained using Comsol’s moving mesh algorithm. Theses comparisons show that the extended finite element method correctly solves the phase change problem with non-constant densities

The LifeV library: engineering mathematics beyond the proof of concept

LifeV is a library for the finite element (FE) solution of partial differential equations in one, two, and three dimensions. It is written in C++ and designed to run on diverse parallel architectures, including cloud and high performance computing facilities. In spite of its academic research nature, meaning a library for the development and testing of new methods, one distinguishing feature of LifeV is its use on real world problems and it is intended to provide a tool for many engineering applications. It has been actually used in computational hemodynamics, including cardiac mechanics and fluid-structure interaction problems, in porous media, ice sheets dynamics for both forward and inverse problems. In this paper we give a short overview of the features of LifeV and its coding paradigms on simple problems. The main focus is on the parallel environment which is mainly driven by domain decomposition methods and based on external libraries such as MPI, the Trilinos project, HDF5 and ParMetis. Dedicated to the memory of Fausto Saleri.Comment: Review of the LifeV Finite Element librar

Space-time unfitted finite element methods for time-dependent problems on moving domains

We propose a space-time scheme that combines an unfitted finite element method in space with a discontinuous Galerkin time discretisation for the accurate numerical approximation of parabolic problems with moving domains or interfaces. We make use of an aggregated finite element space to attain robustness with respect to the cut locations. The aggregation is performed slab-wise to have a tensor product structure of the space-time discrete space, which is required in the numerical analysis. We analyse the proposed algorithm, providing stability, condition number bounds and anisotropic \emph{a priori} error estimates. A set of numerical experiments confirm the theoretical results for a parabolic problem on a moving domain. The method is applied for a mass transfer problem with changing topology

A note on the penalty parameter in Nitsche's method for unfitted boundary value problems

Nitsche's method is a popular approach to implement Dirichlet-type boundary conditions in situations where a strong imposition is either inconvenient or simply not feasible. The method is widely applied in the context of unfitted finite element methods. From the classical (symmetric) Nitsche's method it is well-known that the stabilization parameter in the method has to be chosen sufficiently large to obtain unique solvability of discrete systems. In this short note we discuss an often used strategy to set the stabilization parameter and describe a possible problem that can arise from this. We show that in specific situations error bounds can deteriorate and give examples of computations where Nitsche's method yields large and even diverging discretization errors