[论文解读] Nonlinear Geometric Optics method based multi-scale numerical schemes for highly-oscillatory transport equations

Nicolas Crouseilles, Shi Jin+1 more|arXiv (Cornell University)|May 31, 2016

Advanced Numerical Methods in Computational Mathematics参考文献 23被引用 3

一句话总结

本论文通过采用非线性几何光学假设，将振荡相位嵌入到一个独立变量 τ 中，提出了一种新颖的渐近保持（asymptotic-preserving）数值格式，用于求解高度振荡的输运方程。该方法在波周期内实现了与小参数 ε 无关的均匀一阶精度，无需解析高频振荡，从而实现了网格尺寸和时间步长与 ε 无关的稳定且精确的模拟。数值实验在多个尺度下验证了该方法的有效性。

ABSTRACT

We introduce a new numerical strategy to solve a class of oscillatory transport PDE models which is able to captureaccurately the solutions without numerically resolving the high frequency oscillations {\em in both space and time}.Such PDE models arise in semiclassical modeling of quantum dynamics with band-crossings, and otherhighly oscillatory waves. Our first main idea is to use the nonlinear geometric optics ansatz, which builds theoscillatory phase into an independent variable. We then choose suitable initial data, based on the Chapman-Enskog expansion, for the new model. For a scalar model, we prove that so constructed model will have certain smoothness, and consequently, for a first order approximation scheme we prove uniform error estimates independent of the (possibly small) wave length. The method is extended to systems arising from a semiclassical model for surface hopping, a non-adiabatic quantum dynamic phenomenon. Numerous numerical examples demonstrate that the method has the desired properties.

研究动机与目标

• 开发一种数值方法，能够在不解析小波长 ε 的情况下准确捕捉高度振荡输运方程的解。
• 克服标准格式在高频振荡下因要求 Δx, Δt = O(ε) 而导致的计算不可行性。
• 通过引入几何光学原理，将渐近保持（AP）格式扩展至空间与时间均振荡的问题。
• 通过基于 Chapman-Enskog 展开精心构造的初始条件，确保与 ε 无关的均匀误差估计。
• 将该方法应用于涉及非绝热量子跃迁的半经典表面跳跃模型中的系统。

提出的方法

• 该方法采用非线性几何光学（NGO）假设，将振荡相位转化为独立变量 τ，从而将高频振荡从物理空间和时间中解耦。
• 推导出振幅函数 V 的新方程组，其演化遵循带有相位相关系数的修正输运方程。
• 利用 Chapman-Enskog 类型展开构造初始数据，以确保与原始问题振荡结构一致的光滑性与一致性。
• 采用时间分裂谱方法进行空间离散化，并对输运项与松弛项实现精确时间积分，而 τ 中的快速振荡则通过伪谱方法处理。
• 对于系统，该方法扩展至 2×2 系统以模拟表面跳跃，通过将系统变换为 (x, p, τ) 变量下的形式，并在 τ 上实现精确积分。
• 该格式通过相位变量变换与恰当的初始数据，避免了对 ε 的解析，同时在粗网格（Δx ≫ ε）与固定时间步长下保持精度。

实验结果

研究问题

• RQ1 能否设计一种数值格式，在不解析小波长 ε 的前提下，准确捕捉输运方程的高度振荡解？
• RQ2 非线性几何光学假设是否能实现数值解在 ε → 0 极限下的均匀收敛？
• RQ3 Chapman-Enskog 展开能否用于构造在变换系统中保持光滑性与精度的初始数据？
• RQ4 该方法在涉及非绝热跃迁的系统（如表面跳跃模型）中表现如何？
• RQ5 当 ε 减小至 1/256 时，该方法能否在标准格式计算上变得不可行的情况下，仍保持精度与效率？

主要发现

• 该方法在标量模型中实现了与 ε 无关的均匀一阶误差估计，且误差不随波周期变化，证明了无需 Δx 或 Δt 随 ε 缩放即可实现收敛。
• 当 ε = 1 时，新方法使用 N_x = 32, N_p = 64, N_τ = 8 的参数，其精度与参考解（Δt = 0.05, N_x = 256）相当，而 CPU 时间从 15 秒降低至 1 分钟。
• 当 ε = 1/32 时，新方法在相同粗网格（N_x = 32）下仍保持精度，而参考方法需 Δt = 0.02 与 N_x = 512，CPU 时间增至 75 秒。
• 当 ε = 1/256 时，新方法使用固定参数（Δt = 0.05, N_x = 32）仍能产生精确结果，而参考方法需 Δt = 0.0005 与 N_x = 4096，CPU 时间达 1420 秒。
• 数值结果表明，即使网格粗于空间振荡波长，该方法仍能准确捕捉点态解与密度分布。
• 该方法成功处理了复杂的 2×2 表面跳跃模型，展现出在 ε 变化多个数量级下的鲁棒性与可扩展性。