空气动力学数值方法：格子玻尔兹曼方法(LBM)：空气动力学基础理论.pdfVIP

空气动力学数值方法：格子玻尔兹曼方法(LBM)：空气动力

学基础理论

1空气动力学基础

1.1流体力学基本方程

流体力学是研究流体（液体和气体）的运动和静止状态的学科。在空气动

力学中，我们主要关注气体的流动。流体的运动可以用一组偏微分方程来描述，

这些方程统称为纳维-斯托克斯方程（Navier-Stokesequations）。纳维-斯托克斯

方程是流体力学的基本方程，它包含了连续性方程、动量方程和能量方程。

1.1.1连续性方程和动量方程

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

表示为：

∂

+∇⋅=0

∂

其中，是流体的密度，是流体的速度向量，是时间，∇是梯度算子。

动量方程描述了流体动量的守恒，对于不可压缩流体，可以简化为：

∂2

+∇⋅⊗=−∇+

∂

其中，是流体的压力，是流体的动力粘度，是作用在流体上的外力向量。

1.1.2能量方程和状态方程

能量方程描述了流体能量的守恒，对于不可压缩流体，可以表示为：

∂

+∇⋅+=∇⋅⋅∇+⋅

∂

其中，是流体的总能量，包括内能和动能。

状态方程将流体的物理性质联系起来，对于理想气体，状态方程可以表示

为：

=

其中，是气体常数，是流体的温度。

1

1.2流体的粘性和压缩性

1.2.1流体的粘性

流体的粘性是指流体内部相邻层之间相对运动时产生的内摩擦力。粘性是

流体的一个重要属性，它影响流体的流动状态，如层流和湍流。流体的粘性可

以用动力粘度或运动粘度来描述。

1.2.2流体的压缩性

流体的压缩性是指流体在压力变化下体积发生变化的性质。对于气体，压

缩性是一个显著的特性，因为气体的体积可以显著地随压力和温度的变化而变

化。流体的压缩性可以用压缩系数或速度音速来描述。

1.2.3示例：使用Python求解二维不可压缩流体的纳维-斯托克斯方

程

importnumpyasnp

importmatplotlib.pyplotasplt

#定义网格大小和时间步长

nx,ny=100,100

nt=100

dx=2/(nx-1)

dy=2/(ny-1)

sigma=1.0

omega=1/sigma

rho=np.ones((ny,nx))

u=np.zeros((ny,nx))

v=np.zeros((ny,nx))

#定义边界条件

u[0,:]=0.0

u[-1,:]=0.0

u[:,0]=0.0

u[:,-1]=1.0

v[0,:]=0.0

v[-1,:]=0.0

v[:,0]=0.0

v[:,-1]=0.0

求解纳维斯托克斯方程

#-

2

forninrange(nt):

un=u.copy()