空气动力学基本概念：涡流：涡流与流体动力学方程.docxVIP

PAGE1

PAGE1

空气动力学基本概念：涡流：涡流与流体动力学方程

1空气动力学基础

1.1流体的性质

1.1.1连续性方程

连续性方程描述了流体在流动过程中质量守恒的原理。在空气动力学中，当流体通过一个变化的截面时，流体的质量流量在任何截面都保持不变。这可以数学地表示为：

?

其中，ρ是流体的密度，v是流体的速度向量，??是散度算子，t

1.1.2伯努利方程

伯努利方程是流体动力学中一个重要的能量守恒方程，它描述了在理想流体（无粘性、不可压缩）中，流体的动能、位能和压力能之间的转换关系。伯努利方程可以表示为：

1

其中，ρ是流体的密度，v是流体的速度，g是重力加速度，h是流体的高度，p是流体的压力。这个方程表明，在流体流动过程中，流体的总能量（动能、位能和压力能之和）保持不变。

1.2流体动力学基本方程

1.2.1纳维-斯托克斯方程

纳维-斯托克斯方程是描述粘性流体运动的偏微分方程组，它基于牛顿第二定律，考虑了流体的惯性力、压力梯度力和粘性力。在不可压缩流体中，纳维-斯托克斯方程可以表示为：

?

其中，u,v,w分别是流体在x,y,

1.2.2欧拉方程

欧拉方程是纳维-斯托克斯方程在忽略粘性力时的简化形式，适用于理想流体（无粘性）。在不可压缩流体中，欧拉方程可以表示为：

?

这些方程仅考虑了流体的惯性力和压力梯度力，忽略了粘性力的影响。在空气动力学中，欧拉方程常用于分析高速流动和涡流的形成，因为粘性效应在这些情况下通常可以忽略。

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

下面是一个使用Python和NumPy库求解二维不可压缩流体的纳维-斯托克斯方程的示例。我们将使用一个简单的离散化方法来近似方程。

importnumpyasnp

importmatplotlib.pyplotasplt

#定义网格参数

nx=101

ny=101

nt=100

nit=50

dx=2/(nx-1)

dy=2/(ny-1)

sigma=.1

nu=.1

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

#定义流体的入口速度

u[:,int(.5/dy)]=2

#离散化纳维-斯托克斯方程

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/dx**2+dt/dy**2)*(vn[1:-1,2:]-2*vn[1:-1,1:-1]+vn[1:-1,0:-2]+vn[2:,1:-1]-2*vn[1:-1,1:-1]+vn[0:-2,1:-1])

u[0,:]=0

u[-1,:]=0

v[:,0]=0

v[:,-1]=0

u[:,int(.5/dy)]=2

#应用压力修正

forninran