Implicit discontinuous Galerkin method on agglomerated high-order grids for 3D simulations

High quality of geometry representation is regarded essential for high-order methods to maintain their high-order accuracy. An agglomerated high-order mesh generating method is inves-tigated in combination with discontinuous Galerkin (DG) method for solving the 3D compressible Euler and Navier-Stokes equations. In this method, a fine linear mesh is first generated by standard commercial mesh generation tools. By taking advantage of an agglomeration method, a quadratic high-order mesh is quickly obtained, which is coarse but provides a high-quality geometry represen-tation, thus very suitable for high-order computations. High-order discretizations are performed on the obtained grids with DG method and the discretized system is treated fully implicitly to obtain steady state solutions. Numerical experiments on several flow problems indicate that the agglomer-ated high-order mesh works well with DG method in dealing with flow problems of curved geome-tries. It is also found that with a fully implicit discretized system and a p-sequencing method, the DG method can achieve convergence state within several time steps which shows significant effi-ciency improvements compared to its explicit counterparts.     

降低CC50汽轮机推力瓦运行温度研究

针对4#机组汽轮机出现正推力瓦运行金属温度偏高的状况,在深入调查4#机组汽轮机现状的基础上,通过与现有标准相比较,确定目标为正推力瓦运行金属温度≤95℃;在具体分析其影响因素的基础上,确定了5个主要影响因素,并针对其采取了相应措施,改造结果表明确实使4#机组汽轮机组正推力瓦运行金属温度低于95℃,达到了预期目标,为汽轮机的改造提供了技术指导。     

极坐标系连续常值推力机动分析

连续常值推力是空间飞行常用的轨道机动方式，在空间交会与星际航行使命中具有重要的应用价值。其中，小推力适合于地球轨道航天器交会机动，而切向或周向推力以及较大的正径向推力可用于脱离地球引力场的逃逸飞行，执行星际交会使命。应用常推力作用下的极坐标系质心运动方程，对机动推力的量值没有限制；在航天器交会应用中，对相对距离也无要求。这种方法可直接获得向径与速度等轨道参数随时间或极角（绕地心的转动角）的变化，便于分析轨道转移与逃逸运动，有助于飞行使命与运动轨迹的设计。特别是，若机动转移的初轨为圆轨道，在推力较小、飞行时间不长的情况下，应用无量纲形式运动方程，可获得具有工程应用价值的近似解。文章给出一些有关的结果与应用案例。     

矢量推力飞行控制不确定性系统中鲁棒控制器的设计

采用加权灵敏度函数环路成形的H∞准则,给出了矢量推力喷管飞行控制系统的一种鲁棒控制器设计方法。飞行条件变化时,该对象模型P不稳定,且模型存在动态不确定性,P~(-1)含虚轴极点。在摄动的模型上进行了鲁棒性能的仿真验证。     

喷流流场旋涡激波干扰的数值模拟研究

本文采用隐式CSCM格式,通过数值求解NS方程,对平板上二维超声速主流中横向喷流流场进行了数值模拟,并与实验和其它计算结果进行了分析比较,结果令人满意。     

H—5飞机灵活推力起飞探讨

根据目前训练飞行不挂炸弹,起飞质量较轻,冬季飞行发动机推力较大这一特点,通过理论分析及计算,得出了正常起飞质量下起飞所需转速随大气温度的关系,并通过反复试飞证明这一结论是正确的,且对使用小转速起飞的优点作了分析。     

混合解析/数值方法在含高频振荡点源项的欧拉方程数值计算上的应用

 通过模型方程分析看到了一个重要的现象, 如果源项涉及的时间尺度远小于对流项的时间尺度, 那么基于对流时间尺度作为步长的传统数值方法, 即使源项相对于墩流项是个小量, 也会导致平均尺度上错误的结果。为了克服这种困难, 采用时间分裂方法, 把方程分裂成含对流项部分的偏微分方程( PDE) 和包含源项的常微分方程(ODE) 基础上, PDE 使用传统的数值方法, ODE 用解析的方法求解。该混合方法在数值格式时间步长小于平均流动时间尺度时, 得到正确解, 而与点源所隐含的时间尺度无关。把这个方法应用在含源流动的Euler 方程的计算中, 计算了翼型振荡问题, 取得了理想的结果。     

二元混压超声速进气道三维流动数值分析

按照MacCormack时间分裂方法，应用NND格式、对流迎风矢通量分裂（AUSM）技术和Baldwin-Lomax混合长度代数湍流模型，对贴体坐标下三维雷诺平均N－S方程进行有限差分离散。数值计算了均匀来流和非均匀来流条件下二元混压超声速进气道三维湍流流动态，分析了该进气道主要性能参数随出口反压、来流马赫数、导弹工作高度及飞行攻角的变化特性，为深入研究进气道性能提供了一定辅助手段。     

塞式喷管主喷管角度特性

首先建立了塞式喷管角度特性的研究模型，并运用简化的塞锥设计方法设计塞锥型面，在此基础上，通过连续改变主喷管倾角来研究主喷管倾角对塞式喷管的推力性能的影响，并确定在给定设计条件下的主喷管最佳倾角，还分析研究了底部高度，内膨胀比，总膨胀比，总压和飞行高度与主喷管最佳倾角的关系以及变化规律，针对特定的塞式喷管实验发动机，进行了初步的变角度实验，并与计算结果进行比较和分析。     

高阶精度非线性格式WCNS-E-5在二维流动中的应用研究

本文采用具有5阶精度的加权紧致非线性显式格式(WCNS-E-5)对定常与非定常二维流动进行数值模拟,研究表明该格式对各类间断有很好的分辨捕捉能力,而且对强间断如激波的计算,即使在高马赫数与高雷诺数条件下它仍具有很好的收敛性与可靠的计算结果.此外,WCNS-E-5在粗网格条件下也体现出优越性.类如WCNS-E-5的高精度激波捕捉方法将为以后开展湍流数值模拟工作提供坚实的技术保证.