Local mesh refinement is one of the key steps in the implementations of adaptive finite element methods. This paper presents a parallel algorithm for distributed memory parallel computers for adaptive local refinement...Local mesh refinement is one of the key steps in the implementations of adaptive finite element methods. This paper presents a parallel algorithm for distributed memory parallel computers for adaptive local refinement of tetrahedral meshes using bisection. This algorithm is used in PHG, Parallel Hierarchical Grid Chttp：//lsec. cc. ac. cn/phg/）, a toolbox under active development for parallel adaptive finite element solutions of partial differential equations. The algorithm proposed is characterized by allowing simukaneous refinement of submeshes to arbitrary levels before synchronization between submeshes and without the need of a central coordinator process for managing new vertices. Using the concept of canonical refinement, a simple proof of the independence of the resulting mesh on the mesh partitioning is given, which is useful in better understanding the behaviour of the biseetioning refinement procedure.展开更多
We construct and analyze a family of quadratic finite volume method(FVM)schemes over tetrahedral meshes.In order to prove the stability and the error estimate,we propose the minimum V-angle condition on tetrahedral me...We construct and analyze a family of quadratic finite volume method(FVM)schemes over tetrahedral meshes.In order to prove the stability and the error estimate,we propose the minimum V-angle condition on tetrahedral meshes,and the surface and volume orthogonal conditions on dual meshes.Through the technique of element analysis,the local stability is equivalent to a positive definiteness of a 9 × 9 element matrix,which is difficult to analyze directly or even numerically.With the help of the surface orthogonal condition and the congruent transformation,this element matrix is reduced into a block diagonal matrix,and then we carry out the stability result under the minimum V-angle condition.It is worth mentioning that the minimum V-angle condition of the tetrahedral case is very different from a simple extension of the minimum angle condition for triangular meshes,while it is also convenient to use in practice.Based on the stability,we prove the optimal H^(1) and L^(2) error estimates,respectively,where the orthogonal conditions play an important role in ensuring the optimal L^(2) convergence rate.Numerical experiments are presented to illustrate our theoretical results.展开更多
We extend the weighted essentially non-oscillatory(WENO)schemes on two dimensional triangular meshes developed in[7]to three dimensions,and construct a third order finite volume WENO scheme on three dimensional tetrah...We extend the weighted essentially non-oscillatory(WENO)schemes on two dimensional triangular meshes developed in[7]to three dimensions,and construct a third order finite volume WENO scheme on three dimensional tetrahedral meshes.We use the Lax-Friedrichs monotone flux as building blocks,third order reconstructions made from combinations of linear polynomials which are constructed on diversified small stencils of a tetrahedral mesh,and non-linear weights using smoothness indicators based on the derivatives of these linear polynomials.Numerical examples are given to demonstrate stability and accuracy of the scheme.展开更多
We construct a nonlinear monotone finite volume scheme for threedimensional diffusion equation on tetrahedral meshes.Since it is crucial important to eliminate the vertex unknowns in the construction of the scheme,we ...We construct a nonlinear monotone finite volume scheme for threedimensional diffusion equation on tetrahedral meshes.Since it is crucial important to eliminate the vertex unknowns in the construction of the scheme,we present a new efficient eliminating method.The scheme has only cell-centered unknowns and can deal with discontinuous or tensor diffusion coefficient problems on distorted meshes rigorously.The numerical results illustrate that the resulting scheme can preserve positivity on distorted tetrahedral meshes,and also show that our scheme appears to be approximate second-order accuracy for solution.展开更多
This work introduces a scalable and efficient topological structure for tetrahedral and hexahedral meshes. The design of the data structure aims at maximal flexibility and high performance. It provides a high scalabil...This work introduces a scalable and efficient topological structure for tetrahedral and hexahedral meshes. The design of the data structure aims at maximal flexibility and high performance. It provides a high scalability by using hierarchical representa-tions of topological elements. The proposed data structure is array-based, and it is a compact representation of the half-edge data structure for volume elements and half-face data structure for volumetric meshes. This guarantees constant access time to the neighbors of the topological elements. In addition, an open-source implementation named Open Volumetric Mesh （OVM） of the pro-posed data structure is written in C＋＋ using generic programming concepts.展开更多
For the numerical simulation of compressible flows,normally different mesh sizes are expected in different regions.For example,smaller mesh sizes are required to improve the local numerical resolution in the regions w...For the numerical simulation of compressible flows,normally different mesh sizes are expected in different regions.For example,smaller mesh sizes are required to improve the local numerical resolution in the regions where the physical variables vary violently(for example,near the shock waves or in the boundary layers)and larger elements are expected for the regions where the solution is smooth.h-adaptive mesh has been widely used for complex flows.However,there are two difficulties when employing h-adaptivity for high-order discontinuous Galerkin(DG)methods.First,locally curved elements are required to precisely match the solid boundary,which significantly increases the difficulty to conduct the"refining"and"coarsening"operations since the curved information has to be maintained.Second,h-adaptivity could break the partition balancing,which would significantly affect the efficiency of parallel computing.In this paper,a robust and automatic h-adaptive method is developed for high-order DG methods on locally curved tetrahedral mesh,for which the curved geometries are maintained during the h-adaptivity.Furthermore,the reallocating and rebalancing of the computational loads on parallel clusters are conducted to maintain the parallel efficiency.Numerical results indicate that the introduced h-adaptive method is able to generate more reasonable mesh according to the structure of flow-fields.展开更多
A general framework for the development of high-order compact schemes has been proposed recently.The core steps of the schemes are composed of the following.1).Based on a kinetic model equation,from a generalized init...A general framework for the development of high-order compact schemes has been proposed recently.The core steps of the schemes are composed of the following.1).Based on a kinetic model equation,from a generalized initial distribution of flow variables construct a time-accurate evolution solution of gas distribution function at a cell interface and obtain the corresponding flux function;2).Introduce the WENO-type weighting functions into the high-order time-derivative of the cell interface flux function in the multistage multi-derivative(MSMD)time stepping scheme to cope with the possible impingement of a shock wave on a cell interface within a time step,and update the cell-averaged conservative flow variables inside each control volume;3).Model the time evolution of the gas distribution function on both sides of a cell interface separately,take moments of the inner cell interface gas distribution function to get flow variables,and update the cell-averaged gradients of flow variables inside each control volume;4).Based on the cell-averaged flow variables and their gradients,develop compact initial data reconstruction to get initial condition of flow distributions at the beginning of next time step.A compact gas-kinetic scheme(GKS)up to sixth-order accuracy in space and fourth-order in time has been constructed on 2D unstructured mesh.In this paper,the compact GKS up to fourth-order accuracy on three-dimensional tetrahedral mesh will be further constructed with the focus on the WENO-type initial compact data reconstruction.Nonlinear weights are designed to achieve high-order accuracy for the smooth Navier-Stokes solution and keep super robustness in 3D computation with strong shock interactions.The fourth-order compact GKS uses a large time step with a CFL number 0.6 in the simulations from subsonic to hypersonic flow.A series of test cases are used to validate the scheme.The high-order compact GKS can be used in 3D applications with complex geometry.展开更多
A novel method for boundary constrained tetrahedral mesh generation is proposed based on Advancing Front Technique(AFT)and conforming Delaunay triangulation.Given a triangulated surface mesh,AFT is firstly applied to ...A novel method for boundary constrained tetrahedral mesh generation is proposed based on Advancing Front Technique(AFT)and conforming Delaunay triangulation.Given a triangulated surface mesh,AFT is firstly applied to mesh several layers of elements adjacent to the boundary.The rest of the domain is then meshed by the conforming Delaunay triangulation.The non-conformal interface between two parts of meshes are adjusted.Mesh refinement and mesh optimization are then preformed to obtain a more reasonable-sized mesh with better quality.Robustness and quality of the proposed method is shown.Convergence proof of each stage as well as the whole algorithm is provided.Various numerical examples are included as well as the quality of the meshes.展开更多
Free-moving simulations of airplanes, submarines and other automobiles under extreme and emergency conditions are becoming increasingly important from operational and tactical perspectives. Such simulations are fairly...Free-moving simulations of airplanes, submarines and other automobiles under extreme and emergency conditions are becoming increasingly important from operational and tactical perspectives. Such simulations are fairly challenging due to the extreme unsteady motions and high Re(Reynolds) numbers. The aim of this study is to perform a six-DOF motion simulation of a 6:1prolate spheroid that is falling in a fluid field. Prior to conducting the six-DOF simulation, some verification simulations were performed. First, a laminar flow past an inclined prolate spheroid at a Re number of 1000 and incidence angle of 45. with a tetrahedral mesh was simulated to verify the relevant targeted discrete method for an unstructured mesh. Second, to verify the LES(large eddy simulation) models and dependent parameters for the DDES(delayed detached eddy simulation), a turbulent flow past a sphere was performed at a subcritical Re number of 10000. Third, a steady maneuvering problem about a prolate spheroid pitching up from 0. to 30. incidence at a uniform angular velocity was established based on a dynamic tetrahedral mesh with changing topology and the ALE(arbitrary Lagrangian-Eulerian) method of fluid-structure coupling at a Re number of 4.2 × 10~6.Finally, two six-DOF motions of an inclined 6:1 prolate spheroid at an initial incidence of 45. were simulated at different Re numbers of 10000 and 4.2 × 10~6.展开更多
基金supported by the 973 Program of China 2005CB321702China NSF 10531080.
文摘Local mesh refinement is one of the key steps in the implementations of adaptive finite element methods. This paper presents a parallel algorithm for distributed memory parallel computers for adaptive local refinement of tetrahedral meshes using bisection. This algorithm is used in PHG, Parallel Hierarchical Grid Chttp：//lsec. cc. ac. cn/phg/）, a toolbox under active development for parallel adaptive finite element solutions of partial differential equations. The algorithm proposed is characterized by allowing simukaneous refinement of submeshes to arbitrary levels before synchronization between submeshes and without the need of a central coordinator process for managing new vertices. Using the concept of canonical refinement, a simple proof of the independence of the resulting mesh on the mesh partitioning is given, which is useful in better understanding the behaviour of the biseetioning refinement procedure.
基金supported by National Natural Science Foundation of China(Grant Nos.12071177 and 11701211)the Science Challenge Project(Grant No.TZ2016002)the China Postdoctoral Science Foundation(Grant No.2021M690437)。
文摘We construct and analyze a family of quadratic finite volume method(FVM)schemes over tetrahedral meshes.In order to prove the stability and the error estimate,we propose the minimum V-angle condition on tetrahedral meshes,and the surface and volume orthogonal conditions on dual meshes.Through the technique of element analysis,the local stability is equivalent to a positive definiteness of a 9 × 9 element matrix,which is difficult to analyze directly or even numerically.With the help of the surface orthogonal condition and the congruent transformation,this element matrix is reduced into a block diagonal matrix,and then we carry out the stability result under the minimum V-angle condition.It is worth mentioning that the minimum V-angle condition of the tetrahedral case is very different from a simple extension of the minimum angle condition for triangular meshes,while it is also convenient to use in practice.Based on the stability,we prove the optimal H^(1) and L^(2) error estimates,respectively,where the orthogonal conditions play an important role in ensuring the optimal L^(2) convergence rate.Numerical experiments are presented to illustrate our theoretical results.
基金The research of the second author is supported by NSF grants AST-0506734 and DMS-0510345.
文摘We extend the weighted essentially non-oscillatory(WENO)schemes on two dimensional triangular meshes developed in[7]to three dimensions,and construct a third order finite volume WENO scheme on three dimensional tetrahedral meshes.We use the Lax-Friedrichs monotone flux as building blocks,third order reconstructions made from combinations of linear polynomials which are constructed on diversified small stencils of a tetrahedral mesh,and non-linear weights using smoothness indicators based on the derivatives of these linear polynomials.Numerical examples are given to demonstrate stability and accuracy of the scheme.
基金The authors thank two reviewers for their numerous constructive comments and suggestions that helped to improve the paper significantly.This work is partially supported by NSAF(No.U1430101)the National Natural Science Foundation of China(91330106,11571047,11571048,11401034)+2 种基金China Postdoctoral Science Foundation(20110490328)the natural science foundation of Shandong Province(ZR2012AM019,ZR2013AM023,ZR2014AM013)Independent innovation foundation of Shandong University(2012TS018).
文摘We construct a nonlinear monotone finite volume scheme for threedimensional diffusion equation on tetrahedral meshes.Since it is crucial important to eliminate the vertex unknowns in the construction of the scheme,we present a new efficient eliminating method.The scheme has only cell-centered unknowns and can deal with discontinuous or tensor diffusion coefficient problems on distorted meshes rigorously.The numerical results illustrate that the resulting scheme can preserve positivity on distorted tetrahedral meshes,and also show that our scheme appears to be approximate second-order accuracy for solution.
基金Supported by Fundamental Research Funds for the Central Universities(Nos.2013ZM087,2012zz0062,2012zz0063)Doctoral Fund of Ministry of Education of China(No.20130172120010)
文摘This work introduces a scalable and efficient topological structure for tetrahedral and hexahedral meshes. The design of the data structure aims at maximal flexibility and high performance. It provides a high scalability by using hierarchical representa-tions of topological elements. The proposed data structure is array-based, and it is a compact representation of the half-edge data structure for volume elements and half-face data structure for volumetric meshes. This guarantees constant access time to the neighbors of the topological elements. In addition, an open-source implementation named Open Volumetric Mesh （OVM） of the pro-posed data structure is written in C＋＋ using generic programming concepts.
基金supported by the funding of the Key Laboratory of Aerodynamic Noise Control(No.ANCL20190103)the State Key Laboratory of Aerodynamics(No.SKLA20180102)+1 种基金the Aeronautical Science Foundation of China(Nos.2018ZA52002,2019ZA052011)the National Natural Science Foundation of China(Nos.61672281,61732006)。
文摘For the numerical simulation of compressible flows,normally different mesh sizes are expected in different regions.For example,smaller mesh sizes are required to improve the local numerical resolution in the regions where the physical variables vary violently(for example,near the shock waves or in the boundary layers)and larger elements are expected for the regions where the solution is smooth.h-adaptive mesh has been widely used for complex flows.However,there are two difficulties when employing h-adaptivity for high-order discontinuous Galerkin(DG)methods.First,locally curved elements are required to precisely match the solid boundary,which significantly increases the difficulty to conduct the"refining"and"coarsening"operations since the curved information has to be maintained.Second,h-adaptivity could break the partition balancing,which would significantly affect the efficiency of parallel computing.In this paper,a robust and automatic h-adaptive method is developed for high-order DG methods on locally curved tetrahedral mesh,for which the curved geometries are maintained during the h-adaptivity.Furthermore,the reallocating and rebalancing of the computational loads on parallel clusters are conducted to maintain the parallel efficiency.Numerical results indicate that the introduced h-adaptive method is able to generate more reasonable mesh according to the structure of flow-fields.
基金the National Natural Science Foundation of China(No.12172316)Hong Kong research grant council 16208021 and 16301222CORE as a joint research centre for ocean research between QNLM and HKUST through the project QNLM20SC01-A and QNLM20SC01-E.
文摘A general framework for the development of high-order compact schemes has been proposed recently.The core steps of the schemes are composed of the following.1).Based on a kinetic model equation,from a generalized initial distribution of flow variables construct a time-accurate evolution solution of gas distribution function at a cell interface and obtain the corresponding flux function;2).Introduce the WENO-type weighting functions into the high-order time-derivative of the cell interface flux function in the multistage multi-derivative(MSMD)time stepping scheme to cope with the possible impingement of a shock wave on a cell interface within a time step,and update the cell-averaged conservative flow variables inside each control volume;3).Model the time evolution of the gas distribution function on both sides of a cell interface separately,take moments of the inner cell interface gas distribution function to get flow variables,and update the cell-averaged gradients of flow variables inside each control volume;4).Based on the cell-averaged flow variables and their gradients,develop compact initial data reconstruction to get initial condition of flow distributions at the beginning of next time step.A compact gas-kinetic scheme(GKS)up to sixth-order accuracy in space and fourth-order in time has been constructed on 2D unstructured mesh.In this paper,the compact GKS up to fourth-order accuracy on three-dimensional tetrahedral mesh will be further constructed with the focus on the WENO-type initial compact data reconstruction.Nonlinear weights are designed to achieve high-order accuracy for the smooth Navier-Stokes solution and keep super robustness in 3D computation with strong shock interactions.The fourth-order compact GKS uses a large time step with a CFL number 0.6 in the simulations from subsonic to hypersonic flow.A series of test cases are used to validate the scheme.The high-order compact GKS can be used in 3D applications with complex geometry.
基金Singapore MOE ARC 29/07 T207B2202,MOE RG 59/08 M52110092,NRF 2007IDM-IDM 002-010Natural Science Foundation of China 10971226 and 91130013,973 Program of China 2009CB723800the foundation of State Key Laboratory of Aerodynamics.
文摘A novel method for boundary constrained tetrahedral mesh generation is proposed based on Advancing Front Technique(AFT)and conforming Delaunay triangulation.Given a triangulated surface mesh,AFT is firstly applied to mesh several layers of elements adjacent to the boundary.The rest of the domain is then meshed by the conforming Delaunay triangulation.The non-conformal interface between two parts of meshes are adjusted.Mesh refinement and mesh optimization are then preformed to obtain a more reasonable-sized mesh with better quality.Robustness and quality of the proposed method is shown.Convergence proof of each stage as well as the whole algorithm is provided.Various numerical examples are included as well as the quality of the meshes.
基金supported by the National Natural Science Founation of China(Grant No.11572350)
文摘Free-moving simulations of airplanes, submarines and other automobiles under extreme and emergency conditions are becoming increasingly important from operational and tactical perspectives. Such simulations are fairly challenging due to the extreme unsteady motions and high Re(Reynolds) numbers. The aim of this study is to perform a six-DOF motion simulation of a 6:1prolate spheroid that is falling in a fluid field. Prior to conducting the six-DOF simulation, some verification simulations were performed. First, a laminar flow past an inclined prolate spheroid at a Re number of 1000 and incidence angle of 45. with a tetrahedral mesh was simulated to verify the relevant targeted discrete method for an unstructured mesh. Second, to verify the LES(large eddy simulation) models and dependent parameters for the DDES(delayed detached eddy simulation), a turbulent flow past a sphere was performed at a subcritical Re number of 10000. Third, a steady maneuvering problem about a prolate spheroid pitching up from 0. to 30. incidence at a uniform angular velocity was established based on a dynamic tetrahedral mesh with changing topology and the ALE(arbitrary Lagrangian-Eulerian) method of fluid-structure coupling at a Re number of 4.2 × 10~6.Finally, two six-DOF motions of an inclined 6:1 prolate spheroid at an initial incidence of 45. were simulated at different Re numbers of 10000 and 4.2 × 10~6.
