Derivation of effective kinetic equations describing oscillations in viscoelasticity and in compressible Navier-Stokes

アサナシオス・E・ツァヴァラス（Athanasios E. Tzavaras）による講義ノート「粘性および圧縮性 Navier-Stokes における振動を記述する有効運動方程式の導出」の技術的概要を以下にまとめます。

1. 問題設定 (Problem)

本論文は、双曲 - 放物型偏微分方程式系（hyperbolic-parabolic systems）の解において生じる持続的な振動（persistent oscillations）の数学的性質と、それらを記述するための有効方程式（effective equations）の導出を主題としています。

具体的には、以下の 2 つの力学系を主要な例として扱います：

1. ケルビン・フォイグト型（Kelvin-Voigt type）の粘弾性系：
   • 相転移をモデル化する際に用いられる、二重井戸ポテンシャル（double well potentials）を持つひずみエネルギーを含む系。
   • 非凸な応力 - ひずみ関係により、振動が発生し持続する現象（例：結晶の双晶化）を記述。
2. 圧縮性 Navier-Stokes 方程式（バロトロピック流体）：
   • 非単調な圧力関数（例：Van der Waals 状態方程式）を持つ流体。
   • 密度に振動が生じ、それが時間とともに消滅しないケース。

これらの系において、初期データに振動が含まれている場合、その振動が解にどのように伝播・維持されるか、そしてその巨視的な振る舞いを記述する「均質化（homogenization）」された方程式をどのように導出するかが中心的な課題です。

2. 手法とアプローチ (Methodology)

従来の均質化理論では、ヤング測度（Young measures）が弱収束の極限を記述する標準的な道具として用いられてきました。しかし、ヤング測度はパラメータ化された確率測度の族であり、そのダイナミクス（時間発展）を追跡することが困難であるという欠点があります。

本論文では、スカラー保存則の運動論的定式化（kinetic formulation）のアイデアを、双曲 - 放物型系への拡張に応用することでこの課題を解決します。

• ヤング測度の分布関数としての表現：
  振動する数列 $\{u_\varepsilon\}$ の弱極限を記述するヤング測度 $\nu_{t,x}$ に対し、その分布関数 $F(t, x, \xi) = \nu_{t,x}((-\infty, \xi])$ を導入します。これにより、確率測度という抽象的な対象を、より扱いやすい関数 $F$（0 から 1 の値を取る）として記述できます。
• 運動論的方程式の導出：
  元の非線形偏微分方程式系から出発し、テスト関数を用いた変形や、ヤング測度の性質を利用することで、分布関数 $F$ が満たす輸送型方程式（運動論的方程式）を導出します。この方程式は、巨視的な流れ（速度場や応力場）と結合しており、振動の伝播を記述します。
• 近似解からの極限：
  粘性項を含む近似解（$\varepsilon \to 0$ の極限）の挙動を解析し、その極限において有効な運動論的方程式が成立することを示します。

3. 主要な貢献と結果 (Key Contributions and Results)

A. 振動の存在と持続性の確認

• 線形系における振動：
  線形熱粘弾性系（thermoviscoelasticity）において、散逸行列が特異（singular）である場合、高周波の初期データから生じる振動が、変位や速度ではなく、ひずみ（strain）やひずみ速度、温度勾配などの特定の変数において時間的に持続することを示しました。
• 非線形系における振動の構成：
  • 粘弾性：非単調な応力関数を持つ 1 次元粘弾性棒のモデルに対し、Rankine-Hugoniot 条件を満たす定常衝撃波（standing shocks）を介して、密度（ひずみ）に不連続性を持ちながら速度が連続であるような周期的解を構成しました。
  • 圧縮性 Navier-Stokes：非単調な圧力を持つ場合、1 次元および多次元において、速度場は連続だが密度場が特性曲線（または衝撃波面）をまたいでジャンプする解を構成しました。これらの解は、初期データの振動が時間とともに減衰せず持続することを示しています。

B. 有効運動方程式の導出（均質化）

• 1 次元粘弾性モデル：
  非単調な応力 $\sigma(u)$ を持つ系に対し、ヤング測度の分布関数 $F(t, x, \xi)$ と巨視的な速度 $v$、応力 $S$ を変数とする有効方程式系を導出しました。
  得られた系は以下の形です：
  $$\begin{cases} \partial_t F + \partial_\xi \left( (v_x + \sigma(u) - \sigma(\xi)) F \right) + \sigma'(\xi)F = 0 \\ \partial_t v = \partial_{xx}v + \partial_x (\sigma(u)) \\ S = \sigma(u) + v_x = \int \sigma(\xi) dF_{t,x}(\xi) + v_x \end{cases}$$
  ここで、$F$ の伝播は運動論的方程式で記述され、その係数は巨視的なモーメント（$u, v, S$）に依存します。この系は、不安定領域を横切るデータを与えた際に、相分離（phase separation）を正しく記述することが期待されます。

• 圧縮性 Navier-Stokes 系：
  密度 $\rho$ のみが発振する変数であることに着目し、ヤング測度をその分布関数 $H(t, x, \xi)$（または $F$）で置き換えるアプローチを提案しました。

  • Hillairet [19] の結果を踏まえ、再正規化解（renormalized solutions）の弱極限に対するヤング測度の方程式を厳密に導出しました。
  • 有効粘性フラックス（effective viscous flux）の性質を利用し、分布関数 $H$ が満たす輸送方程式を導出しました。
    $$\partial_t H + \text{div}(uH) + \frac{1}{\lambda+2\mu} \xi \partial_\xi \left( H z \right) - \frac{1}{\lambda+2\mu} \xi p(\xi) \partial_\xi H = 0$$
    （ここで $z$ は有効粘性フラックス、$p$ は圧力）。
    この方程式は、密度の振動の時間発展を巨視的な流れ $u$ と結合して記述するものです。

C. 理論的基盤の強化

• タルタールの定理の別証明：
  スカラー保存則における振動の相殺（cancellation of oscillations）に関するタルタールの定理を、ヤング測度とその分布関数の対応関係を用いて簡潔に証明しました。これにより、振動が持続するのは波速が一定の領域に限られることを示唆しています。

4. 意義と重要性 (Significance)

1. ダイナミクスの記述可能性：
   従来のヤング測度アプローチは、解の存在証明や弱極限の記述には強力ですが、その時間発展（ダイナミクス）を追跡するのは困難でした。本論文で提示された「運動論的方程式（kinetic equation）」のアプローチは、振動の分布関数がどのように時間的に進化するかを記述する微分方程式を提供し、均質化問題における動的な側面を捉えるための新しい枠組みを確立しました。
2. 相転移と多相流への応用：
   非凸なポテンシャルや非単調な状態方程式を持つ系（相転移、多相流など）において、微視的な振動が巨視的な挙動にどう影響するかを記述するモデルを提供しています。特に、不安定領域をまたぐ振動がどのように相分離を起こすかという現象を、運動論的方程式を通じて理解する道筋を示しました。
3. 数学的厳密性と一般性：
   粘性項を持つ双曲 - 放物型系という、物理的に重要なクラスの問題に対して、厳密な数学的導出を行いました。これは、乱流や衝撃波の構造理解など、より複雑な流体現象の解析にも寄与する可能性があります。

総じて、本論文は「振動する解の巨視的記述」という難題に対し、ヤング測度の分布関数と運動論的定式化を融合させることで、有効な微分方程式系を導出する画期的な手法を提示した点に大きな意義があります。