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.展开更多
The space block search technology is used to determine a connected three-dimensional fracture network in polygonal shapes,i.e.,seepage paths.After triangulation on these polygons,a finite element mesh for 3D fracture ...The space block search technology is used to determine a connected three-dimensional fracture network in polygonal shapes,i.e.,seepage paths.After triangulation on these polygons,a finite element mesh for 3D fracture network seepage is obtained.Through introduction of the generalized Darcy's law,conservative equations for both fracture surface and fracture interactions are established.Combined with the boundary condition of Signorini's type,a partial differential equation(PDE) formulation is presented for the whole domain concerned.To solve this problem efficiently,an equivalent variational inequality(VI) formulation is given.With the penalized Heaviside function,a finite element procedure for unconfined seepage problem in 3D fracture network is developed.Through an example in a homogeneous rectangular dam,validity of the algorithm is verified.The analysis of an unconfined seepage problem in a complex fracture network shows that the proposed algorithm is very applicable to complex three-dimensional problems,and is effective in describing some interesting phenomenon usually encountered in practice,such as "preferential flow".展开更多
基金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.
基金supported by the National Natural Science Foundation of China(Grant No.51079110)the National Basic Research Program of China("973"Project)(Grant No.2011CB013506)
文摘The space block search technology is used to determine a connected three-dimensional fracture network in polygonal shapes,i.e.,seepage paths.After triangulation on these polygons,a finite element mesh for 3D fracture network seepage is obtained.Through introduction of the generalized Darcy's law,conservative equations for both fracture surface and fracture interactions are established.Combined with the boundary condition of Signorini's type,a partial differential equation(PDE) formulation is presented for the whole domain concerned.To solve this problem efficiently,an equivalent variational inequality(VI) formulation is given.With the penalized Heaviside function,a finite element procedure for unconfined seepage problem in 3D fracture network is developed.Through an example in a homogeneous rectangular dam,validity of the algorithm is verified.The analysis of an unconfined seepage problem in a complex fracture network shows that the proposed algorithm is very applicable to complex three-dimensional problems,and is effective in describing some interesting phenomenon usually encountered in practice,such as "preferential flow".