空气动力学实验方法：风洞实验与流体动力学数值模拟教程.docxVIP

PAGE1

PAGE1

空气动力学实验方法：风洞实验与流体动力学数值模拟教程

1空气动力学基础

1.1流体力学基本原理

流体力学是研究流体（液体和气体）的运动和静止状态的学科。在空气动力学中，我们主要关注气体的流动特性。流体的运动可以由连续介质假设来描述，即认为流体在微观上由无数连续分布的流体质点组成。

1.1.1连续性方程

连续性方程描述了流体质点的守恒定律，即在任意固定体积内，流体的质量是守恒的。对于不可压缩流体，连续性方程可以简化为：

?

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

1.1.2动量方程

动量方程，即纳维-斯托克斯方程，描述了流体运动中力与加速度的关系。对于不可压缩流体，无粘性流动的简化形式为：

?

?

?

其中，ρ是流体密度，p是流体压力，gx、gy、gz是外力在x、y

1.1.3能量方程

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

?

其中，E是总能量，μ是流体的动力粘度，q是单位体积的热源。

1.2空气动力学中的关键概念

1.2.1雷诺数

雷诺数（Reynoldsnumber）是流体力学中的一个重要无量纲数，用于判断流体流动的性质，是层流还是湍流。雷诺数的计算公式为：

R

其中，ρ是流体密度，v是流体速度，L是特征长度，μ是流体的动力粘度。

1.2.2升力和阻力

在空气动力学中，升力（Lift）和阻力（Drag）是两个关键的力。升力是垂直于流体流动方向的力，而阻力则是沿着流体流动方向的力。这些力的计算通常基于流体动力学方程和物体的几何形状。

1.2.3伯努利定理

伯努利定理（Bernoulli’sprinciple）描述了流体中速度与压力的关系。在流体流动中，流速越快的地方，压力越小；流速越慢的地方，压力越大。伯努利方程可以表示为：

p

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

1.3流体动力学方程组

流体动力学方程组由连续性方程、动量方程和能量方程组成，它们共同描述了流体的运动状态。在实际的空气动力学问题中，这些方程通常需要通过数值方法求解，例如有限差分法、有限体积法或有限元法。

1.3.1数值模拟示例：有限差分法

下面是一个使用Python和NumPy库的简单示例，展示如何使用有限差分法求解二维不可压缩流体的连续性方程和动量方程。这个例子使用了前向差分和中心差分来近似时间导数和空间导数。

importnumpyasnp

importmatplotlib.pyplotasplt

#定义网格参数

nx=101

ny=101

nt=100

dx=2/(nx-1)

dy=2/(ny-1)

nu=0.01

sigma=0.2

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

#计算时间步长

dt=dx*dy/(2*nu*(1/dx**2+1/dy**2))

#定义差分算子

un=np.empty_like(u)

vn=np.empty_like(v)

#时间迭代

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