[論文レビュー] Analytical approximation of Blasius' similarity solution with rigorous error bounds

Ovidiu Costin, S. Tanveer|arXiv (Cornell University)|Mar 6, 2013

Fractional Differential Equations Solutions参考文献 7被引用数 3

ひとこと要約

この論文は、多項式と誤差関数表現を組み合わせたハイブリッド手法を用いて、境界層方程式のブラジウスの相似解に対する厳密な解析的近似を提示する。収縮写像とエネルギー推定を用いて高精度の誤差境界を達成し、f''(0) = 0.469600 ± 0.000022 の、証明可能な収束性とブラジウスを越える非線形常微分方程式への適用可能性を有する。

ABSTRACT

We use a recently developed method \cite{Costinetal}, \cite{Dubrovin} to find accurate analytic approximations with rigorous error bounds for the classic similarity solution of Blasius of the boundary layer equation in fluid mechanics, the two point boundary value problem $f^{\prime \prime \prime} + f f^{\prime \prime} =0$ with $f(0)=f^\prime (0)=0$ and $\lim_{x ightarrow \infty} f^\prime (x) =1$. The approximation is given in terms of a polynomial in $[0, \frac{5}{2}]$ and in terms of the error function in $[\frac{5}{2}, \infty)$. The two representations for the solution in different domains match at $x=\frac{5}{2}$ determining all free parameters in the problem, in particular $f^{\prime \prime} (0) =0.469600 \pm 0.000022 $ at the wall The method can in principle provide approximations to any desired accuracy for this or wide classes of linear or nonlinear differential equations with initial or boundary value conditions. The analysis relies on controlling the errors in the approximation through contraction mapping arguments, using energy bounds for the Green's function of the linearized problem.

研究の動機と目的

• 証明可能な誤差境界を有する高精度な解析的近似を、ブラジウスの相似解に対して開発すること。
• 非線形境界値問題における既存の経験的近似の厳密な誤差制御の欠如を克服すること。
• 明示的な小さいパラメータが存在しない非線形常微分方程式に、初期近似に基づく摂動法を拡張すること。
• 境界条件を伴う線形および非線形微分方程式の広範なクラスに適用可能なフレームワークを確立すること。

提案手法

• 区分的近似の構築：[0, 5/2] で12次多項式、[5/2, ∞) で誤差関数に基づく表現を用いる。
• x = 5/2 における両表現の一致から、f''(0) を含むすべての自由パラメータを決定する。
• 摂動方程式における残差誤差を制御するため、適応されたノルムにおける収縮写像の議論を用いる。
• 線形化問題のグリーン関数に対するエネルギー推定を用いて、誤差の伝播を制御する。
• ラプラス変換技術と漸近的解析を用いて、近似における剰余項を推定する。
• 明示的な積分表現と不等式を用いて、剰余関数 R3 と R4 を境界づけることで、一様な誤差制御を確保する。

実験結果

リサーチクエスチョン

• RQ1 ブラジウス方程式に対して、証明可能な誤差境界を有する高精度な解析的近似を構築することは可能か？
• RQ2 ブラジウス問題に小さいパラメータが存在しない状況において、摂動法を適用するためにどのような対処が可能か？
• RQ3 境界値問題において、局所的な多項式近似とグローバルな特殊関数近似を最適に一致させる方法は何か？
• RQ4 非線形常微分方程式における誤差制御において、詳細なグリーン関数解析を、より柔らかいエネルギーに基づく推定に置き換えることは可能か？
• RQ5 このフレームワークを用いて、ブラジウス問題における f''(0) の最もきつい証明可能な境界は何か？

主な発見

• 本手法により、区分的解析的近似が得られた：[0, 5/2] で12次多項式、[5/2, ∞) で誤差関数表現で、x = 5/2 で一致する。
• 壁面せん断応力は f''(0) = 0.469600 ± 0.000022 と計算され、収縮写像とエネルギー推定から導かれた厳密な誤差境界を有する。
• 剰余項 R3 と R4 は、t ≥ 1.99 のとき R3,m ≤ 0.02057 および R4,m ≤ 0.009042 で有界であり、収束を保証する。
• 本アプローチは、初期値または境界値を伴う広範な非線形および線形微分方程式のクラスに一般化可能である。
• 詳細なグリーン関数解析を、より簡単なエネルギーに基づく推定に置き換えることで、誤差制御フレームワークを単純化した。
• 解の最終的な誤差境界は、t ≥ T のとき 1.02×10⁻⁸ 以内に制御され、高い数値的信頼性を示している。