空气动力学方程：层流和湍流模型：湍流基本概念与特性.pdf
空气动力学方程：层流和湍流模型：湍流基本概念与特性
1空气动力学基础
1.1流体动力学基本方程
流体动力学是空气动力学的基础，它研究流体（液体和气体）的运动规律。
在流体动力学中，有两个基本方程至关重要：连续性方程和动量方程。这些方
程描述了流体在不同条件下的行为，是理解和分析空气动力学现象的关键。
1.1.1连续性方程
连续性方程基于质量守恒原理，即在没有质量源或汇的情况下，流体通过
任意封闭区域的质量流量必须保持恒定。对于不可压缩流体，连续性方程可以
表示为：
∂
+∇⋅=0
∂
∇⋅
其中，是流体的密度，是流体的速度向量，是散度算子。对于不可
压缩流体，密度可以视为常数，因此方程简化为：
∇⋅=0
1.1.2动量方程
动量方程，也称为纳维-斯托克斯方程，描述了流体的动量变化。它基于牛
顿第二定律，即力等于质量乘以加速度。对于不可压缩流体，动量方程可以表
示为：
∂
+⋅∇−2
=∇+∇+
∂
其中，是流体的密度，是流体的速度向量，流体的压力，流体的
∇∇2
动力粘度，是作用在流体上的外力向量，是梯度算子，是拉普拉斯算子。
1.1.3示例：使用Python求解二维不可压缩流体的连续性方程
假设我们有一个二维不可压缩流体的流动，其中速度场由=
,给出。我们可以使用Python和NumPy库来求解连续性方程。
importnumpyasnp
#定义网格尺寸
nx,ny=100,100
x=np.linspace(0,1,nx)
y=np.linspace(0,1,ny)
1
X,Y=np.meshgrid(x,y)
#定义速度场
u=np.sin(2*np.pi*X)*np.cos(2*np.pi*Y)
v=-np.cos(2*np.pi*X)*np.sin(2*np.pi*Y)
#计算散度
div_u=np.gradient(u,axis=0)
div_v=np.gradient(v,axis=1)
divergence=div_u+div_v
#输出结果
print(连续性方程的解（散度）：)
print(divergence)
在这个例子中，我们定义了一个二维网格，并在该网格上创建了一个速度
场。然后，我们使用NumPy的gradient函数来计算速度场的散度，即连续性方
∇⋅=0
程的解。由于我们定义的速度场满足连续性方程（），因此计算出的散
度应该接近于零。
1.2连续性方程与动量方程
连续性方程和动量方程是相互关联的，它们共同描述了流体的运动。连续
性方程确保了流体的质量守恒，而动量方程则描述了流体动量的变化，包括压
力、粘性力和外力的影响。在实际应用中，这两个方程通常需要同时求解，以
获得流体的完整运动状态。
1.2.1示例：使用OpenFOAM求解三维不可压缩流体的纳维-斯托克
斯方程
OpenFOAM是一个开源的CFD（计算流体动力学）软件包，可以用来求解
复杂的流体动力学问题，包括三维不可压缩流体的纳维-斯托克斯方程。下面是
一个使用OpenFOAM求解三维不可压缩流体流动的基本步骤概述：