Modelling merging and fragmentation in multiphase flows with SURFER
Journal of Computational Physics
A level set approach for computing solutions to incompressible two-phase flow
Journal of Computational Physics
A boundary condition capturing method for Poisson's equation on irregular domains
Journal of Computational Physics
Journal of Computational Physics
An adaptive mesh algorithm for evolving surfaces: simulation of drop breakup and coalescence
Journal of Computational Physics
Interfacial dynamics for Stokes flow
Journal of Computational Physics
Boundary integral methods for multicomponent fluids and multiphase materials
Journal of Computational Physics
Level set methods: an overview and some recent results
Journal of Computational Physics
A front-tracking method for the computations of multiphase flow
Journal of Computational Physics
Journal of Computational Physics
A finite element method for three-dimensional unstructured grid smoothing
Journal of Computational Physics
Adaptive tetrahedral meshing in free-surface flow
Journal of Computational Physics
Adaptive unstructured volume remeshing - I: The method
Journal of Computational Physics
Journal of Computational Physics
A moving mesh interface tracking method for 3D incompressible two-phase flows
Journal of Computational Physics
Journal of Computational Physics
Journal of Computational Physics
Hi-index | 31.45 |
In multiphase flows, the length scales of thin regions, such as thin films between nearly touching drops and thin threads formed during the interface pinch-off, are usually several orders of magnitude smaller than the size of the drops. In this paper, a number of extra length criteria for adaptive meshes are developed and implemented in the moving mesh interface tracking method to solve these multiple-length-scale problems with high fidelity. A nominal length scale based on the solutions of Laplace's equations with the unit normal vectors of surfaces as the boundary conditions is proposed for the adaptive mesh refinement in the thin regions. For almost flat interfaces/boundaries which are near to the thin regions, the averaged length of the interior edges sharing the two nodes with the boundary edge is introduced for the mesh adaptation. The averaged length of the interfacial edges is used for the interior elements near the interfaces but outside of the thin regions. For the interior mesh away from the interfaces/boundaries, different averaged length scales based on the initial mesh are employed for the adaptive mesh refining and coarsening. Numerous cases are simulated to demonstrate the capability of the proposed schemes in handling multiple length scales, which include the relaxation and necking of an elongated droplet, droplet-droplet head-on approaching, droplet-wall interactions, and a droplet pair in a shear flow. The smallest length resolved for the thin regions is three orders of magnitude smaller than the largest characteristic length of the problem.