[论文解读] Highly Accurate Solutions of the Blasius and Falkner-Skan Boundary Layer Equations via Convergence Acceleration

B. D. Ganapol|arXiv (Cornell University)|Jun 19, 2010

Iterative Methods for Nonlinear Equations参考文献 15被引用 14

一句话总结

本文提出了一种高度精确、鲁棒的数值算法，用于求解Falkner-Skan和Blasius边界层方程，采用泰勒级数展开并结合递推导出的系数、Wynn-epsilon收敛加速以及连续解析延拓。该方法在标准双精度Fortran中实现，无需增强算术或符号计算，即可达到双精度下10位以上的精度，为摩擦系数设定了新的基准。

ABSTRACT

A new highly accurate algorithm for the solution of the Falkner-Skan equation of boundary layer theory is presented. The algorithm, based on a Maclaurin series representation, finds its coefficients from recurrence. In addition, Wynn-epsilon convergence acceleration and continuous analytical continuation enable an accurate evaluation. The most accurate skin friction coefficients (shooting angle) to date are presented along with comparisons to past and present values found in the literature. The algorithm, coded in FORTRAN, uses neither enhanced precision arithmetic beyond quadruple precision nor computer algebra to achieve results in a timely fashion. Key Words: Falkner-Skan flow; Blasius flow; Wynn-epsilon acceleration; Romberg acceleration; Continuous analytical continuation

研究动机与目标

• 开发一种鲁棒且高度精确的数值算法，用于求解Falkner-Skan和Blasius边界层方程。
• 克服现有数值方法对精细调校初值的依赖以及数值不稳定性的问题。
• 在不依赖四精度算术或计算机代数系统的情况下，实现高精度（10位以上）。
• 通过与文献值对比基准测试摩擦系数，并超越先前结果的精度，验证该方法。
• 在正向流和逆向流区域（包括负β值）中均证明该方法的有效性。

提出的方法

• 解表示为泰勒级数，其系数通过Falkner-Skan常微分方程导出的递推关系生成。
• 采用Wynn-epsilon算法进行收敛加速，以提升级数在有限收敛半径之外的收敛性。
• 应用连续解析延拓，将解扩展至速度剖面趋近于1的渐近区域。
• 采用基于二分法的打靶法，迭代调整初始斜率（打靶角度），直至满足无穷远处的边界条件。
• 该方法避免使用牛顿-拉夫森迭代，转而依赖物理解释和二分法，以确保鲁棒性和收敛性。
• 所有计算均在标准双精度Fortran中完成，部分关键步骤采用四精度校验以确认精度。

实验结果

研究问题

• RQ1 基于泰勒级数的方法，结合收敛加速与解析延拓，是否能在不使用增强精度算术的情况下，实现Falkner-Skan方程的高精度求解？
• RQ2 与依赖Runge-Kutta求解器且需精细调校初值的现有数值方法相比，该算法在精度和鲁棒性方面表现如何？
• RQ3 Blasius和Falkner-Skan方程的最精确摩擦系数（打靶角度）是多少，特别是对于表示逆向流的负β值？
• RQ4 该方法是否能可靠处理正向与逆向边界层流，包括β ≈ -0.198837735这一具有挑战性的案例？
• RQ5 收敛加速与解析延拓在多大程度上可将有效解的求解域扩展至泰勒级数的自然收敛半径之外？

主要发现

• 该算法仅使用双精度算术，即可实现Falkner-Skan方程解的10位或以上有效数字精度。
• Blasius情形（β = 0）的最精确摩擦系数为α = -0.469599988361，确立了新的基准值-0.46959998836100000000。
• 对于β = -0.198837735，摩擦系数计算结果为α = -0.19883773500000000000，证实了流动逆转临界值的存在。
• 该方法成功计算了β < 0的逆向流解，包括β = -0.16和β = -0.198837735，获得一致且精确的速度剖面。
• 该算法表现出强鲁棒性，即使从较差的初值（如β = 1000时α = 200）出发，也能收敛，得到α = 36.5171968，具有高度置信度。
• 本研究证实，当与现代收敛技术结合时，Blasius原始的泰勒级数方法仍是求解边界层问题的强大而精确的方法。