一类逆风ENO线方法的构造及其对Euler方程的应用

本文利用ENO(essentially-non-oscillatory)插值思想,对非线性双曲型守恒律=0提出了一类逆风ENO线方法,使用这种格式其计算量比C.Shu在文[6]中提出的通量分裂方法明显减少,进而同Runge-kutta时间离散相结合,提出了时空高精度差分格式。本文还把这类方法推广到方程组并用于计算Euler方程组的黎曼问题和拟一维喷管问题,得到了比较满意的数值结果。     

一种修正的Osher通量在高阶WCNS中的特性

将一种基于多维修正的Osher通量运用于高阶加权紧致非线性格式(WCNS)中,该修正通量主要在垂直于激波面的界面上增加耗散,能够改善Osher通量的激波捕捉稳定性,同时对边界层和接触间断的分辨率影响非常小。对修正的Osher通量在高阶WCNS中的特性进行研究,通过数值模拟测试了基于Osher通量的WCNS的激波稳定性、热流预测精度、边界层模拟能力、激波边界层干扰模拟能力,并与Steger-Warming通量和Roe通量进行了对比。结果表明修正后的Osher通量比Harten修正的Roe通量具有更好的激波捕捉鲁棒性,而边界层、驻点热流值和激波边界层干扰的模拟则明显优于Steger-Warming通量。上述结果说明了基于修正的Osher通量的高阶WCNS具有较好的激波捕捉特性、热流预测精度和边界层计算能力。     

某类跨大气层飞行器全航程红外辐射特性分析

针对某类跨大气层飞行器飞行过程中需要进出大气层,且飞行速度跨度大等特点,建立了不同空域、不同飞行速度下的机体头部受热模型。当在大气层内飞行时分别使用经验公式和Lees公式计算马赫数小于5和大于5时头部受热情况,当在大气层外飞行时建立了考虑太阳直射、地球反照和地球红外加热情况下的受热模型。利用普朗克定律得到了相同温度下不同波段的红外辐射能量占比情况。最后利用本方法对某类跨大气层飞行器进行了温度和红外辐射计算,得到的计算结果能够大致反映全航程头部温度和红外辐射变化情况。     

利用红外加热笼进行低热流模拟的设计方法研究

在航天器热平衡试验中,散热面区域的低热流密度要求往往难以实现。在以往的试验中,考虑到加热笼覆盖系数对热流密度均匀性的影响,通常选用红外灯阵作为主要的外热流模拟装置。但是红外灯阵体积和热容较大,形状适应性不强,对试验容器空间要求较高。文章利用仿真的手段研究了加热笼的覆盖系数对热流密度均匀性的影响,探讨了利用红外加热笼进行低热流模拟的可行性。     

Upper limit of the total magnetic flux in an active region according to the photometric-magnetic dynamical model

The photometric-magnetic dynamical model handles the evolution of an individual sunspot as an autonomous nonlinear, though integrable, dynamical system. One of its consequences is the prediction of an upper limit of the sunspot areas. This upper limit is analytically expressed by the model parameters, while its calculated value is verified by the observational data. In addition, an upper limit for the magnetic strength inside the sunspot is also predicted, and then, we obtain the following significant result: The upper limit of the total magnetic flux in an active region is found to be of about 7.23 × 10²³ Mx, namely, phenomenologically equal to the magnetic flux concentrated in the totality of the granules of the quiet Sun, having a typical maximum magnetic strength of about 12G. Therefore, the magnetic flux concentrated in an active region cannot exceed the magnetic flux concentrated in the photosphere as a whole.     

CE/SE法模拟爆震波点火及传播过程

运用时空守恒元和解元(CE/SE)法模拟一维爆震波的形成及传播过程,分析产生爆震波的不同机理和DDT(爆燃-爆震转捩)过程的影响因素.采取激波点火和热点点火两种点火方式,分别对应SDT(激波-爆震转捩)和DDT过程,热点点火初始场采用线性温度分布和阶梯性温度分布,整个管道均为常压.用隐式梯形方法实现刚性源项积分.结果表明,SDT和DDT对应于不同的爆震波机理,并且DDT受点火区不同参数的影响,为研究爆震波点火和DDT过程提供了理论依据.     

