[论文解读] Lorentz Transform and Staggered Finite Differences for Advective Acoustics

François Dubois, Eric Duceau+2 more|arXiv (Cornell University)|May 7, 2011

Aerodynamics and Acoustics in Jet Flows参考文献 10被引用 7

一句话总结

本文提出了一种新颖的数值方法，通过应用洛伦兹变换以消除外部流场影响，并将声学方程重新表述为非对流形式，从而求解对流声学问题。采用HaWAY网格上的交错有限差分法，并结合基于洛伦兹不变变换的完美匹配层（PML），该方法实现了稳定且无反射的模拟，当PML厚度从4个单元增加到20个单元时，误差降低了两个数量级。

ABSTRACT

We study acoustic wave propagation in a uniform stationary flow. We develop a method founded on the Lorentz transform and a hypothesis of irrotationality of the acoustic perturbation. After a transformation of the space-time and of the unknown fields, we derive a system of partial differential equations that eliminates the external flow and deals with the classical case of non advective acoustics. A sequel of the analysis is a new set of perfectly matched layers equations in the spirit of the work of Berenger and Collino. The numerical implementation of the previous ideas is presented with the finite differences method HaWAY on cartesian staggered grids. Relevant numerical tests are proposed.

研究动机与目标

• 开发一种稳定且精确的数值方法，用于模拟均匀亚音速流场中的声波传播，适用于航空噪声预测。
• 通过时空坐标与场变量（密度、速度、压力）的洛伦兹变换，从控制方程中消除外部流场的影响。
• 推导一种新型完美匹配层（PML）方程，有效吸收向外传播的声波且无反射。
• 采用HaWAY方法在交错笛卡尔网格上实现该方法，以提升数值稳定性和精度。
• 通过物理与数值实验验证该方法，证明其在不同流速下具有低反射率和有界解。

提出的方法

• 对时空坐标及未知场变量（密度、速度、压力）应用洛伦兹变换，将对流声学系统转化为非对流形式。
• 基于无旋流与等熵流假设，将线性化欧拉方程简化为类似于麦克斯韦方程的标量型系统。
• 通过扩展Berenger与Collino的方法，推导新型PML公式，嵌入零阶阻尼项以确保适定性与波吸收效果。
• 在交错笛卡尔网格上使用HaWAY方法实现数值格式，其中压力与速度分量分别位于单元中心与边长中点。
• 在计算域边界处应用PML层，厚度可变（4至20个单元），确保无虚假反射。
• 采用具有高斯包络与正弦时间调制的点源激励，测试波传播与吸收效率。

实验结果

研究问题

• RQ1 洛伦兹变换能否成功消除声学方程中的外部流场影响，从而将对流声学问题转化为非对流问题？
• RQ2 在变换后的洛伦兹参考系中，如何推导并实现完美匹配层（PML），以确保向外传播声波的有效吸收？
• RQ3 PML厚度对不同流速下压力场中数值反射与误差的影响如何？
• RQ4 HaWAY交错有限差分格式在对流声学模拟中的稳定性与精度表现如何？
• RQ5 该方法在不同流速方向与马赫数下，能否在长时间模拟中保持有界解与低误差？

主要发现

• 洛伦兹变换成功地将外部流场与声学方程解耦，使得在变换后的参考系中可使用标准非对流求解器。
• 由于引入了零阶阻尼项，所推导的PML方程具有适定性，确保了数值稳定性与有效波吸收。
• 当PML厚度从4个单元增至20个单元时，压力场的L2误差降低了约两个数量级，表明吸收精度显著提升。
• 数值实验表明，该方法在长时间模拟中保持稳定，即使马赫数高达0.5，也未出现数值发散现象。
• 该方法在不同流场方向下均表现出低反射水平，包括斜向流场如(u/c₀, v/c₀) = (1/(2√2), 1/(2√2)) 与 (2/√17, 1/(2√17))，验证了其鲁棒性。
• 在大域域（x_max = 150）上计算的参考解表明，小域域结合PML的解在时间窗口t < 300Δt内具有足够精度，验证了PML性能。