准一维拉法尔喷管流解析解

觅°

已于 2024-08-21 11:33:13 修改

阅读量441

点赞数

文章标签：算法 1024程序员节

于 2023-10-24 12:14:44 首次发布

本文链接：https://blog.csdn.net/qq_41817886/article/details/133984148

版权

假设：定常、绝热、无粘、理想气体

已知：总温T0，总压p0，出口静压pe，喷管面积A(x)

求：Ma分布，即Ma(x)

1. 全亚音速工况
根据

$p=\frac{p_0}{(1+\frac{\gamma-1}{2}Ma^2)^{\frac{\gamma}{\gamma-1}}},T=\frac{T_0}{1+\frac{\gamma-1}{2}Ma^2}$

可求出口Mae,Te

$Ma_e=\sqrt{(\frac{p_0}{p_e})^{\frac{\gamma-1}{\gamma}}-1}$

$T_e=\frac{T_0}{1+\frac{\gamma-1}{2}Ma_e^2}$

得

$\rho_e=\frac{p_e}{RT_e}$

获得质量流量qm

$qm=\rho_e A(x_e)V_e$

V=Ma*a，代入，将所有未知量转化为Ma，最终得

$\frac{qm}{A(x)}\frac{\sqrt{T_0}}{p_0}\frac{\sqrt R}{\sqrt \gamma}=Ma(1+\frac{\gamma-1}{2}Ma^2)^{-3}$

对于某一x，左式为已知量，因此上式为一元非线性方程f(Ma)=0，用牛顿迭代法即可求解。

2. 亚音-超音-亚音工况

此时喉部t处为音速，喉部下游存在一个通道正激波，激波下游为亚音。激波前后均为等熵流动，但激波处存在熵增。波前总温等于波后总温，T01=T02=T0，波后总压小于波前总压p02<p01=p0.

首先计算出质量流量

$qm=\frac{p0At}{\sqrt{T0}}\sqrt{\frac{\gamma}{R}(\frac{2}{\gamma +1})^{\frac{\gamma +1}{\gamma -1}}}$

式中At为喉部面积。由等熵关系式

$\frac{p_{e}}{p_{02}}=(1+\frac{\gamma-1}{2}Ma_e^2)^{-\frac{\gamma}{\gamma-1}}$

$(\frac{A_e}{A_{2t}})^2=\frac{1}{Ma_e^2}(\frac{2}{\gamma +1}(1+\frac{\gamma-1}{2}Ma_e^2))^{\frac{\gamma+1}{\gamma-1}}$

A2t为激波后喷管所对应的气动喉部，可得出

$p_{02}A_{2t}=f(Ma_e)$

代入质量流量关系式 $qm=\frac{p0At}{\sqrt{T0}}\sqrt{\frac{\gamma}{R}(\frac{2}{\gamma +1})^{\frac{\gamma +1}{\gamma -1}}}$ 中，即可求出Mae。再由

$p_e=\frac{p_{02}}{(1+\frac{\gamma-1}{2}Ma_e^2)^{\frac{\gamma}{\gamma-1}}}$

计算出激波下游的总压p02。再根据正激波前后总压比关系式计算出波前马赫数Ma1（牛顿迭代法）

$\frac{p_{02}}{p_{01}}=([\frac{\gamma-1}{\gamma+1}+\frac{2}{(\gamma+1)Ma_1^2}]^{\gamma}(\frac{2\gamma}{\gamma+1}Ma_1^2-\frac{\gamma -1}{\gamma +1}))^{-\frac{\gamma}{\gamma -1}}$

那么产生激波的位置所对应的喷管面积为

$A=\frac{qm\sqrt{T0}}{p_{01}Ma_1(1+\frac{\gamma-1}{2}Ma_1^2)^{-3}}$

根据A(x）关系式即可求出产生激波的位置xs。

当x<=xt（喉部所对应的x）时，利用

$(\frac{A(x)}{A_{t}})^2=\frac{1}{Ma^2}(\frac{2}{\gamma +1}(1+\frac{\gamma-1}{2}Ma^2))^{\frac{\gamma+1}{\gamma-1}}$

求Ma(x)，取Ma<=1.的解；

当xs>=x>xt时，利用等熵关系式

$\frac{qm}{A(x)}\frac{\sqrt{T_0}}{p_0}\frac{\sqrt R}{\sqrt \gamma}=Ma(1+\frac{\gamma-1}{2}Ma^2)^{-3}$

求Ma(x)，取Ma>1.的解；

当x>xs时，也利用等熵关系式求Ma(x)，取Ma<1.的解。

3. 亚音-音速-亚音工况

喉部处为音速，其他区域为亚音速，此时可直接利用

$(\frac{A(x)}{A_{t}})^2=\frac{1}{Ma^2}(\frac{2}{\gamma +1}(1+\frac{\gamma-1}{2}Ma^2))^{\frac{\gamma+1}{\gamma-1}}$

关系式求Ma(x)，取Ma<=1的解。

分析解析解，可得到以下性质：

1）在全亚音速工况下，质量流量qm是坐标x和出口压力pe的函数，因此出口条件改变时，质量流量会发生变化；

2）存在激波时，qm是进口总温、总压以及喉部面积的函数，因此出口压力pe变化对质量流量没有影响，但会改变激波位置和强度，而对激波上游流动参数没有干扰。

觅°

关注

0
点赞
踩
2

收藏

觉得还不错? 一键收藏
0
评论
复制链接

分享到 QQ

分享到新浪微博

扫一扫