空气动力学实验方法：流场显示技术：激光多普勒测速技术教程.docxVIP

PAGE1

PAGE1

空气动力学实验方法：流场显示技术：激光多普勒测速技术教程

1空气动力学基础

1.1流体动力学基本概念

流体动力学是研究流体（液体和气体）在静止和运动状态下的行为及其与固体边界相互作用的学科。在空气动力学中，我们主要关注气体的流动特性。流体动力学的基本概念包括：

流体的不可压缩性和可压缩性：在低速流动中，气体的密度变化可以忽略，此时气体被视为不可压缩的。但在高速流动中，如超音速飞行，气体的密度变化显著，必须考虑其可压缩性。

流体的粘性：流体流动时，流体分子之间的相互作用会产生内摩擦力，这种性质称为流体的粘性。

流体的连续性：流体可以被视为连续介质，其物理性质如密度、压力和速度在空间中连续变化。

1.2流场与流线

1.2.1流场

流场是指在空间中流体的物理性质（如速度、压力、温度和密度）随位置和时间变化的分布。流场可以分为：

均匀流场：流体的物理性质在空间中保持不变。

非均匀流场：流体的物理性质随空间位置变化。

1.2.2流线

流线是在某一时刻，流体中各点速度方向的连线。流线的性质包括：

流线不会相交。

流线的密度反映了流速的大小，流线越密集，流速越大。

1.3流体的连续性和动量方程

1.3.1连续性方程

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

?

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

1.3.2动量方程

动量方程描述了流体运动中力与加速度的关系，是牛顿第二定律在流体动力学中的应用。对于不可压缩流体，无粘性流动的动量方程（欧拉方程）可以表示为：

?

其中，u是流体速度矢量，ρ是流体密度，p是流体压力，g是作用在流体上的外力（如重力）。

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

假设我们有一个二维流场，其中流体的速度分量为ux,y,t和v

importnumpyasnp

fromscipy.sparseimportdiags

fromscipy.sparse.linalgimportspsolve

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

nx,ny=100,100

dx,dy=1.0/(nx-1),1.0/(ny-1)

dt=0.01

#初始化速度场

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

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

#定义速度场的边界条件

u[0,:]=1.0#左边界速度为1

u[-1,:]=0.0#右边界速度为0

v[:,0]=0.0#下边界速度为0

v[:,-1]=0.0#上边界速度为0

#定义连续性方程的离散化矩阵

A=diags([-1,1],[-1,1],shape=(nx*ny-1,nx*ny-1)).toarray()

A[0,:nx]=1.0/dx

A[nx-1,:nx]=-1.0/dx

A[::ny,-ny]=1.0/dy

A[-ny,-ny]=-1.0/dy

#定义连续性方程的右侧向量

b=np.zeros(nx*ny-1)

#求解连续性方程

fortinrange(100):

#更新速度场

u[1:-1,1:-1]+=dt*(v[1:-1,2:]-v[1:-1,:-2])/(2*dy)

v[1:-1,1:-1]+=dt*(u[2:,1:-1]-u[:-2,1:-1])/(2*dx)

#计算连续性方程的右侧向量

b[1:-1]=(u[2:,1:-1]-u[:-2,1:-1])/(2*dx)+(v[1:-1,2:]-v[1:-1,:-2])/(2*dy)

#求解连续性方程

p=spsolve(diags([1,-2,1],[-1,0,1],shape=(nx*ny-1,nx*ny-1)),b)

#更新速度场以满足连续性条件

u[1:-1,1:-1]-=dt*(p[1:]-p[:-1])/dx

v[1:-1,1:-1]-=dt*(p[nx:]-p[:-nx])/dy

#可视化速度场

importmatplotlib.pyplotasplt

plt.quiver(u,v)

plt.show()

在这个例子中，我们首先定义了网格大小和时间步长，然后初始化了速度场和边界条件