An extension of DG methods for hyperbolic problems to one-dimensional semi-infinite domains

DOI10.1016/J.AMC.2018.12.057MaRDI QIDQ2009259zbMATH OpenOpenAlexWikidataFDO

Authors Tommaso Benacchio, Luca Bonaventura

Publication date 27 November 2019

Published in Applied Mathematics and Computation (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1803.10997

zbMATH Keywords

discontinuous Galerkin methodsdomain decompositionhyperbolic equationsspectral methodsLaguerre functionsopen boundary conditions

Mathematics Subject Classification ID

Initial-boundary value problems for first-order hyperbolic systems (35L50) Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs (65M60) Spectral, collocation and related methods for initial value and initial-boundary value problems involving PDEs (65M70)

Abstract: We consider spectral discretizations of hyperbolic problems on unbounded domains using Laguerre basis functions. Taking as model problem the scalar advection equation, we perform a comprehensive stability analysis that includes strong collocation formulations, nodal and modal weak formulations, with either inflow or outflow boundary conditions, using either Gauss - Laguerre or Gauss - Laguerre - Radau quadrature nodes and based on either scaled Laguerre functions or scaled Laguerre polynomials. We show that some of these combinations give rise to intrinsically unstable schemes, while the combination of scaled Laguerre functions with Gauss - Laguerre - Radau nodes appears to be stable for both strong and weak formulations. We then show how a modal discretization approach for hyperbolic systems on an unbounded domain can be naturally and seamlessly coupled to a discontinuous finite element discretization on a finite domain. Examples of one dimensional hyperbolic systems are solved with the proposed domain decomposition technique. The errors obtained with the proposed approach are found to be small, enabling the use of the coupled scheme for the simulation of Rayleigh damping layers in the semi-infinite part. Energy errors and reflection ratios of the scheme in absorbing wavetrains and single Gaussian signals show that a small number of nodes in the semi-infinite domain are sufficient to damp the waves. The theoretical insight and numerical results corroborate previous findings by the authors and establish the scaled Laguerre functions-based discretization as a flexible and efficient tool for absorbing layers as well as for the accurate simulation of waves in unbounded regions.

Recommendations

• DG method for solving multiple solutions of one-dimensional semilinear differential equations
• Superconvergence of DG method for one-dimensional singularly perturbed problems
• General DG-methods for highly indefinite Helmholtz problems
• Discontinuous Galerkin methods for solving a hyperbolic inequality
• A Survey of One-Step Splitting Methods for Semi-Discrete First Order Hyperbolic Partial Differential Equations
• A seamless, extended DG approach for advection-diffusion problems on unbounded domains
• Semidiscrete schemes of the finite element method for degenerate hyperbolic equations
• A survey of DG methods for elliptic problems
• Analysis of a DG method for singularly perturbed convection-diffusion problems
• Some properties of the discontinuous Galerkin method for one-dimensional singularly perturbed problems

Cites work

• scientific article; zbMATH DE number 1552918 (Why is no real title available?)
• A coupled Legendre-Laguerre spectral-element method for the Navier-Stokes equations in unbounded domains
• A high-order accurate discontinuous finite element method for the numerical solution of the compressible Navier-Stokes equations
• A high-order super-grid-scale absorbing layer and its application to linear hyperbolic systems
• A perfectly matched layer for the absorption of electromagnetic waves
• A positive finite-difference advection scheme
• A review of infinite element methods for exterior Helmholtz problems
• ASYMPTOTIC DERIVATION OF THE SECTION-AVERAGED SHALLOW WATER EQUATIONS FOR NATURAL RIVER HYDRAULICS
• Absorbing Boundary Conditions for the Numerical Simulation of Waves
• Absorbing boundary conditions: a spectral collocation approach
• An experimental adaptation of Higdon-type non-reflecting boundary conditions to linear first-order systems
• Application of high‐order Higdon non‐reflecting boundary conditions to linear shallow water models
• Approximation of radiation boundary conditions
• High-order non-reflecting boundary conditions for the linearized 2-D Euler equations: no mean flow case
• Mixed finite elements for numerical weather prediction
• Numerical solution of problems on unbounded domains. A review
• On the analysis and construction of perfectly matched layers for the linearized Euler equations
• Perfectly matched layers for convex truncated domains with discontinuous Galerkin time domain simulations
• Pseudospectral method using generalized Laguerre functions for singular problems on unbounded domains
• Some recent advances on spectral methods for unbounded domains
• Stable and Efficient Spectral Methods in Unbounded Domains Using Laguerre Functions
• TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws. III: One-dimensional systems
• Well-posed perfectly matched layers for advective acoustics

Cited in

(2)

• Efficient hyperbolic-parabolic models on multi-dimensional unbounded domains using an extended DG approach
• A seamless, extended DG approach for advection-diffusion problems on unbounded domains

This page was built for publication: An extension of DG methods for hyperbolic problems to one-dimensional semi-infinite domains

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2009259)