High Order Schemes on Three-Dimensional General Polyhedral Meshes — Application to the Level Set Method

2012 ◽  
Vol 12 (1) ◽  
pp. 1-41 ◽  
Author(s):  
Thibault Pringuey ◽  
R. Stewart Cant

AbstractIn this article, we detail the methodology developed to construct arbitrarily high order schemes — linear and WENO — on 3D mixed-element unstructured meshes made up of general convex polyhedral elements. The approach is tailored specifically for the solution of scalar level set equations for application to incompressible two-phase flow problems. The construction of WENO schemes on 3D unstructured meshes is notoriously difficult, as it involves a much higher level of complexity than 2D approaches. This due to the multiplicity of geometrical considerations introduced by the extra dimension, especially on mixed-element meshes. Therefore, we have specifically developed a number of algorithms to handle mixed-element meshes composed of convex polyhedra with convex polygonal faces. The contribution of this work concerns several areas of interest: the formulation of an improved methodology in 3D, the minimisation of computational runtime in the implementation through the maximum use of pre-processing operations, the generation of novel methods to handle complex 3D mixed-element meshes and finally the application of the method to the transport of a scalar level set.

2008 ◽  
Vol 11 (4-6) ◽  
pp. 221-235 ◽  
Author(s):  
S. P. van der Pijl ◽  
A. Segal ◽  
C. Vuik ◽  
P. Wesseling

2019 ◽  
Vol 142 (1) ◽  
Author(s):  
Jianjian Xin ◽  
Fulong Shi ◽  
Qiu Jin ◽  
Lin Ma

Abstract A three-dimensional (3D) gradient-augmented level set (GALS) two-phase flow model with a pretreated reinitialization procedure is developed to simulate violent sloshing in a cuboid tank. Based on a two-dimensional (2D) GALS method, 3D Hermite, and 3D Lagrange polynomial schemes are derived to interpolate the level set function and the velocity field at arbitrary positions over a cell, respectively. A reinitialization procedure is performed on a 3D narrow band to treat the strongly distorted interface and improve computational efficiency. In addition, an identification-correction technique is proposed and incorporated into the reinitialization procedure to treat the tiny droplet which can distort the free surface shape, even lead to computation failure. To validate the accuracy of the present GALS method and the effectiveness of the proposed identification-correction technique, a 3D velocity advection case is first simulated. The present method is validated to have better mass conservation property than the classical level set and original GALS methods. Also, distorted and thin interfaces are well captured on all grid resolutions by the present GALS method. Then, sloshing under coupled surge and sway excitation, sloshing under rotational excitation are simulated. Good agreements are obtained when the present wave and pressure results are compared with the experimental and numerical results. In addition, the highly nonlinear free surface is observed, and the relationship between the excitation frequency and the impulsive pressure is investigated.


2011 ◽  
Vol 314-316 ◽  
pp. 364-368 ◽  
Author(s):  
Jun Zhang ◽  
Jian Xin Zhou ◽  
Ming Yuan Zhang ◽  
Sheng Yong Pang ◽  
Dun Ming Liao ◽  
...  

A three-dimensional incompressible two phase flow model of vertical centrifugal casting is proposed to simulate the fluid flow of mould filling process accurately and effectively. The Projection method is adopted to solve the govern equation of the flow field, and the Level Set method is used to capture the free surface. The mold filling of a complex part with thin-wall is simulated. The numerical result shows that the Projection-Level Set method could simulate centrifugal casting effectively. The present study has a guiding significance to the production of vertical centrifugal casting.


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.


2015 ◽  
Vol 17 (3) ◽  
pp. 615-656 ◽  
Author(s):  
Marc R. J. Charest ◽  
Clinton P. T. Groth ◽  
Pierre Q. Gauthier

AbstractHigh-order discretization techniques offer the potential to significantly reduce the computational costs necessary to obtain accurate predictions when compared to lower-order methods. However, efficient and universally-applicable high-order discretizations remain somewhat illusive, especially for more arbitrary unstructured meshes and for incompressible/low-speed flows. A novel, high-order, central essentially non-oscillatory (CENO), cell-centered, finite-volume scheme is proposed for the solution of the conservation equations of viscous, incompressible flows on three-dimensional unstructured meshes. Similar to finite element methods, coordinate transformations are used to maintain the scheme’s order of accuracy even when dealing with arbitrarily-shaped cells having non-planar faces. The proposed scheme is applied to the pseudo-compressibility formulation of the steady and unsteady Navier-Stokes equations and the resulting discretized equations are solved with a parallel implicit Newton-Krylov algorithm. For unsteady flows, a dual-time stepping approach is adopted and the resulting temporal derivatives are discretized using the family of high-order backward difference formulas (BDF). The proposed finite-volume scheme for fully unstructured mesh is demonstrated to provide both fast and accurate solutions for steady and unsteady viscous flows.


Sign in / Sign up

Export Citation Format

Share Document