储气箱喷管口处速度的计算

bhZ567

于 2024-10-07 11:45:49 发布

阅读量681

点赞数 4

分类专栏： Python 文章标签： python 气体动力学

本文链接：https://blog.csdn.net/2301_76273740/article/details/142737580

版权

Python 专栏收录该内容

3 篇文章 0 订阅

订阅专栏

储气箱内部的空气压力（总压）为 $p_1^*$ ，温度（总温）为 $T_1^*$ ，设总压恢复系数为 $\sigma$ ，出口处为大气压 $p_a$ ，则 $p_2^*=\sigma p_1^*$ 。则出口处 $Ma_2=\sqrt{\frac{2}{k-1}\left( \left( \frac{p_{2}^{*}}{p_a} \right) ^{\frac{k-1}{k}}-1 \right)}$ ，静温 $T_2=\frac{T_{1}^{*}}{1+\frac{k-1}{2}Ma_{2}^{2}}$ ，速度 $v_2=Ma_2c_2=Ma_2\sqrt{kRT_2}$
代码如下

import math
#对于空气
k=1.4
R=287.06

p1_=2.943e5
T1_=288
sigma=1
p2_=p1_*sigma
pa=9.81e4
def p2Ma(p_,p,k):
    return math.sqrt(2/(k-1)*((p_/p)**((k-1)/k)-1))
def Ma2T(T_,k,Ma):
    return T_/(1+(k-1)*Ma**2/2)

Ma2=p2Ma(p2_,pa,k)
T2=Ma2T(T1_,k,Ma2)
v2=Ma2*math.sqrt(k*R*T2)
print(v2)