scholarly journals A geometric mass control approach in level set method to simulate multiphase flows with complex interface topologies, case study: Oblique coalescence of gas bubbles in a liquid

Author(s):  
Amin Hadidi ◽  
Majid Eshagh Nimvari

A geometrical mass control loop is proposed in this research to use in the level set method in order to simulation of multiphase flows with complex topologies of the interface and a case study is investigated using proposed scheme. In this regard oblique interaction and coalescence of bubbles in a liquid is investigated. The level set method is suffering from poor mass conservation in the case of severe changes of interface and complex topologies encountered in a wide range of problems which one of them is oblique coalescence of the bubbles. Despite the use of full re-initialization and reconstruction approach of the level set method as well as application of fine mesh, deviation of mass conservation of the method even becomes 100%. Therefore, simulation of such problems sometimes becomes impossible using this method. So in the geometric mass control loop, mass deviation in each time step is calculated and is compensated in the dispersed phase, which prevents the propagation of mass error entire the simulation. Efficiency of proposed geometrical mass control loop is verified by simulation of oblique interaction and coalescence of gas bubbles in a liquid. The governing equations are continuity and momentum equations which have been discretized using the finite volume method and the SIMPLE algorithm. The results outlined in the present study well agree with the existing experimental and numerical results. Results show that the maximum amount of mass dissipation was less than 4%. Therefore, the level set method with proposed geometric mass control loop could be used properly for the simulation of multiphase flows with sharp and high variations in the interface.

2018 ◽  
Vol 28 (6) ◽  
pp. 1256-1278 ◽  
Author(s):  
Ali Karakus ◽  
Tim Warburton ◽  
Mehmet Haluk Aksel ◽  
Cuneyt Sert

Purpose This study aims to focus on the development of a high-order discontinuous Galerkin method for the solution of unsteady, incompressible, multiphase flows with level set interface formulation. Design/methodology/approach Nodal discontinuous Galerkin discretization is used for incompressible Navier–Stokes, level set advection and reinitialization equations on adaptive unstructured elements. Implicit systems arising from the semi-explicit time discretization of the flow equations are solved with a p-multigrid preconditioned conjugate gradient method, which minimizes the memory requirements and increases overall run-time performance. Computations are localized mostly near the interface location to reduce computational cost without sacrificing the accuracy. Findings The proposed method allows to capture interface topology accurately in simulating wide range of flow regimes with high density/viscosity ratios and offers good mass conservation even in relatively coarse grids, while keeping the simplicity of the level set interface modeling. Efficiency, local high-order accuracy and mass conservation of the method are confirmed through distinct numerical test cases of sloshing, dam break and Rayleigh–Taylor instability. Originality/value A fully discontinuous Galerkin, high-order, adaptive method on unstructured grids is introduced where flow and interface equations are solved in discontinuous space.


Author(s):  
P. A. Beau ◽  
T. Me´nard ◽  
R. Lebas ◽  
A. Berlemont ◽  
S. Tanguy ◽  
...  

The main objective of our work is to develop direct numerical simulation tools for the primary break up of a jet. Results can help to determine closure relation in the ELSA model [1] which is based on a single-phase Eulerian model and on the transport equation for the mean liquid/gas interface density in turbulent flows. DNS simulations are carried out to obtain statistical information in the dense zone of the spray where nearly no experimental data are available. The numerical method should describe the interface motion precisely, handle jump conditions at the interface without artificial smoothing, and respect mass conservation. We develop a 3D code [2], where interface tracking is ensured by Level Set method, Ghost Fluid Method [3] is used to capture accurately sharp discontinuities, and coupling between Level Set and VOF methods is used for mass conservation [4]. Turbulent inflow boundary conditions are generated through correlated random velocities with a prescribed length scale. Specific care has been devoted to improve computing time with MPI parallelization. The numerical methods have been applied to investigate physical processes that are involved in the primary break up of an atomizing jet. The chosen configuration is close as possible of Diesel injection (Diameter D = 0.1 mm, Velocity = 100m/s, Liquid density = 696kg/m3, Gas density = 25kg/m3). Typical results will be presented. From the injector nozzle, the turbulence initiates some perturbations on the liquid surface, that are enhanced by the mean shear between the liquid jet and the surrounding air. The interface becomes very wrinkled and some break-up is initiated. The induced liquid parcels show a wide range of shapes. Statistics are carried out and results will be provided for liquid volume fraction, liquid/gas interface density, and turbulent correlations.


