[论文解读] 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

一句话总结

本文提出了一种结合多项式与误差函数表示的混合方法，对边界层方程的Blasius相似解进行了严格的解析近似。通过压缩映射和能量估计，实现了高精度的误差界，得到 f''(0) = 0.469600 ± 0.000022，且具有可证明的收敛性，该方法可推广至Blasius以外的非线性常微分方程。

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.

研究动机与目标

• 开发一种针对Blasius相似解的高精度解析近似，并附带可证明的误差界。
• 克服现有经验近似在非线性两点边值问题中缺乏严格误差控制的缺陷。
• 将基于初值近似的摄动方法推广至无显式小参数的非线性常微分方程。
• 建立一个适用于广泛线性与非线性微分方程（含初值或边值条件）的通用框架。

提出的方法

• 构造分段近似：在 [0, 5/2] 上使用12次多项式，在 [5/2, ∞) 上使用基于误差函数的表达式。
• 在 x = 5/2 处匹配两部分表达式，以确定所有自由参数，包括 f''(0)。
• 在自适应范数下使用压缩映射论证，控制摄动方程中的残差误差。
• 利用线性化问题的格林函数的能量界，控制误差传播。
• 采用拉普拉斯变换技术和渐近分析，估计近似中余项的大小。
• 通过显式积分表达式和不等式，对余项函数 R3 和 R4 进行有界控制，确保一致的误差控制。

实验结果

研究问题

• RQ1 能否为Blasius方程构造一个具有严格误差界的高精度解析近似？
• RQ2 在Blasius问题中缺乏小参数的情况下，如何应用摄动方法？
• RQ3 在两点边值问题中，如何最优地匹配局部多项式近似与全局特殊函数近似？
• RQ4 能否用基于能量的估计替代详细的格林函数分析，以实现非线性常微分方程的误差控制？
• RQ5 在此框架下，f''(0) 的最紧致可证明误差界是什么？

主要发现

• 该方法生成了分段解析近似：在 [0, 5/2] 上为12阶多项式，在 [5/2, ∞) 上为误差函数表达式，且在 x = 5/2 处实现匹配。
• 壁面剪切应力计算为 f''(0) = 0.469600 ± 0.000022，其严格误差界由压缩映射和能量估计推导得出。
• 余项 R3 和 R4 满足 R3,m ≤ 0.02057 与 R4,m ≤ 0.009042（当 t ≥ 1.99 时），确保了收敛性。
• 该方法可推广至包含初值或边值条件的广泛类线性与非线性微分方程。
• 该方法用更柔和的能量估计替代了详细的格林函数分析，简化了误差控制框架。
• 在 t ≥ T 时，解的最终误差界被控制在 1.02×10⁻⁸ 以内，表现出极高的数值可靠性。