A high-order ghost-point immersed boundary method coupled with auxiliary differential equations for time-domain impedance walls
Abstract
A ghost-point immersed boundary method is coupled with auxiliary differential equations to impose a locally-reacting impedance condition on a wall that does not conform to a Cartesian grid, in the framework of the linearized Euler equations. The ghost-point reconstruction follows the normal-stencil approach: an elliptical least squares cloud of fluid nodes is used to build a high-order polynomial along the local wall normal, which is then extrapolated onto the ghost points by one-dimensional La grange interpolation. The impedance, modeled as a massspringdamper oscillator, is advanced in time through two auxiliary state variables per ghost point, coupled to the reconstruction at every RungeKutta stage. The method is validated against exact and semi-analytical references on a flat immersed wall, a cylindrical wall reflecting an acoustic pulse, and a cylinder scattering a harmonic monopole. Convergence orders ranging from about three on the curved wall to about four and a half on the flat wall are measured, and a cross-check against an exact harmonic reference shows that a persistent offset observed against a semi-analytical reference is a property of that ref erence rather than of the coupled solver. A spectral stability analysis of the linearized semi-discrete operator further shows that the normal-stencil reconstruction reduces, without eliminating, a known instability that occurs when the immersed wall becomes tangent to the grid.
Related
- A ghost-cell immersed boundary method for the shallow water equations on Cartesian grids: mass conservation and wetting and drying over variable topography
- Simultaneous energy and entropy stability of volume-penalisation immersed-boundary methods for the compressible Navier--Stokes equations
- Existence and multiplicity results for third-order ordinary differential equations with derivative dependence subject to the Riemann–Stieltjes integral boundary conditions
- Kolmogorov-Arnold Networks for Free-Boundary Partial Differential Equations
- Boundary value problems for linear differential-algebraic equations: solvability via the Kronecker canonical form
BibTeX
@misc{ouattara2026highorder,
title = {{A high-order ghost-point immersed boundary method coupled with auxiliary differential equations for time-domain impedance walls}},
author = {Abdoulaye Ouattara},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.24785},
eprint = {2609.24785},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.24785},
note = {CC BY-NC-ND 4.0},
}
RIS
TY - UNPB
TI - A high-order ghost-point immersed boundary method coupled with auxiliary differential equations for time-domain impedance walls
AU - Abdoulaye Ouattara
PY - 2026
DA - 2026/09/21
PB - arXiv
DO - 10.48550/arxiv.2609.24785
UR - https://arxiv.org/abs/2609.24785
AB - A ghost-point immersed boundary method is coupled with auxiliary differential equations to impose a locally-reacting impedance condition on a wall that does not conform to a Cartesian grid, in the framework of the linearized Euler equations. The ghost-point reconstruction follows the normal-stencil approach: an elliptical least squares cloud of fluid nodes is used to build a high-order polynomial along the local wall normal, which is then extrapolated onto the ghost points by one-dimensional La grange interpolation. The impedance, modeled as a massspringdamper oscillator, is advanced in time through two auxiliary state variables per ghost point, coupled to the reconstruction at every RungeKutta stage. The method is validated against exact and semi-analytical references on a flat immersed wall, a cylindrical wall reflecting an acoustic pulse, and a cylinder scattering a harmonic monopole. Convergence orders ranging from about three on the curved wall to about four and a half on the flat wall are measured, and a cross-check against an exact harmonic reference shows that a persistent offset observed against a semi-analytical reference is a property of that ref erence rather than of the coupled solver. A spectral stability analysis of the linearized semi-discrete operator further shows that the normal-stencil reconstruction reduces, without eliminating, a known instability that occurs when the immersed wall becomes tangent to the grid.
ER -