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

为了快速预测发动机冷却通道内碳氢燃料在考虑热裂解时的流动传热特性,基于正癸烷热裂解反应机理,建立了一套模拟正癸烷裂解吸热和超临界压力传热现象的快速算法。采用三维物性库查表算法计算裂解反应混合物的热物性,同时对组分输运方程进行简化,简化后只需求解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方程,给出了单段翼型和多段翼型的算例,并与实验结果进行了比较和分析.     

O形圈密封载荷衰减研究

对两种硅橡胶材料O形圈的载荷衰减规律进行了试验研究,结果发现,O形圈的载荷衰减规律与O形圈材料的性能有关;温度对O形圈的衰减规律也有很大影响,温度越高,O形圈的载荷衰减速度越快.在试验的基础上对O形圈载荷衰减的理论模型进行了研究,采用无穷多个Maxwell模型并联的方式,找到了一个能较好地表示O形圈载荷衰减规律的模型.      

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

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

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

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

薄壁管轴压胀形关键工艺参数及成形极限

利用粘性介质传压,对薄壁管轴压胀形进行了研究;基于塑性理论和分散性失稳理论,推导了轴对称胀形内压力分布规律,得出了内压加载区间和轴向加载关系;建立了基于载荷控制的加载函数,提出了送料控制方法,优化了加载模式;基于集中性失稳理论,计算了胀形压力极限,建立了均匀变形极限和极限应力比的概念,形成了零件可成形性的基本判据;分析了管材成形过程中应力变化特点,揭示了管材胀形区由拉-压应力状态迅速发展到双向受拉的机理.通过实验并利用有限元方法对薄壁管轴压胀形过程进行了模拟计算和研究,验证了理论分析结果.      

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

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