空气动力学方程：纳维-斯托克斯方程与湍流详解.pdf

空气动力学方程：纳维-斯托克斯方程与湍流详解

1空气动力学基础

1.1流体的连续性方程

流体的连续性方程是空气动力学中描述流体质量守恒的基本方程。在不可

压缩流体中，连续性方程可以简化为：

∂∂∂

++=0

∂∂∂

其中，、、分别是流体在、、方向的速度分量。这个方程表明，在

任意固定体积内，流体的质量不会随时间改变，即流入的质量等于流出的质量。

1.1.1示例

=2=−2

假设我们有一个二维流场，其中速度分量为和。我们可以验

证连续性方程是否成立。

importsympy

#定义符号变量

x,y=sympy.symbols(xy)

#定义速度分量

u=2*x

v=-2*y

#计算连续性方程的左侧

continuity_eq=sympy.diff(u,x)+sympy.diff(v,y)

#输出结果

print(continuity_eq)

运行上述代码，输出结果为0，这表明连续性方程在给定的流场中成立。

1.2伯努利方程与流体动力学

伯努利方程描述了流体在无粘性、不可压缩、稳定流动条件下的能量守恒。

方程表达为：

1

2

++ℎ=常数

2

其中，是流体的压力，是流体的密度，是流体的速度，是重力加速度，

ℎ是流体的高度。这个方程表明，在流体流动过程中，压力、动能和势能之间可

1

以相互转换，但总能量保持不变。

1.2.1示例

假设我们有两个高度不同的点，点A和点B，其中点A的高度为0，压力

为100kPa，速度为0m/s；点B的高度为10m，速度为10m/s。我们可以计算点

B的压力。

#定义常数

rho=1.225#空气密度，单位：kg/m^3

g=9.81#重力加速度，单位：m/s^2

p_A=100e3#点A的压力，单位：Pa

v_A=0#点A的速度，单位：m/s

h_A=0#点A的高度，单位：m

v_B=10#点B的速度，单位：m/s

h_B=10#点B的高度，单位：m

#计算点B的压力

p_B=p_A+0.5*rho*(v_A**2-v_B**2)+rho*g*(h_A-h_B)

#输出结果

print(f点B的压力为：{p_B:.2f}Pa)

运行上述代码，我们可以得到点B的压力值。

1.3流体的粘性与牛顿流体模型

流体的粘性是流体抵抗变形的能力，它影响流体的流动特性。牛顿流体模

型假设流体的剪切应力与速度梯度成正比，即：

∂

=

∂

∂

其中，是剪切应力，是动力粘度，是速度梯度。这个模型适用于许多

∂

常见的流体，如水和空气。

1.3.1示例

−5

假设我们有一层流体，其动力