一种利用正交流线坐标系求解跨音速管流的新方法

本文以正交流线坐标系为工具,采用位流模型,导出了正交流线坐标系下平面或轴对称流动的控制方程组,利用流函数与位函数的存在性,进而导出了两个首次积分。设计了一套求解匀熵管流的分凝迭代方案,模拟了管内无激波的平面或轴对称喷管的跨音速流动。数值实验表明该方法能够处理小喉道曲率半径平面或轴对称喷管跨音速流动,与已有的数值计算结果及实验结果符合较好。     

轴对称跨音速流中无激波超音速区的设计

本文直接从音速线出发计算了一个旋成体跨音速绕流超音速泡中的特征线网。提出了从音速线开拓到一条超音速线的方法,计算了无旋轴对称流的特征线网,发现其特点是不能跨越对称轴线,并对其原因作了初步的讨论。本文的计算结果可应用于物面修形。     

数值模拟高超音速化学非平衡流的隐式格式

针对高超音速化学非平衡流提出了两种隐式格式.一种是近似因子分解算法对角化形式在隐式TVD(全变差递减)格式上的推广,另一种可视为隐式TVD格式的LU分解形式.数值试验表明,本文给出的算法保持了TVD格式的高分辨率特性,同时计算量大大减少,是高超音速化学非平衡流数值模拟的有效方法.     

火焰筒头部流场研究之二——二元轴对称等温廻流区计算

本文将火焰筒头部的等温流场,用二元轴对称不可压定常流来处理并用数值方法解出初步结果。基本方法是用涡量流函数法,双方程粘度模型,采用交替方向隐含格式差分求解,在边界及迭代初值提出自己的看法,初步计算结果与[10]之实验定性一致。     

用复位表示的一般二维热弹性问题

本文采用位移函数统一地导出既有体力,也有温差的平面热弹性问题的复位表示式,既使在无体力,无温差的特殊情形下;其公式形式也与Kolosov-Muskhelishvili公式的形式不同,比Kolosov-Muskhelishvili 公式的形式要整齐,推导方法更简明。对偶的也用统一的应力函数法和推广Stevenson法(既不用位移函数,也不用应力函数)分别导出有体力,有温差的Kolosov-Muskhelishvili形式的公式。最后,得出结论,用统一的位移函数是最好的一种方法:推导方法简明,公式形式整齐。     

进气旋流对推力喷管性影响的研究

介绍了模拟近代大涵道比发动机排气状况的旋流模型,和在喷管落压比1.2到6.5范围内,进气旋流对轴对称喷管,矩形收-扩喷管,楔形二元喷管,和单边膨胀斜板形喷管等性能和射流流态影响等的一些实验结果。实验数据指出:进气旋流可使喷管出口射流柱中的旋涡加强,射流与外界气流间的混合区扩大,中心核心流区缩小,但它对轴对称喷管出口流态的影响很小。进气旋流同时对喷管推力系数和流量系数均起不利影响,增大进气旋流角度会使喷管推力系数和流量都下降,特别是当旋流角超过某一“临界”值后二者降低很急剧,并且二元喷管降低的程度比轴对称喷管更严重,因之在具体应用情况中,保留多大发动机的旋流,要针对具体情况综合考虑各种因素来选定。     

以流面坐标系任意特征线计算超音速绕圆—椭圆三角翼的流动

本文以流面坐标系沿任一条特征线的物理量变化关系式对音速前缘圆—椭圆三角翼的绕流作了计算。其在流面坐标系沿任一条特征线的物理量变化关系式是根据文献[1]所推导得的公式。计算得的机翼表面压强分布与文献[3]、[4]及文献[5]、[6]作了比较;并与文献[7]、[8]的实验结果也进行了比较。结果表明,本文Δζ~3=0.19635的计算结果与文献[3]、[4]、[5]、[6]的计算结果符合较好,而Δζ~3=0.087266的有向实验靠近的趋势。     

俯仰阻尼导数分量的CFD数值模拟

提出一种直接求解直接阻尼导数的方法,该方法不仅适用于轴对称外形,也适用于非轴对称外形.数值模拟飞行器的非定常强迫沉浮运动和强迫角振动,并在Etkin的非定常气动力模型基础上,辨识得到飞行器的洗流时差导数和俯仰阻尼导数,研究飞行器俯仰通道各阻尼导数的数值计算方法.对弹道外形和基本带翼导弹标模外形及Hyflex升力体外形进行研究分析表明,将阻尼导数分量相加得到的俯仰阻尼导数与直接求解强迫角振动得到的俯仰阻尼导数与试验结果吻合很好,各阻尼导数分量随质心位置的变化趋势也与理论预测相符.对于带翼飞行器,超声速条件下,洗流时差导数在俯仰阻尼导数中占主导作用.      

高超音速化学非平衡流的隐式格式迭代解法

针对高超音速化学非平衡流计算中的隐式TVD格式给出了一种迭代解法.该方法是在原有的Harten-Yee格式基础上,在每一时间步对无粘部分进行迭代运算,目的是为了加速收敛.借助对模型方程的收敛性分析,证明该迭代算法是收敛的.二维数值实验表明该方法对于高马赫数时的化学非平衡流能明显提高隐式格式的计算效率,其计算结果是满意的.     

高超音速钝头体绕流反应速率与壁面催化问题

