基于动网格的喉栓式推力可调喷管内流场数值模拟

以FLUENT软件为工具,利用动网格技术,建立了喉栓式推力可调喷管内流场动态特性分析模型。通过网格的合并与分割,较好地解决了喉栓调节运动所导致的计算区域的瞬变问题,并分析了喉栓调节运动速度对喷管轴对称二维内流场动态特性的影响。分析结果表明,随着喉栓的调节运动,喷管内流场动态下的压强建立与稳态下的压强建立相比存在着明显的延迟,且随着喉栓调节运动速度的增大,延迟现象越明显。当喉栓完全进入喷管几何喉部位置时,喷管推力达到最大值。     

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

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

火箭发动机塞式喷管流场的数值模拟研究

本文建立了计算塞式喷管流的物理数学模型，通过求解采用Ｋ－ε紊流模型的二维Ｎ－Ｓ方程组，发展了相应的数值计算方法，对在不同环境压强下某塞式喷管的流场进行了数值模拟。数值模拟的结果表明：塞式喷管扩张段的膨胀过程能够自动适应环境压强的变化；环境压强由高变低时，回流我由开式结构变成闭式结构，形成闭式结构时，塞锥底部压强近似于常数，受环境压强影响不大；塞式喷管的塞锥长度减小到一定程度后，塞锥长度对流场结构影     

复合结构喷管热强度有限元分析

本文采用有限元法对复合结构喷管进行了二维瞬态温度场和应力场计算。在温度场计算中采用了网络变换法对考虑烧蚀动边界并有效地消除了大温度梯度情况下温度振荡现象，算例表明，本文提供的分析方法简单可靠，易于工程实践采用。     

喷管内横流流场的数值模拟

采用显式的原始ＭａｃＣｏｒｍａｃｋ时间分裂格式和改进的ＴＶＤ ＭａｃＣｏｒｍａｃｋ格式分别求解了二维轴对称喷管内横流流场的Ｎ－Ｓ方，得到了与理论分析结果相近似的流场结果图谱，并对以上两方法得到的结果进行了比较，结果表明：ＴＶＤ的方法具有更高的精度。文中还介绍了一种基于ＴＴＭ方法的改进的网格生成方法，用于生成贴壁性和正交性都比较好的网格，实现了完全由边界网点分布来确定内部网格点的分布。     

RBCC发动机性能分析模型改进方法研究

对已有RBCC发动机性能分析模型的改进方法进行了研究。采用最小吉布斯自由能法,建立了一次火箭燃烧室热力计算模型,并采用冻结流假设进行喷管热力计算;针对以H2O2(L)/JP-10(L)为推进剂的一次火箭进行了热力计算,并与CEA和CHEMKIN等化学热力学软件计算结果进行了对比校验,验证了所建模型计算结果的准确性。建立了采用CFD软件求解二维Euler方程的后体性能分析方法,并基于CFD软件提供的接口函数和PYTHON脚本技术设计了后体自动性能分析流程,自主实现了后体流场分析过程的几何造型、网格划分、边界定义、CFD求解器设置、CFD方程求解及性能参数计算。     

固体火箭发动机喷管效率计算

采用半经验法计算了固体火箭发动机喷管的效率，即用计算流场的方法确字喷管二维两相流损失和边界层损失，用SPP经验法预示了喷管的化学动力学损失，喷管烧蚀损失和喷管潜入损失。利用该方法对几个实际固体火箭发动机喷管效率进行了计算，计算结果与实际结果比较符合，精度偏差在1％之内。     

一种新的界面映射推进方法及其在气动弹性力学中的应用

非线性气动弹性研究中,涉及到非线性的流体动力学(CFD)和非线性的结构动力学(CSD )耦合问题,在耦合界面上不仅要进行两场之间信息的传递,而且要有匹配的时间推进格式 。针对计算流体动力学（CFD）和计算结构动力学（CSD）的耦合计算方法,基于边界能 量守恒,发展了一种新的界面处理方法,包括：① 界面映射方法,该方法基于耦合边界上 局部的网格信息,将两场之间载荷信息和位移信息的转换放在同一个映射矩阵中来处理,克 服了占用大量CPU时间和内存的需求;② 界面推进方法,该法基于松耦合计算流程,通过引 入半步交错推进,达到了二阶时间精度,提高了计算精度和效率。最后将该界面处理方法应 用于柔性大展弦比机翼的气动弹性计算和AGARD445.6机翼的动响应分析中。结果表明该方法 能够高效、高精度地处理不同物理场之间的数据交换和时间推进,并具有处理复杂几何体非 线性气动弹性问题的能力。     

组合式数值网格生成

本文针对流场特点,采用组合式数值网格,即流场内部采用球对称贴体网格、外部采用轴对称网格,对该种结构网格的分布和两组网格在交界处的结合进行了研究.通过泊松方程把物理网格变换成计算网格,并对网格分布的控制进行了研究.提出并采用了滑移边界,保证了边界上网格的正交性.研究方法和结论可以应用于类似问题的数值分析中.     

空间目标碰撞概率计算方法研究

空间目标碰撞概率的计算是航天器进行空间碎片预警和规避机动的基础,提出了一种基于压缩空间和坐标旋转的碰撞概率计算方法。定义了相遇坐标系和积分计算坐标系,通过坐标转换关系将位置误差协方差投影到相遇坐标系中,旋转相遇坐标系后得到积分计算坐标系,将计算碰撞概率的问题转化为2维概率密度函数在圆域内的积分问题。通过压缩空间的方法,将不等方差概率密度函数在圆域内的积分化为等方差概率密度函数在椭圆域内的积分,然后通过极坐标变换将二重积分化为一重积分形式的概率表达式。对一次空间碎片碰撞的实例计算得到碰撞概率为10 -3量级,表明该方法是正确有效的。     

