Hybrid Fourier pseudospectral/discontinuous Galerkin time-domain method for wave propagation

Research output: Contribution to journalArticleAcademicpeer-review

5 Citations (Scopus)
25 Downloads (Pure)

Abstract

The Fourier Pseudospectral time-domain (Fourier PSTD) method was shown to be an efficient way of modelling acoustic propagation problems as described by the linearized Euler equations (LEE), but is limited to real-valued frequency independent boundary conditions and predominantly staircase-like boundary shapes. This paper presents a hybrid approach to solve the LEE, coupling Fourier PSTD with a nodal Discontinuous Galerkin (DG) method. DG exhibits almost no restrictions with respect to geometrical complexity or boundary conditions. The aim of this novel method is to allow the computation of complex geometries and to be a step towards the implementation of frequency dependent boundary conditions by using the benefits of DG at the boundaries, while keeping the efficient Fourier PSTD in the bulk of the domain. The hybridization approach is based on conformal meshes to avoid spatial interpolation of the DG solutions when transferring values from DG to Fourier PSTD, while the data transfer from Fourier PSTD to DG is done utilizing spectral interpolation of the Fourier PSTD solutions. The accuracy of the hybrid approach is presented for one- and two-dimensional acoustic problems and the main sources of error are investigated. It is concluded that the hybrid methodology does not introduce significant errors compared to the Fourier PSTD stand-alone solver. An example of a cylinder scattering problem is presented and accurate results have been obtained when using the proposed approach. Finally, no instabilities were found during long-time calculation using the current hybrid methodology on a two-dimensional domain.
Original languageEnglish
Pages (from-to)416-432
Number of pages17
JournalJournal of Computational Physics
Volume348
Issue number1 November 2017
DOIs
Publication statusPublished - 31 Jul 2017

Fingerprint

Wave propagation
wave propagation
Euler equations
Boundary conditions
Interpolation
Acoustics
Galerkin methods
Data transfer
Scattering
boundary conditions
interpolation
Geometry
methodology
acoustic propagation
stairways
Galerkin method
mesh
constrictions
acoustics
geometry

Keywords

  • Hybrid time-domain method
  • Linearized Euler equations
  • Fourier pseudospectral
  • Discontinuous Galerkin
  • Acoustic propagation problems
  • Complex geometries

Cite this

@article{359e3c28d1704b93b93144d2c4e9aef4,
title = "Hybrid Fourier pseudospectral/discontinuous Galerkin time-domain method for wave propagation",
abstract = "The Fourier Pseudospectral time-domain (Fourier PSTD) method was shown to be an efficient way of modelling acoustic propagation problems as described by the linearized Euler equations (LEE), but is limited to real-valued frequency independent boundary conditions and predominantly staircase-like boundary shapes. This paper presents a hybrid approach to solve the LEE, coupling Fourier PSTD with a nodal Discontinuous Galerkin (DG) method. DG exhibits almost no restrictions with respect to geometrical complexity or boundary conditions. The aim of this novel method is to allow the computation of complex geometries and to be a step towards the implementation of frequency dependent boundary conditions by using the benefits of DG at the boundaries, while keeping the efficient Fourier PSTD in the bulk of the domain. The hybridization approach is based on conformal meshes to avoid spatial interpolation of the DG solutions when transferring values from DG to Fourier PSTD, while the data transfer from Fourier PSTD to DG is done utilizing spectral interpolation of the Fourier PSTD solutions. The accuracy of the hybrid approach is presented for one- and two-dimensional acoustic problems and the main sources of error are investigated. It is concluded that the hybrid methodology does not introduce significant errors compared to the Fourier PSTD stand-alone solver. An example of a cylinder scattering problem is presented and accurate results have been obtained when using the proposed approach. Finally, no instabilities were found during long-time calculation using the current hybrid methodology on a two-dimensional domain.",
keywords = "Hybrid time-domain method, Linearized Euler equations, Fourier pseudospectral, Discontinuous Galerkin, Acoustic propagation problems, Complex geometries",
author = "{Pag{\'a}n Mu{\~n}oz}, R. and M.C.J. Hornikx",
year = "2017",
month = "7",
day = "31",
doi = "10.1016/j.jcp.2017.07.046",
language = "English",
volume = "348",
pages = "416--432",
journal = "Journal of Computational Physics",
issn = "0021-9991",
publisher = "Academic Press Inc.",
number = "1 November 2017",

}

