We propose a high-order discontinuous Galerkin scheme for nonlinear acoustic waves on polytopic meshes. To model sound propagation with losses through homogeneous media, we use Westervelt’s nonlinear wave equation with strong damping. Challenges in the numerical analysis lie in handling the nonlinearity in the model, which involves the derivatives in time of the acoustic velocity potential, and in preventing the equation from degenerating. We rely in our approach on the Banach fixed-point theorem combined with a stability and convergence analysis of a linear wave equation with a variable coefficient in front of the second time derivative. By doing so, we derive an a priori error estimate for Westervelt’s equation in a suitable energy norm for the polynomial degree p ≥ 2. Numerical experiments carried out in two-dimensional settings illustrate the theoretical convergence results. In addition, we demonstrate efficiency of the method in a three- dimensional domain with varying medium parameters, where we use the discontinuous Galerkin approach in a hybrid way.

A high-order discontinuous Galerkin method for nonlinear sound waves

Antonietti, Paola F.;Mazzieri, Ilario;
2020

Abstract

We propose a high-order discontinuous Galerkin scheme for nonlinear acoustic waves on polytopic meshes. To model sound propagation with losses through homogeneous media, we use Westervelt’s nonlinear wave equation with strong damping. Challenges in the numerical analysis lie in handling the nonlinearity in the model, which involves the derivatives in time of the acoustic velocity potential, and in preventing the equation from degenerating. We rely in our approach on the Banach fixed-point theorem combined with a stability and convergence analysis of a linear wave equation with a variable coefficient in front of the second time derivative. By doing so, we derive an a priori error estimate for Westervelt’s equation in a suitable energy norm for the polynomial degree p ≥ 2. Numerical experiments carried out in two-dimensional settings illustrate the theoretical convergence results. In addition, we demonstrate efficiency of the method in a three- dimensional domain with varying medium parameters, where we use the discontinuous Galerkin approach in a hybrid way.
File in questo prodotto:
File Dimensione Formato  
2020_Antonietti_Mazzieri_MUHR_NIKOLIC_WOHLMUTH_JCP.pdf

Accesso riservato

Descrizione: Articolo principale
: Publisher’s version
Dimensione 1.66 MB
Formato Adobe PDF
1.66 MB 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: http://hdl.handle.net/11311/1136030
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 3
  • ???jsp.display-item.citation.isi??? 4
social impact