We
have developed a cut finite element method for one-phase Stefan
problems, with applications in laser manufacturing. The geometry of the
workpiece is represented implicitly via a level set function. Material
above the melting/vaporisation temperature is represented by a
fictitious gas phase. The moving interface between the workpiece and the
fictitious gas phase may cut arbitrarily through the elements of the
finite element mesh, which remains fixed throughout the simulation,
thereby circumventing the need for cumbersome remeshing operations. The
primal/dual formulation of the linear one-phase Stefan problem is recast
into a primal non-linear formulation using a Nitsche-type approach,
which avoids the difficulty of constructing inf-sup stable primal/dual
pairs. Through the careful derivation of stabilisation terms, we show
that the proposed Stefan-Signorini-Nitsche CutFEM method remains stable
independently of the cut location. In addition, we obtain optimal
convergence with respect to space and time refinement. Several 2D and 3D
examples are proposed, highlighting the robustness and flexibility of
the algorithm, together with its relevance to the field of
micro-manufacturing.
Simulations by S. Claus, S. Bigot and P. Kerfriden in the FEniCS library CutFEM.
S. Claus, S. Bigot and P. Kerfriden, CutFEM Method for Stefan--Signorini Problems with Application in Pulsed Laser Ablation, SIAM J. Sci. Comput., 40(5), 2018
Funding: Sêr Cymru National Research Network
No comments:
Post a Comment