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

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

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

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

考虑燃料雾化的气液两相连续旋转爆轰数值模拟

为研究气液两相连续旋转爆轰发动机燃烧室内爆轰波的传播特性,以汽油为燃料,富氧空气为氧化剂,建立了欧拉-拉格朗日模型进行二维数值仿真,其中气相方程采用时空守恒元与求解元方法求解,液相方程采用标准四阶龙格库塔法求解。在两相旋转爆轰模型还考虑了液滴雾化破碎过程。计算结果表明：起爆后形成的初始爆轰波经过初始燃料填充区域后逐渐衰减,随后入口附近新生成的压力波经过一系列发展形成了自持稳定传播的旋转爆轰波;旋转爆轰波的传播模态受燃料与氧化剂的喷注压力和氧化剂填充比影响,在不同工况下旋转爆轰波呈现出4种传播模态,即稳定单波模态、稳定双波模态、不稳定双波模态和不稳定单波模态;在双波模态工况下,燃烧室内初始只形成1个爆轰波,后由入口附近局部爆炸产生的压力波发展为新的爆轰波,转化为双波模态后爆轰波的强度略有下降,但燃烧室整体推力更加平稳。     

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

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

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

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

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

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

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

为了分析乘员在没有良好通风的压力舱内停留时,呼出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"现象,该类熵修正的数值粘性较小,不足以使膨胀过程正确求解;直接采用 δ 值代替特征值的特征值修正方法效果好于传统的特征值修正方法.     

激光功率与底面状态对选区激光熔化球化的影响

为研究激光功率与底面状态对选区激光熔化熔池流动的影响,基于离散单元法建立了选区激光熔化铺粉模型,采用粒径分布与实验相符的马氏体时效钢粉末分别铺展到平坦底面和增材底面上,将计算获得的粉末分布导入到基于有限体积法建立的选区激光熔化熔池计算流体力学模型中,研究激光功率和基板底面粗糙度对熔池流动和熔道表面形貌的影响。采用激光单道扫描实验验证铺粉模型和选区激光熔化模型。结果表明：随着激光功率的降低,单位长度的球化数量增加;由于增材底面使熔池润湿性变差,同时又对熔池流动行为产生扰动,使得增材粗糙底面上熔道的球化数量增加。选区激光熔化铺粉模拟及激光单道扫描模拟结果与实验结果吻合较好。本研究可为选区激光熔化工艺中工艺参数的选择提供理论指导。