空气动力学实验方法：流动显示技术：空气动力学基础理论.docxVIP

PAGE1

PAGE1

空气动力学实验方法：流动显示技术：空气动力学基础理论

1空气动力学基础

1.1流体力学基本概念

流体力学是研究流体（液体和气体）的运动和静止状态的学科。在空气动力学中，我们主要关注气体的特性及其对物体的作用。流体的基本特性包括：

密度（ρ）:单位体积的流体质量。

压力（P）:流体垂直作用于单位面积上的力。

速度（v）:流体在某一点的运动速度。

温度（T）:流体的热状态，影响其密度和压力。

流体的运动状态可以通过流线、迹线和流管来描述。流线表示在某一时刻流体的运动方向；迹线是流体中某一质点随时间的运动轨迹；流管则是由流线构成的管状区域，其流体流量保持不变。

1.2伯努利定理与连续性方程

1.2.1伯努利定理

伯努利定理描述了在理想流体（无粘性、不可压缩）中，流体速度增加时，其压力或势能会减少。数学表达式为：

P

其中，P是压力，ρ是流体密度，v是流体速度，g是重力加速度，h是高度。

1.2.2连续性方程

连续性方程基于质量守恒原理，表示在流体流动中，流过任意截面的流体质量是恒定的。对于不可压缩流体，连续性方程简化为：

ρ

其中，A是流体流过的截面积。

1.3流体的粘性和可压缩性

1.3.1粘性

流体的粘性是指流体内部相邻层之间相对运动时的内摩擦力。粘性影响流体的流动状态，特别是在边界层中，它会导致流体速度梯度的形成。

1.3.2可压缩性

流体的可压缩性是指流体在压力变化下体积的变化。对于气体，其可压缩性远大于液体。在高速流动中，气体的可压缩性变得显著，需要使用更复杂的方程来描述其行为。

1.4流体动力学方程组

流体动力学方程组包括：

连续性方程：描述流体质量守恒。

动量方程：基于牛顿第二定律，描述流体动量守恒。

能量方程：描述流体能量守恒。

在理想流体中，这些方程可以简化为欧拉方程。但在实际应用中，流体的粘性和可压缩性必须考虑，因此使用纳维-斯托克斯方程。

1.4.1纳维-斯托克斯方程

纳维-斯托克斯方程是描述粘性流体运动的偏微分方程。在三维空间中，方程组可以表示为：

?

其中，v是流体速度向量，μ是动力粘度，f是外部力向量，E是总能量。

1.4.2欧拉方程

在忽略流体粘性的情况下，纳维-斯托克斯方程简化为欧拉方程：

?

1.4.3解方程组

解流体动力学方程组通常需要数值方法，如有限差分法、有限体积法或有限元法。以下是一个使用Python和NumPy库解决二维不可压缩流体流动的简单示例：

importnumpyasnp

importmatplotlib.pyplotasplt

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

nx,ny=100,100

nt=100

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))

#初始化压力场

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

#定义边界条件

u[0,:]=0

u[-1,:]=0

v[:,0]=0

v[:,-1]=0

#更新速度场和压力场

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])-