Scattered wavefield in the stochastic homogenization regime

Josselin Garnier,Laure Giovangigli,Quentin Goepfert,Pierre Millien

2023-09-14

Abstract:In the context of providing a mathematical framework for the propagation of ultrasound waves in a random multiscale medium, we consider the scattering of classical waves (modeled by a divergence form scalar Helmholtz equation) by a bounded object with a random composite micro-structure embedded in an unbounded homogeneous background medium. Using quantitative stochastic homogenization techniques, we provide asymptotic expansions of the scattered field in the background medium with respect to a scaling parameter describing the spatial random oscillations of the micro-structure. Introducing a boundary layer corrector to compensate the breakdown of stationarity assumptions at the boundary of the scattering medium, we prove quantitative $L^2$- and $H^1$- error estimates for the asymptotic first-order expansion. The theoretical results are supported by numerical experiments.

Analysis of PDEs

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is to provide a mathematical framework for describing the propagation of ultrasonic waves in random multi - scale media. Specifically, the author focuses on the problem of scattering of classical waves (modeled by the scalar Helmholtz equation in divergence form) by a finite object with a random composite microstructure in an unbounded homogeneous background medium. By using the quantitative stochastic assimilation technique, the author provides an asymptotic expansion of the scattering field with respect to the scaling parameter that describes the spatial random oscillation of the microstructure. To compensate for the failure of the stationary assumption at the boundary of the scattering medium, a boundary - layer corrector is introduced, and the L2 - and H1 - error estimates of the first - order asymptotic expansion are proved. The theoretical results are supported by numerical experiments.

### Main contributions
1. **Mathematical model**: A mathematical model is proposed to describe the acoustic properties of soft tissues and the quantitative asymptotic expansion of the scattering field with respect to the scale of acoustic heterogeneity in the medium.
2. **Boundary - layer corrector**: A boundary - layer corrector is introduced to enforce the transmission condition on the object boundary, and the L2 - and H1 - norm convergence rates are proved.
3. **Asymptotic expansion of the scattering field**: Using the Lippman - Schwinger equation and the results on commutator fluctuations, a quantitative first - order expansion of the scattering wave outside the object is derived.
4. **Numerical verification**: Numerical illustrations of the solution of the multi - scale problem, its corrector, and the first - order expansion are presented, and the numerical convergence rates are calculated to support the theoretical results.

### Key techniques and methods
- **Stochastic assimilation technique**: Used to deal with partial differential equations with rapidly oscillating coefficients, especially in the stochastic case.
- **Boundary - layer corrector**: Used to compensate for the failure of the stationary assumption at the boundary.
- **Lippman - Schwinger equation**: Used to describe the integral representation of the scattering field.
- **Numerical simulation**: To verify the validity of the theoretical results.

### Mathematical formulas
- **Helmholtz equation**:
\[
-\nabla\cdot(a^{\varepsilon}\nabla u^{\varepsilon})-k^{2}n^{\varepsilon}u^{\varepsilon} = 0\quad\text{in }\mathbb{R}^{d}
\]
- **Asymptotic expansion**:
\[
u^{\varepsilon}(x)\approx u_{0}(x)+\varepsilon u_{1,\varepsilon}(x)
\]
- **Boundary - layer corrector**:
\[
v^{\varepsilon}\text{ satisfies }\begin{cases}
-\Delta v^{\varepsilon}-k^{2}n_{0}v^{\varepsilon}=0&\text{in }B_{R}\setminus D,\\
-\nabla\cdot(a^{\varepsilon}\nabla v^{\varepsilon})-k^{2}n^{\varepsilon}v^{\varepsilon}=0&\text{in }D,\\
v^{-}_{\varepsilon}-v^{+}_{\varepsilon}=\varepsilon u_{1,\varepsilon}&\text{on }\partial D,\\
\nabla v^{-}_{\varepsilon}\cdot\nu - a^{\varepsilon}\nabla v^{+}_{\varepsilon}\cdot\nu=\varepsilon\sum_{i = 1}^{d}(\nabla\cdot(\sigma^{\varepsilon}_{i}\partial_{i}u_{0})+\cdot\nu - k^{2}\varepsilon(\beta^{\varepsilon}u_{0})+\cdot\nu)&\text{on }\partial D,\\
\nu
\end{cases}
\]