A hybrid approach to solve the high-frequency Helmholtz equation with source singularity in smooth heterogeneous media

Fang, Jun; Qian, Jianliang; Zepeda-Núñez, Leonardo; Zhao, Hongkai

doi:10.1016/j.jcp.2018.03.011

Cited by 13 publications

(11 citation statements)

References 78 publications

(166 reference statements)

Supporting

Mentioning

Contrasting

Order By: Relevance

“…Regarding (ii): the first such bounds were proved for the exterior Helmholtz equation in [8] (for A variable and n = 1) and [9], (for n variable and A = I). Recent such bounds were proved in [59,38,32,79] (and for problems posed on bounded domains with impedance boundary conditions [23,Chapter 2], [6,11,67,78,38,39]), the renewed interest due to growing interest in the numerical analysis of Helmholtz equation with variable coefficients [24,6,11,67,30,33,25,26,70,39,32,71,50,36,51].…”

Section: Discussion Of the Novelty Of The Bound In Theorem 16mentioning

confidence: 99%

Scattering by finely-layered obstacles: frequency-explicit bounds and homogenization

Chaumont-Frelet¹,

Spence²

2021

Preprint

View full text Add to dashboard Cite

We consider the scalar Helmholtz equation with variable, discontinuous coefficients, modelling transmission of acoustic waves through an anisotropic penetrable obstacle. We first prove a well-posedness result and a frequency-explicit bound on the solution operator, with both valid for sufficiently-large frequency and for a class of coefficients that satisfy certain monotonicity conditions in one spatial direction, and are only assumed to be bounded (i.e., L ∞ ) in the other spatial directions. This class of coefficients therefore includes coefficients modelling transmission by penetrable obstacles with a (potentially large) number of layers (in 2-d) or fibres (in 3-d). Importantly, the frequency-explicit bound holds uniformly for all coefficients in this class; this uniformity allows us to consider highly-oscillatory coefficients and study the limiting behaviour when the period of oscillations goes to zero. In particular, we bound the H 1 error committed by the first-order bulk correction to the homogenized transmission problem, with this bound explicit in both the period of oscillations of the coefficients and the frequency of the Helmholtz equation; to our knowledge, this is the first homogenization result for the Helmholtz equation that is explicit in these two quantities and valid without the assumption that the frequency is small.

show abstract

Section: Discussion Of the Novelty Of The Bound In Theorem 16mentioning

confidence: 99%

Scattering by finely-layered obstacles: frequency-explicit bounds and homogenization

Chaumont-Frelet¹,

Spence²

2021

Preprint

View full text Add to dashboard Cite

show abstract

