空气动力学基本概念：流体力学基础：流体动力学中的控制体积法.pdfVIP

空气动力学基本概念：流体力学基础：流体动力学中的控

制体积法

1流体力学概述

1.1流体的性质

流体，包括液体和气体，具有不同于固体的特性。流体的性质主要包括：

连续性：流体可以被视为连续介质，其物理性质（如密度、压力、

速度）在空间中连续变化。

压缩性：气体可以被压缩，而液体在一般情况下被认为是不可压

缩的。

粘性：流体内部存在内摩擦力，称为粘性，它影响流体的流动状

态。

惯性：流体具有质量，因此在加速或减速时会表现出惯性。

表面张力：流体表面存在一种力，使表面尽可能缩小，这种力称

为表面张力。

1.2流体动力学的基本方程

流体动力学的基本方程是描述流体运动的数学模型，主要包括：

1.2.1连续性方程

连续性方程描述了流体质量的守恒。对于不可压缩流体，连续性方程可以

简化为：

∂∂∂

++=0

∂∂∂

其中，、、分别是流体在、、方向的速度分量。

1.2.2动量方程

动量方程，即纳维-斯托克斯方程，描述了流体动量的守恒。对于不可压缩

流体，无粘性流动（理想流体）的动量方程简化为欧拉方程：

∂1

+⋅∇=−∇+

∂

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

量。

1

1.2.3能量方程

能量方程描述了流体能量的守恒，包括动能和内能。对于不可压缩流体，

能量方程可以表示为：

∂

+∇⋅=∇⋅−+⋅

∂

其中，是单位质量的总能量，是热传导矢量。

1.2.4示例：使用Python求解一维不可压缩流体的连续性方程

假设我们有一个简单的管道流动问题，流体在管道中沿轴流动，我们可以

通过数值方法求解连续性方程。

importnumpyasnp

importmatplotlib.pyplotasplt

#定义网格

nx=101

dx=2/(nx-1)

nt=25

dt=0.025

c=1

#初始化速度分布

u=np.ones(nx)

u[int(.5/dx):int(1/dx+1)]=2

#定义边界条件

u[0]=1.0

u[-1]=1.0

#时间迭代

forninrange(nt):

un=u.copy()

u[1:-1]=un[1:-1]-c*dt/dx*(un[1:-1]-un[:-2])

#绘制结果

plt.plot(np.linspace(0,2,nx),u)

plt.show()

代码解释：-我们首先定义了网格参数，包括网格点数、网格间距、时间

0.512

步长等。-初始化速度分布，其中在到的区域，流体速度为，其余区域

1

为。-通过时间迭代，使用显式差分格式更新速度分布。-最后，使用

matplotlib绘制速度分布随时间变化的结果。

此代码示例展示了如何使用Py