空气动力学实验方法：热电偶测量与空气动力学基础理论.docxVIP

PAGE1

PAGE1

空气动力学实验方法：热电偶测量与空气动力学基础理论

1空气动力学基础理论

1.1流体力学基本概念

流体力学是研究流体（液体和气体）在静止和运动状态下的力学性质的学科。在空气动力学中，我们主要关注气体的行为，尤其是空气。流体的基本概念包括：

密度（ρ）:单位体积流体的质量。

压力（P）:流体垂直作用于单位面积上的力。

速度（V）:流体在某一点的运动速度。

温度（T）:流体的热状态，影响其密度和压力。

粘度（μ）:流体流动时内摩擦力的度量。

1.2空气动力学中的流体特性

空气动力学中的流体特性包括：

可压缩性:空气的密度随压力和温度变化而变化。

粘性:空气流动时，分子间的相互作用产生摩擦力。

热传导:空气中的热量传递。

湍流与层流:描述流体流动状态的两种类型。

1.2.1湍流与层流的区别

层流:流体流动平滑，各层之间互不干扰。

湍流:流体流动混乱，存在大量涡旋和混合。

1.3流体动力学方程

流体动力学方程是描述流体运动的基本方程，包括：

连续性方程:描述流体质量守恒的方程。

动量方程:描述流体动量守恒的方程。

能量方程:描述流体内能和动能守恒的方程。

1.3.1连续性方程

对于不可压缩流体，连续性方程简化为：

?

其中，ρ是流体密度，v是流体速度矢量，t是时间。

1.3.2动量方程

动量方程（纳维-斯托克斯方程）描述了流体的运动状态：

ρ

其中，μ是流体的动态粘度，f是作用在流体上的外力。

1.3.3能量方程

能量方程描述了流体的能量守恒：

ρ

其中，e是流体的单位质量能量。

1.4边界层理论

边界层理论描述了流体在固体表面附近的行为，包括：

边界层厚度:流体与固体表面接触时，速度从固体表面的零逐渐增加到自由流速度的距离。

边界层分离:当流体遇到物体表面的逆压梯度时，边界层可能分离，形成涡旋。

边界层的形成和分离对物体的阻力和升力有重要影响。

1.5示例：计算流体动力学（CFD）中的湍流模型

在计算流体动力学中，湍流模型用于模拟湍流行为。其中，最常用的模型之一是k-ε模型。

1.5.1k-ε模型

k-ε模型基于湍流动能（k）和湍流耗散率（ε）的方程。方程如下：

?

?

其中，μt是湍流粘度，σk和σε是湍流模型的常数，Pk是湍流动能的产生项，C

1.5.2代码示例

以下是一个使用Python和SciPy库来求解k-ε模型的简化示例。请注意，实际应用中需要更复杂的网格和边界条件处理。

importnumpyasnp

fromscipy.sparseimportdiags

fromscipy.sparse.linalgimportspsolve

#定义网格参数

nx=100

ny=100

dx=1.0/(nx-1)

dy=1.0/(ny-1)

#定义湍流模型参数

rho=1.225#空气密度

mu=1.7894e-5#空气动力粘度

mu_t=np.zeros((nx,ny))#湍流粘度

sigma_k=1.0

sigma_e=1.3

C1=1.44

C2=1.92

#初始化湍流动能和耗散率

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

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

#定义湍流动能和耗散率的产生项

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

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

#定义速度场

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

#构建离散方程的矩阵

A=diags([-1,2,-1],[-1,0,1],shape=(nx-2,nx-2)).toarray()/dx**2

B=diags([-1,2,-1],[-1,0,1],shape=(ny-2,ny-2)).toarray()/dy**2

#求解k-ε模型

foriinrange(1,nx-1):

forjinrange(1,ny-1):

#计算湍流粘度

mu_t[i,j]=0.07*(k[i,j]/epsilon[i,j])

#计算湍流动能和耗散率的产生项

P_k[i,j]=mu_t[i,j]*(np.gradient(v[i,j,0])**2+np.gradient(v[i,j,1])**2)

P_e[i,j]=C1*rho*epsilon[i,j]