三维磁流体动力学管道流动加减速控制数值研究

阐述了磁流体动力学(MHD)控制流场作用机理.在小磁雷诺数条件下,运用数值模拟方法,对不同外加电磁场条件下三维MHD管道流动的流场情况进行研究,得出不同磁场、电场等MHD参数作用下MHD加速器的性能.对MHD加速器的应用前景进行了展望.计算结果显示,电场取6000V/m,磁场为0.92T 情况下,流场的加速性能可达13.46%,并且可以通过减小磁场或增大电场进一步提高加速性能;在仅添加0.92T磁场条件下,速度减速可达16.35%.     

六类高超声速激波-激波干扰的数值模拟研究

基于粘性全N-S方程，利用有限体积法进行数值离散，对Ⅰ到Ⅵ六种类型的激波一激波干扰现象进行了数值模拟，其中采用了二阶Harten-TVD格式，LU隐式算法以及Baldwin-Lomax代数湍流模型。每种类型都得到了清晰的干扰流场结构，其中Ⅳ型干扰的模拟结果与实验结果进行了比较，符合较好。数值模拟结果的对比指出，由入射位置差异引起的不同类型激波-激波干扰结构有很大不同，对壁面压强分布的影响也存在较大差异。     

三维高超声速底部喷流干扰流场数值模拟与试验研究

通过数值模拟和风洞试验两种手段对来流马赫数M∞=4、喷流压比Pj/P∞=156.8、不同迎角下的三维高超声速底部喷流干扰流场进行了研究.研究结果表明超声速底部喷流干扰流场结构复杂,有、无喷流时底部流场有很大不同,对气动力系数影响显著;在大喷流压比情况下,喷流干扰使导弹纵向气动力系数下降、压心前移.最后,对数值模拟与风洞试验在结果上的差异进行了分析.     

超燃冲压发动机唇口气动热计算研究与分析

采用数值模拟方法,对超燃冲压发动机进气道唇口的气动热进行了计算和分析.依据唇口与前体激波所处相对位置的不同,分三种情形进行了计算研究.结果表明,三种情形下,进气道前缘驻点热流均大于头部.通过算例验证,证明该方法和结论的正确性.研究表明,激波干扰以及前体激波压缩后超声速来流使唇口热流增大,热防护和结构设计时应当予以充分考虑.      

大攻角/舵偏角条件下超声速导弹轴向力工程计算方法

以超声速导弹为研究对象,计算了较大攻角和舵偏角共同作用下导弹的压差阻力、粘性摩擦力和底部阻力等轴向力.计算零攻角状态的轴向力,再根据实验数据对攻角和舵偏角的影响作非线性修正.分析结果表明:该法适用的攻角与舵偏角范围较大,且计算精度满足工程要求.     

超声速底部喷流干扰流场数值模拟

数值模拟了不同马赫数,不同喷流压比下的轴对称超声速底部喷流干扰流场,采用LU隐式算法进行数值求解并引入了Baldwin-Lomax代数湍流模型.采用分区网格将弹身与底部区域合为一个整体进行计算,得到了清晰的流场结构和弹体表面及底部的压力分布,并与试验结果进行了比较.数值模拟结果表明超声速底部喷流干扰流场结构复杂,有、无喷流时底部流场有很大不同, 喷流对底压分布有明显影响,进而对轴向力系数影响显著.     

低压比煤油涡轮数值计算与优化设计

小流量煤油涡轮泵可用于膨胀循环超燃冲压发动机燃料供应系统,针对特定工况提出了超临界/裂解态煤油基低压比涡轮的数值计算方法和优化设计策略。根据液体火箭发动机中典型的涡轮设计方法获得了低压比煤油涡轮的设计方案,采用湍流模拟方法结合煤油的多组分代理模型对25kr/min转速下的涡轮内部超临界态流动进行数值计算,发现设计方案的轴功率超过所需轴功率的120%,不利于涡轮泵系统在设计点工况下的稳定运转。取涡轮轴功率大于所需轴功率为约束条件,选择涡轮结构尺寸为设计变量,以两个目标量(优化方案的轴功率和效率相对于设计方案的变化率)的加权函数值最大为目标,基于响应面模型和多岛遗传算法开展渐进优化,优化过程中采用i SIGHT平台集成了3维参数化建模和流场仿真等C++程序和软件以实现数值计算自动化。利用试验设计方法建立样本数据库,并进行了涡轮轴功率和效率关于设计变量的灵敏度分析,发现二者成合作关系;所得涡轮优化方案的两个目标量分别下降16.5%和2.9%,以较低的效率损失为代价实现了轴功率的良好配合。     

预处理方法在化学非平衡流场数值模拟中的应用

将预处理方法发展应用于化学非平衡流场的数值模拟中。采用有限差分,LU-SGS隐式方法对二维可压化学非平衡控制方程进行了耦合求解,对化学源项和无粘项隐式处理,粘性项则采用显式。通过化学反应剪切层和激波点火的算例表明,所采用的数值方法能够有效求解化学非平衡流场,证明了本算法的可行性,为下一步的工程应用奠定了基础。     

高超声速数值激波失稳的网格依赖性

高超声速强激波的稳定捕捉仍极具挑战性。目前工程计算中普遍应用的数值格式在模拟多维强激波时通常会遭遇明显的激波不稳定现象,且数值格式的激波稳定性对计算网格表现出严重的依赖性。基于矩阵稳定性分析法,对比了具有不同耗散性质的数值格式稳定捕捉激波的能力,分析了空间二阶精度格式的激波稳定性及限制器对激波稳定性的影响。在此基础上,重点探究了计算网格对激波稳定性的影响规律,研究了网格长宽比和畸变角度对激波稳定性的影响。结果显示,在激波附近通过增大网格长宽比或改变网格畸变角度可有效改善激波捕捉的稳定性;相比于增大网格长宽比,改变网格畸变角度提升激波捕捉稳定性的效果更加明显。在此基础上,结合数值耗散分析对高马赫数下数值激波失稳现象的网格依赖性机制进行了探讨。     

气动热数值模拟中的网格相关性及收敛

针对气动热数值模拟中的网格相关性和收敛问题,采用计算流体力学(CFD)计算方法,以典型的钝头为算例,进行了数值模拟研究。研究结果表明:网格是气动热数值模拟中的关键因素,壁面附近法向网格间距最为敏感,热流结果随网格间距不同会出现数倍乃至数量级的差异,在得出正确的壁面压力结果的网格上不能保证得到正确的热流结果。分析认为壁面附近客观存在的巨大温度梯度是网格要求敏感的原因。计算表明:气动热计算的收敛过程比压力和流场收敛慢得多;在压力和流场完全收敛时,热流与最终的迭代收敛值偏差达20%以上;以方程残值的下降量级和压力收敛作为热流收敛的判别标准是不合适的。在气动热计算中,应直接观察热流数据的收敛,确保得到真正的收敛解。