Hybrid Fourier pseudospectral/discontinuous Galerkin time-domain method for wave propagation. / Pagán Muñoz, R.; Hornikx, M.C.J.

In: Journal of Computational Physics, Vol. 348, No. 1 November 2017, 31.07.2017, p. 416-432.

Research output: Contribution to journalArticleAcademicpeer-review

TY - JOUR

T1 - Hybrid Fourier pseudospectral/discontinuous Galerkin time-domain method for wave propagation

AU - Pagán Muñoz, R.

AU - Hornikx, M.C.J.

PY - 2017/7/31

Y1 - 2017/7/31

N2 - The Fourier Pseudospectral time-domain (Fourier PSTD) method was shown to be an efficient way of modelling acoustic propagation problems as described by the linearized Euler equations (LEE), but is limited to real-valued frequency independent boundary conditions and predominantly staircase-like boundary shapes. This paper presents a hybrid approach to solve the LEE, coupling Fourier PSTD with a nodal Discontinuous Galerkin (DG) method. DG exhibits almost no restrictions with respect to geometrical complexity or boundary conditions. The aim of this novel method is to allow the computation of complex geometries and to be a step towards the implementation of frequency dependent boundary conditions by using the benefits of DG at the boundaries, while keeping the efficient Fourier PSTD in the bulk of the domain. The hybridization approach is based on conformal meshes to avoid spatial interpolation of the DG solutions when transferring values from DG to Fourier PSTD, while the data transfer from Fourier PSTD to DG is done utilizing spectral interpolation of the Fourier PSTD solutions. The accuracy of the hybrid approach is presented for one- and two-dimensional acoustic problems and the main sources of error are investigated. It is concluded that the hybrid methodology does not introduce significant errors compared to the Fourier PSTD stand-alone solver. An example of a cylinder scattering problem is presented and accurate results have been obtained when using the proposed approach. Finally, no instabilities were found during long-time calculation using the current hybrid methodology on a two-dimensional domain.

AB - The Fourier Pseudospectral time-domain (Fourier PSTD) method was shown to be an efficient way of modelling acoustic propagation problems as described by the linearized Euler equations (LEE), but is limited to real-valued frequency independent boundary conditions and predominantly staircase-like boundary shapes. This paper presents a hybrid approach to solve the LEE, coupling Fourier PSTD with a nodal Discontinuous Galerkin (DG) method. DG exhibits almost no restrictions with respect to geometrical complexity or boundary conditions. The aim of this novel method is to allow the computation of complex geometries and to be a step towards the implementation of frequency dependent boundary conditions by using the benefits of DG at the boundaries, while keeping the efficient Fourier PSTD in the bulk of the domain. The hybridization approach is based on conformal meshes to avoid spatial interpolation of the DG solutions when transferring values from DG to Fourier PSTD, while the data transfer from Fourier PSTD to DG is done utilizing spectral interpolation of the Fourier PSTD solutions. The accuracy of the hybrid approach is presented for one- and two-dimensional acoustic problems and the main sources of error are investigated. It is concluded that the hybrid methodology does not introduce significant errors compared to the Fourier PSTD stand-alone solver. An example of a cylinder scattering problem is presented and accurate results have been obtained when using the proposed approach. Finally, no instabilities were found during long-time calculation using the current hybrid methodology on a two-dimensional domain.

KW - Hybrid time-domain method

KW - Linearized Euler equations

KW - Fourier pseudospectral

KW - Discontinuous Galerkin

KW - Acoustic propagation problems

KW - Complex geometries

U2 - 10.1016/j.jcp.2017.07.046

DO - 10.1016/j.jcp.2017.07.046

M3 - Article

VL - 348

SP - 416

EP - 432

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

IS - 1 November 2017

ER -