Efficient two-phase mass-conserving level set method for simulation of incompressible turbulent free surface flows with large density ratio

2016 ◽  
Vol 136 ◽  
pp. 212-227 ◽  
Author(s):  
J.M. Cubos-Ramírez ◽  
J. Ramírez-Cruz ◽  
M. Salinas-Vázquez ◽  
W. Vicente-Rodríguez ◽  
E. Martinez-Espinosa ◽  
...  
2012 ◽  
Vol 11 (1) ◽  
pp. 215-248 ◽  
Author(s):  
Xin Lv ◽  
Qingping Zou ◽  
D.E. Reeve ◽  
Yong Zhao

AbstractWe present a three dimensional preconditioned implicit free-surface capture scheme on tetrahedral grids. The current scheme improves our recently reported method [10] in several aspects. Specifically, we modified the original eigensystem by applying a preconditioning matrix so that the new eigensystem is virtually independent of density ratio, which is typically large for practical two-phase problems. Further, we replaced the explicit multi-stage Runge-Kutta method by a fully implicit Euler integration scheme for the Navier-Stokes (NS) solver and the Volume of Fluids (VOF) equation is now solved with a second order Crank-Nicolson implicit scheme to reduce the numerical diffusion effect. The preconditioned restarted Generalized Minimal RESidual method (GMRES) is then employed to solve the resulting linear system. The validation studies show that with these modifications, the method has improved stability and accuracy when dealing with large density ratio two-phase problems.


2005 ◽  
Vol 19 (28n29) ◽  
pp. 1743-1746
Author(s):  
HANBIN GU ◽  
YANBAO LI ◽  
PENGZHI LIN

A three-dimension numerical wave model (3DWAVE) has been developed to simulate free surface flows. The model solves Navier-Stokes equations for two-phase flows of air and water. The level set method is employed to track water surfaces. The model is tested for water sloshing in a 3-D confined tank. The relative error in the mass and total energy computation is less than 1%. Excellent agreements between numerical results and analytical solution are obtained for free surface calculation. The nonlinearity in the 3-D fluid sloshing is analyzed. These have laid a foundation on research of breaking waves.


Author(s):  
Doug Enright ◽  
Duc Nguyen ◽  
Frederic Gibou ◽  
Ron Fedkiw

In this paper, we present an enhanced resolution capturing method for topologically complex two and three dimensional incompressible free surface flows. The method is based upon the level set method of Osher and Sethian to represent the interface combined with two recent advances in the treatment of the interface, a second order accurate discretization of the Dirichlet pressure boundary condition at the free surface (2002, J. Comput. Phys.176, 205) and the use of massless marker particles to enhance the resolution of the interface through the use of the particle level set method (2002, J. Comput. Phys., 183, 83). Use of these methods allow for the accurate movement of the interface while at the same time preserving the mass of the liquid, even on coarse computational grids. Also, these methods complement the level set method in its ability to handle changes in interface topology in a robust manner. Surface tension effects can be easily included in our method. The method is presented in three spatial dimensions, with numerical examples in both two and three spatial dimensions.


2018 ◽  
Vol 01 (02) ◽  
pp. 1840001 ◽  
Author(s):  
Sauro Manenti

The standard weakly compressible Smoothed Particle Hydrodynamics (WCSPH) is successfully applied to multi-phase problems involving fluids with similar densities, but when density ratio increases at some order of magnitude, serious instability phenomena occur at the interface. Several remedies have been proposed based on numerical correctives that deviate from standard formulation, increasing the algorithm complexity and, sometimes, the computational cost. In this study, the standard SPH has been adapted to treat free-surface multi-phase flows with a large density ratio through a modified form of the governing equations which is based on the specific volume (i.e. the inverse of particle volume) instead of density: the former is continuous across the fluid interface while the latter is not and generates numerical instability. Interface sharpness is assured without cohesion forces; kernel truncation at the interface is avoided. The model, relatively simple to implement, is tested by simulating two-phase dam breaking for two configurations: kinematic and dynamic features are compared with reference data showing good agreement despite the reduced computational effort.


Sign in / Sign up

Export Citation Format

Share Document