Improved Conservative Phase Field Method for Contact Line Motion With Solidification

Phase field method is appealing due to its implicit capture of the interface. It mainly contains two types, one being the Allen–Cahn and the other being the Cahn-Hilliard. The latter conserves mass, but is fourth order, not as amenable to numerical methods as the former. Recently, a conservative Allen–Cahn phase field method has been put forward, improving the boundedness of the phase field. The new method conserves mass well and is only second order. This paper couples the conservative Allen–Cahn phase field model and the Navier–Stokes equation and the heat balance equation. Interfacial tension is expressed in a potential form. The model is discretized using a finite difference method on a half-staggered grid, and the velocity–pressure coupling is decoupled using an explicit projection method. The heat balance equation is solved using an iterative method. The model was tested against the classic static bubble case and was applied to impact in various conditions. Moreover, the model is parallelized using shared memory parallelism, OpenMP.