Combined DG-CG finite element method for the Westervelt equation

Fuente: arXiv
Salvato in:
Dettagli Bibliografici
Autori principali: Gómez, Sergio, Nikolić, Vanja
Natura: Preprint
Pubblicazione: 2024
Soggetti:
Accesso online:
Tags: Aggiungi Tag
Nessun Tag, puoi essere il primo ad aggiungerne!!
_version_ 1866911011022307328
author Gómez, Sergio
Nikolić, Vanja
author_facet Gómez, Sergio
Nikolić, Vanja
contents We propose and analyze a space-time finite element method for Westervelt's quasilinear model of ultrasound waves in second-order formulation. The method combines conforming finite element spatial discretizations with a discontinuous-continuous Galerkin time stepping. Its analysis is challenged by the fact that standard Galerkin testing approaches for wave problems do not allow for bounding the discrete energy at all times. By means of redesigned energy arguments for a linearized problem combined with Banach's fixed-point argument, we show the well-posedness of the scheme, \emph{a priori} error estimates, and robustness with respect to the strong damping parameter $δ$. Moreover, the scheme preserves the asymptotic preserving property of the continuous problem; more precisely, we prove that the discrete solutions corresponding to $δ>0$ converge, in the singular vanishing dissipation limit, to the solution of the discrete inviscid problem. We use several numerical experiments in $(2 + 1)$ dimensions to validate our theoretical results.
format Preprint
id arxiv_https___arxiv_org_abs_2412_09095
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle Combined DG-CG finite element method for the Westervelt equation
Gómez, Sergio
Nikolić, Vanja
Numerical Analysis
Analysis of PDEs
65M12, 65M60, 35L72
We propose and analyze a space-time finite element method for Westervelt's quasilinear model of ultrasound waves in second-order formulation. The method combines conforming finite element spatial discretizations with a discontinuous-continuous Galerkin time stepping. Its analysis is challenged by the fact that standard Galerkin testing approaches for wave problems do not allow for bounding the discrete energy at all times. By means of redesigned energy arguments for a linearized problem combined with Banach's fixed-point argument, we show the well-posedness of the scheme, \emph{a priori} error estimates, and robustness with respect to the strong damping parameter $δ$. Moreover, the scheme preserves the asymptotic preserving property of the continuous problem; more precisely, we prove that the discrete solutions corresponding to $δ>0$ converge, in the singular vanishing dissipation limit, to the solution of the discrete inviscid problem. We use several numerical experiments in $(2 + 1)$ dimensions to validate our theoretical results.
title Combined DG-CG finite element method for the Westervelt equation
topic Numerical Analysis
Analysis of PDEs
65M12, 65M60, 35L72
url https://arxiv.org/abs/2412.09095