宽粒径范围颗粒对气相膨胀的影响

利用双流体模型对固体火箭发动机喷管内气固两相流进行计算,研究了两相流中颗粒相与气相的相互作用,分析了颗粒粒径及颗粒质量分数对两相流中气相在喷管中膨胀的影响。结果表明,颗粒粒径造成两相流损失的大小是颗粒比表面积、气相湍流动能及颗粒在喷管内的滞留时间等因素共同作用的结果,两相流损失随着颗粒粒径的增大先增大、后减小,同时随着颗粒质量分数增大而增大;JPL喷管中粒径为2μm的颗粒造成的两相流损失最大,此粒径下,当颗粒质量分数从10％增大到40％时,两相流损失从10％增大到26％。     

固体发动机喷管中气体—颗粒流的数值解研究

本文介绍了一种分析固体发动机喷管中的气体—颗粒流的数值解程序.由轨线模型导出的控制方程,采用时间相关的MacCormack显格式求解;颗粒的速度、轨线、温度用Lagrangian方法计算;受湍流影响的颗粒扩散速度,用网格解颗粒数密度方程得出.本文结果给出了气体—颗粒的耦合程度以及不同粒度和质量分数的颗粒所引起的各种差别.     

固体火箭发动机喷管两相流动计算

采用MacCormark格式与特征线法联合的方法求解了固体火箭发动机喷管两相流场。在跨音速段,气相控制方程用基于MacCormark格式的时间相关法求解,粒子方程采用跟踪粒子轨迹的特征线法求解;在超音速段,气相控制方程用MacCormark两步显格式空间步进求解,粒子方程仍采用跟踪粒子轨迹的特征线法求解,两相充分偶合。在控制方程中考虑了三氧化二铝颗粒的相变。最后对JPL喷管进行了计算,并讨论了不同质量分数对流场的影响。     

固体火箭发动机燃烧室冷态突扩气固多相流动的数值模拟

基于欧拉-欧拉双流体方法,建立固体火箭发动机燃烧室轴对称突扩气固多相流动的统一二阶矩两相湍流模型,采用k-ε双方程模拟气相湍流流动和颗粒相湍流雷诺应力kpg方程描述大尺度颗粒湍流运动、小尺度颗粒各向异性弥散特性以及气固相间作用。数值模拟了突扩流动中的颗粒轴向、径向平均速度和脉动速度以及颗粒-气相脉动速度关联轴向分量的分布特性,计算结果与试验吻合较好。研究结果表明:颗粒轴向速度脉动大于径向速度脉动,轴向平均速度梯度影响脉动速度峰值的发生位置和大小;两相脉动速度关联轴向分量大约是径向分量的2倍。     

非结构网格高超声速MHD流场逆风格式数值模拟

应用逆风通量分裂格式在非结构混合网格上对二维高超声速理想磁流体绕钝头体的流场进行了数值模拟。控制方程为Euler方程耦合Maxwell方程的理想MHD方程,空间离散采用AUSM格式,时间推进采用显式5步Runge-Kutta格式。引用双曲型散度清除技术加强·Ｂ＝0的条件。计算模型为二维钝头体,在高超声速来流条件下,对不同磁感应强度的均匀分布磁场干扰下的流场进行了计算,得到了较满意的结果,并与有限的参考文献进行了对比。结果表明本文发展的方法可用于高超声速MHD方程的求解。     

空间重力环境对稠密气固两相流动特性影响的数值模拟

文章基于稠密气体分子运动论和颗粒动理学,建立气固流化反应器的稠密气固两相流动的欧拉-欧拉双流体数学模型,考虑了颗粒与颗粒之间离散介质特性,大涡模拟求解气相湍流流动,数值模拟反应器内稠密气固两相流体流动行为,得到了不同重力条件下的颗粒浓度、气相和颗粒相速度等计算结果。结果表明：空间重力环境是影响稠密气固两相流体动力特性的重要因素,颗粒脉动强度将随着重力条件的减小而增大。     

