火箭发动机气体-颗粒两相流双流体模型研究

探讨了两相流双流体模型中颗粒相压强、两相界面压强以及两相界面速度等若干关键问题,从理论上推导了双流体模型的控制方程,并讨论了颗粒相控制方程的双曲性.把LU隐式时间格式和Van-Leer矢通量分裂格式应用到两相流双流体模型中,并对一维等截面传热摩擦管道和轴对称JPL喷管中气固两相流进行了计算,计算结果与文献结果吻合良好,所建数学模型正确,所用数值计算方法可靠.     

空间高阶差分近似对三维ADI-FDTD数值色散的影响

研究采用空间高阶差分近似的三维交替方向隐式时域有限差分法(ADI-FDTD,Alternating-Direction Implicit Finite-Difference Time-Domain)数值色散问题,首先推导了采用空间高阶差分近似的数值色散迭代公式,分别对最小相速方向和最大相速方向进行二阶、四阶、六阶、十阶空间差分数值色散误差的数值计算和比较,并分析数值色散对空间差分阶数的依赖关系,最后发现在均匀网格中采用四阶空间差分能得到较好的效果.     

基于自适应扩展卡尔曼滤波的载波跟踪算法

精确的载波相位测量是精密测距中一个很重要的研究点。针对传统扩展卡尔曼滤波(EKF)的固定设计在先验信息不充分和动态变化环境中存在的不足,提出了一种基于自适应扩展卡尔曼滤波(AEKF)的载波跟踪算法。该算法通过实时监测滤波器新息或残差的动态变化,以修正状态噪声方差和观测噪声方差,进而调整滤波器增益,控制状态预测值和观测值在滤波结果中的权重。理论分析和仿真结果表明,本算法充分利用了观测信号的统计特性,克服了传统扩展卡尔曼滤波算法的不足,能够获得更好的载波跟踪性能。     

各向异性Chaboche黏塑性本构方程隐式应力积分算法

针对各向异性Chaboche黏塑性本构方程,利用有效参数控制法推导了该方程的隐式积分算法,该算法把所有未知量都用同一个控制参数表示,并形成控制方程,从而把本构方程应力积分的问题转化成一个非线性控制方程的求解问题.该积分算法程序,与通用有限元软件Marc结合,通过单个单元对单向拉伸、循环、松弛的响应的实验数据进行了数值模拟,通过与显式算法对比,该隐式算法的稳定性和收敛性有所提高.构件的对比分析也表明隐式算法比显式算法具有较高的效率.      

隐式耦合求解燃烧室内喷雾两相流场

在任意非结构网格下,采用开源程序对某型涡喷发动机环形燃烧室内的喷雾燃烧过程进行数值模拟.首先选用k-ε双方程模型、EBU(eddy break up)-Arrhenius紊流燃烧模型、颗粒随机轨道模型等模拟流场中的物理过程.然后用隐式耦合算法求解气相流动、拉格朗日方法跟踪液滴相在流动区域的运动,并考虑气相与液滴相的完全双向耦合作用.仿真实验得到的 燃烧室内温度分布、马赫数分布、速度分布、压力分布均较合理,温度分布与文献中的基本一致,这表明隐式耦合方法能够较好地模拟燃烧室内流场,反映其内部流动特点.      

激光冲击强化提高主动连杆振动疲劳性能

针对某型发动机放气活门主动连杆转接R处的疲劳断裂问题,进行了主动连杆激光冲击强化工艺研究;完成了原型件和激光冲击强化、喷丸及加大R3种处理工艺后主动连杆振动疲劳对比试验.结果表明,在317MPa的应力水平下,激光冲击强化、喷丸和加大R后主动连杆疲劳寿命分别为原型件的5.24倍、3.89倍及1.36倍;与其它2种工艺相比,激光冲击强化具有强化效果好、工艺稳定性高、易于实现等优点,是截至目前减少主动连杆疲劳断裂故障的最佳处理工艺.主动连杆激光冲击强化后表面粗糙度较小,产生深度约1.4mm的高数值残余压应力层,表层晶粒细化至纳米晶并伴随高密度位错.三者共同作用可以有效缓解转接R处的应力集中,抑制疲劳裂纹在转接R处的萌生和扩展,是主动连杆振动疲劳寿命提高的主要原因.     

3D-Cf / Mg 复合材料的热残余应力研究

采用随动硬化有限元模型,研究了三维碳纤维增强镁基复合材料(3D-Cf/Mg)基体热残余应力的大小、分布以及不同工艺处理对热残余应力的影响。同时,对复合材料进行了高、低温处理,利用XRD定性分析了处理前后复合材料基体热残余应力的变化。计算结果表明:经过-196℃低温处理后,复合材料基体的平均Mises热残余应力由169.06 MPa减小至55.29 MPa;高温处理后,基体平均热残余应力几乎不变。该结果与实验结果吻合,证明了低温处理能明显降低复合材料基体的热残余应力。     

旋翼跨音速非定常黏性绕流的高效CFD模拟方法

为提高旋翼跨音速黏性绕流CFD模拟的效率,建立了一套基于隐式LU-SGS算法和OpenMP并行策略的旋翼非定常流场高效数值求解方法。首先,基于二维剖面翼型的Poisson方程求解和网格插值、翻折方法生成绕桨叶的O-C-O型贴体正交网格,并采用高效的"扰动衍射"(Disturbance diffraction method,DDM)挖洞方法、并行化″Inverse map″(IM)的贡献单元搜索方法,构建了旋翼运动嵌套网格系统。在此基础上,以耦合S-A湍流模型的非定常RANS方程为主控方程,对流通量采用高精度的Roe-MUSCL格式进行离散,时间推进采用隐式LU-SGS双时间方法,同时采用OpenMP并行策略加速旋翼流场求解。最后,运用所建立的方法分别对悬停ONERA 7A旋翼、前飞Caradonna-Tung(C-T)旋翼无升力状态及SA349/2旋翼有升力状态的气动特性及涡尾迹特征进行了模拟,通过与试验值对比验证了文中方法在旋翼非定常流场CFD模拟中的有效性和高效性。     

基于并联贮箱结构的卫星推进剂剩余量测量方法

星上推进系统推进剂剩余量测量是卫星在轨管理的重要工作,事关卫星剩余寿命估计和离轨时机的选择,对于提高卫星效率具有重大的意义．中国现有的推进剂剩余量计算方法只适用于传统的双贮箱结构,不适用于桁架式卫星平台的多贮箱并联结构,为此必须开发新的计算方法,同时提高计算精度．比较多种推进剂剩余量测量方法,重点论述两种可用于并联贮箱结构的测量方法．     

碳酚醛-铝板中二维X射线热击波数值模拟

文章利用平面应变正交各向异性材料率相关性的弹塑性本构模型和各向同性材料理想弹塑性本构模型,采用有限元方法编写程序,对X射线辐照碳酚醛-铝板时诱导的热击波进行数值模拟,讨论了热击波的传播规律。结果表明,在1 keV的 X射线辐照下,X射线穿透极浅,材料表面发生气化现象,应力峰值较大;在3 keV的X射线辐照下,X射线穿透较深,热击波呈现双波结构,材料的层裂更明显。