基于中心差分法的Blasius边界层流动数值分析及其分数阶导数拓展研究

《Next Nanotechnology》：Numerical analysis of Blasius boundary layer flow over a flat plate using the central difference method

　　

本研究探讨了不可压缩Navier-Stokes方程的数学表述，并结合平板诱导的边界层流动进行分析。首先通过变换推导出描述空间分数阶导数的Blasius方程。研究人员给出了边界层区域内Blasius方程的显式解，并通过MATLAB软件递增地呈现数值结果。数值结果采用打靶法（shooting method）获得，该方法结合了移位Grünwald公式（shifted Grünwald formula）和经典RK-Five（Runge-Kutta五阶）方法作为迭代方案。与采用适当小步长的其他方法相比，这一新解为验证所有用于求解设计领域微分方程的方法的精度水平提供了准确的参考基准。因此，导数近似方法为解决与Blasius方程具有相似难度的其他微分方程提供了一种新颖、简便的途径。根据基本流廓线研究，在相同外部流动条件下，黏弹性流体（viscoelastic fluid）在表面边界层中的流动速度比黏性流体（viscous fluid）更快。随着参数值的减小，黏弹性流体更靠近平板运动。

论文解读文章

一、研究背景与意义

船舶水动力学是流体力学在船舶设计中的直接应用，它分析船舶如何改变两种连续介质共享的表面，其中一种介质在水中变化，另一种在空气中变化。因此，在运动过程中，水对浸没船体提供水动力阻力，而空气对上层建筑提供空气动力阻力。工业中常见的流体普遍具有黏弹性（viscoelastic）性质，即材料同时表现出弹性和黏性特性。Pelton等人指出液体的黏弹性特性会导致基于实验数据和Navier-Stokes方程的数值预测之间存在显著差异。Chakraborty和Sader讨论了修正的线性黏弹性流动，该流动在低频率下表现为流体性质，在高频率下表现为固体性质。Meral等人提出了在流固耦合问题中经常表现非牛顿性质和黏弹性性质的流体。Alqahtani等人研究了通过线性收缩或膨胀表面的纳米流体的稳态非对齐驻点流动，考虑了温度和速度滑移。边界层假设是数学模型建立的基础。

分数阶导数是经典导数的推广，允许对任意实数阶定义积分和导数。目前已提出多种不同类型的分数阶导数，包括Riemann-Liouville分数阶导数、Hadamard分数阶导数、Caputo导数和Grünwald-Letnikov导数等。Bagley和Torvik为时间分数阶导数的适用性提供了理论依据。Wharmby和Bagley考虑了受黏弹性材料谱弛豫模式影响的分数阶导数。Jannelli提出了非线性两点上Blasius流动的Caputo导数。Huang和Shen解决了任意有界域上分数阶Laplacian的数值逼近问题，他们使用了非结构化简化网格和均匀网格覆盖底层域，并基于均匀网格提供有限差分逼近。近期，Bas等人介绍了Riemann-Liouville分数阶导数和积分。Elnady等人提出了新的广义局部分数阶导数，涵盖了现有的几种分数阶导数，该新导数在保持整数阶微积分标准特征的同时扩展了分数阶上下文定理。Tolami等人研究了发散通道流动中壁面注入对浓度廓线的分数阶导数影响。

本研究的主要创新点在于采用中心有限差分方法对Blasius边界层问题进行数值研究，并将其可能应用于分数阶导数。与传统主要使用标准打靶技术或Runge-Kutta求解器获得经典Blasius解的方法不同，本研究聚焦于控制非线性方程的离散化、边界条件的实施以及数值精度和收敛性的评估。经典Blasius解验证了所提数值程序的可靠性。该研究还考察了移位Gründwald近似在分数阶导数中的可能应用，为研究经典边界层模型的非局部数学拓展提供了一个框架。因此，主要贡献是对所选离散化技术及其在流体力学边界层计算中的适用性进行了系统的数值评估。

二、关键技术方法

研究人员采用了以下主要关键技术方法：

1. 1.
   Prandtl边界层方程推导：从二维稳态不可压缩Navier-Stokes方程出发，利用边界层假设简化得到Prandtl边界层方程组，包括连续性方程和动量方程。
2. 2.
   相似变换与Blasius方程：引入流函数和相似变量，将偏微分方程转化为常微分形式的Blasius方程 f? + ff″ = 0。
3. 3.
   中心差分法离散化：采用中心有限差分格式对Blasius方程进行空间离散化，将非线性常微分方程转化为代数方程组。
4. 4.
   打靶法与RK-Five方法：将边界值问题转化为初值问题，结合打靶法和经典五阶Runge-Kutta方法进行迭代求解。
5. 5.
   移位Grünwald公式：针对分数阶导数情形，采用移位Grünwald公式对分数阶导数进行数值逼近，实现分数阶Blasius方程的数值求解。

三、研究结果

2.1 问题公式化（边界区域方程）

研究人员考察了半无限厚平板上二维稳态黏弹性流体流动。沿(x, y)笛卡尔坐标方向的速度分量为 u* = (u, v)，密度为 ρ，压力为 p。假设平板周围流动为平面流动，z方向速度分量设为零。给出了二维稳态Prandtl边界层方程，包括连续性方程 ?_x*u + ?_y*v* = 0 和动量方程 u?_x*u + v?_y*u = (1/ρ)[-?_x*p + ?_y*τ_x*y]。

数值方法与求解策略

研究人员通过相似变换将边界层偏微分方程转化为Blasius常微分方程，然后采用中心差分法对该非线性方程进行离散化处理。对于经典整数阶Blasius方程，通过打靶法结合RK-Five方法获得了数值解，并与已知的经典Blasius解进行对比验证了方法的可靠性。在此基础上，研究人员进一步将方法扩展到分数阶导数情形，采用移位Grünwald公式近似分数阶导数，建立了分数阶Blasius方程的数值求解框架。

结果分析与讨论

根据基本流廓线研究，在相同外部流动条件下，黏弹性流体在表面边界层中的流动速度比黏性流体更快。随着参数值的减小，黏弹性流体更靠近平板运动。数值结果表明，所提出的中心差分法结合打靶法的方案能够有效求解Blasius方程，并且对于分数阶导数情形同样适用。与传统的纯打靶法或Runge-Kutta方法相比，该方法在处理非线性边界层问题时展现出良好的精度和稳定性。通过对比不同步长下的数值结果，验证了算法的收敛性。

四、讨论与结论

讨论部分：本研究的核心价值在于提出了一种结合中心差分法和打靶法的数值框架，不仅适用于经典Blasius方程的求解，还为分数阶边界层问题的研究提供了可行途径。与现有文献中主要依赖纯打靶技术或Runge-Kutta求解器的方法相比，本文的离散化策略在数值精度和收敛性评估方面更加系统化。研究结果证实了黏弹性流体在边界层内的流动特性与黏性流体存在本质差异，这对理解复杂流体的近壁面行为具有重要意义。该数值方法为解决具有类似非线性特征的微分方程提供了新的思路，特别是在工程设计中需要高精度参考解的场合。

研究结论：本研究通过中心有限差分方法对Blasius边界层问题进行了数值研究，并将其成功拓展至分数阶导数领域。经典Blasius解验证了所提数值程序的可靠性。该研究提供了一个研究经典边界层模型非局部数学拓展的框架。主要贡献是对所选离散化技术及其在流体力学边界层计算中的适用性进行了系统的数值评估。根据基本流廓线研究，在相同外部流动条件下，黏弹性流体在表面边界层中的流动速度比黏性流体更快。随着参数值的减小，黏弹性流体更靠近平板运动。