“…Can other techniques be combined with transmission conditions as a remedy? We emphasize the recent progress of some techniques for solving wave problems: for example, coarse spaces (Conen, Dolean, Krause and Nataf 2014, Bonazzoli, Dolean, Graham, Spence and Tournier 2018, Bootland, Dolean, Jolivet and Tournier 2021a or deflation (Dwarka and Vuik 2020, Dwarka, Tielen, Möller and Vuik 2021, Bootland, Dwarka, Jolivet, Dolean and Vuik 2021b, time domain solvers (Grote and Tang 2019, Grote, Nataf, Tang and Tournier 2020, Appelo, Garcia and Runborg 2020, Stolk 2021, absorption or shifted-Laplace (Gander, Graham and Spence 2015, Graham, Spence and Vainikko 2017, Graham et al 2020, Hocking and Greif 2021, H-matrix (Beams, Gillman and Hewett 2020, Lorca, Beams, Beecroft and Gillman 2021, Liu, Ghysels, Claus and Li 2021, Bonev and Hesthaven 2022, DPG (Petrides and Demkowicz 2021) and high-frequency asymptotic (Lu, Qian and Burridge 2016, Fang, Qian, Zepeda-Núñez and Zhao 2018, Jacobs and Luo 2021. There are plenty of questions still to be answered in the future.…”

Section: Double Sweep Schwarz Methods For the Layered Medium Wave Pro...mentioning

confidence: 99%

Schwarz methods by domain truncation

Gander¹,

Hui²

2022

Acta Numerica

View full text Add to dashboard Cite

Schwarz methods use a decomposition of the computational domain into subdomains and need to impose boundary conditions on the subdomain boundaries. In domain truncation one restricts the unbounded domain to a bounded computational domain and must also put boundary conditions on the computational domain boundaries. In both fields there are vast bodies of literature and research is very active and ongoing. It turns out to be fruitful to think of the domain decomposition in Schwarz methods as a truncation of the domain onto subdomains. Seminal precursors of this fundamental idea are papers by Hagstrom, Tewarson and Jazcilevich (1988), Després (1990) and Lions (1990). The first truly optimal Schwarz method that converges in a finite number of steps was proposed by Nataf (1993), and used precisely transparent boundary conditions as transmission conditions between subdomains. Approximating these transparent boundary conditions for fast convergence of Schwarz methods led to the development of optimized Schwarz methods – a name that has become common for Schwarz methods based on domain truncation. Compared to classical Schwarz methods, which use simple Dirichlet transmission conditions and have been successfully used in a wide range of applications, optimized Schwarz methods are much less well understood, mainly due to their more sophisticated transmission conditions.A key application of Schwarz methods with such sophisticated transmission conditions turned out to be time-harmonic wave propagation problems, because classical Schwarz methods simply do not work in this case. The past decade has given us many new Schwarz methods based on domain truncation. One review from an algorithmic perspective (Gander and Zhang 2019) showed the equivalence of many of these new methods to optimized Schwarz methods. The analysis of optimized Schwarz methods, however, is lagging behind their algorithmic development. The general abstract Schwarz framework cannot be used for the analysis of these methods, and thus there are many open theoretical questions about their convergence. Just as for practical multigrid methods, Fourier analysis has been instrumental for understanding the convergence of optimized Schwarz methods and for tuning their transmission conditions. Similar to local Fourier mode analysis in multigrid, the unbounded two-subdomain case is used as a model for Fourier analysis of optimized Schwarz methods due to its simplicity. Many aspects of the actual situation, e.g. boundary conditions of the original problem and the number of subdomains, were thus neglected in the unbounded two-subdomain analysis. While this gave important insight, new phenomena beyond the unbounded two-subdomain models were discovered.This present situation is the motivation for our survey: to give a comprehensive review and precise exploration of convergence behaviours of optimized Schwarz methods based on Fourier analysis, taking into account the original boundary conditions, many-subdomain decompositions and layered media. We consider as our model problem the operator $-\Delta + \eta $ in the diffusive case $\eta>0$ (screened Laplace equation) or the oscillatory case $\eta <0$ (Helmholtz equation), in order to show the fundamental difference in behaviour of Schwarz solvers for these problems. The transmission conditions we study include the lowest-order absorbing conditions (Robin), and also more advanced perfectly matched layers (PMLs), both developed first for domain truncation. Our intensive work over the last two years on this review has led to several new results presented here for the first time: in the bounded two-subdomain analysis for the Helmholtz equation, we see strong influence of the original boundary conditions imposed on the global problem on the convergence factor of the Schwarz methods, and the asymptotic convergence factors with small overlap can differ from the unbounded two-subdomain analysis. In the many-subdomain analysis, we find the scaling with the number of subdomains, e.g. when the subdomain size is fixed, robust convergence of the double-sweep Schwarz method for the free-space wave problem, either with fixed overlap and zeroth-order Taylor conditions or with a logarithmically growing PML, and we find that Schwarz methods with PMLs work like smoothers that converge faster for higher Fourier frequencies; in particular, for the free-space wave problem, plane waves (in the error) passing through interfaces at a right angle converge more slowly. In addition to our main focus on analysis in Sections 2 and 3, we start in Section 1 with an expository historical introduction to Schwarz methods, and in Section 4 we give a brief interpretation of the recently proposed optimal Schwarz methods for decompositions with cross-points from the viewpoint of transmission conditions. We conclude in Section 5 with a summary of open research problems. In Appendix A we provide a Matlab program for a block LU form of an optimal Schwarz method with cross-points, and in Appendix B we give the Maple program for the two-subdomain Fourier analysis.

show abstract

“…However, in room acoustics, since the speed of sound is typically a perturbation about a constant value due to small fluctuations in temperature, our simple model here should be reasonably accurate. There are other ways of modeling the amplitude near a point source that have been developed for use with high-order sweeping methods [107,48]. In Section 1.6, we develop reflection and diffraction BCs for the amplitude.…”

Section: The Point Source Amplitudementioning

confidence: 99%

Numerical geometric acoustics

Potter

Cameron

Duraiswami

2020

The Journal of the Acoustical Society of America

View full text Add to dashboard Cite

Sound propagation in air is accurately described by a small perturbation of the ambient pressure away from a quiescent state. This is the realm of linear acoustics, where the propagation of a time-harmonic wave can be modeled using the Helmholtz equation. When the wavelength is small relative to the size of a scattering obstacle, techniques from geometric optics are applicable. Geometric methods such as raytracing are often used for computational room acoustics simulations in situations where the geometry of the built environment is sufficiently complicated. At the same time, the high-frequency approximation of the Helmholtz equation is described by two partial differential equations: the eikonal equation, whose solution gives the first arrival time of a geometric acoustics/optics wavefront as a field; and a transport equation, the solution of which describes the amplitude of that wavefield. Phenomena related to high-frequency acoustic diffraction are frequently omitted from these models because of their complexity. These phenomena can be modeled using a high-frequency diffraction theory, such as the uniform theory of diffraction. Despite their shortcomings, geometric methods for room acoustics provide a useful trade-off between realism and computational efficiency.

show abstract

A hybrid approach to solve the high-frequency Helmholtz equation with source singularity in smooth heterogeneous media

Cited by 13 publications

References 78 publications

Scattering by finely-layered obstacles: frequency-explicit bounds and homogenization

Scattering by finely-layered obstacles: frequency-explicit bounds and homogenization

Schwarz methods by domain truncation

Numerical geometric acoustics

Contact Info

Product

Resources

About