A 3D Unsplit Forward/Backward Volume-of-Fluid Approach and Coupling to the Level Set Method

Chenadec, Vincent Le; Pitsch, Heinz

doi:10.1016/j.jcp.2012.07.019

Cited by 61 publications

(37 citation statements)

References 43 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…As an example, on the sequential simulation of the Richtmyer-Meshkov instability [19] with an unstructured tetrahedral mesh of 250K cells, our VOF solver represents 22% of the computing costs. A similar percentage was reported by Le Chenadec and Pitsch [20], on the solution of a diesel jet with a Cartesian mesh of 256×256×1152 cells. Anyway, beyond the percentage obtained for any particular simulation, it is a certainty that, in the high performance computing context, the cost of the VOF calculations will become more and more important while the algorithmic solutions adopted disregard parallel performance issues.…”

Section: Introductionsupporting

confidence: 83%

Parallel load balancing strategy for Volume-of-Fluid methods on 3-D unstructured meshes

Jofre

Borrell

Lehmkuhl

et al. 2015

Journal of Computational Physics

View full text Add to dashboard Cite

UPCommonsPortal del coneixement obert de la UPC http://upcommons.upc.edu/e-prints Abstract: Volume-of-Fluid (VOF) is one of the methods of choice to reproduce the interface motion in the simulation of multi-fluid flows. One of its main strengths is its accuracy in capturing sharp interface geometries, although requiring for it a number of geometric calculations. Under these circumstances, achieving parallel performance on current supercomputers is a must. The main obstacle for the parallelization is that the computing costs are concentrated only in the discrete elements that lie on the interface between fluids. Consequently, if the interface is not homogeneously distributed throughout the domain, standard domain decomposition (DD) strategies lead to imbalanced workload distributions. In this paper, we present a new parallelization strategy for general unstructured VOF solvers, based on a dynamic load balancing process complementary to the underlying DD. Its parallel efficiency has been analyzed and compared to the DD one using up to 1024 CPU-cores on an Intel SandyBridge based supercomputer. The results obtained on the solution of several artificially generated test cases show a speedup of up to 12x with respect to the standard DD, depending on the interface size, the initial distribution and the number of parallel processes engaged. Moreover, the new parallelization strategy presented is of general purpose, therefore, it could be used to parallelize any VOF solver without requiring changes on the coupled flow solver. Finally, note that although designed for the VOF method, our approach could be easily adapted to other interface-capturing methods, such as the Level-Set, which may present similar workload imbalances. SUMMARY OF CHANGESAccordingly to the comments raised by the reviewer, the changes on the manuscript are:1. As suggested by Reviewer #1: (1) Eq. 3 is valid only for divergence-free velocity fields, hence, this condition has been introduced in the text; (2) the limits of the sum in the numerator of Eqs. 9 -11 were incorrect, therefore, they have been corrected to start at 0 and end at P-1; (3) the description of Fig. 10 has been improved.2. The reviewer observed that Alg. 4 did not scale in terms of memory, since vector WI gathered information from all processes and, thus, grew with the size of the problem and number of CPU-cores. This may had become a bottleneck on larger scale tests. Therefore, we have developed a new version of the algorithm overcoming such limitations. The new version is based on non-blocking point-to-point communications.The new algorithm has substituted the previous one on the manuscript, and all tests for the LB strategy have been repeated and the results updated. Since Alg. 4 represents only 5% of the algorithm time, its upgrading has provided rather small improvements on the results. Even so, the important point is that the memory bottleneck has been eliminated.3. Furthermore, we have corrected some minor errors and we have tried to improve the wording across all the manuscrip...

show abstract

Section: Introductionsupporting

confidence: 83%

Parallel load balancing strategy for Volume-of-Fluid methods on 3-D unstructured meshes

Jofre

Borrell

Lehmkuhl

et al. 2015

Journal of Computational Physics

View full text Add to dashboard Cite

show abstract

