Invariants Preserving Time-Implicit Local Discontinuous Galerkin Schemes for High-Order Nonlinear Wave Equations

Abstract

This paper presents a novel approach to preserving invariants in the time-implicit numerical discretization of high-order nonlinear wave equations. Many high-order nonlinear wave equations have an infinite number of conserved quantities. Designing time-implicit numerical schemes that preserve many conserved quantities simultaneously is challenging. The proposed method utilizes Lagrange multipliers to reformulate the local discontinuous Galerkin (LDG) discretization as a conservative discretization. Combining an implicit spectral deferred correction method for time discretization can achieve a high-order scheme in both temporal and spatial dimensions. We use the Korteweg-de Vries equations as an example to illustrate the implementation of the algorithm. The algorithm can be easily generalized to other nonlinear wave problems with conserved quantities. Numerical examples for various high-order nonlinear wave equations demonstrate the effectiveness of the proposed methods in preserving invariants while maintaining the high accuracy.

Keywords

Local discontinuous Galerkin (LDG) method / Conservative and dissipative schemes / Korteweg-de Vries type equations / Spectral deferred correction (SDC) method / Lagrange multiplier / 65M60 / 35Q53

References

[1] Bihlo A, Valiquette F. Symmetry-preserving finite element schemes: an introductory investigation. SIAM J. Sci. Comput., 2019, 41(5): A3300-A3325

[2] Chen, Y., Dong, B., Pereira, R.: A new conservative discontinuous Galerkin method via implicit penalization for the generalized Korteweg-de Vries equation. SIAM J. Numer. Anal. 60(6), 3078–3098 (2022)

[3] Cockburn, B., Hou, S., Shu, C.-W.: The Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws. IV. The multidimensional case. Math. Comput. 54(190), 545–581 (1990)

[4] Cockburn B, Lin SY, Shu C-W. TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws. III. One-dimensional systems. J. Comput. Phys., 1989, 84(1): 90-113

[5] Cockburn, B., Shu, C.-W.: TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws. II. General framework. Math. Comput. 52(186), 411–435 (1989)

[6] Cockburn B, Shu C-W. The local discontinuous Galerkin method for time-dependent convection-diffusion systems. SIAM J. Numer. Anal., 1998, 35(6): 2440-2463

[7] Cockburn, B., Shu, C.-W.: The Runge-Kutta discontinuous Galerkin method for conservation laws. V. Multidimensional systems. J. Comput. Phys 141(2), 199–224 (1998)

[8] Cui, Y., Mao, D.: Numerical method satisfying the first two conservation laws for the Korteweg-de Vries equation. J. Comput. Phys. 227(1), 376–399 (2007)

[9] Deuflhard, P.: Newton methods for nonlinear problems: affine invariance and adaptive algorithms. In: Springer Series in Computational Mathematics, vol. 35. Springer, Heidelberg (2011)

[10] Dey B, Khare A, Nagaraja Kumar C. Stationary solitons of the fifth order KdV-type. Equations and their stabilization. Phys. Lett. A, 1996, 223(6): 449-452

[11] Dutt A, Greengard L, Rokhlin V. Spectral deferred correction methods for ordinary differential equations. BIT Numer. Math., 2000, 40(2): 241-266

[12] Guo, R., Xia, Y., Xu, Y.: Semi-implicit spectral deferred correction methods for highly nonlinear partial differential equations. J. Comput. Phys. 338, 269–284 (2017)

[13] Guo, R., Xu, Y.: Fast solver for the local discontinuous Galerkin discretization of the KdV type equations. Commun. Comput. Phys. 17(2), 424–457 (2015)

[14] Guo, R., Xu, Y.: Semi-implicit spectral deferred correction method based on the invariant energy quadratization approach for phase field problems. Commun. Comput. Phys. 26(1), 87–113 (2019)

[15] Hairer, E., Lubich, C., Wanner, G.: Geometric numerical integration. In: Springer Series in Computational Mathematics, vol. 31, 2nd edn. Springer, Berlin (2006)

[16] Karaduman G, Yang M, Li R-C. A least squares approach for saddle point problems. Jpn. J. Ind. Appl. Math., 2023, 40(1): 95-107

[17] Karpman VI. Stabilization of soliton instabilities by higher order dispersion: KdV-type equations. Phys. Lett. A, 1996, 210(1): 77-84

[18] Kim P. Invariantization of numerical schemes using moving frames. BIT Numer. Math., 2007, 47(3): 525-546

