NUMERICAL SIMULATION OF RAYLEIGH–TAYLOR INSTABILITY BASED ON AN IMPROVED PARTICLE LEVEL SET METHOD

2014 ◽  
Vol 11 (04) ◽  
pp. 1350094 ◽  
Author(s):  
HUI TIAN ◽  
GUOJUN LI ◽  
XIONGWEN ZHANG

An improved particle correction procedure for particle level set method is proposed and applied to the simulation of Rayleigh–Taylor instability (RTI) of the incompressible two-phase immiscible fluids. In the proposed method, an improved particle correction method is developed to deal with all the relative positions between escaped particles and cell corners, which can reduce the disturbance arising in the distance function correction process due to the non-normal direction movement of escaped particles. The improved method is validated through accurately capturing the moving interface of the Zalesak's disk. Furthermore, coupled with the projection method for solving the Navier–Stokes equations, the time-dependent evolution of the RTI interface over a wide range of Reynolds numbers, Atwood numbers and Weber numbers are numerically investigated. A good agreement between the present results and the existing analytical solutions is obtained and the accuracy of the proposed method is further verified. Moreover, the effects of control parameters including viscosity, density ratio, and surface tension coefficient on the evolution of RTI are analyzed in detail, and a critical Weber number for the development of RTI is found.

2009 ◽  
Vol 622 ◽  
pp. 115-134 ◽  
Author(s):  
ANTONIO CELANI ◽  
ANDREA MAZZINO ◽  
PAOLO MURATORE-GINANNESCHI ◽  
LARA VOZELLA

The Rayleigh–Taylor instability of two immiscible fluids in the limit of small Atwood numbers is studied by means of a phase-field description. In this method, the sharp fluid interface is replaced by a thin, yet finite, transition layer where the interfacial forces vary smoothly. This is achieved by introducing an order parameter (the phase-field) continuously varying across the interfacial layers and uniform in the bulk region. The phase-field model obeys a Cahn–Hilliard equation and is two-way coupled to the standard Navier–Stokes equations. Starting from this system of equations we have first performed a linear analysis from which we have analytically rederived the known gravity–capillary dispersion relation in the limit of vanishing mixing energy density and capillary width. We have performed numerical simulations and identified a region of parameters in which the known properties of the linear phase (both stable and unstable) are reproduced in a very accurate way. This has been done both in the case of negligible viscosity and in the case of non-zero viscosity. In the latter situation, only upper and lower bounds for the perturbation growth rate are known. Finally, we have also investigated the weakly nonlinear stage of the perturbation evolution and identified a regime characterized by a constant terminal velocity of bubbles/spikes. The measured value of the terminal velocity is in agreement with available theoretical prediction. The phase-field approach thus appears to be a valuable technique for the dynamical description of the stages where hydrodynamic turbulence and wave-turbulence come into play.


2005 ◽  
Vol 83 (6-7) ◽  
pp. 479-490 ◽  
Author(s):  
Douglas Enright ◽  
Frank Losasso ◽  
Ronald Fedkiw

2004 ◽  
Vol 126 (4) ◽  
pp. 578-585 ◽  
Author(s):  
Hiroyuki Takahira ◽  
Tomonori Horiuchi ◽  
Sanjoy Banerjee

For the present study, we developed a three-dimensional numerical method based on the level set method that is applicable to two-phase systems with high-density ratio. The present solver for the Navier-Stokes equations was based on the projection method with a non-staggered grid. We improved the treatment of the convection terms and the interpolation method that was used to obtain the intermediate volume flux defined on the cell faces. We also improved the solver for the pressure Poisson equations and the reinitialization procedure of the level set function. It was shown that the present solver worked very well even for a density ratio of the two fluids of 1:1000. We simulated the coalescence of two rising bubbles under gravity, and a gas bubble bursting at a free surface to evaluate mass conservation for the present method. It was also shown that the volume conservation (i.e., mass conservation) of bubbles was very good even after bubble coalescence.


PAMM ◽  
2004 ◽  
Vol 4 (1) ◽  
pp. 506-507 ◽  
Author(s):  
Daniel Gaudlitz ◽  
Nikolaus A. Adams

Author(s):  
Daniel Gaudlitz ◽  
Nikolaus A. Adams

The hybrid particle-level-set method was proposed by Enright et al. (2002) for improving the efficiency and accuracy of the original level-set method. In this paper, we examine the discretization of the surface-tension term and of the reinitialization procedure with respect to their effect on the mass-conservation properties of the overall method. We apply the hybrid particle-level-set method for the computation of two-phase flows and assess its performance with different parameter choices in comparison with the original level-set method.


Sign in / Sign up

Export Citation Format

Share Document