[論文レビュー] Div-Curl Problems and Stream Functions in 3D Lipschitz Domains

Matthias Kirchhart, Erick Schulz|arXiv (Cornell University)|Jan 1, 2020

Advanced Numerical Methods in Computational Mathematics参考文献 1被引用数 2

ひとこと要約

本稿では、3次元リプシッツ領域における発散なし速度場 U ∈ L²(Ω) を、与えられた渦度 F = curl U から回復するための適切に定式化されたストリーム関数の構成を提案する。解は、U = curl A を満たすベクトルポテンシャル A ∈ H¹(Ω) を用いて得られ、適切な関数解析的条件下で存在性と一意性が確立される。

ABSTRACT

We consider the problem of recovering the divergence-free velocity field $\mathbf{U}\in\mathbf{L}^2(\Omega)$ of a given vorticity $\mathbf{F}=\mathrm{curl}\,\mathbf{U}$ on a bounded Lipschitz domain $\Omega\subset \mathbb{R}^3$. To that end, we solve the 'div-curl problem' for a given $\mathbf{F}\in\bigl[ \mathbf{H}_0(\mathrm{curl};\Omega)\bigr]'$. The solution is given in terms of a vector potential (or stream function) $\mathbf{A}\in\mathbf{H}^1(\Omega)$ such that $\mathbf{U}=\mathrm{curl}\,{\mathbf{A}}$. After discussing existence and uniqueness of solutions and associated vector potentials, we propose a well-posed construction for the stream function. A numerical example of the construction is presented at the end.

研究の動機と目的

• 有界な3次元リプシッツ領域 Ω において、与えられた渦度 F = curl U から発散なし速度場 U ∈ L²(Ω) を再構成する逆問題に対処すること。
• F ∈ [H₀(curl; Ω)]′ に対して、発散-回転問題の解の存在性と一意性を確立すること。
• U = curl A を満たす適切に定式化されたベクトルポテンシャル A ∈ H¹(Ω) を構成し、解が安定的かつ物理的に意味のあるものであることを保証すること。
• 数値的に取り扱いやすい方法でストリーム関数 A を構築するための構成的枠組みを提供し、数値的例を提示すること。

提案手法

• 一般化済み渦度場 F を取り扱うために、[H₀(curl; Ω)]′ の双対空間に発散-回転問題を定式化すること。
• ソボレフ空間における関数解析的道具を用いて、与えられた F に対して L²(Ω) 内の解 U の存在性と一意性を証明すること。
• U = curl A を満たすベクトルポテンシャル A ∈ H¹(Ω) を構成し、U が構成上発散なしであることを保証すること。
• A に対して適切に定式化された変分問題を確立し、H¹(Ω) のノルムにおける安定性と収束性を保証すること。
• ガラーキン型アプローチまたは有限要素法を用いて数値計算を実施し、数値的例により検証すること。

実験結果

リサーチクエスチョン

• RQ1 3次元リプシッツ領域において、発散なし速度場 U ∈ L²(Ω) は、その渦度 F = curl U から一意に再構成可能か？
• RQ2 F ∈ [H₀(curl; Ω)]′ に対して、U = curl A を満たすベクトルポテンシャル A ∈ H¹(Ω) の存在性と一意性を保証する条件は何か？
• RQ3 不規則な3次元領域において、ストリーム関数 A の安定的かつ適切に定式化された構成はどのように達成可能か？
• RQ4 一般化渦度を扱う発散-回転問題の定式化において、双対空間 [H₀(curl; Ω)]′ の役割は何か？
• RQ5 ストリーム関数 A は、保証された収束性と精度を満たすように、どのように数値的に計算可能か？

主な発見

• F ∈ [H₀(curl; Ω)]′ に対して、ストリーム関数 A ∈ H¹(Ω) の適切に定式化された構成が確立され、U = curl A が L²(Ω) 内で一意解であることが保証される。
• リプシッツ領域上のソボレフ空間における関数解析的手法を用いて、解 U の存在性と一意性が証明される。
• ベクトルポテンシャル A は必要な正則性と境界条件を満たし、U が発散なしでかつ L²(Ω) に属することを保証する。
• 本手法は安定的かつ数値的に実装可能な枠組みを提供し、構成プロセスの数値的例により実証されている。
• 本手法は、古典的なストリーム関数理論を、リプシッツ領域を含む境界がやや不規則な3次元領域へと拡張し、流体力学および偏微分方程式理論における応用範囲を広げている。