含缺陷固体发动机两相流场计算方法研究

通过对国内外含缺陷固体发动机性能预示方法的分析比较，提出采用“N－S方程＋粒子仿真”方法来计算两相流动的新途径，并对某轴对称两相流动的流场进行仿真计算，计算结果表明，该计算方法可行，计算模型合理，能够较为准确地反映两相流中颗粒相对气体流动的影响。     

基于粘附功的复合推进剂AP/基体界面损伤宏细观仿真

为研究复合推进剂AP/基体界面脱湿机理,建立了基于粘聚力界面模型的双尺度有限元损伤分析平台.为确定界面模型的输入参数,采用点滴法和Washburm毛细管上升法,分别测得基体与AP颗粒的接触角,再经Young's方程计算得到界面的粘附功.通过对不同体积分数下简化配方推进剂试件的拉伸试验结果与计算结果的对比显示,两者具有较...     

无喷管固体火箭发动机点火瞬变过程内流场计算

本文提出了计算无喷管固体火箭发动机压力建立过程的 P(x,t)模型,它的控制方程是一组一维非定常两相非平衡流和一组一维非定常两相非平衡流动力学方程,该方程采用 MacCormack 显示差分格式求解.本文还建立了在跨音速和超音速气流流动下的侵蚀燃烧模型,该模型适用于无喷管固体火箭发动机.利用本文的模型可精确预示无喷管固体火箭发动机点火瞬变过程的内弹道性能,并可研究无喷管固体火箭发动机的内流场变化规律.     

复合固体推进剂脱湿模型的研究

研究了复合固体推进剂颗粒界面的脱湿机理.运用弹性-粘弹性对应性原理建立了复合固体推进剂粘弹性脱湿两相球颗粒分析模型 ,利用能量守恒得到了脱湿时间的控制方程,并导出了以界面粘结能表示的临界脱湿应力的简洁表达式.计算结果表明,颗粒的体积分数、粒径和加载速率对脱湿时间有显著的影响.     

潜入式摆动喷管两相内流场数值模拟

气相采用三维薄层近似N-S方程，固相采用拉格朗日坐标系下颗粒轨道模型，计算方法采用显式ENO差分格式和分区多块网格技术，建立了固体火箭发动机内流场模拟软件。对轴对称JPL喷管进行模拟，计算结果与文献进行比较，符合较好。推广到实际三维装药型面时，发现颗粒数目增加导致计算效率降低。采用软件中的气相计算模块对潜入式摆动喷管内流场进行数值模拟，比较了喷管摆动角在发动机喉道附近引起流动变化，研究了非对称流动对喷管内表面压强的影响。     

超声速固体粒子侵蚀的欧拉模型

利用计算流体力学（CFD）技术,引入粒子容积分数的概念,基于欧拉两相流模型,计算超声速飞行器侵蚀表面粒子速度和容积分数分布,模拟侵蚀后的三维表面外形。以超声速计算结果为基础,针对现有拉格朗日法计算侵蚀后退率的局限性,推导出适用于欧拉模型的计算公式。采用本文的数值模拟方法,模拟侵蚀后飞行器的外形,并与实验结果进行对比,验证了本文方法的正确性。计算和分析了不同环境参数对侵蚀量的影响,结果表明,马赫数和粒子容积分数对侵蚀量有较大影响,而粒子容积分数不变的情况下,粒子直径对侵蚀量影响不大。     

耦合凝相-气相动力学机理的二硝酰胺铵燃烧模型

