流过泰氟隆烧蚀表面的化学非平衡粘性激波层数值研究

本文对有泰氟隆烧蚀的化学非平衡粘性激波层流场进行了数值模拟。应用空间推进与总体迭代相结合的求解方法，采用１９组元、２９种反应的空气—泰氟隆化学反应系统，对化学非平衡粘性激波层流场与泰氟隆壁面烧蚀传热过程进行耦合求解，研究了平衡催化壁和非催化壁条件对粘性激波层烧蚀流场的影响。     

三维化学非平衡粘性激波层流场数值模拟

本文通过引进主流贴体坐标 ,将粘性激波层概念推广到三维流场 ,得到了三维粘性激波层方程 ;应用空间推进与总体迭代相结合的求解方法 ,对三维化学非平衡粘性激波层流场进行了数值模拟     

飞船轴对称热化学非平衡流场数值求解

本文从轴对称热化学非平衡Ｎ－Ｓ方程出发,利用时间相关法,采用双温度、七组元反应气体模型,利用隐式ＮＮＤ有限差分格式和时间预处理技术数值求解了ＦＩＲＥＩ飞船热化学非平衡流场,得到了较为准确的结果,并分析了热力非平衡对流场的影响。     

AUSMPW+格式在高超声速热化学非平衡流数值模拟中的应用

将AUSMPW+格式应用到高超声速热化学非平衡流场的数值模拟中。为提高精度,采用了三阶MUSCL插值方法。与LU SGS方法结合,提高了单步计算效率和收敛性。采用热化学非平衡十一组元气体模型求解了非定常轴对称Navier Stokes方程组,得到了收敛结果。数值模拟结果与文献结果进行了对比,并在弹道靶中进行了钢质圆球的实验验证。计算结果与文献、实验的对比说明,AUSMPW+格式可以在热化学非平衡流的数值模拟中精确地捕捉到强弓形激波,得到准确的空气动力系数。     

小钝头锥体全目标流场与光电特性一体化数值模拟

从简化的或全ＮＳ方程出发,采用化学非平衡和热辐射非平衡模型,数值模拟了有烧蚀产物引射条件下小钝头锥体全目标流场与光电特性,旨在为防护设计提供气动物理环境数据。计算结果表明,本文流场模拟数据与实验数据基本吻合,壁面烧蚀产物引射和壁面催化特性对流场电子数密度均有量级影响。     

飞船再入舱三维化学非平衡流数值模拟

对飞船再入舱高超声速化学非平衡三维流场进行数值模拟。控制方程为含化学反应源项的全Ｎａｖｉｅｒ－ Ｓｔｏｋｅｓ 方程组, 化学模型为高温空气的七组元模型, 其组元成分为Ｎ２ 、Ｏ２、ＮＯ、Ｎ、Ｏ、ＮＯ＋ 、ｅ－ 。差分格式采用张涵信院士提出的ＮＮＤ激波捕捉格式并用全耦合方法处理流动方程及化学反应方程, 对化学反应源项采用全隐式处理, 运用时间预处理技术以加快收敛。同时, 也对完全气体模型进行计算, 以分析真实气体效应。与实验数据对比, 说明本文结果是可信的。     

化学非平衡尾迹流场的数值计算

高超音速化学非平衡尾迹流场的分析及计算是再入空气动力学的一个重要课题。本文采用轴线极限方程,成功地消除了控制方程在轴线上的奇性;采用二阶精度的中心差隐格式,求解了抛物化的Ｎ－Ｓ方程,得到了２０ｋｍ及６０ｋｍ再入时尾迹流场的数值解。计算结果表明,在远尾迹场,对电子数密度的主要影响是Ｏ_２的吸附反应。     

一种新型的计算化学非平衡流动的解耦方法

给出一种计算化学非平衡流动的解耦方法。利用该方法从薄层近似N S方程出发 ,采用NND2M差分格式 ,数值模拟了H2 /O2 燃烧的超音速绕流场 ;计算表明这种方法可以有效地解决刚性问题 ,而且与现有的隐式算法相比 ,具有公式推导简便 ,占用内存少 ,计算效率高等优点 ,特别适合于多组分和多反应的情况     

高超声速乘波飞行器三维流场的并行数值模拟

针对高超声速乘波飞行器三维绕流流场,在基于LINUX+MPI系统的分布式并行计算平台上,并行求解了三维雷诺平均的N-S方程。并行数值方法采用的是有限体积方法(FVM)、OC-TVD差分格式、B-L代数湍流模型及流场分区的并行方法。计算结果表明,所采用的并行数值模拟方法能够求解包含强激波的流场,激波穿越区域边界时无断层、错位等通量不守恒的现象。并行计算效率高,8个处理机计算时的并行加速比达到了6 8。     

加权紧致非线性格式在热化学非平衡流数值模拟中的应用

高精度格式因其在较低的计算资源消耗下仍能保持足够的精度而在精细流动结构的模拟中具有独特优势。将显式5阶加权紧致非线性格式(WCNS-E-5)引入二维高温非平衡流计算,进行自由流速度5.0~5.7 km/s的高焓风洞和高空条件下的圆柱非平衡流场数值模拟,得到正确的流场参数分布、壁面压力和壁面热流等结果与实验值吻合较好。计算结果表明,WCNS-E-5格式较2阶的MUSCL格式具有更好的热流网格收敛特性。分析并针对算例条件初步解决了高阶格式在刚性化学反应流模拟中的鲁棒性问题,为进一步采用高阶格式开展流场结构更加复杂的热化学非平衡流数值研究奠定了良好基础。     

