[論文レビュー] Slow motion for compressible isentropic Navier--Stokes equations

Corrado Mascia, Marta Strani|arXiv (Cornell University)|Mar 22, 2013

Navier-Stokes equation solutions参考文献 11被引用数 4

ひとこと要約

本稿では、有界区間上の圧縮性断熱ナビエ＝ストークス流れにおける粘性遷移層の遅い運動を記述する微分方程式を、近似定常状態の多様体とスペクトル射影を用いて導出する。主な結果は、層の速度が粘性εに対して指数的に小さくなることを示しており、衝撃型解における境界起因のずれを定量的に明らかにする。

ABSTRACT

We consider the compressible Navier-Stokes equations for isentropic dynamics with real viscosity on a bounded interval. In the case of boundary data defining an admissible shock wave for the corresponding unviscous hyperbolic system, we determine a scalar differential equation describing the motion of the internal transition layer. In particular, for small viscosity, the velocity of the motion is exponentially small. The approach is based on the construction of a one-parameter manifold of approximate solutions and on an appropriate projection of the evolution of the complete Navier-Stokes system towards such manifold.

研究の動機と目的

• 有界区間上における圧縮性断熱ナビエ＝ストークス方程式における内部遷移層の遅い運動を記述すること。
• 粘性εが小さいときの境界条件が、粘性衝撃プロファイルの力学に与える影響を定量すること。
• 1パラメータの近似解の多様体を用いて、遷移層の位置ξ(t)に関するスカラー常微分方程式を導出すること。
• 層の運動がεに関して指数的に遅いこと、すなわち準安定的挙動と整合することを確立すること。

提案手法

• x=ξにおける連続性条件を満たすように、(−ℓ,ξ)および(ξ,ℓ)における正確な解をマッチングすることにより、近似定常状態W(·,ξ)からなる1パラメータ多様体を構築する。
• 各W(·,ξ)のまわりでナビエ＝ストークス系を線形化し、演算子𝒪_ξを定義する。このとき、その第1固有値λ₁(ξ)が実数かつ単純であるものとする。
• 随伴演算子𝒪_ξ*の第1固有関数(φ(·,ξ),ψ(·,ξ))へのスペクトル射影を用いて、ξ(t)の時間発展方程式を導出する。
• κ₊(ξ)−κ₋(ξ)を含むdξ/dtの明示的表現を導出し、εが小さい場合の漸近的積分による近似を用いる。
• u₊およびu₋の近傍でF(u,vₐ)および∂ᵤG(u,vₐ)を近似し、dξ/dtの式に指数的減衰・増大項を導入する。
• 最終的に、式(37)としてdξ/dtの式を、2つの指数的に小さい項の和として得る。各項はξに依存する重みおよび第1固有値λ₁(ξ)で調制されている。

実験結果

リサーチクエスチョン

• RQ1 境界データは、粘性圧縮性ナビエ＝ストークス系における内部遷移層の遅いずれをどのように誘導するか？
• RQ2 粘性ε→0のとき、層の運動速度の明確な漸近的形は何か？
• RQ3 遷移層の力学は、近似解の多様体上のスカラーODEで記述可能か？
• RQ4 線形化演算子𝒪_ξの第1固有値λ₁(ξ)の符号と大きさは何かを決定する要因は何か？また、安定性にどのように影響するか？
• RQ5 この形式的アプローチは、非断熱的または半直線設定へどの程度拡張可能か？

主な発見

• 遷移層の運動は、境界効果による遅いずれを記述するスカラーODE (37) に従う。
• 層の速度dξ/dtはεに関して指数的に小さく、主要な振るまいはe^{−C/ε}と指数的に減衰する2つの項が支配的である。
• dξ/dtへの主な寄与は、勾配がやや小さい境界付近の指数的減衰項に起因し、その率は∂ᵤF⁺/∂ᵤG⁺によって決定される。
• dξ/dtの式には、(u₊−uₐ)および(uₐ−u₋)を含む係数が含まれており、衝撃構造の非対称性を反映している。
• 線形化演算子𝒪_ξの第1固有値λ₁(ξ)は運動の振幅を制御し、その符号は負であると予想され、安定性を示唆する。
• 導出された方程式は形式的に妥当で柔軟性を有しており、非断熱的ナビエ＝ストークスを含む他の双曲型–放物型系への応用が期待される。