空气动力学仿真技术：直接数值模拟(DNS)：高精度数值方法在DNS中的应用.pdf

空气动力学仿真技术：直接数值模拟(DNS)：高精度数值方

法在DNS中的应用

1空气动力学基础理论

1.1流体力学基本方程

流体力学是研究流体（液体和气体）的运动和静止状态的学科，其基本方

程是纳维-斯托克斯方程（Navier-Stokesequations），描述了流体的动量守恒、

质量守恒和能量守恒。在空气动力学中，我们主要关注的是不可压缩流体的纳

维-斯托克斯方程，其形式如下：

1.1.1质量守恒方程（连续性方程）

∂

+∇⋅=0

∂

其中，是流体的密度，是流体的速度向量，是时间。

1.1.2动量守恒方程

∂

⋅⊗−⋅

+∇=∇+∇+

∂

其中，流体的压力，应力张量，是重力加速度向量。

1.1.3能量守恒方程

∂

+∇⋅=∇⋅⋅+⋅

∂

其中，是流体的总能量，包括内能和动能。

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

importnumpyasnp

importmatplotlib.pyplotasplt

#定义网格参数

nx=101

ny=101

nt=100

nit=50

1

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

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

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

dt/(2*rho*