出现一个只有通过超算才会知道的误区:
不导入压缩性系数的话,整个公式将会完全报废。
因此在钱五师询问意见后,徐云立刻提出了自己的看法——如果钱五师不问,徐云就会主动开口。
而在徐云身边。
钱五师闻言也点了点头:
“正合我意。”
于是很快。
钱五师便计算起了背压比。
所谓背压比。
指的喷嘴出口静压力与喷嘴上游滞止压力之比,不过在设计方案中指的是锥流场与气体的耦合比。
当锥流场刚好达到临界条件时。
外部气体达到音速,同时气体质量流量达到最大值,此时的背压比即称为最大背压比。
这个概念有点类似后世的bpr,不过释义上更接近下游。
接着很快。
徐云也估量了一番自己的右手状态。
今天他的右手还没用过,负载为0,因此他便也拿起笔和纸协助写了起来。
众所周知。
如果激波为正激波,且不考虑激波厚度,那么激波控制体的形状就会很对称:
你比划个剪刀的手势,然后指尖向下。
这就是激波控制体的图示了。
而控制体cv基本方程,则由三个连续方程组成:
dΦdt=ddt∫v?(r,t)dv=??t∫v?(r,t)dv+∮s?(r,t)u?nda
Δn=(?iiσpdt+?iiiσpdt))t+Δt?(?iiσppdt)t
liΔt→0(?iσpdt)t+ΔtΔt=??σ?v→?da→=?σpvs?αda(起点这排版将就着看吧)
其中t为时间;
fx为控制体内流体的受力在x轴上的分量;
v为流体速度失量;
a为控制体表面面积失量;
v为控制体体积。
同时考虑气体稳定流动,再假设速度、能量在激波截面上是均匀的。
便有∫csv·da=ca。
随后徐云把截面态联立在了一起,准备继续推导下去。
然而半分钟后。
徐云忽然眉头一皱,嘴里啧了一声,轻轻摇了头:
“不行,要是这样拟合的话,就没法继续计算了.....”
结果话音刚落。
徐云的耳边忽然传来了一道声音:
“韩立同志,为什么没法继续计算?”
“?”
徐云顿时一怔,顺势朝发声者看去。
转过头后。
发现数算小组的那位被叫做什么“大于”的圆脸中年人,不知何时已经来到了自己身边。
徐云见状扫了眼正在低头计算的钱五师,压低声音解释道:
“大于同志,这不是很明显吗?”
“激波后的温度高于激波产生前,压力间断性地急剧上升,扩散段的方程显然是算不出来的。”
说罢。
徐云便摇了摇头,准备试着思考另一种方法。
然而令他有些意外的是。
圆脸中年人闻言后没有再说话,而是同样低头拿着笔和纸写了起来。
徐云见状也不再说什么,继续做起了思考。
过了大概三四分钟。
中年人忽然将算纸递到了徐云面前,说道:
“韩立同志,你看看这个。”
徐云这会儿还处在思路断档期,被人反复打搅,心中多少还是有些想法的。
反感谈不上。
但不耐烦肯定有点儿。
毕竟这可是后世的2023年都已经形成定式的准公理,在徐云看来没太多讨论的必要。
不过出于对这个时代先辈的敬重,徐云还是决定先帮忙这位同志找出问题,给他简单的上上一课。
结果在看到算纸内容的第一时间。
徐云便顿时童孔一缩。
只见此时此刻。
算纸上赫然写着一段推导:
【已知pd/pu=(k+1)a2u/2+(k?1)a2u】
【以及y=pd/pu√[2kx2?(k?1)k+1]^1/2】
【对以上二方程进行联立,建立二维柱坐标下的可压缩粘性气体的连续性方程、n-s方程、能量方程和气体状态方程】
【通过变式可知,截面态会在扩散段后半段中逐渐增大,引入气体边界层影响后可得最终式......】
【∑f→cv=??t?cvv→pdb+?→)da......】
“??????”