A discontinuous Galerkin method and its error estimate for nonlinear fourth-order wave equations
Introduction
In recent years, many numerical methods have been defined and analyzed for the wave equations [1], [2], [3], [4], [5], [6], [7], [8], [9]. The nonlinear fourth-order wave equations arise commonly from the studies of vibration of beams and thin plates [10]. In this paper, we are interested in the numerical methods for a class of nonlinear fourth-order wave equations [11], [12], [13], [14], [15], [16], [17], [18], [19], with the initial conditions For the sake of simplicity, we only consider the periodic boundary condition. The solution , is a real-valued function, and the initial conditions and are assumed to be as smooth as necessary. Levandosky [16] proved that our problem (1.1) admits a unique local solution for nonlinearity which satisfies where for and for , denotes the critical exponent for the embedding of into , .
There are many numerical methods proposed in the literature for solving the fourth-order equations [20], [21], [22], [23], [24], [25], [26]. In [20], Achouri designed a second-order conservative finite difference scheme for the two-dimensional fourth-order nonlinear wave equation. The mixed finite elements for the fourth-order wave equations also have been studied by He et al. in [25]. They considered mixed finite element method with explicit and implicit discretization in time and derived the optimal error estimate in the norm. In [22], [23], Baccouch applied the local discontinuous Galerkin (LDG) method for the fourth-order Euler–Bernoulli partial differential equation (PDE) in one dimension, including superconvergence analysis and a posterior error estimate.
We consider an ultraweak-local discontinuous Galerkin (UWLDG) method introduced in [27] for (1.1). The DG method is a class of finite element methods using completely discontinuous basis functions. The first DG method was introduced in 1973 by Reed and Hill [28] in the framework of neutron transport. It was later developed for time-dependent nonlinear hyperbolic conservation laws, coupled with the Runge–Kutta time discretization, by Cockburn et al. [29], [30], [31]. Since then, the DG method has been intensively studied and successfully applied to various problems in a wide range of applications due to its flexibility with meshing, its compactness and its high parallel efficiency. The UWLDG method is a discontinuous Galerkin method designed for PDEs with high order spatial derivatives, which combines the advantage of the LDG method and the ultra-weak DG (UWDG) method. The idea of the LDG method [7], [32], [33], [34] is to rewrite the equations with higher order spatial derivatives into a first order system, then apply the DG method to this system and design suitable numerical fluxes to ensure stability. The UWDG method [35] is based on repeated integration by parts to move all spatial derivatives to the test function in the weak formulation, and to ensure stability by carefully choosing numerical fluxes. In our method, at first, we rewrite Eq. (1.1) as a second-order system. Then we repeat the application of integration by parts, and choose suitable numerical fluxes to ensure stability. Compared to the LDG method, we introduce fewer auxiliary variables, thereby reducing memory and computational costs. Compared to the UWDG method, we do not need any internal penalty terms to ensure stability.
We define the energy where and . For Eq. (1.1) is a constant. Therefore, we would like to design a numerical method that conserves the energy . Energy conserving DG methods for wave equations have been developed in [36], [37], [38], [39], [40]. Recently, Chou et al. [36], [40] developed an optimal energy-conserving local DG method for multi-dimensional second-order wave equation in heterogeneous media. Later, Fu and Shu [37] proposed an optimal energy conserving DG method for linear symmetric hyperbolic systems on general unstructured meshes. They proved a priori optimal error estimates for the semi-discrete scheme in one dimension, and also in multi-dimensions for Cartesian meshes when using tensor-product polynomials. They also proposed an energy-conserving ultra-weak DG method for the generalized Korteweg–de Vries (KdV) equations in one dimension [38], and proved its optimal error estimate. In this work, we design an optimally convergent energy-conserving method for the nonlinear fourth-order equations. We choose the alternating fluxes, and prove that the energy is conserved both in one-dimensional and two-dimensional cases. We also prove the optimal error estimates in the -norm for the solution itself as well as for the auxiliary variables.
The organization of the paper is as follows. In Section 2, we introduce some notations and the UWLDG method. In Section 3, the energy conserving property of our scheme will be discussed. In Section 4, we will introduce some projections and give the optimal error estimates in the -norm for one-dimensional and two-dimensional cases. Time discretization will be shown in Section 5. The theoretical results are confirmed numerically in Section 6. In Section 7, we give some concluding remarks.
Section snippets
Notations
Let us introduce some notations. Throughout this paper, we adopt standard notations for the Sobolev spaces such as on the subdomain equipped with the norm . If , we omit the index ; and if , we set , ; and we use to denote the norm in .
Let denote a tessellation of with shape-regular elements , and the union of the boundary faces of elements , denoted as . We denote the diameter of by , and set
Energy conservation
In this section, we will demonstrate that the UWLDG scheme (2.3)–(2.4) conserves the discrete energy. Experience shows that the scheme conserving the discrete energy can often behave better, especially in long time simulation.
Theorem 3.1 The energy is conserved by the semi-discrete UWLDG method (2.3)–(2.4), with numerical fluxes (2.5)–(2.8) for all time.
Proof Without loss of generality, we choose the flux (2.5). In Eq. (2.3), we take the test function to be :
Error estimates
We study the optimal error estimates for the UWLDG method defined in (2.3)–(2.4) for Eq. (1.1). In Section 4.1, we introduce some projections and inequalities that will be used in our proof. In Section 4.2, we give the error estimate in the norm.
Time discretization
In this section, we consider the fully discrete method of scheme (2.3)–(2.4). We use the UWLDG method for the spacial discretization, it can be of high order accuracy. Therefore, we also would like to introduce an explicit, energy conserving, high order time stepping method. As in LDG method, the auxiliary variable in our method could be solved in terms of in an element-by-element fashion. After eliminating , we can get a linear second-order ordinary differential system as follows:
Numerical examples
In this section, we present numerical examples to verify our theoretical convergence properties of the UWLDG method.
Example 6.1 First example, we consider the linear fourth-order equations in one-dimension with the periodic boundary condition. and the initial conditions The exact solution of the problem is We implement the DG method (2.3)–(2.4) with the alternating fluxes (2.5) and use the time discrete method (5.1). The errors
Concluding remarks
In this paper, we have developed an UWLDG method for a class of nonlinear fourth-order wave equations. The UWLDG methods combine the LDG and UWDG methods for solving time-dependent PDEs with high order spatial derivatives. The numerical fluxes have been carefully chosen to make our scheme energy conserving. We have proved the optimal error estimate in the -norm for the solution itself as well as for the auxiliary variables approximating its derivatives in the semi-discrete method, and have
References (42)
- et al.
A discontinuous Galerkin method for the wave equation
Comput. Methods Appl. Mech. Engrg.
(2011) A local discontinuous Galerkin method for the second-order wave equation
Comput. Methods Appl. Mech. Engrg.
(2012)- et al.
Energy conservation issues in the numerical solution of the semilinear wave equation
Appl. Math. Comput.
(2015) - et al.
Numerical resolution of discontinuous Galerkin methods for time dependent wave equations
Comput. Methods Appl. Mech. Engrg.
(2011) - et al.
Dynamics of transversely vibrating beams using four engineering theories
J. Sound Vib.
(1999) Existence and uniqueness results for the bending of an elastic beam equation at resonance
J. Math. Anal. Appl.
(1988)- et al.
Global existence and global nonexistence of solutions of the Cauchy problem for a nonlinearly damped wave equation
J. Math. Anal. Appl.
(1998) Decay estimates for fourth order wave equations
J. Differential Equations
(1998)Conservative finite difference scheme for the nonlinear fourth-order wave equation
Appl. Math. Comput.
(2019)- et al.
An efficient method for computing eigenelements of Sturm–Liouville fourth-order boundary value problems
Appl. Math. Comput.
(2006)
Analysis of mixed finite element methods for fourth-order wave equations
Comput. Math. Appl.
The Runge–Kutta discontinuous Galerkin method for conservation laws. V. Multidimensional systems
J. Comput. Phys.
Optimal energy conserving local discontinuous Galerkin methods for second-order wave equation in heterogeneous media
J. Comput. Phys.
Optimal energy-conserving discontinuous Galerkin methods for linear symmetric hyperbolic systems
J. Comput. Phys.
An energy-conserving ultra-weak discontinuous Galerkin method for the generalized Korteweg–De Vries equation
J. Comput. Appl. Math.
Optimal discontinuous Galerkin methods for wave propagation
SIAM J. Numer. Anal.
Optimal discontinuous Galerkin methods for the acoustic wave equation in higher dimensions
SIAM J. Numer. Anal.
Discontinuous Galerkin finite element method for the wave equation
SIAM J. Numer. Anal.
Optimal error estimates of the semi-discrete local discontinuous Galerkin methods for high order wave equations
SIAM J. Numer. Anal.
A conservative discontinuous Galerkin method for nonlinear electromagnetic Schroedinger equations
SIAM J. Sci. Comput.
Nonlinear wave interactions in nonlinear nonintegrable systems
Stud. Appl. Math.
Cited by (6)
Efficient high-order physical property-preserving difference methods for nonlinear fourth-order wave equation with damping
2023, Computers and Mathematics with ApplicationsAn energy-conserving finite element method for nonlinear fourth-order wave equations
2023, Applied Numerical MathematicsDiscretization-Robust Safety Barrier of Partial Differential Equation
2023, 2023 11th International Conference on Control, Mechatronics and Automation, ICCMA 2023Superconvergence Analysis of the Ultra-Weak Local Discontinuous Galerkin Method for One Dimensional Linear Fifth Order Equations
2021, Journal of Scientific Computing
- 1
Research supported by China Scholarship Council.
- 2
Research supported by Science Challenge Project, China TZZT2019-A2.3, National Numerical Windtunnel, China Project NNW2019ZT4-B08, NSFC, China grants 11722112.
- 3
Research supported by NSF grant DMS-1719410.