[論文レビュー] 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セルに増加させた際、誤差が最大2桁低下するなど、安定で反射のないシミュレーションが実現された。

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格子に適用し、数値的安定性と精度を向上させる。
• 物理的および数値的実験を通じた手法の検証により、さまざまな流れ速度下で反射が低減され、解が有界であることを示す。

提案手法

• 空間時間座標および未知数（密度、速度、圧力）にローレンツ変換を適用し、移流的音響系を非移流的形に変換する。
• 渇きのない流れと断熱的流れの仮定を用いて、線形化されたエーラー方程式をマクスウェル方程式に類似したスカラー型系に簡略化する。
• ベレンジャーおよびコリーノの手法を拡張し、ゼロ次微分項を組み込むことで、新たなPML定式化を導出。これにより、適切な定式化と波の吸収が保証される。
• ステアガード有限差分法を用いて、圧力と速度成分がそれぞれセル中心およびエッジに配置された、ステアガード型の直交格子上で数値スキームを実装する。
• 計算領域の端縁にPML層を適用し、変動する厚さ（4～20セル）を設定することで、不自然な反射を防ぐ境界条件を実装する。
• ガウス包絡と正弦波時間変調を用いた点源励起を用い、波の伝播および吸収効率を評価する。

実験結果

リサーチクエスチョン

• RQ1 ローレンツ変換を用いて、音響方程式から外部流れを完全に分離し、移流的音響問題を非移流的問題に変換できるか？
• RQ2 変換されたローレンツ座標系において、完全吸収層（PML）をどのように導出し、実装するか。これにより、出射音波の効果的吸収が保証されるか？
• RQ3 PMLの厚さが、さまざまな流れ速度下における数値的反射および圧力場の誤差に与える影響は何か？
• RQ4 HaWAYステアガード有限差分スキームは、移流的音響シミュレーションにおいて、安定性と精度の面でどの程度の性能を示すか？
• RQ5 異なる流れ方向およびマッハ数を想定した長時間シミュレーションにおいて、解が有界であり、誤差が低く保たれるか？

主な発見

• ローレンツ変換により、外部流れが音響方程式から明確に分離され、変換された座標系では標準的な非移流的ソルバーが利用可能になる。
• ゼロ次微分項を含むPML方程式は、適切な定式化が保証され、数値的安定性と波の効果的吸収を実現する。
• PMLの厚さを4セルから20セルに増加させた場合、圧力場におけるL2誤差は約2桁低下し、吸収精度の著しい向上が確認された。
• 数値実験では、マッハ数0.5までに達する場合でさえ、長時間にわたり数値的爆発の兆しなしに安定したシミュレーションが維持された。
• さまざまな流れ方向、特に( u/c₀, v/c₀ ) = (1/(2√2), 1/(2√2)) および (2/√17, 1/(2√17)) の斜め流れに対しても、反射が低く抑えられ、高いロバスト性が確認された。
• 大域的な領域（x_max = 150）で計算された基準解は、t < 300Δt の時間窓内において、PMLを用いた小域的領域解が高精度であることを裏付け、PMLの性能を検証した。