无喷管火箭发动机侵蚀燃速辨识

本文提出用无喷管火箭发动机终止燃烧后测得的燃层厚度,辅之以测得的p(x,t)曲线,来辨识无喷管火箭发动机工作条件下推进剂侵蚀燃速规律的一种方法。该方法克服了以往用P-t曲线间接辨识燃速方法中存在的问题,避免了复杂的内弹道计算,大大缩短了计算机时,实践表明,这种方法辨识结果稳定,用其结果计算的燃层厚度与实验值符合较好。     

矢通量分裂法在喷管超声速喷流计算中的应用

研究了一种适用于喷管超声速喷流场的数值模拟方法。采用矢通量分裂法对喷管超声速无粘喷流场进行了数值模拟。将矢通量分裂法应用到轴对称Euler方程,并对发动机喷管尾部超声速射流进行了数值计算。计算结果与实验纹影图所反映的流动特性相吻合,与其他高精度数值格式的计算结果相比,矢通量分裂法的计算结果在轴线反射点附近具有很高分辨率,表明该方法对激波具有较好的间断捕捉能力,不会产生伪振荡或抹平现象。     

固体火箭发动机燃烧室喷管统一流场计算

采用ＳＩＭＰＬＥ算法的可压缩形式求解固体火箭发动机燃烧室喷管内的全速度统一流场;在非交错的配置网格基础上用有限体积法离散Ｎ－Ｓ方程,采用强隐式算法（ＳＩＰ）求解线性代数方程组;采用边界标志法实现复杂几何形状下的多区域统一计算;采用喷管外型逼近法保障计算稳定性;通过与同类计算结果的对比表明：上述算法的组合能够成功地获取固体火箭发动机内流场的整体结构。     

用于定常二维跨声速喷管流动分析的N—S解程序

本文介绍了用于分析平面和轴对称定常跨声速喷管流动的有效的Navier—Stokes解算法。本法是利用跳格数值算法在物理上取的坐标中求解守衡型非定常可压缩N—S方程。分析了几种定常跨声速粘性喷管层流,给出了计算结果,并与实验结果进行比较,证实了数值方法的精度和效率。     

无喷管助推器非定常动边界内流场数值模拟

无喷管助推器依靠药柱的扩张面实现对气流的加速,燃面的线性增加及喉面的指数增加,使其内流场呈现明显的非定常性;侵蚀燃烧、压强的不均匀分布及药柱变形,使无喷管助推器的工作过程十分复杂。通过对FLUENT软件进行二次开发,建立了无喷管助推器的非定常动边界内流场数值模型,考虑了推进剂的点火过程和侵蚀燃烧,可获得各个时刻的燃面位置及内流场分布,从而可对发动机的内弹道性能进行较准确地预示。通过对某无喷管助推器的工作过程进行模拟,内弹道计算曲线与试验曲线吻合得较好。     

图像型火灾探测技术中的火灾定位研究

针对图像型火灾探测技术中的定位问题,提出一种在摄像机内外参数未知的条件下解决火灾二维定位的方法。即首先根据定标板角点的世界坐标和图像坐标计算透视投影矩阵,并在采用图像处理技术获取火灾图像的中心坐标的基础之上,通过求解线性系统获取火灾的实际位置。     

叶-盘系统振动特性有限元分析程序实现

引入空间八节点等参数单元,对叶-盘结构进行了有限元网格划分.利用自行编制的有限元循环对称结构振动分析程序对叶-盘结构(泵叶轮及涡轮叶轮)进行了振动特性分析,并与ANSYS软件计算结果进行了对比.经过对比,本程序的计算结果与ANSYS软件的计算结果吻合良好,从而验证了本计算程序具有良好的工程计算精度.     

大椭圆轨道卫星姿态基准坐标系与运动学研究

针对大椭圆轨道卫星在轨运行中特殊光照条件导致的能源供给、热控设计等工程研制难点,兼顾光照时间和遥感任务,提出将基于太阳矢量与地面物点矢量的日地观测坐标系作为卫星姿态控制的基准坐标系。研究了相应的卫星姿态运动学,针对日地观测坐标系的角速度矢量在物理层面难以写出解析表达式的问题,提出一种解析计算方法,给出了算法流程,推导了日地观测坐标系相对惯性坐标系的角速度矢量在惯性坐标系中的计算公式。与数值差分法的仿真对比结果表明:该解析算法正确可行,其曲线特性和计算精度优于数值差分法。研究对卫星研制具有理论意义和应用价值。     

旋转发动机的喷管流场计算

采用时间相关法计算了旋转发动机的喷管气相流场,内部点使用 MacCormack 二步显格式,对入口边界和固壁边界点的参数计算,除采用了物理边界条件外,还采用了参考平面上的特征方程,计算结果表明:旋转效应沿着喷管轴向逐渐减弱,当转速不很高时,旋转对气相流场参数分布影响不大,可近似按不旋转计算;当转速很高时,应考虑旋转效应。     

大型火箭发动机喷管内流场数值分析

我对某大型火箭发动机的喷管的喷管结构，采用显式的TVD Mac-Cormack格式，结合两方程k- ε紊流模型求解喷管内流场的二维轴对称N－S方程，获得详细的喷管内流场结构分布，其结果同理论分析的结果非常相似，从验证了用数值模拟方法取代部分实际试验的可行性及可信性。