“…More recently, Le Chenadec and Pitsch [127] propose an Hybrid Lagrangian-Eulerian Method for multiphase flows (HyLEM). The method solves two-phase fluid flows with a combination of VOF and an ALE-like step: the mesh is deformed in a Lagrangian way for advecting the volume fraction, and then it is remapped to the target (original) mesh by using a strategy based on the LS function.…”

Section: Fixed Discretization Methodsmentioning

confidence: 99%

Numerical Modeling and Experimental Validation of Free Surface Flow Problems

Cruchaga

Battaglia

Storti

et al. 2014

Arch Computat Methods Eng

View full text Add to dashboard Cite

In this paper we present a summary of numerical methods for solving free surface and two fluid flow problems. We will focus the attention on level set formulations extensively used in the context of the finite element method. In particular, numerical developments to achieve accurate solutions are described. Specific topics of the algorithms, like mass preservation and interface redefinition, are evaluated. To illustrate these aspects, numerical strategies previously developed are applied to the solution of a sloshing and a water column collapse problems. To assess the capabilities of these techniques, the numerical results are compared against each other and with experimental data. Considering these aspects, the present work is aimed to outline some well reported aspects of level set-like formulations.Fil: Cruchaga, Marcela. Universidad de Santiago de Chile; ChileFil: Battaglia, Laura. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Santa Fe. Centro de Investigaciones En Metodos Computacionales. Universidad Nacional del Litoral. Centro de Investigaciones En Metodos Computacionales; ArgentinaFil: Storti, Mario Alberto. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Santa Fe. Centro de Investigaciones En Metodos Computacionales. Universidad Nacional del Litoral. Centro de Investigaciones En Metodos Computacionales; ArgentinaFil: D'elia, Jorge. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Santa Fe. Centro de Investigaciones En Metodos Computacionales. Universidad Nacional del Litoral. Centro de Investigaciones En Metodos Computacionales; Argentin

show abstract

“…The convex join operation is necessary to reduce the convex set to a more optimal decomposition of the node volume. Note also that geometric VOF advection algorithms (see for example [9,[31][32][33]) also generate polyhedra of known convex decompositions, obviating the need for an expensive and complicated convexdecomposition algorithm. To test the accuracy of the EHF method, it is also necessary to design a second-order accurate method to initialize F from a signed-distance function.…”

Section: Geometric Toolsmentioning

confidence: 99%

Accurate interface normal and curvature estimates on three-dimensional unstructured non-convex polyhedral meshes

Ivey

Moin

2015

Journal of Computational Physics

View full text Add to dashboard Cite

This paper presents a framework for extending the height-function technique for the calculation of interface normals and curvatures to unstructured non-convex polyhedral meshes with application to the piecewise-linear interface calculation volume-of-fluid method. The methodology is developed with reference to a collocated node-based finitevolume two-phase flow solver that utilizes the median-dual mesh, requiring a set of data structures and algorithms for non-convex polyhedral operations: truncation of a polyhedron by a plane, intersection of two polyhedra, joining of two convex polyhedra, volume enforcement of a polyhedron by a plane, and volume fraction initialization by a signed-distance function. By leveraging these geometric tools, a geometric interpolation strategy for embedding structured height-function stencils in unstructured meshes is developed. The embedded height-function technique is tested on surfaces with known interface normals and curvatures, namely cylinder, sphere, and ellipsoid. Tests are performed on the median duals of a uniform cartesian mesh, a wedge mesh, and a tetrahedral mesh, and comparisons are made with conventional methods. Across the tests, the embedded height-function technique outperforms contemporary methods and its accuracy approaches the accuracy that the traditional height-function technique exemplifies on uniform cartesian meshes.

show abstract

A 3D Unsplit Forward/Backward Volume-of-Fluid Approach and Coupling to the Level Set Method

Cited by 61 publications

References 43 publications

Parallel load balancing strategy for Volume-of-Fluid methods on 3-D unstructured meshes

Parallel load balancing strategy for Volume-of-Fluid methods on 3-D unstructured meshes

Numerical Modeling and Experimental Validation of Free Surface Flow Problems

Accurate interface normal and curvature estimates on three-dimensional unstructured non-convex polyhedral meshes

Contact Info

Product

Resources

About