Isogeometric simulation of acoustic radiation

Victoria Hernández Mederos, Eduardo Moreno Hernández, Jorge Estrada Sarlabous, Isidro A. Abelló Ugalde and Domenico Lahaye

Mathematics and Computers in Simulation (MATCOM), 2024, vol. 225, issue C, 914-938

Abstract: In this paper we discuss the numerical solution of the Helmholtz equation with mixed boundary conditions on a 2D physical domain Ω. The so called radiation problem depends on the constant wavenumber k, that in some medical applications can be of order of thousands. For these values of k the classical Finite Element Method (FEM) faces up several numerical difficulties. To mitigate these limitations we apply the Isogeometric Analysis (IgA) to compute the approximated solution uh. Main steps of IgA are discussed and specific proposals for their fulfillment are addressed, with focus on some aspects not covered in available publications. In particular, we introduce a low distortion quadratic NURBS parametrization of Ω that represents exactly its boundary and contributes to the accuracy of uh. Our approach is non-isoparametric since uh is a bicubic tensor product polynomial B-spline function on Ω. This allows to improve the numerical solution refining the approximation space and keeping the coarser parametrization of the domain. Moreover, we discuss the role of the number of degrees of freedom in the directions perpendicular and longitudinal to wave front and its relationship with the noise and the shift in amplitude and phase of uh. The linear system derived from IgA discretization of the radiation problem is solved using GMRES and we show through experiments that the incomplete factorization of the Complex Shifted Laplacian provides a very good preconditioner. To solve the radiation problem, we have implemented IgA approach using the open source package GeoPDEs. A comparison with FEM is included, to provide evidence that IgA approach is superior since it is able to reduce significantly the pollution error, especially for high values of k, producing additionally smoother solutions which depend on less degrees of freedom.

Keywords: Isogeometric analysis; Helmholtz equation; Radiation problem (search for similar items in EconPapers)
Date: 2024
References: View references in EconPapers View complete reference list from CitEc
Citations:

Downloads: (external link)
http://www.sciencedirect.com/science/article/pii/S0378475423004184
Full text for ScienceDirect subscribers only

Related works:
This item may be available elsewhere in EconPapers: Search for items with the same title.

Export reference: BibTeX RIS (EndNote, ProCite, RefMan) HTML/Text

Persistent link: https://EconPapers.repec.org/RePEc:eee:matcom:v:225:y:2024:i:c:p:914-938

DOI: 10.1016/j.matcom.2023.09.018

Access Statistics for this article

Mathematics and Computers in Simulation (MATCOM) is currently edited by Robert Beauwens

More articles in Mathematics and Computers in Simulation (MATCOM) from Elsevier
Bibliographic data for series maintained by Catherine Liu ().