ADI方法求解完全跨声速非定常小扰动方程

本文修改了Whitlow的ADI格式,以便求解另一变形的完全跨声速小扰动方程。本文给出的NACA64A006和NACA0012翼型的非定常跨声速流动的计算结果与其它数值结果和实验结果的比较表明本文格式是合理的且便于实际应用的。     

多重网格法求二元跨声速流在一种半自适应网格上的欧拉解

本文采用MacCormack两步显式格式,用有限体积法求解了二元跨声速欧拉流。推导了物面边界条件,采用了特征远场边界条件及远场环量修正。利用保角变换方法生成O型贴体网格,并修改得到了一种在激波处局部加密的半自适应贴体网格。采用多重网格及焓阻尼加速收敛技术计算了NACA0012翼型的跨声速气动特性,得到了十分满意的收敛过程和计算结果。     

超声速X形鸭翼-弹身组合体涡迹发展

应用蒸汽屏方法研究超声速X形鸭翼-弹身组合体涡迹发展。观察了起源于鸭翼后缘的四个翼涡在横截面上形成的“蛙跃”和上反角二翼涡与弹身一对对称脱体涡形成的“混合式蛙跃”现象。在临近蛙跃距离时,有不稳定特性发生。文中还给出了细长体理论计算的涡迹路径跟实验数据比较,结果表明:如果各个旋涡的初始位置和相对强度适合,这种数学模型可计算导弹上的各个旋涡路径,二者存在的差别,可能是由于计算未能模拟涡面和涡量耗散的缘故。为了有助于理解导弹的气动特性,用少量的离散涡计算涡迹路径,作为工程估算是适宜的。     

风对建筑绕流流动的数值方法

本文利用k-ε湍流模型及SIMPLE方法,在二维条件下,就风对建筑绕流流动进行了数值模拟,建筑及地面的边界条件采用以低Re数k-ε湍流模型中导出的壁函数处理近壁湍流。自由边界采用外推格式。对于计算域中的建筑物则利用控制粘性系数的方法处理。并对计算结果和风洞实验进行了比较。 为了研究建筑物之间的气流及压力分布,本文还对并列两栋建筑的不同情况进行了数值计算和分析。     

化学非平衡边界层绕流计算

本文介绍用化学非平衡边界层方程组,数值计算高超声速钝头体绕流流场、组元浓度和电子密度分布情况。计算中考虑了电离空气的七种主要成分以及相互间的化学反应,比较了变熵、普朗特数及Lewis数的影响,讨论了完全催化壁和非催化壁两种情况,计算结果与文献及试验进行了比较,结果令人满意。     

超声速燃烧室性能一维数值模拟

由气体动力学基本定律，推导出燃烧室性能分析的一种一元流计算方法。该方法考虑了面积变化、质量添加、两相流、多组分效应、化学反应、壁面摩擦、壁面散热等因素的影响。与试验数据的对比分析表明，该方法能够准确描述超声速燃烧室两相参数分布，从而为燃烧室方案设计阶段的快速性能评估和设计优化提供了一种有效的手段。     

跨音速大迎角Euler方程数值分析

应用有限体积法离散三维Euler方程,选用三步显式格式进行时间推进求解;并通过当地时间步长、残值光顺、焓修正等方法加速收敛。对一维波动方程的von Neumann稳定性分析说明这种三步格式的最大Courant数是2。对三角翼跨音速大迎角流动的数值分析表明文中的三步格式及焓修正方法优于人们经常采用的四步、五步格式以及焓阻尼技术。     

飞机电网零序阻抗的数值计算

从麦克斯韦方程出发,推导了飞机电网零序阻抗参数的计算公式。对几种不同敷设方式载流导线的零序阻抗进行了计算,其结果与标准值较为接近。这种数值计算方法比传统的解析近似估算方法更为合理和准确。     

高阶面元法计算多段翼型之位流

本文给出了一种求解多段翼型位流的高阶面元法。面元和奇点的分布分别采用了二阶和三阶的形式。由于高阶面元法较为精确地反映了绕流体的几何形状和奇点表面分布,本文给出了比低阶面元法更为精确的解,并且消除了低阶面元法计算中经常出现的翼型后缘处解的波动。     

室火轰燃的大涡模拟

采用一套修正的适用于低马赫数浮力流的N-S方程组来描述火灾场,运用Smagorinsky亚格子模型的大涡模拟方法对室内池火的轰燃发展过程进行了数值模拟,把模拟与实验结果进行了对比和分析。结果表明,该方法能较好地预测轰燃。同时模拟结果显示火焰区向房间后部偏转稍大,近壁处温度预测的误差大一些,随着亚格子模型的不断完善和计算能力的提高,用该技术来预测轰燃现象会更加准确。