[論文レビュー] Infinite-thin shock layer solutions for stationary compressible conical flows and numerical results via Fourier spectral method

Aifang Qu, Xueying Su+1 more|arXiv (Cornell University)|Jul 4, 2023

Computational Fluid Dynamics and AerodynamicsEngineering被引用数 3

ひとこと要約

本稿は、静的圧縮性円錐流れにおける無限に薄い衝撃層解を記述する高特異的非線形常微分方程式（ODE）を解くための数値フレームワークを提案する。特に、チャプリンガスおよび迎角を有する超音速流れを対象としている。主な貢献は、円錐表面に重み付きデルタ関数を用いたRandon測度解を構築し、一般化されたニュートン＝ブースマン圧力則を導出し、圧力分布および粒子軌道の物理的期待と整合する高精度な数値解を得ることにある。

ABSTRACT

We consider the problem of uniform steady supersonic Euler flows passing a straight conical body with attack angles, and study Radon measure solutions describing the infinite-thin shock layers, particularly for the Chaplygin gas and limiting hypersonic flows. As a byproduct, we obtain the generalized Newton-Busemann pressure laws. To construct the Radon measure solutions containing weighted Dirac measures supported on the edge of the cone on the 2-sphere, we derive some highly singular and non-linear ordinary differential equations (ODE). A numerical algorithm based on the combination of Fourier spectral method and Newton's method is developed to solve the physically desired nonnegative and periodic solutions of the ODE. The numerical simulations for different attack angles exhibit proper theoretical properties and excellent accuracy, thus would be useful for engineering of hypersonic aerodynamics.

研究の動機と目的

• 迎角を有する定常超音速円錐流れにおける無限に薄い衝撃層をRandon測度解を用いてモデル化すること。
• 円錐表面におけるデルタ関数の重みを記述する特異的非線形ODE系を導出・解釈すること。
• 特異的衝撃層モデルから生じる非負かつ周期的なODEの解を高精度に計算するための堅牢な数値アルゴリズムを開発すること。
• 圧力分布および粒子軌道の数値シミュレーションを通じて、解の物理的整合性を検証すること。
• 迎角がゼロでない場合のニュートン＝ブースマン圧力則を一般化すること。

提案手法

• 円錐体と迎角を有する一様な超音速流れを想定し、球座標系で問題を定式化すること。
• 無限に薄い衝撃層を表すデルタ関数の重みを記述する特異的非線形ODE系を導出すること。
• 周期的かつ非負のODEを効率的に離散化・解釈するため、Fourierスペクトル法を適用すること。
• Fourierスペクトル法とニュートン法を組み合わせることで、非線形系の収束性と精度を向上させること。
• ODEの解を用いて速度成分 $ u^t, w, w_\rho $ を再構成し、円錐表面における粒子軌道を解析すること。
• 圧力および流れ変数の単調性・対称性といった物理的整合性を示す数値シミュレーションにより、結果を検証すること。

実験結果

リサーチクエスチョン

• RQ1 迎角がゼロでない圧縮性円錐流れに対して、Randon測度解を用いて無限に薄い衝撃層解をどのように構築できるか？
• RQ2 任意の迎角におけるチャプリンガスに対する一般化されたニュートン＝ブースマン圧力則は何か？
• RQ3 衝撃層の重みを記述する高特異的非線形ODEを、高い精度と非負制約を満たしてどのように数値的に解けるか？
• RQ4 数値解において、圧力分布および流れ変数が示す物理的性質（単調性・対称性）は何か？
• RQ5 円錐表面における粒子軌道はどのように変化し、迎角の変化に伴う流れの挙動を何を示唆するか？

主な発見

• 数値法は特異的ODEの非負かつ周期的な解を高精度に計算でき、物理的に整合する無限に薄い衝撃層解の存在を裏付けた。
• 迎角 $ \alpha_0 $ が増加するにつれ、最大圧力 $ \max W_C $ は増加し、最小圧力 $ \min W_C $ は減少する。これは風上側の圧力が高くなり、背風側が低くなるという物理的直観と整合する。
• 圧力関数 $ W_C $ の最大値および最小値が、理論的曲線 $ \sin^2(\theta_0 + \alpha_0) $ および $ \sin^2(\theta_0 - \alpha_0) $ と正確に一致し、導出された一般化されたニュートン＝ブースマン則の妥当性が検証された。
• 速度成分 $ u^t, w, w_\rho $ は $ \phi = 0, \pm\pi $ で期待される特異性を示すが、これは解析的期待である $ C^1 $ 正則性とは対照的であり、数値的アーチファクトとされる。
• 円錐表面における粒子軌道は $ \phi = 0 $ 側に収束し、背風側への流れの蓄積を示しており、$ x^1 $-軸に関して対称的である。
• 数値結果は単調な傾向を示しており、$ f(\phi) $ は $ \alpha_0 $ が増加するにつれて増加し、$ w_\rho $ は $ \phi = \pm\pi $ 近辺で増加し、$ \phi = 0 $ 近辺で減少する。これは粒子衝突密度の変化を反映している。