反应速率与壁面催化是反应流中两个重要问题。本文以高超音速化学反应钝头体绕流为例对这两个问题做了初步的分析和讨论。反应速率问题实质上是反应速率常数实验数据的误差及其影响,本文根据实验数据的分散度大小以样本数据组和等效速率常数不确定性因子这两个概念建立反应速率常数分析数据库,通过对数据误差和计算结果的分析与比较总结出高超音速钝体绕流高温特性与化学反应特性计算结果依赖于反应速率数据库的某些规律。壁面催化问题是固壁对与之接触的气体化学反应的催化影响,通过引入一个新的壁面催化率参数来考虑壁面催化性质对高超音速钝体绕流计算结果的影响。由于引入和强化了上述概念,给化学反应速率与壁面催化问题的分析与研究提供了有力的手段,因而提高了数值模拟结果的可信度。     

超音速轴对称有旋流特征线法的计算程序

编写了一个ＦＯＲＴＲＡＮ７７程序，计算了Ｍａ０为２．０￣５．０的超速轴对称有旋流Ｃａｕｃｈｙ问题的一些特征线网，但推进方向是从初值线在物面，这样，可以对物形进行反设计，在设计科波机时，可以选用这种流场中的流面作为其升力面的组成部分。     

双三角翼外翼前缘钝化对涡流特性的影响

通过风洞试验，利用７孔探头对７５°／４５°和７５°／６０°两种典型双三角翼进行了空间流场测量，研究了双三角翼前缘形状对大迎角涡流特性的影响．试验结果表明，双三角翼外翼前缘钝化使双涡态的双三角翼内、外翼涡互相靠拢、增加干扰，当出现合并涡态时，使合并涡涡核向内下方偏移；外翼前缘钝化使外翼涡或合并涡的Ｃｐ０、Ｖｘ较尖前缘时为高，最终使外翼涡或合并涡推迟破裂．     

卫星瞬变热流红外模拟技术研究

文章对卫星瞬变热流红外灯模拟的基本要求、系统组成及各个分系统的研究作了扼要地叙述,并给出了5路不同瞬变热流红外灯模拟试验结果。     

引射喷管的气动性能和传热计算

对某型飞机引射喷管装三种机尾罩时的推力性能和机尾罩温度进行了计算，主流和次流的特性分别用特征线法和一维方法计算，为了计算机尾罩的温度，建立了一种简化的一维模型，计算结果与试飞实测数据吻合良好，引射喷管装三种机尾罩的比较表明：切短的机尾罩可以在推力下降不太多的同时改善机尾罩温度。     

引射混合器内流场的数值分析

计算了以波瓣为主喷管的三维引射掺混流场；分析了包括复杂型面和引射二次流在内的边界条件处理；再现了流向涡的产生，计算结果与实验结果作了对比分析，吻合程度较好。     

毛细管束中层流粘性振荡流强化传热的若干问题

本文是作者从事振荡管流换热研究的总结。研究结果表明:流体在管内作轴向振荡运动时,可以极大地强化轴向的传热。本文给出了计及管壁效应的振荡管流换热的精确解析解,公布了变Pr数的实验结果。在此基础上,讨论了影响振荡管流轴向换热的6个无因次参数的物理含义,并指出其中的某些无因次参数在物理本质上就是傅里叶数的倒数的平方根,这将有利于理解振荡管流换热这种非定常换热过程。由本文的精确解析解,经过简化,可以得到国内外已公布的近似解。实验数据和解析解预言的基本规律符合很好。     

失重飞机搭载气 /液两相流实验研究

利用俄罗斯IL-76MDK失重飞机对失重条件下方形截面(12mm×12mm)管道内气/液两相流动进行了实验研究,观测到失重条件下气/液两相泡状流、弹状流和环状流等三种主要流型,同时得到了气、液两相介质的流量、温度和实验段内压力等结果,分析了失重条件下方形截面管道内气/液两相流型转换的条件.      

非定常集中涡破裂形态及其对扰动的影响

本文中利用三角翼前缘涡襟翼的强制振动来推迟涡破裂。流场显示实验结果表明,在有强制扰动时,集中涡显示出两种不同形式的破裂,可分成六种不同的破裂形态。非定常集中涡的一个破裂过程涉及几种破裂形态之间的转换。涡襟翼振动产生的非定常效应对涡破裂的影响与三角翼后掠角有密切的关系,随后掠角增大非定常效应的影响变小。就实验迎角范围而言(α<35°),对50°后掠角三角翼,涡襟翼振动可显著推迟涡破裂;但对70°后掠角三角翼,振动却能促进涡的破裂。涡襟翼振动改变了集中涡的性态,使集中涡趋向于和前缘平行。振动对涡破裂的作用部分地与这种效应有关。     

对连续流式电泳简化数学模型的研究

用有限分析法对连续流式电泳的电流体力学过程进行了三维数值模拟。数值结果表明,Ｐｏｌｅｚｈａｅｖ（［１,２］）等人提出的电泳简化数学模型在处理重力条件下有回流问题时,所得结果与三维准确方程的计算结果严重不符。因此该模型不能用于重力有回流问题。为了克服上述缺陷,在解析分析的基础上,文章提出了一个新的电泳问题简化数学模型。数值计算显示,该模型所得结果与三维准确方程的计算结果十分吻合。