A monotonicity preserving conservative sharp interface flow solver for high density ratio two-phase flows

Chenadec, Vincent Le; Pitsch, Heinz

doi:10.1016/j.jcp.2013.04.027

Cited by 57 publications

(37 citation statements)

References 32 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…Furthermore, the present approach assumes that the velocity of the two phases match at the interface, which is valid for incompressible viscous fluids without phase change (Tryggvason et al, 2011;Rudman, 1998;Le Chenadec and Pitsch, 2013).…”

Section: Governing Equationsmentioning

confidence: 99%

Multiscale simulation of atomization with small droplets represented by a Lagrangian point-particle model

Ling

Zaleski

Scardovelli

2015

International Journal of Multiphase Flow

130

107

View full text Add to dashboard Cite

a b s t r a c tModeling and simulation of atomization is challenging due to the existence of a wide range of length scales. This multiscale nature of atomization introduces a fundamental challenge to numerical simulation. A pathway to comprehensive modeling is still to be found. The present study proposes a multiscale multiphase flow model for atomization simulations, where the large-scale interfaces are resolved by the Volume-of-Fluid (VOF) method and the small droplets by the Lagrangian point-particle (LPP) model. Particular attention is focused on the momentum coupling between LPP and resolved flow and the conversion between droplets represented by VOF and LPP. A series of multiphase flow problems are considered to validate the model. The results obtained by a number of simulations are compared against direct numerical simulation (DNS) results and experimental data. In particular, the model is applied to simulate the gas-assisted atomization experiment, and the numerical results are compared to the experimental measurements for a quantitative validation.

show abstract

Section: Governing Equationsmentioning

confidence: 99%

Multiscale simulation of atomization with small droplets represented by a Lagrangian point-particle model

Ling

Zaleski

Scardovelli

2015

International Journal of Multiphase Flow

130

107

View full text Add to dashboard Cite

show abstract

“…Consistent coupling between the VOF method and momentum transport is essential for robust simulations of realistic flows [1]. Le Chenadec and Pitsch [20] used a geometric semi-Lagrangian methodology [16] to transport both the VOF and momentum to improve consistency. Their approach transports a computational cell backwards or forwards in time during a time-step and an update is computed using geometric operations.…”

Section: Introductionmentioning

confidence: 99%

“…2. In Le Chenadec and Pitsch's scheme [20] the updates for liquid volume fraction (mass) and momentum are computed at different locations and are not consistent. This inconsistency manifests as a conservation error.…”

Section: Introductionmentioning

confidence: 99%

“…The ambiguity is in the definition of the velocity or scalar which can be taken as a constant within each computational cell or more accurately computed with a higher-order reconstruction. Following the ideas of Le Chenadec and Pitsch [20], we use a second-order, conservative method with a linear reconstruction that does not alter the magnitude of the conserved quantity within each cell. For example, the x component of velocity is reconstructed within a u-momentum cell using…”

mentioning

confidence: 99%

See 1 more Smart Citation

A mass and momentum conserving unsplit semi-Lagrangian framework for simulating multiphase flows

Owkes

Desjardins

2017

Journal of Computational Physics

View full text Add to dashboard Cite

In this work, we present a computational methodology for convection and advection that handles discontinuities with second order accuracy and maintains conservation to machine precision. This method can transport a variety of discontinuous quantities and is used in the context of an incompressible gas-liquid flow to transport the phase interface, momentum, and scalars. The proposed method provides a modification to the three-dimensional, unsplit, second-order semi-Lagrangian flux method of Owkes & Desjardins (JCP, 2014). The modification adds a refined grid that provides consistent fluxes of mass and momentum on a staggered grid and discrete conservation of mass and momentum, even with flows with large density ratios. Additionally, the refined grid doubles the resolution of the interface without significantly increasing the computational cost over previous non-conservative schemes. This is possible due to a novel partitioning of the semi-Lagrangian fluxes into a small number of simplices. The proposed scheme is tested using canonical verification tests, rising bubbles, and an atomizing liquid jet.

show abstract

“…Several variants of the CLSVOF approach can be found in literature, e.g. Le Chenadec and Pitsch [15], Menard et al [16] and Wang et al. [17].…”

mentioning

confidence: 99%

Accurate and efficient surface reconstruction from volume fraction data on general meshes

Scheufler

Roenby

2019

Journal of Computational Physics

105

View full text Add to dashboard Cite

Simulations involving free surfaces and fluid interfaces are important in many areas of engineering. There is, however, still a need for improved simulation methods. Recently, a new efficient geometric VOF method called isoAdvector for general polyhedral meshes was published. We investigate the interface reconstruction step of isoAdvector, and demonstrate that especially for unstructured meshes the applied isosurface based approach can lead to noisy interface orientations. We then introduce a novel computational interface reconstruction scheme based on calculation of a reconstructed distance function (RDF). By iterating over the RDF calculation and interface reconstruction, we obtain second order convergence of both the interface normal and position within cells even with a strict L ∞ error norm. In 2D this is verified with reconstruction of a circle on Cartesian meshes and on unstructured triangular and polygonal prism meshes. In 3D the second order convergence is verified with reconstruction of a sphere on Cartesian meshes and on unstructured tetrahedral and polyhedral meshes. The new scheme is combined with the interface advection step of the isoAdvector algorithm. Significantly reduced absolute advection errors are obtained, and for CFL number 0.2 and below we demonstrate second order convergence on all the mentioned mesh types in 2D and 3D. The implementation of the proposed interface reconstruction schemes is straightforward and the computational cost is significantly reduced compared to contemporary methods. The schemes are implemented as an extension to the Computational Fluid Dynamics (CFD) Open Source software package, OpenFOAM R . The extension module and all test cases presented in this paper are released as open source. and are still used today e.g. in OpenFOAM and Fluent R , but are often not sufficiently accurate. Models assuming an infinitesimal interface thickness comprise of two steps: Interface reconstruction and interface advection. The most common reconstruction method is the piecewise linear interface reconstruction method (PLIC), where the interface within a cell is represented by a plane cutting the cell into two subcells. The resulting surface lacks C 0 continuity. Several approaches for the computation of the normal vector of the plane can be found in literature. In the widely used method by Youngs [3], the normal vector is calculated from the gradient of the volume fraction field. This method is easy to implement, also on three dimensional unstructured meshes, but lacks accuracy. In the least squares volume-of-fluid interface reconstruction algorithm (LVIRA) [4] an interface plane is projected into the neighbouring cells. The plane is used to calculate a volume fraction. The error of the calculated volume fraction in the neighbouring cells is minimized by changing the orientation of the plane. This delivers accurate results on all mesh types in three dimensions. However, the method needs a two dimensional minimizer in three dimensions which complicates the implementation and makes the m...

show abstract

A monotonicity preserving conservative sharp interface flow solver for high density ratio two-phase flows

Cited by 57 publications

References 32 publications

Multiscale simulation of atomization with small droplets represented by a Lagrangian point-particle model

Multiscale simulation of atomization with small droplets represented by a Lagrangian point-particle model

A mass and momentum conserving unsplit semi-Lagrangian framework for simulating multiphase flows

Accurate and efficient surface reconstruction from volume fraction data on general meshes

Contact Info

Product

Resources

About