空气动力学实验方法：流动显示技术与流体动力学数值模拟教程.docxVIP

PAGE1

PAGE1

空气动力学实验方法：流动显示技术与流体动力学数值模拟教程

1空气动力学基础

1.1流体力学基本方程

流体力学基本方程是描述流体运动的数学模型，主要包括连续性方程、动量方程和能量方程。这些方程基于质量守恒、动量守恒和能量守恒的原理，是空气动力学分析和计算的基础。

1.1.1连续性方程

连续性方程描述了流体质量的守恒。在不可压缩流体中，连续性方程简化为：

?

其中，u、v、w分别是流体在x、y、z方向的速度分量。

1.1.2动量方程

动量方程，即纳维-斯托克斯方程，描述了流体动量的守恒。在简化的情况下，二维不可压缩流体的动量方程可以表示为：

?

?

其中，ρ是流体密度，p是压力，ν是动力粘度。

1.1.3能量方程

能量方程描述了流体能量的守恒，包括动能和内能。在不可压缩流体中，能量方程可以简化为：

?

其中，T是温度，α是热扩散率，μ是动力粘度。

1.2空气动力学中的流体特性

空气动力学研究中，流体特性对流动行为有重要影响。主要特性包括：

密度：空气的密度随温度和压力变化，影响流体的流动状态。

粘度：空气的粘度决定了流动的层流或湍流状态。

压缩性：在高速流动中，空气的压缩性效应显著，需考虑其对流动的影响。

热导率：热导率影响流体的热传递，对热流体流动分析至关重要。

1.3流动类型与分类

空气动力学中的流动类型多样，主要分类包括：

层流与湍流：层流流动中，流体层间互不干扰；湍流流动中，流体层间存在复杂的相互作用。

亚音速与超音速：亚音速流动中，流体速度低于音速；超音速流动中，流体速度超过音速，产生激波。

定常与非定常：定常流动中，流体参数不随时间变化；非定常流动中，流体参数随时间变化。

不可压缩与可压缩：不可压缩流动中，流体密度被视为常数；可压缩流动中，流体密度随压力和温度变化。

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

importnumpyasnp

importmatplotlib.pyplotasplt

#定义网格

nx=101

ny=101

nt=100

nit=50

dx=2/(nx-1)

dy=2/(ny-1)

sigma=.1

nu=.05

dt=sigma*dx*dy/nu

#初始化速度场

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

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

#边界条件

u[0,:]=0

u[-1,:]=0

v[:,0]=0

v[:,-1]=0

#内部边界条件

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

#求解纳维-斯托克斯方程

forninrange(nt):

un=u.copy()

vn=v.copy()

u[1:-1,1:-1]=(un[1:-1,1:-1]-

un[1:-1,1:-1]*dt/dx*(un[1:-1,1:-1]-un[1:-1,0:-2])-

vn[1:-1,1:-1]*dt/dy*(un[1:-1,1:-1]-un[0:-2,1:-1])-

dt/(2*rho*dx)*(p[1:-1,2:]-p[1:-1,0:-2])+

nu*(dt/dx**2+dt/dy**2)*(un[1:-1,2:]-2*un[1:-1,1:-1]+un[1:-1,0:-2]+un[2:,1:-1]-2*un[1:-1,1:-1]+un[0:-2,1:-1]))

v[1:-1,1:-1]=(vn[1:-1,1:-1]-

un[1:-1,1:-1]*dt/dx*(vn[1:-1,1:-1]-vn[1:-1,0:-2])-

vn[1:-1,1:-1]*dt/dy*(vn[1:-1,1:-1]-vn[0:-2,1:-1])-

dt/(2*rho*dy)*(p[2:,1:-1]-p[0:-2,1:-1])+

nu*(dt/