In this paper,a DG(Discontinuous Galerkin)method which has been widely employed in CFD(Computational Fluid Dynamics)is used to solve the twodimensional time-domain Maxwell’s equations for complex geometries on unstru...In this paper,a DG(Discontinuous Galerkin)method which has been widely employed in CFD(Computational Fluid Dynamics)is used to solve the twodimensional time-domain Maxwell’s equations for complex geometries on unstructured mesh.The element interfaces on solid boundary are treated in both curved way and straight way.Numerical tests are performed for both benchmark problems and complex cases with varying orders on a series of grids,where the high-order convergence in accuracy can be observed.Both the curved and the straight solid boundary implementation can give accurate RCS(Radar Cross-Section)results with sufficiently small mesh size,but the curved solid boundary implementation can significantly improve the accuracy when using relatively large mesh size.More importantly,this CFD-based high-order DG method for the Maxwell’s equations is very suitable for complex geometries.展开更多
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.展开更多
基金National Natural Science Foundation of China:No.11272152 and Aeronautical Science Foundation of China:No.20101552018.
文摘In this paper,a DG(Discontinuous Galerkin)method which has been widely employed in CFD(Computational Fluid Dynamics)is used to solve the twodimensional time-domain Maxwell’s equations for complex geometries on unstructured mesh.The element interfaces on solid boundary are treated in both curved way and straight way.Numerical tests are performed for both benchmark problems and complex cases with varying orders on a series of grids,where the high-order convergence in accuracy can be observed.Both the curved and the straight solid boundary implementation can give accurate RCS(Radar Cross-Section)results with sufficiently small mesh size,but the curved solid boundary implementation can significantly improve the accuracy when using relatively large mesh size.More importantly,this CFD-based high-order DG method for the Maxwell’s equations is very suitable for complex geometries.
基金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.