[論文レビュー] On the Riccati dynamics of 2D Euler-Poisson equations with attractive forcing

Yongki Lee|arXiv (Cornell University)|Mar 2, 2022

Navier-Stokes equation solutions参考文献 25被引用数 1

ひとこと要約

本稿では、速度勾配テンソルのリッカティ型ダイナミクスを分析することにより、引力的力が作用する2次元Euler-Poisson系に対して、グローバルな滑らかな解の存在を確立する。非局所的で爆発を助長する項（密度のリーマン変換に関連）の増大が指数関数的レートを超えない場合、解がグローバルに正則のままであることを証明する。3次元補助系と比較原理を用いて、初期データが広大で明示的な配置集合に属する場合の不変解空間を同定する。

ABSTRACT

The Euler-Poisson (EP) system describes the dynamic behavior of many important physical flows. In this work, a Riccati system that governs two-dimensional EP equations is studied. The evolution of divergence is governed by the Riccati type equation with several nonlinear/nonlocal terms. Among these, the vorticity accelerates divergence while others further amplify the blow-up behavior of a flow. The growth of these blow-up amplifying terms are related to the Riesz transform of density, which lacks a uniform bound makes it difficult to study global solutions of the multi-dimensional EP system. We show that the Riccati system can afford to have global solutions, as long as the growth rate of blow-up amplifying terms is not higher than exponential, and admits global smooth solutions for a large set of initial configurations. To show this, we construct an auxiliary system in 3D space and find an invariant space of the system, then comparison with the original 2D system is performed. Some numerical examples are also presented.

研究の動機と目的

• 非局所的力が作用する多次元Euler-Poisson系のグローバル正則性に関する未解決問題を解消すること。
• 爆発を助長する非局所項が存在する中でも、グローバル滑らかな解が存在する初期配置の広大な集合を同定すること。
• 速度勾配テンソル ∇u のリッカティ型進化とその成分（発散 d、渦度 ω、非等方性 η, ξ）を分析すること。
• 密度のリーマン変換が一様有界でないという課題を克服するため、3次元補助系を構築すること。
• 元の2次元系と補助的3次元系との間の比較原理を確立し、グローバル正則性を導出すること。

提案手法

• 2次元Euler-Poisson方程式から、速度勾配テンソル M = ∇u に対する閉形式のリッカティ系を導出する。
• M をスカラー成分に分解：発散 d = tr(M)、渦度 ω = M21 − M12、非等方性 η = M11 − M22、ξ = M12 + M21。
• d と ρ の挙動をモデル化する3次元補助常微分方程式系 (a(t), b(t)) を構築し、動力学 ȧ = −b a、ḃ = −½b² − eᵗa² − a + 1 を満たす。
• 補助系の解が有界かつグローバルに定義される不変領域 ΩTMB ⊂ ℝ² を特定する。
• 比較原理を適用：初期データ (ρ₀, d₀) が Ω ⊂ ΩTMB の部分集合に属する場合、d(t) は有界のまま保たれ、ρ(t) は減少するため、グローバル正則性が保証される。
• 発散方程式における非局所項を制御するため、−eᵗ ≤ A(t) ≤ ½(ω₀/ρ₀)² の境界を用いる。

実験結果

リサーチクエスチョン

• RQ1 引力的力が作用する2次元Euler-Poisson系が、初期データのどのような条件下でグローバル滑らかな解を持つのか？
• RQ2 非局所項（特に密度のリーマン変換）は、発散 d の爆発挙動にどのように影響を与えるか？
• RQ3 リーマン変換に一様有界性がないという問題を克服し、グローバル存在を証明できるか？
• RQ4 2次元リッカティ系の本質的ダイナミクスを捉える比較系が存在するか？
• RQ5 渦度と非等方性は、有限時間爆発を抑制するか、あるいは加速するか？

主な発見

• 初期データが不変領域 Ω ⊂ ΩTMB に属する広大で明示的な集合に属する場合、引力的力が作用する2次元Euler-Poisson系に対してグローバル滑らかな解が存在する。
• ρ(t) が一様に有界であれば、d(t) はすべての t ≥ 0 で上界から有界のまま保たれる。これは ρ₀ < ½ かつ初期条件が Ω に属する場合に成立する。
• 補助的3次元系 (a(t), b(t)) は、すべての t ≥ 0 で 0 < a(t) ≤ ½ および −½ ≤ b(t) ≤ max{|b₀|, √2} を満たすグローバル解を有する。
• 比較原理により、d(0) > b(0) かつ ρ(0) < a(0) ならば、すべての t > 0 で d(t) > b(t) かつ ρ(t) < a(t) が成り立ち、有界性が保たれる。
• 数値シミュレーションにより、|ρ(t)|∞ および |∆⁻¹(ρ − c)|∞ が非増加または穏やかに増加することが確認され、非局所項の穏やかな増大を支持する。
• 非局所項 A(t) = R₁₁[ρ − c] + R₂₂[ρ − c] は −eᵗ ≤ A(t) ≤ ½(ω₀/ρ₀)² を満たし、その増大が指数関数的以下に制御される。これによりグローバル存在が可能になる。