同轴双剪切气-气喷嘴数值模拟

通过求解 k-ε 湍流模型的Navier-Stokes 方程组,对以气氢/气氧为推进剂的同轴双剪切喷嘴燃烧室进行数值模拟研究.研究结果表明:与传统的同轴剪切喷嘴相比,双剪切气-气喷嘴使推进剂有两个剪切燃烧面,且出口尺寸变化不大;双剪切喷嘴中心氢与氢总质量流量的比例是双剪切喷嘴的关键设计参数,当比例值为0.3时,能充分发挥双剪切喷嘴两个燃烧面的优势,使双剪切喷嘴能在大流量工况下实现高的燃烧效率.     

三维Euler方程自适应八叉树结构各向异性直角网格算法

采用各向异性直角网格对相应的三维复杂流场进行了模拟 ,提出了物面三角化方法及网格与物面的相交判断过程。用以中心差分为基础的 Jameson的有限体积法对 Euler方程进行求解。计算结果表明 ,各向异性直角网格算法对相应流场进行模拟 ,计算时间较短而效果较好     

Euler computation of the nearfield wake vortex of an aircraft in take-off configuration

Numerical simulations based on the three-dimensional Euler equations are used to investigate the predictive capability of an Euler code for calculations of the nearfield wake of a narrow-body airliner wind tunnel model in take-off configuration up to a half span behind the wing tip trailing edge. Simulation results on both structured and unstructured grids are presented. The results on the block-structured grid were obtained within the scope of the EU-project EUROWAKE. The simulation quality of the vortex formation and spatial development is analysed by comparison to wind tunnel measurements of the spanwise lift distribution available from the EU-project DUPRIN II and to experimental PIV data available from the EU-project EUROWAKE.     

带副翼偏转的三角翼自由滚转运动数值模拟

 通过耦合求解非定常Euler/Navier-Stokes方程和单自由度滚转运动方程,对带副翼偏转的65°后掠角尖前缘三角翼WI1-SLE自由滚转运动进行了研究,Navier-Stokes方程的求解采用基于Spalart-Allmaras湍流模型的脱体涡模拟(DES)。在多块结构网格上,应用基于弧长的无限插值理论(TFI)生成变形网格,实现副翼偏转,而三角翼的滚转运动则通过网格的整体旋转实现。结果表明:Euler方程和DES方法均准确地模拟出了三角翼在滚转运动过程中存在的3个平衡位置。出现平衡位置的原因分别是:①流动对称性;②机翼左侧发生涡破裂的分离涡与右侧分离涡相互平衡使得滚转力矩为0,并且平衡位置仅与三角翼两侧涡强的差有关;③副翼偏转和左右机翼不对称分离涡涡强差产生的滚转力矩相互平衡。此外,滚转运动对副翼偏角幅值很敏感,幅值的微小改变会影响最终的平衡位置和向平衡位置运动的路径。     

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

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

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

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

一组新的四元数轨道要素建模方法

针对传统以欧拉角为参数的轨道要素的奇异性、不确定性以及计算效率低等问题,提出了一种新的四元数轨道要素。建立了新轨道要素与经典轨道要素,以及新轨道要素与惯性系下位置、速度的相互转换关系,推导了基于新轨道要素的高斯摄动方程。以J₂项摄动下的轨道推演为例进行仿真验证,结果表明新轨道要素不仅在圆轨道与赤道平面轨道处不存在奇异性和不确定性,而且由于新轨道要素不涉及三角函数运算,新高斯摄动方程积分效率和计算精度明显提高。     

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

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

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

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

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

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