空气动力学方程：层流和湍流模型：湍流基本概念与特性.pdfVIP

空气动力学方程：层流和湍流模型：湍流基本概念与特性

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求解三维不可压缩流体流动的基本步骤概述：