PPubMed14 Jan 2015Cited 39×

Boundary regularized integral equation formulation of the Helmholtz equation in acoustics BRIEF of the Helmholtz equation

Qiang Sun, Evert Klaseboer, Boo-Cheong Khoo, Derek Y. C. Chan

Abstract

A boundary integral formulation for the solution of the Helmholtz equation is developed in which all traditional singular behaviour in the boundary integrals is removed analytically. The numerical precision of this approach is illustrated with calculation of the pressure field owing to radiating bodies in acoustic wave problems. This method facilitates the use of higher order surface elements to represent boundaries, resulting in a significant reduction in the problem size with improved precision. Problems with extreme geometric aspect ratios can also be handled without diminished precision. When combined with the CHIEF method, uniqueness of the solution of the exterior acoustic problem is assured without the need to solve hypersingular integrals.

A figure from Boundary regularized integral equation formulation of the Helmholtz equation in acoustics BRIEF of the Helmholtz equation — fig. from the paper

§ The Valyu brief

The paper, in five minutes

Read from the full text via Valyu Contents

A new singularity‑free boundary integral method (BRIEF) solves the Helmholtz equation with higher accuracy and fewer degrees of freedom.

What’s new

i.Derives the Boundary Regularized Integral Equation Formulation (BRIEF) that analytically removes all singularities from the Helmholtz boundary integral equation.
ii.Shows that BRIEF allows the use of higher‑order surface elements (linear, quadratic, spline) without computing solid angles, reducing degrees of freedom and computational time.
iii.Combines BRIEF with the CHIEF method to guarantee uniqueness without hypersingular integrals.
iv.Provides numerical evidence of error reductions of 1–3 orders of magnitude and condition‑number improvements up to a factor of 60 versus conventional BIE.
v.Demonstrates robustness for problems with extreme geometric aspect ratios and accurate field‑point evaluation near boundaries.

How they did it

The authors construct auxiliary functions g and f that satisfy the Helmholtz equation and specific boundary conditions, then form an auxiliary field ψ to subtract from the conventional boundary integral equation, eliminating singular kernels. This yields a singularity‑free integral equation (BRIEF) which is solved numerically using constant, linear and quadratic surface elements on acoustic test cases. Results are compared against conventional boundary integral methods and analytical solutions.

What they found

Maximum relative error of BRIEF Quadratic is 1–3 orders of magnitude lower than CBIM, achieving the same precision with about one‑tenth the degrees of freedom (Fig. 1b).
Condition number at kR/π = 1.0038 drops from 6000 (constant) to 81 (linear) and 100 (quadratic), a reduction of roughly 60× (Table 1).
BRIEF Quadratic maintains low error until |kR‑π| < 10⁻⁴, while constant elements misplace the resonant wavenumber (Fig. 2a).
For two closely spaced spheres (h/R = 0.001), BRIEF Linear produces smooth pressure fields and avoids the roughness seen in CBIM Linear (Figs. 3, 4).

Caveats

✕The BRIEF‑CHIEF combination may break down in certain configurations.
✕At discrete resonant wavenumbers the solution is non‑unique, a limitation shared with conventional boundary integral methods.
✕Using constant surface elements within BRIEF can locate the resonant wavenumber incorrectly.

§ The Valyu brief

Reading the full paper and taking notes. This takes a few seconds…

§ Ask this paper

Ask a question about this paper

Valyu reads the full text and answers from what the paper actually says.

§ Related work in other archives

Boundary regularized integral equation formulation of the Helmholtz equation in acoustics BRIEF of the Helmholtz equation

The paper, in five minutes

Ask a question about this paper

Similar papers in all four archives

Similar papers in all four archives

The paper, in five minutes