[19] Korteweg DJ, de Vries G. On the change of form of long waves advancing in a rectangular canal, and on a new type of long stationary waves. Philos. Mag. (5), 1895, 39(240): 422-443

[20] Kuznetsov EA, Rubenchik AM, Zakharov VE. Soliton stability in plasmas and hydrodynamics. Phys. Rep., 1986, 142(3): 103-165

[21] Levy D, Shu C-W, Yan J. Local discontinuous Galerkin methods for nonlinear dispersive equations. J. Comput. Phys., 2004, 196(2): 751-772

[22] Li, D., Li, X., Zhang, Z.: Linearly implicit and high-order energy-preserving relaxation schemes for highly oscillatory Hamiltonian systems. J. Comput. Phys. 477, 111925 (2023)

[23] Liu, H., Yi, N.: A Hamiltonian preserving discontinuous Galerkin method for the generalized Korteweg-de Vries equation. J. Comput. Phys. 321, 776–796 (2016)

[24] McLachlan RI, Quispel GRW, Robidoux N. Geometric integration using discrete gradients. R. Soc. Lond. Philos. Trans. Ser. A Math. Phys. Eng. Sci., 1999, 357(1754): 1021-1045

[25] Miura, R.M., Gardner, C.S., Kruskal, M.D.: Korteweg-de Vries equation and generalizations. II. Existence of conservation laws and constants of motion. J. Math. Phys. 9, 1204–1209 (1968)

[26] Montoison A, Orban D. GPMR: an iterative method for unsymmetric partitioned linear systems. SIAM J. Matrix Anal. Appl., 2023, 44(1): 293-311

[27] Nouri, F.Z., Sloan, D.M.: A comparison of Fourier pseudospectral methods for the solution of the Korteweg-de Vries equation. J. Comput. Phys. 83(2), 324–344 (1989)

[28] Nowak, U., Weimann, L.: A family of Newton codes for systems of highly nonlinear equations. Technical Report TR-91-10, ZIB, Takustr. 7, 14195 Berlin (1992)

[29] Quispel GRW, McLaren DI. A new class of energy-preserving numerical integration methods. J. Phys. A, 2008, 41(4): 045206

[30] Reed, W.H., Hill, T.R.: Triangular mesh methods for the neutron transport equation. Technical Report, Los Alamos Scientific Lab., N. Mex., USA (1973)

[31] Rosenau P, Hyman JM. Compactons: solitons with finite wavelength. Phys. Rev. Lett., 1993, 70: 564-567

[32] Simo JC, Tarnow N, Wong KK. Exact energy-momentum conserving algorithms and symplectic schemes for nonlinear dynamics. Comput. Methods Appl. Mech. Eng., 1992, 100(1): 63-116

[33] van der Vegt JJW, Xia Y, Xu Y. Positivity preserving limiters for time-implicit higher order accurate discontinuous Galerkin discretizations. SIAM J. Sci. Comput., 2019, 41(3): A2037-A2063

[34] Xia Y, Xu Y, Shu C-W. Efficient time discretization for local discontinuous Galerkin methods. Discrete Contin. Dyn. Syst. Ser. B, 2007, 8(3): 677-693

[35] Xu, Y., Shu, C.-W.: Local discontinuous Galerkin methods for three classes of nonlinear wave equations. J. Comput. Math. 22, 250–274 (2004)

[36] Xu Y, Shu C-W. Error estimates of the semi-discrete local discontinuous Galerkin method for nonlinear convection-diffusion and KdV equations. Comput. Methods Appl. Mech. Eng., 2007, 196(37/38/39/40): 3805-3822

[37] Yin X, Cao W. A class of efficient Hamiltonian conservative spectral methods for Korteweg-de Vries equations. J. Sci. Comput., 2023, 94(1): 10-29

[38] Zabusky NJ, Kruskal MD. Interaction of solitons in a collisionless plasma and the recurrence of initial states. Phys. Rev. Lett., 1965, 15: 240-243

[39] Zhang Q, Xia Y. Conservative and dissipative local discontinuous Galerkin methods for Korteweg-de Vries type equations. Commun. Comput. Phys., 2019, 25(2): 532-563

[40] Zhao PF, Qin MZ. Multisymplectic geometry and multisymplectic Preissmann scheme for the KdV equation. J. Phys. A, 2000, 33(18): 3613-3626

Funding

National Natural Science Foundation of China(12071455)

RIGHTS & PERMISSIONS

Shanghai University