skip to main content

Title: Central discontinuous Galerkin methods on overlapping meshes for wave equations
In this paper, we study the central discontinuous Galerkin (DG) method on overlapping meshes for second order wave equations. We consider the first order hyperbolic system, which is equivalent to the second order scalar equation, and construct the corresponding central DG scheme. We then provide the stability analysis and the optimal error estimates for the proposed central DG scheme for one- and multi-dimensional cases with piecewise P k elements. The optimal error estimates are valid for uniform Cartesian meshes and polynomials of arbitrary degree k  ≥ 0. In particular, we adopt the techniques in Liu et al . ( SIAM J. Numer. Anal. 56 (2018) 520–541; ESAIM: M2AN 54 (2020) 705–726) and obtain the local projection that is crucial in deriving the optimal order of convergence. The construction of the projection here is more challenging since the unknowns are highly coupled in the proposed scheme. Dispersion analysis is performed on the proposed scheme for one dimensional problems, indicating that the numerical solution with P 1 elements reaches its minimum with a suitable parameter in the dissipation term. Several numerical examples including accuracy tests and long time simulation are presented to validate the theoretical results.
Authors:
; ; ;
Award ID(s):
2010107 1719410
Publication Date:
NSF-PAR ID:
10225997
Journal Name:
ESAIM: Mathematical Modelling and Numerical Analysis
Volume:
55
Issue:
1
Page Range or eLocation-ID:
329 to 356
ISSN:
0764-583X
Sponsoring Org:
National Science Foundation
More Like this
  1. In this paper, we construct, analyze, and numerically validate conservative discontinuous Galerkin (DG) schemes for approximating the Schr\"{o}dinger-Poisson equation. The proposed schemes all satisfy both mass and energy conservation. For the semi-discrete DG scheme optimal $L^2$ error estimates are obtained. Efficient iterative solvers are also constructed to solve the second order implicit time discretization. A number of numerical tests are presented to demonstrate the method’s accuracy and robustness, confirming that both mass and energy are well preserved over long time simulations.
  2. In this paper, we propose a local discontinuous Galerkin (LDG) method for nonlinear and possibly degenerate parabolic stochastic partial differential equations, which is a high-order numerical scheme. It extends the discontinuous Galerkin (DG) method for purely hyperbolic equations to parabolic equations and shares with the DG method its advantage and flexibility. We prove the L 2 -stability of the numerical scheme for fully nonlinear equations. Optimal error estimates ( O ( h (k+1) )) for smooth solutions of semi-linear stochastic equations is shown if polynomials of degree k are used. We use an explicit derivative-free order 1.5 time discretization scheme to solve the matrix-valued stochastic ordinary differential equations derived from the spatial discretization. Numerical examples are given to display the performance of the LDG method.
  3. In this paper, we study the optimal error estimates of the classical discontinuous Galerkin method for time-dependent 2-D hyperbolic equations using P k elements on uniform Cartesian meshes, and prove that the error in the L 2 norm achieves optimal ( k  + 1)th order convergence when upwind fluxes are used. For the linear constant coefficient case, the results hold true for arbitrary piecewise polynomials of degree k  ≥ 0. For variable coefficient and nonlinear cases, we give the proof for piecewise polynomials of degree k  = 0, 1, 2, 3 and k  = 2, 3, respectively, under the condition that the wind direction does not change. The theoretical results are verified by numerical examples.
  4. In this paper, we study superconvergence properties of the ultraweak-local discontinuous Galerkin (UWLDG) method in Tao et al. [To appear in Math. Comput. DOI: https://doi.org/10.1090/mcom/3562 (2020).] for an one-dimensional linear fourth-order equation. With special initial discretizations, we prove the numerical solution of the semi-discrete UWLDG scheme superconverges to a special projection of the exact solution. The order of this superconvergence is proved to be k + min(3, k ) when piecewise ℙ k polynomials with k ≥ 2 are used. We also prove a 2 k -th order superconvergence rate for the cell averages and for the function values and derivatives of the UWLDG approximation at cell boundaries. Moreover, we prove superconvergence of ( k + 2)-th and ( k + 1)-th order of the function values and the first order derivatives of the UWLDG solution at a class of special quadrature points, respectively. Our proof is valid for arbitrary non-uniform regular meshes and for arbitrary k ≥ 2. Numerical experiments verify that all theoretical findings are sharp.
  5. Purpose The purpose of this paper is as follows: to significantly reduce the computation time (by a factor of 1,000 and more) compared to known numerical techniques for real-world problems with complex interfaces; and to simplify the solution by using trivial unfitted Cartesian meshes (no need in complicated mesh generators for complex geometry). Design/methodology/approach This study extends the recently developed optimal local truncation error method (OLTEM) for the Poisson equation with constant coefficients to a much more general case of discontinuous coefficients that can be applied to domains with different material properties (e.g. different inclusions, multi-material structural components, etc.). This study develops OLTEM using compact 9-point and 25-point stencils that are similar to those for linear and quadratic finite elements. In contrast to finite elements and other known numerical techniques for interface problems with conformed and unfitted meshes, OLTEM with 9-point and 25-point stencils and unfitted Cartesian meshes provides the 3-rd and 11-th order of accuracy for irregular interfaces, respectively; i.e. a huge increase in accuracy by eight orders for the new 'quadratic' elements compared to known techniques at similar computational costs. There are no unknowns on interfaces between different materials; the structure of the global discrete system is themore »same for homogeneous and heterogeneous materials (the difference in the values of the stencil coefficients). The calculation of the unknown stencil coefficients is based on the minimization of the local truncation error of the stencil equations and yields the optimal order of accuracy of OLTEM at a given stencil width. The numerical results with irregular interfaces show that at the same number of degrees of freedom, OLTEM with the 9-points stencils is even more accurate than the 4-th order finite elements; OLTEM with the 25-points stencils is much more accurate than the 7-th order finite elements with much wider stencils and conformed meshes. Findings The significant increase in accuracy for OLTEM by one order for 'linear' elements and by 8 orders for 'quadratic' elements compared to that for known techniques. This will lead to a huge reduction in the computation time for the problems with complex irregular interfaces. The use of trivial unfitted Cartesian meshes significantly simplifies the solution and reduces the time for the data preparation (no need in complicated mesh generators for complex geometry). Originality/value It has been never seen in the literature such a huge increase in accuracy for the proposed technique compared to existing methods. Due to a high accuracy, the proposed technique will allow the direct solution of multiscale problems without the scale separation.« less