A nodal discontinuous Galerkin formulation based on Lagrange polynomials basis is used to simulate the acoustic wave propagation. Its dispersion and dissipation properties for the advection equation are investigated b...A nodal discontinuous Galerkin formulation based on Lagrange polynomials basis is used to simulate the acoustic wave propagation. Its dispersion and dissipation properties for the advection equation are investigated by utilizing an eigenvalue analysis. Two test problems of wave propagation with initial disturbance consisting of a Gaussian profile or rectangular pulse are performed. And the performance of the schemes in short,intermediate,and long waves is evaluated. Moreover,numerical results between the nodal discontinuous Galerkin method and finite difference type schemes are compared,which indicate that the numerical solution obtained using nodal discontinuous Galerkin method with a pure central flux has obviously high frequency oscillations for initial disturbance consisting of a rectangular pulse,which is the same as those obtained using finite difference type schemes without artificial selective damping. When an upwind flux is adopted,spurious waves are eliminated effectively except for the location of discontinuities. When a limiter is used,the spurious short waves are almost completely removed. Therefore,the quality of the computed solution has improved.展开更多
In this paper,Nodal discontinuous Galerkin method is presented to approxi-mate Time-domain Lorentz model equations in meta-materials.The upwind flux is cho-sen in spatial discrete scheme.Low-storage five-stage fourth-...In this paper,Nodal discontinuous Galerkin method is presented to approxi-mate Time-domain Lorentz model equations in meta-materials.The upwind flux is cho-sen in spatial discrete scheme.Low-storage five-stage fourth-order explicit Runge-Kutta method is employed in time discrete scheme.An error estimate of accuracy O(τ^(4)+h^(n))is proved under the L^(2)-norm,specially O(τ^(4)+h^(n+1))can be obtained.Numerical exper-iments for transverse electric(TE)case and transverse magnetic(TM)case are demon-strated to verify the stability and the efficiency of the method in low and higher wave frequency.展开更多
In this short note we present a derivation of the Spectral Difference Scheme from a Discontinuous Galerkin(DG)discretization of a nonlinear conservation law.This allows interpretation of the Spectral Difference Scheme...In this short note we present a derivation of the Spectral Difference Scheme from a Discontinuous Galerkin(DG)discretization of a nonlinear conservation law.This allows interpretation of the Spectral Difference Scheme as a particular discretization under the quadrature-free nodal DG paradigm.Moreover,it enables identification of the key differences between the Spectral Difference Scheme and standard nodal DG schemes.展开更多
The present study develops implicit physical domain-based discontinuous Galerkin(DG)methods for efficient scale-resolving simulations on mixed-curved meshes.Implicit methods are essential to handle stiff systems in ma...The present study develops implicit physical domain-based discontinuous Galerkin(DG)methods for efficient scale-resolving simulations on mixed-curved meshes.Implicit methods are essential to handle stiff systems in many scale-resolving simulations of interests in computational science and engineering.The physical domain-based DGmethod can achieve high-order accuracy using the optimal bases set and preserve the required accuracy on non-affinemeshes.When using the quadraturebased DG method,these advantages are overshadowed by severe computational costs on mixed-curved meshes,making implicit scale-resolving simulations unaffordable.To address this issue,the quadrature-free direct reconstruction method(DRM)is extended to the implicit DG method.In this approach,the generalized reconstruction approximates non-linear flux functions directly in the physical domain,making the computing-intensive numerical integrations precomputable at a preprocessing step.The DRM operator is applied to the residual computation in the matrix-free method.The DRM operator can be further extended to the system matrix computation for the matrix-explicit Krylov subspace method and preconditioning.Finally,the A-stable Rosenbrock-type Runge–Kutta methods are adopted to achieve high-order accuracy in time.Extensive verification and validation from the manufactured solution to implicit large eddy simulations are conducted.The computed results confirm that the proposed method significantly improves computational efficiency compared to the quadrature-based method while accurately resolving detailed unsteady flow features that are hardly captured by scale-modeled simulations.展开更多
采用数值通量的方式进行场量交互是时域非连续伽略金(Discontinuous Galerkin Time Domain,DGTD)算法区别于时域有限元(Finite Element Time Domain,FETD)方法的主要方面.从二维TM情形弱解方程出发,讨论了当前三角形单元和相邻单元进行...采用数值通量的方式进行场量交互是时域非连续伽略金(Discontinuous Galerkin Time Domain,DGTD)算法区别于时域有限元(Finite Element Time Domain,FETD)方法的主要方面.从二维TM情形弱解方程出发,讨论了当前三角形单元和相邻单元进行场量交互时数值通量物理意义和不同形式.结合数值通量和弱解方程得到了DGTD算法的迭代计算式.给出了线元辐射和双线元干涉的数值算例,算例结果表明了文中方法的正确性.展开更多
基金Supported by the National Natural Science Foundation of China(51106099,50976072)the Leading Academic Discipline Project of Shanghai Municipal Education Commission(J50501)
文摘A nodal discontinuous Galerkin formulation based on Lagrange polynomials basis is used to simulate the acoustic wave propagation. Its dispersion and dissipation properties for the advection equation are investigated by utilizing an eigenvalue analysis. Two test problems of wave propagation with initial disturbance consisting of a Gaussian profile or rectangular pulse are performed. And the performance of the schemes in short,intermediate,and long waves is evaluated. Moreover,numerical results between the nodal discontinuous Galerkin method and finite difference type schemes are compared,which indicate that the numerical solution obtained using nodal discontinuous Galerkin method with a pure central flux has obviously high frequency oscillations for initial disturbance consisting of a rectangular pulse,which is the same as those obtained using finite difference type schemes without artificial selective damping. When an upwind flux is adopted,spurious waves are eliminated effectively except for the location of discontinuities. When a limiter is used,the spurious short waves are almost completely removed. Therefore,the quality of the computed solution has improved.
基金supported by NSFC.China(NOs.11201501,11571389)the Program for Innovation Research in Central University of Finance and Economics+1 种基金The second author is Supported by NSFC.China(Grant Nos.11471296,11101384)the third author is supported in part by Defense Industrial Technology Development Program(B1520133015).
文摘In this paper,Nodal discontinuous Galerkin method is presented to approxi-mate Time-domain Lorentz model equations in meta-materials.The upwind flux is cho-sen in spatial discrete scheme.Low-storage five-stage fourth-order explicit Runge-Kutta method is employed in time discrete scheme.An error estimate of accuracy O(τ^(4)+h^(n))is proved under the L^(2)-norm,specially O(τ^(4)+h^(n+1))can be obtained.Numerical exper-iments for transverse electric(TE)case and transverse magnetic(TM)case are demon-strated to verify the stability and the efficiency of the method in low and higher wave frequency.
基金Financial support from the Deutsche Forschungsgemeinschaft(German Research Association)through grant GSC 111the Air Force Office of Scientific Research,Air Force Materiel Command,USAF,under grant number FA8655-08-1-3060,is gratefully acknowledged。
文摘In this short note we present a derivation of the Spectral Difference Scheme from a Discontinuous Galerkin(DG)discretization of a nonlinear conservation law.This allows interpretation of the Spectral Difference Scheme as a particular discretization under the quadrature-free nodal DG paradigm.Moreover,it enables identification of the key differences between the Spectral Difference Scheme and standard nodal DG schemes.
基金the financial support provided by the Defense Acquisition Program Administration(DAPA)under Grant UD200046CD(Data-driven Flow Modeling Research Laboratory)the Korea Research Institute for defense Technology planning and advancement(KRIT)under Grant KRIT-CT-22-030(Reusable Unmanned Space Vehicle Research Center,2023)supported by the program of the National Research Foundation of Korea(NRF-2021R1A2C2008348).
文摘The present study develops implicit physical domain-based discontinuous Galerkin(DG)methods for efficient scale-resolving simulations on mixed-curved meshes.Implicit methods are essential to handle stiff systems in many scale-resolving simulations of interests in computational science and engineering.The physical domain-based DGmethod can achieve high-order accuracy using the optimal bases set and preserve the required accuracy on non-affinemeshes.When using the quadraturebased DG method,these advantages are overshadowed by severe computational costs on mixed-curved meshes,making implicit scale-resolving simulations unaffordable.To address this issue,the quadrature-free direct reconstruction method(DRM)is extended to the implicit DG method.In this approach,the generalized reconstruction approximates non-linear flux functions directly in the physical domain,making the computing-intensive numerical integrations precomputable at a preprocessing step.The DRM operator is applied to the residual computation in the matrix-free method.The DRM operator can be further extended to the system matrix computation for the matrix-explicit Krylov subspace method and preconditioning.Finally,the A-stable Rosenbrock-type Runge–Kutta methods are adopted to achieve high-order accuracy in time.Extensive verification and validation from the manufactured solution to implicit large eddy simulations are conducted.The computed results confirm that the proposed method significantly improves computational efficiency compared to the quadrature-based method while accurately resolving detailed unsteady flow features that are hardly captured by scale-modeled simulations.
文摘采用数值通量的方式进行场量交互是时域非连续伽略金(Discontinuous Galerkin Time Domain,DGTD)算法区别于时域有限元(Finite Element Time Domain,FETD)方法的主要方面.从二维TM情形弱解方程出发,讨论了当前三角形单元和相邻单元进行场量交互时数值通量物理意义和不同形式.结合数值通量和弱解方程得到了DGTD算法的迭代计算式.给出了线元辐射和双线元干涉的数值算例,算例结果表明了文中方法的正确性.
