空气动力学方程：层流和湍流模型：直接数值模拟(DNS)技术教程.pdfVIP

空气动力学方程：层流和湍流模型：直接数值模拟(DNS)技

术教程

1空气动力学基础

1.1流体动力学基本概念

流体动力学是研究流体（液体和气体）在静止和运动状态下的行为的学科。

在空气动力学中，我们主要关注气体，尤其是空气。流体动力学的基本概念包

括：

流体的连续性：流体可以被视为连续介质，没有离散的颗粒。

流体的不可压缩性：在低速流动中，空气的密度可以认为是常数，

这种假设简化了方程。

流体的粘性：流体内部存在摩擦力，影响流体的流动特性。

边界层：流体与固体表面接触时，由于粘性作用，流体速度从固

体表面的零逐渐增加到自由流速度，形成边界层。

层流与湍流：层流是流体平滑流动的状态，而湍流则是流体流动

中存在大量随机涡旋的状态。

1.2连续性方程解析

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

简化为：

∂∂∂

++=0

∂∂∂

其中，、、分别是流体在、、方向的速度分量。

1.2.1示例代码

importnumpyasnp

#定义网格尺寸

nx,ny,nz=100,100,100

dx,dy,dz=1.0,1.0,1.0

#初始化速度场

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

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

w=np.zeros((nx,ny,nz))

1

#假设在某个时间步，速度场为：

u[50,50,50]=1.0

v[50,50,50]=2.0

w[50,50,50]=3.0

#计算连续性方程的左侧

continuity=(np.gradient(u,dx)[0]+np.gradient(v,dy)[1]+np.gradient(w,dz)[2])

#输出连续性方程的值

print(连续性方程的值：,continuity[50,50,50])

1.3动量方程和能量方程

动量方程描述了流体运动中力与加速度的关系，是牛顿第二定律在流体动

力学中的应用。对于不可压缩流体，动量方程可以表示为：

222

∂∂∂∂1∂∂∂∂

+++=−+++

∂∂∂∂∂∂∂

能量方程描述了流体中能量的守恒，包括动能、位能和内能的转换。

1.3.1示例代码

importnumpyasnp

#定义网格尺寸和时间步

nx,ny,nz=100,100,100

dx,dy,dz=1.0,1.0,1.0

dt=0.01

#初始化速度场和压力场

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

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

w=np.zeros((nx,ny,nz))

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

#假设在某个时间步，速度场和压力场为：

u[50,50,50]=1.0

v[50,50,50]=2.0

w[50,50,50]=3.0

p[50,50,50]=100.0

#定义流体密度和粘性系数

rho=1.225#空气密度，单位：kg/m^3

nu=1.5e-5#空气动力粘度，单位：m^2/s

2

#计算动量方程的左侧

momentum_left=(np.