Abstract
We present a new class of discontinuous Galerkin methods for the space discretization of the time-dependent Maxwell equations whose main feature is the use of time derivatives and/or time integrals in the stabilization part of their numerical traces. These numerical traces are chosen in such a way that the resulting semidiscrete schemes exactly conserve a discrete version of the energy. We introduce four model ways of achieving this and show that, when using the mid-point rule to march in time, the fully discrete schemes also conserve the discrete energy. Moreover, we propose a new three-step technique to devise fully discrete schemes of arbitrary order of accuracy which conserve the energy in time. The first step consists in transforming the semidiscrete scheme into a Hamiltonian dynamical system. The second step consists in applying a symplectic time-marching method to this dynamical system in order to guarantee that the resulting fully discrete method conserves the discrete energy in time. The third and last step consists in reversing the above-mentioned transformation to rewrite the fully discrete scheme in terms of the original variables.
Dedicated to the memory of Francisco-Javier Sayas
Funding source: National Science Foundation
Award Identifier / Grant number: 1912646
Funding source: Fondo Nacional de Desarrollo Científico y Tecnológico
Award Identifier / Grant number: 11180284
Funding statement: B. Cockburn was supported in part by the NSF through grant DMS-1912646. M. A. Sánchez was supported by FONDECYT Iniciación grant n. 11180284.
A Fully Discrete Methods: Second Choice of DG Traces
A.1 ESPRK Methods
We have that
for all (v,w)∈Vh×Wh .
As a consequence, each step of the new fully discrete DG scheme reads
for all (v,w)∈Vh×Wh , and where Fn,ih=Fnh-Δt∑i-1j=1˜bjHn,jh .
We still have to provide the definition of (E0h,F0h) . Since typically the initial value given is E0 , we must define (E0h,F0h)∈Vh×Wh at once in terms of the initial data F0 . This is done in a way similar to the previous case. Let us note that such definition must be such that F0h∈Wh satisfies the equation
A.2 Symplectic DIRK Methods
By similar arguments used in Section 4.5, we obtain the symplectic DIRK time-marching scheme for the 𝑯-𝑭 formulation as follows:
where Rn,iG is defined by (A.1).
We next rewrite the scheme in terms of the 𝑬-𝑯 formulation. First, we multiply the second equation by (aiiΔt)-1 and obtain
Observing that
we obtain
where
Note that
which implies that
Thus, the DIRK scheme for the 𝑬-𝑯 formulation reads as follows.
Perform the following computations for i=1,…,s :
Compute the right-hand-side terms for G=Eh,Hh,Fh,J ,
(A.1) Rn,iG=GnaiiΔt+i-1∑j=1aijaii(Gn,jajjΔt-Rn,jG),where Jn,i=J(tn,i) and tn,i=tn+ciΔt .
Compute the solution (En,ih,Hn,ih) of the system
1aiiΔt(μHn,ih,r)Th+(En,ih,∇×r)Th+⟨n׈Ψ(En,ih,-Fn,ih),r⟩∂Th=(μRn,iH,r)Th,1aiiΔt(ϵEn,ih,v)Th-(Hn,ih,∇×v)Th-⟨n׈Φ(Hn,ih,1aiiΔtEn,ih),v⟩∂Th=Θn,i,with Fn,ih=aiiΔt(Rn,iF-Hn,ih) and
Θn,i:=(ϵRn,iE,v)Th-⟨n׈Φ(0,Rn,iE),v⟩∂Th-(J(tn,i)aiiΔt,v)Th+(Rn,iJ,v)Th.
Advance the solution Gn+1=Gn+Δt∑si=1bikiG for G=Eh,Fh,Hh , where
kiG=Gn,i-GnaiiΔt-i-1∑j=1aijaiikjG.
B Fully Discrete Methods: Third Choice of DG Traces
B.1 ESPRK Methods
Based on the ESPRK scheme (4.1) and the semidiscrete Hamiltonian system (5.2), we obtain the fully discrete ESPRK scheme for the DG methods with the third choice of the numerical traces,
B.2 Symplectic DIRK Methods
Here, we apply the symplectic DIRK scheme summarized in Section 4.5 to the semidiscrete system (5.2). Then the full-discrete scheme can be summarized in the following steps.
Perform the following computations for i=1,…,s :
Compute the right-hand-side terms for G=Eh,Hh ,
Rn,iG=GnaiiΔt+i-1∑j=1aijaii(Gn,jajjΔt-Rn,jG).Recall that Jn,i=J(tn,i) , gn,iE=gE(tn,i) and tn,i=tn+ciΔt . Compute the solution (En,ih,Hn,ih) ,
1aiiΔt((ϵEn,ih,v)Th+⟨C11⟦En,ih⟧,⟦v⟧⟩Fh)=(Hn,ih,∇×v)Th+⟨n×({{Hn,ih}}+C¯T12⟦Hn,ih⟧),v⟩∂Th∖Γ+⟨n×Hn,ih,v⟩Γ-(Jn,i,v)Th+(ϵRn,iE,v)Th+⟨C11⟦Rn,iE⟧,⟦v⟧⟩Fh,1aiiΔt((μHn,ih,r)Th+⟨C22⟦Hn,ih⟧,⟦r⟧⟩F0h)=-(En,ih,∇×r)Th-⟨{{En,ih}}+C¯12⟦En,ih⟧,r⟩∂Th∖Γ-⟨gn,iE,r⟩Γ+(μRn,iH,r)Th+⟨C22⟦Rn,iH⟧,⟦r⟧⟩Fh.
Advance the solution Gn+1=Gn+Δt∑si=1bikiG for G=Eh,Hh , where
kiG=Gn,i-GnaiiΔt-i-1∑j=1aijaiikjG.
C Fully Discrete Methods: Fourth Choice of DG Traces
C.1 ESPRK Methods
Based on the ESPRK scheme (4.1) and the semidiscrete Hamiltonian system (6.1), we obtain the fully discrete ESPRK scheme for the DG methods with the fourth choice of the numerical traces,
C.2 Symplectic DIRK Methods
Here, we apply the symplectic DIRK scheme summarized in Section 4.5 to the semidiscrete system (6.1). Then the full-discrete scheme can be summarized in the following steps.
Perform the following computations for i=1,…,s :
Compute the right-hand-side terms for G=Eh,Fh,Hh,Ah ,
Rn,iG=GnaiiΔt+i-1∑j=1aijaii(Gn,jajjΔt-Rn,jG).Recall that Jn,i=J(tn,i) , gn,iE=gE(tn,i) and tn,i=tn+ciΔt . Compute the solution (En,ih,Fn,ih,Hn,ih,An,ih) ,
1aiiΔt(ϵEn,ih,v)Th=(Hn,ih,∇×v)Th+⟨n×({{Hn,ih}}+C¯T12⟦Hn,ih⟧),v⟩∂Th∖Γ+⟨n×Hn,ih,v⟩Γ-(Jn,i,v)Th-⟨C11⟦An,ih⟧,⟦v⟧⟩Fh+(ϵRn,iE,v)Th,1aiiΔtFn,ih=Hn,ih+Rn,iF on∂Th∖Γ,1aiiΔt(μHn,ih,r)Th=-(En,ih,∇×r)Th-⟨n×({{En,ih}}+C¯12⟦En,ih⟧),r⟩∂Th∖Γ-⟨gn,iE,r⟩Γ-⟨C22⟦Fn,ih⟧,⟦r⟧⟩F0h+(μRn,iH,r)Th,1aiiΔtAn,ih=En,ih+Rn,iA on∂Th,(An,ih)ext=gA(tn,i)×n onΓ.
Advance the solution Gn+1=Gn+Δt∑si=1bikiG for G=Eh,Fh,Hh,Ah , where
kiG=Gn,i-GnaiiΔt-i-1∑j=1aijaiikjG.
References
[1] M. Ainsworth, Dispersive and dissipative behaviour of high order discontinuousGalerkin finite element methods, J. Comput. Phys. 198 (2004), no. 1, 106–130. 10.1016/j.jcp.2004.01.004Search in Google Scholar
[2] T. S. Brown, S. Du, H. Eruslu and F.-J. Sayas, Analysis of models for viscoelastic wave propagation, Appl. Math. Nonlinear Sci. 3 (2018), no. 1, 55–96. 10.21042/AMNS.2018.1.00006Search in Google Scholar
[3] B. Cockburn, The pursuit of a dream, Francisco Javier Sayas and the HDG methods, SeMA J. 79 (2022), no. 1, 37–56. 10.1007/s40324-021-00273-ySearch in Google Scholar
[4] B. Cockburn, Z. Fu, A. Hungria, L. Ji, M. A. Sánchez and F.-J. Sayas, Stormer–Numerov HDG methods for acoustic waves, J. Sci. Comput. 75 (2018), no. 2, 597–624. 10.1007/s10915-017-0547-zSearch in Google Scholar
[5] B. Cockburn, J. Gopalakrishnan and R. Lazarov, Unified hybridization of discontinuous Galerkin, mixed, and continuous Galerkin methods for second order elliptic problems, SIAM J. Numer. Anal. 47 (2009), no. 2, 1319–1365. 10.1137/070706616Search in Google Scholar
[6] B. Cockburn, N. C. Nguyen and J. Peraire, HDG methods for hyperbolic problems, Handbook of Numerical Methods for Hyperbolic Problems, Handb. Numer. Anal. 17, Elsevier/North-Holland, Amsterdam (2016), 173–197. 10.1016/bs.hna.2016.07.001Search in Google Scholar
[7] B. Cockburn and V. Quenneville-Bélair, Uniform-in-time superconvergence of the HDG methods for the acoustic wave equation, Math. Comp. 83 (2014), no. 285, 65–85. 10.1090/S0025-5718-2013-02743-3Search in Google Scholar
[8] S. Du and F.-J. Sayas, A unified error analysis of hybridizable discontinuous Galerkin methods for the static Maxwell equations, SIAM J. Numer. Anal. 58 (2020), no. 2, 1367–1391. 10.1137/19M1290966Search in Google Scholar
[9] G. Fu and C.-W. Shu, Optimal energy-conserving discontinuous Galerkin methods for linear symmetric hyperbolic systems, J. Comput. Phys. 394 (2019), 329–363. 10.1016/j.jcp.2019.05.050Search in Google Scholar
[10] Z. Fu, L. F. Gatica and F.-J. Sayas, Algorithm 949: MATLAB tools for HDG in three dimensions, ACM Trans. Math. Software 41 (2015), no. 3, Article ID 20. 10.1145/2658992Search in Google Scholar
[11] J. Gopalakrishnan, M. Solano and F. Vargas, Dispersion analysis of HDG methods, J. Sci. Comput. 77 (2018), no. 3, 1703–1735. 10.1007/s10915-018-0781-zSearch in Google Scholar
[12] A. Hungria, D. Prada and F. J. Sayas, HDG methods for elastodynamics, Comput. Math. Appl. 74 (2017), 2671–2690. 10.1016/j.camwa.2017.08.016Search in Google Scholar
[13] N. C. Nguyen, J. Peraire and B. Cockburn, High-order implicit hybridizable discontinuous Galerkin methods for acoustics and elastodynamics, J. Comput. Phys. 230 (2011), no. 10, 3695–3718. 10.1016/j.jcp.2011.01.035Search in Google Scholar
[14] N. C. Nguyen, J. Peraire and B. Cockburn, Hybridizable discontinuous Galerkin methods for the time-harmonic Maxwell’s equations, J. Comput. Phys. 230 (2011), no. 19, 7151–7175. 10.1016/j.jcp.2011.05.018Search in Google Scholar
[15] M. A. Sánchez, C. Ciuca, N. C. Nguyen, J. Peraire and B. Cockburn, Symplectic Hamiltonian HDG methods for wave propagation phenomena, J. Comput. Phys. 350 (2017), 951–973. 10.1016/j.jcp.2017.09.010Search in Google Scholar
[16] M. A. Sánchez, B. Cockburn, N.-C. Nguyen and J. Peraire, Symplectic Hamiltonian finite element methods for linear elastodynamics, Comput. Methods Appl. Mech. Engrg. 381 (2021), Paper No. 113843. 10.1016/j.cma.2021.113843Search in Google Scholar
[17] J. M. Sanz-Serna, Symplectic Runge–Kutta and related methods: Recent results, Phys. D 60 (1992), 293–302. 10.1016/0167-2789(92)90245-ISearch in Google Scholar
[18] F.-J. Sayas, T. S. Brown and M. E. Hassell, Variational Techniques for Elliptic Partial Differential Equations: Theoretical Tools and Advanced Applications, CRC Press, Boca Raton, 2019. 10.1201/9780429507069Search in Google Scholar
[19] M. Stanglmeier, N. C. Nguyen, J. Peraire and B. Cockburn, An explicit hybridizable discontinuous Galerkin method for the acoustic wave equation, Comput. Methods Appl. Mech. Engrg. 300 (2016), 748–769. 10.1016/j.cma.2015.12.003Search in Google Scholar
© 2022 Walter de Gruyter GmbH, Berlin/Boston
Articles in the same Issue
- Frontmatter
- Numerical Analysis & No Regrets. Special Issue Dedicated to the Memory of Francisco Javier Sayas (1968–2019)
- Implicit/Explicit, BEM/FEM Coupled Scheme for Acoustic Waves with the Wave Equation in the Second Order Formulation
- Discontinuous Galerkin Methods with Time-Operators in Their Numerical Traces for Time-Dependent Electromagnetics
- A Discontinuous Galerkin Method for the Stationary Boussinesq System
- Korn’s Inequality and Eigenproblems for the Lamé Operator
- Error Estimates for FE-BE Coupling of Scattering of Waves in the Time Domain
- An Efficient Numerical Method for Time Domain Electromagnetic Wave Propagation in Co-axial Cables
- The Time Domain Linear Sampling Method for Determining the Shape of Multiple Scatterers Using Electromagnetic Waves
- A Boundary Integral Formulation and a Topological Energy-Based Method for an Inverse 3D Multiple Scattering Problem with Sound-Soft, Sound-Hard, Penetrable, and Absorbing Objects
- Afternote to “Coupling at a Distance”: Convergence Analysis and A Priori Error Estimates
- Spurious Resonances in Coupled Domain-Boundary Variational Formulations of Transmission Problems in Electromagnetism and Acoustics
- Spectral, Tensor and Domain Decomposition Methods for Fractional PDEs
Articles in the same Issue
- Frontmatter
- Numerical Analysis & No Regrets. Special Issue Dedicated to the Memory of Francisco Javier Sayas (1968–2019)
- Implicit/Explicit, BEM/FEM Coupled Scheme for Acoustic Waves with the Wave Equation in the Second Order Formulation
- Discontinuous Galerkin Methods with Time-Operators in Their Numerical Traces for Time-Dependent Electromagnetics
- A Discontinuous Galerkin Method for the Stationary Boussinesq System
- Korn’s Inequality and Eigenproblems for the Lamé Operator
- Error Estimates for FE-BE Coupling of Scattering of Waves in the Time Domain
- An Efficient Numerical Method for Time Domain Electromagnetic Wave Propagation in Co-axial Cables
- The Time Domain Linear Sampling Method for Determining the Shape of Multiple Scatterers Using Electromagnetic Waves
- A Boundary Integral Formulation and a Topological Energy-Based Method for an Inverse 3D Multiple Scattering Problem with Sound-Soft, Sound-Hard, Penetrable, and Absorbing Objects
- Afternote to “Coupling at a Distance”: Convergence Analysis and A Priori Error Estimates
- Spurious Resonances in Coupled Domain-Boundary Variational Formulations of Transmission Problems in Electromagnetism and Acoustics
- Spectral, Tensor and Domain Decomposition Methods for Fractional PDEs