[论文解读] Exact solution of a differential problem in analytical fluid dynamics: use of Airy's functions

Gianluca Argentini|ArXiv.org|Jun 28, 2006

Fluid Dynamics and Turbulent Flows参考文献 7被引用 4

一句话总结

本文通过沿流线将控制方程约化为Riccati型常微分方程，提出了一种二维定常可压缩Navier-Stokes边值问题的精确解析解。利用Airy函数的变换，解被表示为Ai(t)和Bi(t)及其导数的有理组合，从而在恒定运动粘度和压力-密度耦合的特定本构假设下，给出了封闭形式的速度分布。

ABSTRACT

Treating a boundary value problem in analytical fluid dynamics, translation of 2D steady Navier-Stokes equations to ordinary differential form leads to a second order equation of Riccati type. In the case of a compressible fluid with constant kinematic viscosity along streamlines, it is possible to find an exact solution of the differential problem by rational combination of Airy's functions and their derivatives.

研究动机与目标

• 推导二维定常可压缩Navier-Stokes边值问题沿流线的精确解析解。
• 通过流线参数化与速度分解，将Navier-Stokes方程约化为常微分方程。
• 在恒定运动粘度及压力-密度耦合的特定本构关系假设下，利用Airy函数求解所得的Riccati型方程。
• 通过初始与终值速度条件确定三个积分常数，建立边值解。

提出的方法

• 使用标量参数s对流线进行参数化，假设v₁ ≠ 0且φ₁(s) = s以简化问题。
• 利用链式法则与速度分解u₁ = v₁(φ(s))，将Navier-Stokes方程转化为沿流线的常微分方程。
• 应用Riccati变换u₁ = -2ν (ż/z)，将非线性二阶常微分方程转化为Airy型的线性二阶常微分方程。
• 利用Airy函数Ai(t)与Bi(t)求解变换后的方程，其中t为涉及粘度与压力-密度参数的s的线性函数。
• 重构速度解为u₁ = 2ν (-a)^{1/3} (c₁Ai’(t) + c₂Bi’(t)) / (c₁Ai(t) + c₂Bi(t))，以确保在a < 0时的物理解释有效性。
• 根据边界条件u₁(0) = u₁₀, ẋ₁(0) = ẋ₁₀, 且u₁(L) = u₁L，确定常数c₁, c₂, 与c。

实验结果

研究问题

• RQ1 能否利用特殊函数为二维定常可压缩Navier-Stokes方程沿流线推导出精确解析解？
• RQ2 当运动粘度及压力梯度与密度之比沿流线恒定时，速度分布的形式是什么？
• RQ3 如何利用Airy函数精确求解由Navier-Stokes方程约化得到的Riccati型常微分方程？
• RQ4 在何种物理解释条件下，所推导的解保持为实数且具有物理意义？
• RQ5 如何通过三重边界条件确定积分常数以满足三点边值问题？

主要发现

• 沿流线的Navier-Stokes方程的精确解可表示为Airy函数及其导数的有理函数形式。
• 当参数a = (1/(2ν²))((q̇/ρ) - f₁)为负时，解有效，这在压力减小（q̇ < 0）与正体积力（f₁ > 0）的合理假设下成立。
• 速度分量u₁仅通过u₂分量依赖于流线，后者通过dφ₂/ds对u₁进行缩放。
• 具有三个条件（u₁(0), ẋ₁(0), u₁(L)）的边值问题通过代数方程组唯一确定三个常数c, c₁, 与c₂。
• 当流线为直线且流体不可压缩时，解退化为恒定速度场，与经典流体力学一致。
• 该方法为具有恒定运动粘度及特定压力-密度本构关系的一类可压缩、定常、二维流动提供了封闭形式的解析解。