[论文解读] Frozen Gaussian approximation for high frequency wave propagation in periodic media

Ricardo Delgadillo, Jianfeng Lu+1 more|arXiv (Cornell University)|Apr 30, 2015

Advanced Mathematical Modeling in Engineering参考文献 23被引用 3

一句话总结

该论文通过结合布洛赫分解与相空间渐近分析，提出了一种用于周期介质中高频波传播的冻结高斯近似（FGA）方法。该方法在计算成本与小参数ε无关的情况下实现了对真实解的收敛，克服了传统方法的网格分辨率限制，同时避免了高斯光束方法固有的光束展宽问题。

ABSTRACT

Propagation of high-frequency wave in periodic media is a challenging problem due to the existence of multiscale characterized by short wavelength, small lattice constant and large physical domain size. Conventional computational methods lead to extremely expensive costs, especially in high dimensions. In this paper, based on Bloch decomposition and asymptotic analysis in the phase space, we derive the frozen Gaussian approximation for high-frequency wave propagation in periodic media and establish its converge to the true solution. The formulation leads to efficient numerical algorithms, which are presented in a companion paper [Delgadillo, Lu and Yang, arXiv:1509.05552].

研究动机与目标

• 解决由于短波长和小晶格常数导致的周期介质中高频波传播模拟的计算不可行性。
• 克服传统数值方法在高维情况下的网格分辨率限制。
• 将冻结高斯近似（FGA）扩展至周期介质，同时在焦散区域保持效率与精度。
• 严格建立FGA方法对具有周期性和光滑势的薛定谔方程真实解的收敛性。

提出的方法

• 利用布洛赫分解与渐近分析，为具有周期势的薛定谔方程推导出相空间积分表示。
• 在相空间中将FGA解构造为一系列冻结高斯函数的叠加，每个函数具有固定宽度，并通过哈密顿流演化。
• 通过ε的幂级数展开推导振幅和相位函数的演化方程。
• 利用布洛赫算子的谱控制解分量的正则性与有界性。
• 应用傅里叶积分算子理论与能量估计，证明FGA近似的收敛性。
• 利用贝塞尔不等式与布洛赫能带中的谱隙，建立解分量的统一有界性。

实验结果

研究问题

• RQ1 能否将冻结高斯近似扩展至具有周期势的周期介质中的高频波传播？
• RQ2 在周期系统中存在焦散时，FGA方法如何保持精度并避免光束展宽？
• RQ3 在半经典极限ε → 0下，FGA近似对真实解的收敛速率如何？
• RQ4 能否使计算成本与ε无关，同时在周期介质中解析短波长？

主要发现

• FGA方法在半经典区域实现了对薛定谔方程真实解的收敛，误差由近似解残差乘以ε⁻¹控制。
• 该方法在焦散处仍保持稳定与精确，避免了几何光学与eikonal方程方法的失效。
• 由于布洛赫能带中的谱隙，解分量（包括振幅与相位）在时间和空间上保持统一有界。
• FGA的公式化导出了计算成本与ε无关的数值算法，支持高维中的高效模拟。
• 收敛性通过能量估计与布洛赫算子逆的有界性建立，依赖于谱隙条件。