2018 ◽  
Author(s):  
Inderpreet Kaur ◽  
Anton Butenko ◽  
Gianni Pagnini

Abstract. Fire-spotting is often responsible for a dangerous flare up in the wildfire and causes secondary ignitions isolated from the primary fire zone leading to perilous situations. In this paper a complete physical parametrisation of fire-spotting is presented within a formulation aimed to include random processes into operational fire spread models. This formulation can be implemented into existing operational models as a post-processing scheme at each time step, without calling for any major changes in the original framework. In particular, the efficacy of this formulation has already been shown for wildfire simulators based on an Eulerian moving interface method, namely the Level Set Method (LSM) that forms the baseline of the operational software WRF-SFIRE, and for wildfire simulators based on a Lagrangian front tracking technique, namely the Discrete Event System Specification (DEVS) that forms the baseline of the operational software FOREFIRE. The simple and computationally less expensive parametrisation includes the important parameters necessary for describing the landing behavior of the firebrands. The results from different simulations with a simple model based on the LSM highlight the response of the parametrisation to varying fire intensities, wind conditions and different firebrand radii. The contribution of the firebrands towards increasing the fire perimeter varies according to different concurrent conditions and the simulation results prove to be in agreement with the physical processes. Among the many rigorous approaches available in literature to model the firebrand transport and distribution, the approach presented here proves to be simple yet versatile for application to operational fire spread models.


2019 ◽  
Vol 2019 ◽  
pp. 1-12
Author(s):  
Saeed Mesgari ◽  
Mehrdad Bazazzadeh ◽  
Alireza Mostofizadeh

This study deals with the application of optimization in Finocyl grain design with ballistic objective functions using a genetic algorithm. The classical sampling method is used for space filling; a level-set method is used for simulating the evaluation of a burning surface of the propellant grain. An algorithm is developed beside the level-set code that prepares the initial grain configuration using a computer-aided design (CAD) to export generated models to the level-set code. The lumped method is used to perform internal ballistic analysis. A meta-model is used to surrogate the level-set method in an optimization design loop. Finally, a case study is done to verify the proposed algorithm. Observed results show that the grain design method reduced the design time significantly, and this algorithm can be used in designing any grain type.


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.


2011 ◽  
Vol 4 (1) ◽  
pp. 497-545 ◽  
Author(s):  
J. Mandel ◽  
J. D. Beezley ◽  
A. K. Kochanski

Abstract. We describe the physical model, numerical algorithms, and software structure of WRF-Fire. WRF-Fire consists of a fire-spread model, implemented by the level-set method, coupled with the Weather Research and Forecasting model. In every time step, the fire model inputs the surface wind, which drives the fire, and outputs the heat flux from the fire into the atmosphere, which in turn influences the atmosphere. The level-set method allows submesh representation of the burning region and flexible implementation of various kinds of ignition. WRF-Fire is distributed as a part of WRF and it uses the WRF parallel infrastructure for parallel computing.


2018 ◽  
Vol 30 (4) ◽  
pp. 040908 ◽  
Author(s):  
H.-Z. Yuan ◽  
C. Shu ◽  
Y. Wang ◽  
S. Shu

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.


Author(s):  
Ruquan Liang ◽  
Satoru Komori

We present a numerical strategy for a propagating interface in multiphase flows using a level set method combined with a local mesh adaptative technique. We use the level set method to move the propagating interface in multiphase flows. We also use the local mesh adaptative technique to increase the grid resolution at regions near the propagating interface and additionally at the regions near points of high curvature with a minimum of additional expense. For illustration, we apply the adaptive coupled level set method to a collection of bubbles moving under passive transport. Good agreement has been obtained in the comparision of the numerical results for the collection of bubbles using an adaptative grid with those using a single grid. We also apply the adaptive coupled level set method to a droplet falling on a step where it is important to accurately model the effect of surface tension force and the motion of the free-surface, and the numerical results agree very closely with available data.


Sign in / Sign up

Export Citation Format

Share Document