[論文レビュー] On the existence of self-similar converging shocks for arbitrary equation of state

Zachary M. Boyd, Scott D. Ramsey+1 more|arXiv (Cornell University)|Jul 12, 2017

Computational Fluid Dynamics and Aerodynamics参考文献 29被引用数 7

ひとこと要約

本稿は、任意の状態方程式における自己相似収束衝撃波解の存在に必要な熱力学的条件——具体的には断熱体積弾性率 $ K_S $ が $ K_S = p f(\rho) $ の形をとること——を確立する。非線形固有値問題への還元が示され、定常衝撃速度（$ \lambda = 0 $）は解が非有界となるため、非理想物質に対する標準的な Guderley フレームワークを無効化する。

ABSTRACT

We extend Guderley's problem of finding a self-similar scaling solution for a converging cylindrical or spherical shock wave from the ideal gas case to the case of flows with an arbitrary equation of state closure model, giving necessary conditions for the existence of a solution. The necessary condition is a thermodynamic one, namely that the adiabatic bulk modulus, $K_S$, of the fluid be of the form $pf(ρ)$ where $p$ is pressure, $ρ$ is mass density, and $f$ is an arbitrary function. Although this condition has appeared in the literature before, here we give a more rigorous and extensive treatment. Of particular interest is our novel analysis of the governing ordinary differential equations (ODEs), which shows that, in general, the Guderley problem is always an eigenvalue problem. The need for an eigenvalue arises from basic shock stability principles -- an interesting connection to the existing literature on the relationship between self-similarity of the second kind and stability. We also investigate a special case, usually neglected by previous authors, where assuming constant shock velocity yields a reduction to ODEs for every material, but those ODEs never have a bounded, differentiable solution. This theoretical work is motivated by the need for more realistic test problems in the verification of inviscid compressible flow codes that simulate flows in a variety of non-ideal gas materials.

研究の動機と目的

• 理想気体から任意の状態方程式（EOS）への Guderley の自己相似収束衝撃波問題の拡張を通じて、非理想物質のより現実的なコード検証を可能にする。
• 圧縮流における有界で微分可能な自己相似解の存在に必要な熱力学的条件——具体的には $ K_S = p f(\rho) $ ——を特定する。
• 得られた常微分方程式（ODE）を非線形固有値問題として厳密に解析し、自己相似解における安定性とスケーリングの役割を明確にする。
• 定常衝撃速度（$ \lambda = 0 $）の特殊ケースを検討し、ODE が単純化されるが、有界な解が得られないことを示す。
• 非理想・圧縮流における流体動力学コード検証のための正確な解析的テスト問題を構築する基盤を提供する。

提案手法

• 任意の状態方程式における圧縮流のオイラー方程式の自己相似還元を、空間的および時間的スケーリングに基づく相似変数を用いて導出する。
• 自己相似解の存在に不可欠な熱力学的量として断熱体積弾性率 $ K_S $ を特定し、$ K_S = p f(\rho) $ が解の存在に必要であることを示す。
• 支配方程式である偏微分方程式（PDE）を相似変数 $ \xi = r / |t|^\lambda $ における常微分方程式（ODE）系に還元し、$ \lambda $ がスケーリング特性を決定する。
• 有界で微分可能な解の存在を解析し、衝撃波の安定性制約により問題が本質的に非線形固有値問題であることが判明する。
• 衝撃波先端における ODE の特異性を解析し、解が有界であるためには固有値問題として解かなければならないことを示す。
• 定常衝撃速度（$ \lambda = 0 $）のケースを検討し、ODE が単純化されるが、分母が消える一方分子が非ゼロのため、有界で微分可能な解が存在しないことを示す。

実験結果

リサーチクエスチョン

• RQ1 任意の流体における自己相似収束衝撃波解の存在に必要な状態方程式の熱力学的条件は何か？
• RQ2 一般状態方程式下での圧縮流のオイラー方程式の自己相似還元はどのように振る舞い、衝撃波安定性から生じる制約は何か？
• RQ3 非理想物質における Guderley 問題がなぜ非線形固有値問題に還元されるのか、その物理的根拠は何か？
• RQ4 $ \lambda = 0 $（定常衝撃速度）のケースは、Guderley フレームワーク内で有界で微分可能な解を生じるか？もし不能ならば、その理由は何か？
• RQ5 定常衝撃速度のケースが失敗することは、非理想物質における流体動力学コード検証のための解析的テスト問題の構築にどのような影響を及えるか？

主な発見

• 任意の状態方程式における自己相似収束衝撃波の存在に必要な・十分な熱力学的条件は、$ K_S = p f(\rho) $ であり、ここで $ f $ は密度の任意関数である。
• 任意の状態方程式に対する Guderley 問題は、常に衝撃波安定性と解の有界性の相互作用により、非線形固有値問題に還元される。
• 標準的な自己相似スケーリング（$ \lambda \neq 0 $）では、固有値条件が満たされない限り ODE は適切に定式化されず、有界で微分可能な解が保証されない。
• $ \lambda = 0 $ のケース（定常衝撃速度）は ODE 系が単純化されるが、分母が消える一方で分子が非ゼロのため、有界で微分可能な解は存在しない。
• $ \lambda = 0 $ のケースは、衝撃速度が定常でなければならないという条件を満たすが、有界性条件に反するため、有効な Guderley 解を生成できない。
• 本解析は、特に複雑または非理想状態方程式を有する非理想物質における流体動力学コード検証のための正確な解析的テスト問題を構築するための厳密な基盤を提供する。