基于成分的燃烧室出口燃气温度计算方法

针对燃烧室不同温度下燃气组成的不同,提出建立不同燃气模型以推算燃气温度的两种计算方法。基于物质守恒、化学平衡和能量守恒原理建立温度计算方程,根据能量守恒方程的不同表达形式给出平均比定压热容法和焓值法,同时给出两种求解非线性方程组的数学解法:三变量迭代法和牛顿法,并采用VB(Visual Basic)编制燃气分析测试程序实现温度算法。列举计算实例比较两种算法的区别,结果显示:随着油气比的增加,两种计算方法得到的燃气温度相差越大,考虑热离解的高温算法结果更小,原因在于燃气中离解组分的摩尔分数随着温度升高而增加,燃烧实际放热量减小,导致燃气温度减小,所以当燃气温度较高时建议采用高温算法更符合实际燃烧情况。      

基于传热反问题方法的N_2O/C_2H_4预混推进剂燃烧室热流测量研究

测量液体火箭发动机的热载荷是获取燃烧室内部信息的重要方法。为了获取N_2O/C_2H_4预混推进剂燃烧室内壁的热载荷,建立了液体火箭发动机的热流计算的反问题方法,该方法基于对燃烧室壁面温度场的直接求解,通过对轴向多个位置测量温度的反演计算得到燃烧室内壁热流和温度。研究表明:应用文中建立的传热反问题方法能够较为准确地获得热流随时间及空间的分布;热电偶的位置对计算准确性有明显的影响,与理论深度偏差在0.2mm以内的随机深度偏差可导致超过4%热流反演误差;N_2O/C_2H_4预混推进剂燃烧室热流及温度沿轴向逐渐降低,表明燃烧室内的反应释热过程主要在燃烧室头部附近发生。     

模块化永磁直线游标电机等效磁网络及漏磁分析

在模块化永磁直线游标电机的初始设计及优化设计中,使用有限元方法(FEM)计算漏磁系数时需要耗费大量的建模和计算时间。针对这一问题,建立该种电机的等效磁网络模型,在建模期间,考虑气隙中的边缘效应,并计算永磁体组与定子齿在不同相对位置时的等效磁导,得到该种电机在一个周期内不同相对位置时的漏磁磁导和漏磁系数的数学解析表达式。采用FEM验证了该数学解析表达式的准确性。研究结果为模块化永磁直线游标电机的设计提供了理论依据。     

偏心撞击对撞击式喷嘴雾化特性的影响

为研究偏心撞击对撞击式喷嘴雾化特性的影响,建立了求解自燃推进剂冷态射流撞击雾化过程的数值模拟方案,计算了不同偏心度条件下的射流撞击雾化过程。采用树形自适应加密算法直接求解不可压Navier-Stokes方程组,由分段线性的流体体积(VOF)方法对流体界面进行捕捉。结果表明偏心撞击会导致雾场发生偏转,当无量纲偏心度E为1/8时,雾场偏转角度约为9.2°,应控制加工偏差小于该值。随着偏心度的增大,液膜的偏转角度增大,理论推导得到的液膜偏转角度要小于数值计算得到的液膜偏转角度。正心撞击时燃料与氧化剂流强峰值接近,雾场的流强分布呈单峰分布。当发生偏心撞击时,由于燃料与氧化剂部分射流未参与撞击导致流强峰值出现交错,雾场的流强分布呈双峰分布,混合比的空间分布发生较大改变。正心撞击时撞击点下游液滴的速度分布近似呈轴对称分布,而偏心撞击之后的速度分布则呈中心对称分布。偏心撞击导致的射流动量损失使得雾化性能变差,当无量纲偏心度E为1/8时,一甲基肼(MMH)的Sauter平均直径增大约4.8%,四氧化二氮(NTO)的Sauter平均直径增大约5.8%。