超高速撞击中靶后碎片云的内外边界方程研究

研究超高速撞击防护靶形成的碎片云的形状模型对航天器防护结构设计具有重要意义。文章用数值仿真的方法研究靶后碎片云的内外边界模型:将碎片云根据其材料来源分为靶板碎片云和弹丸碎片云,分别根据各碎片云的分布特点建立碎片云形状的内外边界方程;在双纽线边界方程的基础上,修正为由双纽线方程和圆弧方程组成的边界方程。对比整体碎片云的适用率表明,由双纽线方程和圆弧方程组成的边界方程比单一的双纽线边界方程描述碎片云实际分布更为贴合。     

中空纤维膜机载制氮装置的数学建模分析

建立了考虑浓差极化现象的微分方程数学模型,并用正交配置法求解,对中空纤维膜机载制氮装置(OBIGGS)进行了分析,并对部分状态点的计算结果进行了实验验证.结果表明:富氮气体(NEA)中氧气的质量分数随着进气温度的升高而降低,在达到最小值后又呈上升趋势;在进气与排气压力差保持不变的情况下,随着中空纤维膜排气压力的下降,中空纤维膜的富氮气体质量流量逐渐增加;随着中空纤维膜丝长度的增加,丝内气体的压降和富氮气体质量流量均有所增加;中空纤维膜空气分离制得的富氮气体质量流量越大,则所需进气的质量流量越大,且富氮气体氧气质量分数越高.      

基于三维非结构网格的连续伴随优化方法

基于非结构网格,采用Euler方程、连续伴随方法以及自由形面变形技术(Free-form deformation,FFD)结合序列二次规划(Sequential quadratic programming,SQP)优化算法,开发出了三维飞行器的气动外形优化代码。应用Fortran 90语言,动态存储数组和指针,节省了内存。采用Roe形式的二阶迎风格式和LU-SGS隐式时间离散格式求解连续伴随方程,同时采用一种新的参数化坐标求解方法,使得FFD技术更加鲁棒。和有限差分方法求得的梯度进行比较,验证了连续伴随方法所求梯度的正确性,实现了ONERA M6机翼的减阻优化,证明本文开发的优化代码能够实现三维气动外形的减阻优化。     

DLR-F6翼身组合体的高阶精度数值模拟

基于雷诺平均Navier-Stokes(RANS)方程和结构网格技术,采用五阶空间离散精度的加权紧致非线性格式(WCNS)和剪切应力输运(SST)两方程湍流模型,开展了DLR-F6翼身组合体的高阶精度数值模拟研究。主要目的是确认WCNS模拟跨声速典型运输机构型的能力。采用粗、中、细3套网格开展了网格收敛性研究,从气动特性、压力分布、表面流态等方面研究了网格密度对DLR-F6翼身组合体气动特性的影响;采用中等网格开展了来流迎角对气动特性的影响研究。通过与试验数据、CFL3D软件和TRIP软件计算结果的对比,表明网格密度主要影响激波位置和压差阻力系数,同时对翼身结合部分离区大小有一定影响;采用高阶精度计算方法显著提高了气动力系数的模拟精度,力矩系数数值模拟结果与试验的差异有待进一步分析。     

基于风洞试验的非定常气动力微分方程建模方法

为准确描述大攻角非定常运动中非定常气动力的气动特性,研究了以微分方程为基础的非定常气动建模方法。在某战斗机大振幅偏航、滚转单自由度及偏航-滚转耦合运动风洞试验的基础上,对偏航力矩、滚转力矩建立微分方程模型,分析研究不同频率、同一侧滑角情况下,气动力的迟滞时间;然后对高频运动下的气动力进行拟合,进而建立高频运动下的气动力微分方程模型。研究表明:模型可精确预测不同机动下非定常气动特性;建模方法具有较强的工程可行性。     

基于傅立叶平均法下的连续小推力动力学分析

通过将小推力展开为偏近点角的傅立叶级数,并对高斯摄动方程在一个轨道周期上的平均,将原方程的推力转化为仅由14个傅立叶系数表示的控制变量。仿真计算表明,平均化后的高斯方程使计算量与牛顿积分相比显著减少,且对小推力而言有足够的精度。对利用平均化后的高斯方程计算轨道根数时产生误差的原因进行了研究,并进一步分析小推力的范围和小推力近似表达式对上述误差的影响,为今后小推力下非开普勒轨道动力学分析提供了理论依据和参数。     

轨道要素奇异问题的改进四元数方法

运用四元数方法可以在一定范围内解决描绘飞行器轨道运动时轨道要素的奇异问题。但四元数固有的双值性使得在对运动方程进行积分时,其正负选取很困难。为了解决这一问题,采用了改进约束方程的四元数方法,并用该方法描述了拉格朗日行星摄动方程,然后研究了地球扁率摄动对地球同步卫星轨道的影响。仿真结果表明：当偏心率小于1时,四元数可以很好地解决轨道要素奇异性问题。与改进的春分点轨道要素相比,四元数的方法有着更加明确的物理意义和几何意义,用四元数表示的运动变量方程的计算更为简单,积分计算效率更高,而且其计算误差也能达到精度要求。     

非结构网格下可压缩黏性流体的隐式算法

采用牛顿线化和Roe近似Jacobians矩阵隐式算法,实现基于非结构网格的三维定常黏性流动问题的求解.分别对平板、跨声涡轮叶片Mark-Ⅱ和CW-22这三种具有代表性的流动进行数值模拟,并详细分析了流动细节,模拟结果与理论以及实验值吻合良好,准确预测了流场内激波、二次流等重要的流动现象.验证了基于非结构网格的可压缩流体隐式计算方法对于黏性流动问题的计算是正确的.      

二维翼型绕流的CE/SE方法

采用时空守恒元和解元CE/SE(space-time Conservation Element and Solution Element method)法,求解二维Euler方程,开展了翼型绕流的无粘数值模拟研究.用非敏感克朗数计算格式消除克朗数过小引起的数值耗散对解的污染,结合当地时间步长法,解决网格不均匀引起的当地克朗数变化跨度大的问题.对NACA0012翼型的无激波流场进行了二维数值模拟,并与AGARD算例做了对比.结果表明:CE/SE方法的计算结果与AGARD结果吻合得很好,为该数值计算方法对翼型绕流数值模拟的进一步应用奠定了基础.      

CE／SE并行计算研究及在Euler方程中的应用

时空守恒元和解元方法(CE/SE)是一种高精度的数值求解方法.采用改进的CE/SE方法和并行分区技术,针对非结构网格,发展了一套求解二维Euler方程的并行程序.对NACA0012翼型的亚跨声速流动和多段翼型复杂绕流进行了分区并行求解,多区并行计算的结果与实验的结果吻合较好,表明了本文研究方法的正确性和有效性.同时文中的分区方法能实现各处理器之间的负载均衡,有效地节省计算时间.