Boundary integral equation based solvers for elastic transmission problems in complex media
Abstract
We present a boundary integral equation (BIE) framework for two-dimensional time-harmonic elastic transmission problems. Using the Helmholtz decomposition, the elastodynamic displacement fields are expressed in terms of scalar pressure and shear potentials, so that only the four classical Helmholtz boundary integral operators are needed. The key step is a rewriting of the Navier traction operator in the local normal--tangential frame of the interface, in which second-order normal derivatives are eliminated in favour of tangential derivatives and curvature terms. This leads to a one-parameter family of indirect combined field integral equations, depending on a coupling constant α≠ 0, which we prove to be uniquely solvable. The principal part of the resulting operator is, however, defective, and the discretized systems are severely ill-conditioned. Using the pseudodifferential calculus of the Helmholtz BIOs, we construct explicit left and right analytical preconditioners and prove that the preconditioned operator is a compact perturbation of the identity. Numerical experiments based on high-order Nyström discretizations confirm that the preconditioner clusters the spectrum: the number of GMRES iterations becomes independent of the discretization and is reduced by one order of magnitude or more, also in the high-frequency regime and for non-convex scatterers.
Related
- Hierarchical low-rank compression exploiting symmetric boundary integral equations in transmission scattering
- Volterra Integral Reduction for Boundary Diffusion Problems
- From multi-layered problems to multiple two-layered problems: a novel frequency-time hybrid multiple-scattering integral equation solver
- A 2+1-dimensional extended Dym equation: moving boundary problems solvable via Painlevé II symmetry reduction
- An Extended 2+1-Dimensional Gardner Equation: On Moving Boundary Problems Solvable via Ermakov-Painlevé II Symmetry Reduction
BibTeX
@misc{domınguez2026boundary,
title = {{Boundary integral equation based solvers for elastic transmission problems in complex media}},
author = {V. Dom\'ınguez and C. Turc},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.11836},
eprint = {2610.11836},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.11836},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Boundary integral equation based solvers for elastic transmission problems in complex media
AU - V. Dom\'ınguez
AU - C. Turc
PY - 2026
DA - 2026/10/08
PB - arXiv
DO - 10.48550/arxiv.2610.11836
UR - https://arxiv.org/abs/2610.11836
AB - We present a boundary integral equation (BIE) framework for two-dimensional time-harmonic elastic transmission problems. Using the Helmholtz decomposition, the elastodynamic displacement fields are expressed in terms of scalar pressure and shear potentials, so that only the four classical Helmholtz boundary integral operators are needed. The key step is a rewriting of the Navier traction operator in the local normal--tangential frame of the interface, in which second-order normal derivatives are eliminated in favour of tangential derivatives and curvature terms. This leads to a one-parameter family of indirect combined field integral equations, depending on a coupling constant α≠ 0, which we prove to be uniquely solvable. The principal part of the resulting operator is, however, defective, and the discretized systems are severely ill-conditioned. Using the pseudodifferential calculus of the Helmholtz BIOs, we construct explicit left and right analytical preconditioners and prove that the preconditioned operator is a compact perturbation of the identity. Numerical experiments based on high-order Nyström discretizations confirm that the preconditioner clusters the spectrum: the number of GMRES iterations becomes independent of the discretization and is reduced by one order of magnitude or more, also in the high-frequency regime and for non-convex scatterers.
ER -