Convergence of frozen Gaussian approximation for high frequency wave propagation

Jianfeng Lu,Xu Yang

DOI: https://doi.org/10.48550/arXiv.1012.5055

2010-12-23

Abstract:The frozen Gaussian approximation provides a highly efficient computational method for high frequency wave propagation. The derivation of the method is based on asymptotic analysis. In this paper, for general linear strictly hyperbolic system, we establish the rigorous convergence result for frozen Gaussian approximation. As a byproduct, higher order frozen Gaussian approximation is developed.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the numerical calculation difficulties in high - frequency wave propagation. Specifically, the direct numerical discretization method needs to select a grid size comparable to or smaller than the wavelength to obtain accurate solutions when dealing with high - frequency wave propagation. However, since the size of the domain is usually large, this leads to very high computational costs. To solve this problem, the author introduced an efficient computational method - the Frozen Gaussian Approximation (FGA).

### Main contributions of the Frozen Gaussian Approximation

1. **Efficient computation**: The Frozen Gaussian Approximation provides an efficient computational tool for solving high - frequency wave equations.
2. **Convergence proof**: This paper established strict convergence results for the Frozen Gaussian Approximation and proved its effectiveness for general linear strictly hyperbolic systems.
3. **Higher - order approximation**: As a by - product, this paper also developed a higher - order Frozen Gaussian Approximation method.

### Summary of mathematical expressions

- **Initial condition**: Let $\{u^\varepsilon_0\}_{\varepsilon > 0}$ be a family of asymptotically high - frequency initial conditions, satisfying $\|u^\varepsilon_0\|_{L^2(\mathbb{R}^3)} \leq M$.
- **Propagator**: Denote $P_t$ as the propagator of the first - order strictly hyperbolic system, and $P^\varepsilon_{t,K}$ as the propagator of the $K$-th order Frozen Gaussian Approximation.
- **Main theorem**: For the above - mentioned initial conditions, the following convergence result holds:
\[
\|P_t u^\varepsilon_0 - P^\varepsilon_{t,K} u^\varepsilon_0\|_{L^2(\mathbb{R}^3)} \lesssim \varepsilon^K M
\]

### Advantages of the method

- **Overcoming the limitations of geometric optics and Gaussian beam methods**:
- The geometric optics method fails near caustic points.
- The Gaussian beam method has a decreasing accuracy when the beam spreads.

The Frozen Gaussian Approximation avoids these problems through phase - plane asymptotic analysis and provides broader applicability and higher computational efficiency.

### Conclusion

This paper, through strict mathematical derivations, proved the effectiveness and convergence of the Frozen Gaussian Approximation in high - frequency wave propagation problems, providing a solid theoretical basis for practical applications.