超临界压力下碳氢燃料裂解与流动传热模拟的快速算法

为了快速预测发动机冷却通道内碳氢燃料在考虑热裂解时的流动传热特性,基于正癸烷热裂解反应机理,建立了一套模拟正癸烷裂解吸热和超临界压力传热现象的快速算法。采用三维物性库查表算法计算裂解反应混合物的热物性,同时对组分输运方程进行简化,简化后只需求解1个组分输运方程。通过与现有的实验和数值结果进行比较,检验了快速算法的计算效率和可靠性。结果表明,快速算法与求解全组分输运方程的算法精度相当,但计算效率提升了约20倍。最后采用该算法对三维矩形截面冷却通道内的超临界压力正癸烷裂解与传热过程进行数值模拟,进一步考察了本文快速算法的计算精度及其工程应用价值。     

剪切稀化非牛顿射流撞击液膜破碎直接数值模拟

液态射流撞击是液体火箭推进系统中广泛采用的一种燃料雾化方法,其破碎特征直接影响燃料最终的掺混及燃烧效率。采用直接数值模拟（DNS）工具,研究了低雷诺数（Re_l=41）和中等韦伯数（We_l=163）条件下剪切稀化非牛顿射流撞击液膜破碎的问题,着重分析了对角液膜的三维结构、破碎特征和非牛顿特性等。研究结果表明：在所研究的射流参数下,该非牛顿撞击液膜破碎属于Open Rim类别,破碎过程具有三维特性并伴随液丝与边缘的融合、液丝向液滴的转变等时域流动特征。液体的总表面积随时间不断增长,但单位表面积随液膜破碎的发生而下降,液膜扩张半角随时间逐渐增加并趋于恒定值43°,而后部液膜的长度几乎不随时间发生变化。此外,撞击液膜表现出明显的剪切稀化特性,液体内部最低黏性系数仅为零剪切黏性系数的1/5。     

高超声速飞行器宽速域翼型/机翼设计与分析

对于现代高超声速飞行器的设计而言,除了需要保证高超声速的性能外,还必须兼顾满足工程需求的亚跨超声速特性。首先,采用雷诺平均Navier-Stokes（RANS）方程流动求解器,结合基于Kriging模型的代理优化算法,开展了高超声速飞行器宽速域翼型的优化设计研究,设计出了一种下表面具有双"S"形特征的新翼型。综合性能评估结果表明,该翼型相比于常规的高超声速翼型,在跨声速和高超声速下具有更加优良的气动特性;其跨声速状态下的升阻比达到78.9,高超声速状态下的升阻比达到5.94,能够实现宽速域内良好的综合气动性能。其次,开展了仿德国"桑格尔号"（SANGER）空天飞机运载机机翼的气动特性研究,对配置宽速域翼型与常规高超声速翼型的机翼进行了气动力特性综合对比分析。结果表明,配置新翼型的机翼在宽速域范围内整体气动性能更优,说明所设计的宽速域翼型在三维机翼上也具有一定的实用价值。     

二维非结构网格的生成及其Euler方程解

 采用推进阵面法生成三角形非结构网格,用多重网格技术求解泊桑方程以取得结构背景网格.提出了一种高效的边界剖分法,并用它来生成初始阵面.另外,还提出了一种新颖的数据结构——数组链接表,用它来存贮网格数据.用有限体积法在非结构网格中求解了二维可压缩Euler方程,给出了单段翼型和多段翼型的算例,并与实验结果进行了比较和分析.     

塞式喷管三维流场的数值模拟

从曲线坐标下的三维平均雷诺的N-S方程出发,并用k-ε两方程湍流模型封闭方程组,采用二阶精度的精度无波动、无自由参数的耗散差分格式(NND格式),发展了模拟塞式喷管三维流场的数值程序,并对线性塞式喷管和瓦状塞式喷管在不同落压比下的三维流场进行了数值模拟研究.结果较好地预示了塞式喷管的三维流场特征,并且显示瓦状塞式喷管除了圆柱形结构比平板更能承载外,其塞锥压强分布可以进一步改善其受力情况.     

滚珠自转对高速球轴承环下润滑内部流动的影响