为了研究二硝酰胺铵ADN(Ammonium Dinitramide)固体推进剂燃烧物理化学过程并预测其燃 烧特性,建立一个耦合凝相-气相动力学机理的ADN燃烧模型。该模型基于凝相与气相的总连 续方程、组元连续方程、能量守恒方程及有限速率化学动力学方程而建立,并引入多组元系 统状态方程封闭方程组。模型中包含34种组元,1个固相（凝相）ADN分解总化学反应和165 个气相细节（基元）化学反应,并使用以温度函数表示的物性参数进行计算。应用气相燃烧 模型 对0.6 MPa下ADN燃烧火焰温度、组元摩尔浓度分布进行预测;应用耦合凝相-气相的燃烧模 型对0.2 MPa～36 MPa压强区域内柱状端燃ADN推进剂燃速、燃烧表面温度进行预测,计算结 果与文献报道试验数据较吻合。说明该燃烧模型能够较准确描述ADN气相燃烧波结构和ADN固 体推进剂燃速特性。
     

两相跨音速喷管流动

对气相和颗粒控制方程分别应用时间相关法和特征线法,完成了两相跨音速喷管流场的计算。这是一个两相完全耦合的数值解方法。计算表明,该计算方法是成功的。     

铝冰发动机内流场的数值计算

为了使用数值模拟的方法计算铝冰发动机的性能,用颗粒表面反应模型和气相反应模型模拟铝颗粒在铝冰发动机燃烧室中与水蒸气的燃烧过程,用欧拉-拉格朗日方法计算颗粒沿轨迹的参数,分析了数值模拟的结果,并进行了相同尺寸的铝冰发动机实验,把数值模拟结果与实验结果进行了比较。数值计算得到的燃烧室稳态工作压强约为9.38 MPa,与实验结果接近,燃烧室平均温度为2950.65 K,相比热力计算得到的推进剂燃烧温度略低。通过对铝冰发动机的内流场数值计算,得到了与实验相符合的结果,验证了数值计算模型的有效性。     

过载状态下固体火箭发动机燃烧室内二相流动数值模拟

对过载状态下固体火箭发动机燃烧室内二相流动进行了数值研究，建立了过载状态下二维平面颗粒轨道模型，应用二阶迎风格式构造了有限体积计算框架，对有过载和无过载多种工况燃烧室内流动进行了计算，得出过载对颗粒相和气相的某些影响，结果对于发动机工作过程研究及发动机热结构设计具有一定指导意义。     

粘弹塑性界面的断裂特性

研究张开型粘弹塑性界面断裂。用傅立叶正、余弦变换及逐段定积分变换方法将边值问题的控制方程化为奇异积分方程组。解方程后计算了裂纹尖端塑性区尺寸及裂纹尖端张开位移(COD：crack—tip opening displacement)，并给出了能量释放率算式。结果表明，裂纹尖端塑性区尺寸和COD均随两种材料的最小屈服极限的增加而减小；随时间的增大，COD先增长后衰减，最后渐近地逼近于定值。     

微型脉冲推力器点火启动过程计算与点火药量选择

建立了微型脉冲推力器点火启动模型,模型中考虑了两相流动和颗粒对推进剂表面的冲击传热。根据数值计算结果,给出了不同产物颗粒含量下的点火药量选择参考范围。该参考范围对于产物颗粒含量小于50%的点火药是适用的;由于过多的颗粒含量将对装药表面产生强烈的热力冲蚀破坏作用,产物颗粒含量超过50%的点火药不适合微型火箭发动机。     

基于SPH方法的幂律流体突缩型管路流动过程仿真

应用光滑粒子流体动力学(SPH)方法进行了幂律型流体的突缩型管路流动过程仿真。推导了SPH形式的流体动力学控制方程,应用罚方法施加边界条件,应用人工应力消除拉伸不稳定现象,应用XSPH方法规范粒子秩序。提出了幂律型本构关系的SPH求解方法,推导了剪切速率及粘性项的计算公式。应用Poiseuille流算例对本文的幂律型本构求解方法进行了验证,获取了幂律型流体在突缩管中的流动特性,并与水的流动特性进行了对比分析。讨论了流动特性的成因。