[论文解读] 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 = p f(\rho) $ 的形式。本文表明，古德雷問題可歸結為非線性特徵值問題，並指出常數衝擊速度（$ \lambda = 0 $）會導致解無界，從而使標準古德雷框架在非理想材料中失效。

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），以實現對非理想材料更真實的代碼驗證。
• 識別出自相似解存在所需的熱力學條件——具體而言，$ K_S = p f(\rho) $——以確保壓縮流中解的有界性與可微性。
• 嚴謹分析所得常微分方程（ODE）系統作為非線性特徵值問題，釐清穩定性與自相似解中尺度變換的作用。
• 探討常數衝擊速度（$ \lambda = 0 $）的特殊情況，表明即使ODE簡化，仍無法產生有界解。
• 為構建非理想壓縮流中水動力代碼驗證的精確解析測試問題提供基礎。

提出的方法

• 基於徑向與時間尺度的相似變數，推導任意狀態方程下壓縮流歐拉方程的自相似簡化形式。
• 識別絕熱體積模量 $ K_S $ 為關鍵熱力學量，並證明 $ K_S = p f(\rho) $ 為自相似解存在的必要條件。
• 將控制偏微分方程（PDE）簡化為相似變數 $ \xi = r / |t|^\lambda $ 的常微分方程（ODE）系統，其中 $ \lambda $ 決定尺度行為。
• 分析ODE系統以獲得有界且可微的解，揭示由於衝擊穩定性約束，該問題本質上為非線性特徵值問題。
• 對ODE在衝擊前緣處進行奇點分析，表明除非將系統作為特徵值問題求解，否則會產生無界解。
• 考慮 $ \lambda = 0 $ 情形（常數衝擊速度），證明即使ODE簡化，仍無法獲得有界且可微的解，原因在於分母趨於零而分子非零。

实验结果

研究问题

• RQ1 在任意流體中，自相似收斂衝擊波解存在的狀態方程需滿足何種熱力學條件？
• RQ2 在一般狀態方程下，壓縮流歐拉方程的自相似簡化形式如何表現？衝擊穩定性會帶來何種約束？
• RQ3 為何非理想材料的古德雷問題會歸結為非線性特徵值問題？其物理基礎為何？
• RQ4 在古德雷框架中，$ \lambda = 0 $ 情形（常數衝擊速度）能否產生有界且可微的解？若不能，原因為何？
• RQ5 $ \lambda = 0 $ 情形的失敗對構建非理想材料水動力代碼驗證中解析測試問題有何影響？

主要发现

• 在任意狀態方程下，自相似收斂衝擊波解存在的必要且充分熱力學條件為 $ K_S = p f(\rho) $，其中 $ f $ 為密度的任意函數。
• 對於任意狀態方程，古德雷問題始終可歸結為非線性特徵值問題，這是衝擊穩定性與解有界性相互作用的結果。
• 標準自相似尺度（$ \lambda \neq 0 $）下，僅當滿足特徵值條件時，ODE系統才具有良好的適定性，從而確保有界且可微的解。
• $ \lambda = 0 $ 情形（常數衝擊速度）雖使ODE系統簡化，但無法產生有界且可微的解，原因在於ODE的分母趨於零而分子非零。
• $ \lambda = 0 $ 情形無法產生有效的古德雷解，因為衝擊速度無法為常數而不違反有界性條件，儘管它使ODE簡化。
• 本分析為構建非理想材料（特別是具有複雜或非理想狀態方程者）水動力代碼驗證中的精確解析測試問題提供了嚴謹基礎。