三维高超音速干扰流场数值模拟与实验验证

采用Osher Chakraverthy的TVD格式、Baldwin Lomax湍流模型和LU SSOR隐式方法求解了完全NS方程 ,数值模拟了三维高超音速绕流与横向喷流干扰流场。并在高超音速炮风洞中开展了喷流实验研究 ,对该喷流流场的数值计算结果进行实验验证     

预处理方法在低速粘性流动中的应用

用时间相关法、LU SGS隐式求解经过预处理后的二维可压Navier Stokes方程 ,数值模拟了低速情况下不同马赫数、不同雷诺数的驱动空腔流动。为了使时间相关法能应用于计算低速流动及加快收敛 ,在控制方程的时间导数项前乘了一个预处理矩阵 ,实现了对方程的预处理。数值模拟结果表明 ,预处理确实加快了计算的收敛性 ,而定常解却没有改变。所得结果与可能的计算结果进行了比较     

全隐LU-SGS算法在高超声速热化学非平衡流刚性问题中的应用

在高超声速热化学非平衡流动计算中,当地气体能量松弛时间、化学反应特征时间与流动时间推进步长量级差异过大会带来严重数值刚性问题,且在高雷诺数条件下,壁面、拐角等强干扰区网格加密使得该问题加剧,导致初始最大CFL数极小,收敛速度缓慢.原始LU-SGS算法仅考虑化学反应源项和对流项的隐式处理,通过推导黏性项Jacobian矩...     

二元混压式高超声速进气道流场数值模拟

使用有限差分法,TVD二阶迎风格式,Baldwin-Lomax湍流模型,对二元混压式高超声速进气道设计状态和非设计状态的粘性流场和性能进行了数值模拟,并分析了出口反压及隔离段长度对进气道流场的影响。结果表明隔离段长高比对进气道出口参数影响较大;隔离段出口反压不能太高,否则会引起进气道不起动。     

再入飞行器高超声速三维粘性流场的数值模拟

用数值方法研究大钝头倒锥体再入飞行器高超车速完全气体有攻角绕流的流场特性。控制方程为完全气体三维Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,数值方法采用激波捕捉的ＮＮＤ隐式算法,计算区域为含底部的全流场。文中给出了流场特性和气功力系数,并与实验结果进行了对比。     

涵道风扇外形参数对气动特性的影响

以某涵道风扇为原型,从理论上分析涵道扩张角对涵道风扇气动特性的影响。运用滑移网格模型,采用三维不可压黏性Navier-Stokes控制方程,利用SST k-ω湍流模型,计算两叶桨气动特性,并与试验结果对比,验证该方法的可行性。分别计算涵道风扇在悬停状态下,3000~8500 r/min转速范围内,涵道唇口外形、扩张角和涵道高度对气动特性的影响,并对流场进行分析。椭圆形唇口的涵道风扇总拉力系数小,气动效率低;当涵道扩张角在8.2°附近时,功率系数相对最小,随着扩张角增大,在桨盘下方靠近涵道壁面附近出现气流分离;涵道拉力系数对涵道风扇高度的变化敏感度低,随着高度增加功率系数略有下降。     

高超音速粘性流场数值计算中差分格式的对比分析

本文通过三维高超音速粘性流场的求解对三种常用的差分格式： Ｈａｒｔｅｎ 的 Ｔ Ｖ Ｄ 格式、修正后的 Ｏｓｈｅｒ Ｃｈａｋｒａｖａｒｔｈｙ 的 Ｔ Ｖ Ｄ 格式以及 Ｊａｎｅｓｏｎ 的当地极值递减 （ Ｌ Ｅ Ｄ） 格式对激波的分辨率、边界层计算的准确性和格式的计算效率等三个方面进行了较细致的对比分析， 为三维高超音速复杂流场的数值模拟打下基础。     

高空滑移条件对化学非平衡粘性流动的影响

本文用数值方法计算了不同高度滑移边界条件对钝体绕流流场特性的影响。采用化学非平衡粘性激波层控制方程,边界条件含壁面滑移、激波滑移和无滑移条件。计算结果表明,考虑滑移条件后,激波层较厚,激波处的电子数密度为一有限值。文中给出了计算结果,并与其它文献的结果作了比较。     

涡流阀流量调节特性数值模拟

用数值模拟研究了涡流阀的流量调节特性。在来流总压恒定的条件下,用时间相关法求解了雷诺平均的Ｎａｖｉｅｒ—Ｓｔｏｋｅｓ方程。计算了无控制流和有控制流两种情况,得到了流动参数在流场中的分布和控制流量与喷管流量的关系。计算表明,涡流阀有良好的流量调节特性。     

瞬变流二维理论模型及其应用

提出了求解弯管、孔板、三通等管件中的瞬变流动的二维理论模型和数值计算方法。对瞬变流动的二维控制方程进行坐标变换,得到计算坐标平面上的控制方程,并给出了计算平面和物理平面网格点一一对应关系方程和数值求解方法。建立了计算平面控制方程的边界条件,给出了差分求解方法,并提供了计算实例。