A Eulerian fixed grid approach to simulate the dynamics of two-phase interfaces in the presence of surface tension forces is presented. This level set/vortex sheet method consists of a simplified system of equations that contain individual source terms describing the relevant physical processes at the phase interface explicitly. Hence, this approach provides a framework that will allow for a si...

Phase-field models of solidification with convection often assume the existence of a single (mixture) velocity at any location inside the diffuse interface, and the phase-field, φ, is advected by this mixture velocity. In this paper, the advection of the phase-field is examined for a one-dimensional normal flow to a solidification front induced by a density difference between the solid and liqu...

A novel diffuse interface model is presented for the direct numerical simulation of microstructure evolution in solidification processes involving convection in the liquid phase. The solidification front is treated as a moving interface in the diffuse approximation as known from phase-field theories. The no-slip condition between the melt and the solid is realized via a drag resistivity in the ...

The numerical approximation of liquid vapor flows within the compressible regime is a challenging task because complex physical effects at the phase interfaces govern the global flow behavior. We develop a sharp interface approach which treats the phase boundary like a shock wave discontinuity and takes capillarity effects into account. The approach relies on the solution of Riemann problems ac...

A two-phase diffuse interface model previously developed by the authors is used to simulate the buoyancy-driven flow and Rayleigh-Taylor instability of fluid layers inside a Hele-Shaw cell. The model assumes that the two phases coexist inside the diffuse interface with different properties and velocities. A separate momentum equation is used to calculate the slip velocity between the two phases...

We propose a phase field model that approximates its limiting sharp interface model (free boundary problem) up to second order in interface thickness. A broad range of double-well potentials can be utilized so long as the dynamical coefficient in the phase equation is adjusted appropriately. This model thereby assures that computation with particular value of interface thickness ε, will differ ...