We consider spectral discretizations of hyperbolic problems on unbounded domains us- ing 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 and based on either scaled Laguerre functions or scaled Laguerre polynomials. We show that some of these combinations give rise to intrinsically unstable discretizations, while the combination of scaled Laguerre functions with Gauss–Laguerre–Radau quadrature 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. An example of one dimensional hyperbolic system is solved with the proposed domain decomposition tech- nique. 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 modes 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 dis- cretization as a flexible and efficient tool for absorbing layers as well as for the accurate simulation of waves in unbounded regions.

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

T. Benacchio;L. Bonaventura
2019-01-01

Abstract

We consider spectral discretizations of hyperbolic problems on unbounded domains us- ing 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 and based on either scaled Laguerre functions or scaled Laguerre polynomials. We show that some of these combinations give rise to intrinsically unstable discretizations, while the combination of scaled Laguerre functions with Gauss–Laguerre–Radau quadrature 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. An example of one dimensional hyperbolic system is solved with the proposed domain decomposition tech- nique. 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 modes 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 dis- cretization as a flexible and efficient tool for absorbing layers as well as for the accurate simulation of waves in unbounded regions.
2019
File in questo prodotto:
File Dimensione Formato  
benacchio_bonaventura_2019.pdf

Accesso riservato

Descrizione: Articolo principale
: Publisher’s version
Dimensione 866.81 kB
Formato Adobe PDF
866.81 kB Adobe PDF   Visualizza/Apri

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11311/1073607
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 2
  • ???jsp.display-item.citation.isi??? 2
social impact