高速球轴承正常工作时,处于高速旋转状态,此时油相和气相因重力和离心力的作用在轴承环间剧烈运动。为了更准确地分析轴承环间的两相流动,采取VOF模型进行内部流场的模拟,采用多重旋转坐标系描述部件运动。建立球轴承环下润滑计算模型,分析了考虑滚珠自转因素下轴承内部的流动,并在此基础上探究了转速及供油量对轴承工作状态的影响。结果表明,对比整体模型与滚珠自转模型,发现滚珠自转使得轴承内部油体积分数增大,同时也使得滑油穿透间隙达到外环的能力增加;在考虑滚珠自转情况下,转速的不断增大,使得轴承内部的油相体积分数不断减小,在低转速情况下滚珠自转对流体运动影响较为明显,在高转速情况下公转速度对流体运动起到主导作用,滚珠自转对流体运动影响减弱;供油量不断的增大,使得滚珠自转模型内部的油相体积分数也在不断增大,而且滚珠自转运动会加强滑油在轴承内部的分布。     

高速可压缩来流与凹腔驻涡相互作用及其物质输运特性研究

为了深入认知凹腔驻涡与高速可压缩来流之间的相互作用过程以及来流马赫数对上述动力学过程的内在影响规律,对长深比为1.0的单凹腔基础构型在常温、常压、Ma=0.3~0.7的高速可压缩来流下进行冷态非定常数值模拟,并基于拉格朗日视角研究来流与凹腔驻涡相互作用及其物质输运特性。结果表明：凹腔上方剪切层内旋涡结构仍然具有类似自由剪切层的周期性演化过程,主流与凹腔通过三次掺混进行物质输运;随着主流马赫数的增加,旋涡运动方式及凹腔-主流物质掺混机制都发生了本质上的改变,剪切层内旋涡和凹腔主涡涡量均呈现数量级式增长,剪切层旋涡结构脱落频率近似线性增长并逐渐卷吸至凹腔中部,进入凹腔的主流流体占比降低,但扩散程度和速度加快;此外,源于凹腔的流体流出唇口线后经剪切层的夹带作用回流占比增大,降低了凹腔-主流的整体质量交换水平;驻留参数整体衰减,高马赫数时仅为低速工况的1/2,显著降低了回流区对于流体驻留作用的贡献。     

空间站乘员睡眠区二氧化碳聚集现象

为了分析乘员在没有良好通风的压力舱内停留时,呼出CO₂的聚集过程并评估其可能导致的医学危害,建立了乘员在特定空间站睡眠舱的计算流体动力学模型,对乘员在睡眠区的呼吸过程进行了非稳态模拟.模拟中乘员的CO₂呼出浓度和肺通气量根据生理实验随CO₂吸入浓度变化.通过模拟给出了睡眠区典型位置处CO₂分压随时间的变化规律.最后提出了无通风条件下乘员在睡眠区停留的最大容许时间,睡眠区CO₂浓度监测仪的布置原则和类似情况下为保障乘员生命安全可采取的一些措施.      

气动热CFD计算的格式效应研究

通过对N-S方程的求解计算了二维圆柱绕流及三维绕钝双锥高超声速流动过程中的气动热问题.对于求解方程,在空间离散的格式上,选用了3种上风格式——Roe的FDS,Van Leer的FVS和AUSM+格式以及一种中心差分格式,将各种格式所得的热流结果与实验结果进行了比较.对于AUSM+格式,还应用了minmod limiter 和 双 minmod limiter 两种限制器,并对结果进行了对比.      

Roe格式中不同类型熵修正性能分析

针对Roe格式中较为流行的3类熵修正:Muller型、Harten-Yee型和Harten-Hyman型熵修正,从理论分析入手, 辅以Euler方程为控制方程的激波管问题,前台阶流动和运动激波的双马赫反射3个数值实验,对不同熵修正的性能做了深入研究,得出如下结论: Muller型和Harten-Yee型熵修正方法作用于激波和非物理的膨胀激波2种情况:激波情况下引入的数值耗散,能根本改善"Carbuncle"现象,膨胀激波时其数值粘性也足以使膨胀激波得以耗散;而类Harten-Hyman型熵修正只对膨胀过程起修正作用,对激波情况无效,不能改善"Carbuncle"现象,该类熵修正的数值粘性较小,不足以使膨胀过程正确求解;直接采用 δ 值代替特征值的特征值修正方法效果好于传统的特征值修正方法.