空气动力学方程：连续性方程与湍流模型教程.pdf
空气动力学方程：连续性方程与湍流模型教程
1空气动力学方程：连续性方程与湍流模型
1.1基础概念
1.1.1空气动力学简介
空气动力学是流体力学的一个分支，主要研究气体（尤其是空气）在物体
周围流动时所产生的力和运动效应。它在航空、汽车设计、风力发电等领域有
着广泛的应用。空气动力学的核心在于理解和预测流体的流动特性，这涉及到
流体动力学方程的求解，其中连续性方程和动量方程是基础。
1.1.2连续性方程的物理意义
连续性方程描述了流体在流动过程中质量守恒的原理。对于不可压缩流体，
连续性方程可以表示为：
∂
+∇⋅=0
∂
∇⋅
其中，是流体的密度，是流体的速度向量，是散度算子。对于不可
压缩流体，密度是常数，因此方程简化为：
∇⋅=0
这表明在任何时刻，流体通过任意闭合表面的净流入量必须等于净流出量，
从而保证了流体的质量守恒。
1.1.3湍流模型的基本概念
湍流是流体流动的一种复杂状态，其特征是流体速度的随机波动和能量的
非线性传递。湍流模型是用来简化湍流流动的数学描述，以便于数值模拟和工
程应用。常见的湍流模型包括：
雷诺应力模型（RSM）：基于雷诺平均纳维-斯托克斯方程，考虑
了湍流应力对流动的影响。
k-ε模型：基于湍动能和湍动能耗散率的方程组，是工业应用
中最常见的湍流模型之一。
k-ω模型：与k-ε模型类似，但使用了涡旋频率来代替耗散率，
在边界层和近壁面流动中表现更佳。
湍流模型的核心是通过引入额外的方程来描述湍流的统计特性，从而将复
杂的湍流流动简化为可计算的形式。
1
1.2连续性方程与湍流模型的结合
在处理湍流流动时，连续性方程与湍流模型的结合是至关重要的。例如，
在k-ε模型中，连续性方程和动量方程被雷诺平均，得到雷诺平均纳维-斯托克
斯方程（RANS），同时引入了湍流应力项。湍流模型的方程组（k方程和ε方程）
则用来计算湍流应力项。
1.2.1示例：k-ε模型的数值求解
下面是一个使用Python和SciPy库来求解k-ε模型的简化示例。假设我们
有一个简单的二维湍流流动问题，其中k和ε的分布需要通过数值方法求解。
importnumpyasnp
fromscipy.sparseimportdiags
fromscipy.sparse.linalgimportspsolve
#定义网格参数
nx,ny=100,100
dx,dy=1.0/nx,1.0/ny
x=np.linspace(0,1,nx)
y=np.linspace(0,1,ny)
X,Y=np.meshgrid(x,y)
#定义湍流模型参数
k=np.zeros((nx,ny))
epsilon=np.zeros((nx,ny))
nu=0.01#动力粘度
Cmu=0.09#湍流模型常数
sigma_k=1.0#k方程的Prandtl数
sigma_epsilon=1.3#ε方程的Prandtl数
#定义边界条件
左边界值
k[:,0]=1.0#k
左边界值
epsilon[:,0]=0.1#ε
#构建差分矩阵
A=diags([-1,2,-1],[-1,0,1],shape=(nx-2,nx-2)).toarray()/