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

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

低雷诺数翼型和跨声速翼型边界层计算

本文利用边界层动量积分方程和平均流动能积分方程兼容计算了翼面的层流和湍流边界层流动。文中采用一种e~9型转捩判别公式预测翼面存在层流分离气泡情况的转捩位置。并引入对剪应力系数C_r的滞后方程以体现湍流中雷诺应力的历程效应。为避免在翼面流动分离时边界层方程的奇异,采用反解法计算。算例表明,计算与测量结果吻合良好。     

粘性超声速绕体流动的有限谱法计算

首先推导出粘性可压缩流体动力学基本方程组在广义坐标下矢通量分裂的一般形式，并将有限谱ENO格式用于求解超声速流动问题。通过典型算例，验证了这一格式的精度和稳定性。然后进行了圆柱和机翼超声速绕流的数值模拟，得到了与实验和其他计算相符的结果。这表明有限谱法不单具有精度高，应用灵活等多方面的优点，而且突破了传统谱方法全域性的局限，可以成功地应用到一般形状流场可压缩流动的数值计算中。     

轴向非均匀凹槽叶顶的实验与数值研究

为探究不同轴向非均匀凹槽间隙控制泄漏流动的效果和机理,采用实验与数值结合的研究方法,对两种间隙条件下的叶顶凹槽结合渐缩型间隙、均匀型间隙和渐扩型间隙方案对涡轮叶顶泄漏流动的控制效果进行研究。结果表明:减小泄漏流流量与控制叶栅总压损失之间没有直接联系,渐扩型间隙增大了26.7%的泄漏流流量,但在小间隙和大间隙条件下分别减小了2.44%和3.53%的总压损失;渐扩间隙减小总压损失,是通过有效减小通道涡在节距向和展向的尺度,并在一定程度上减小泄漏涡在节距向的尺度实现的;渐扩间隙减小通道涡和泄漏涡的尺度,其原理在于渐扩间隙的布置增强了凹槽内的径向流动,使压力面再附线更靠近叶片,吸力面泄漏涡分离线位置更靠后,从而抑制了泄漏流动。      

乘波飞行器气动力、热特性的数值模拟研究

采用多块分区网格以及并行计算技术对给定乘波构形的高超声速飞行器进行了数值模 拟,分析了飞行器前缘小半径钝化对飞行器气动性能的影响,计算了前缘钝化后飞行器表面 的热流分布状况。结果表明,前缘钝化对飞行器的升力影响不大,对阻力和升阻比的影响较 大。对于曲率半径为1cm的钝化前缘,与原尖前缘飞行器相比,其升力降低了0.78％,阻 力增加6.96％,升阻比下降7.21％。前缘钝化后,乘波飞行器仍具有较好的气动性能, 飞行器前机身可为发动机提供比较均匀的气流,飞行器整体仍基本保持了乘波的状态,热流 主要集中在飞行器前缘上。为了有效防热,需要采用主动式冷却技术。     

基于人工粘性的间断Galerkin有限元方法

发展了一种基于人工粘性的间断Galerkin有限元方法,作出改进以增强适应性,并向非均匀网格推广.选取了典型算例对方法进行验证.一维激波管算例表明,改进的方法在保持计算结果高分辨率的同时,能够更好地抑制非物理振荡.分析得知,当人工粘性方法用于二阶DG格式时,所得计算结果的数值耗散较大,而当格式精度大于二阶时,采用人工粘性方法所得的结果的分辨率较高.通过计算圆柱高超声速粘性绕流,将三阶DG格式与三阶MUSCL格式和五阶WENO格式的结果进行对比,结果表明,该人工粘性方法对于高超声速流动计算也具有一定的优势.     

高速粘性内流的高分辨率高精度迎风型杂交格式

本文在有限体积离散和ＬＵ分解的基础上，构造出一个新的隐式迎风型杂交格式，并用于求解定常流动的稳态解。     

跨音速压气机级三维紊流流场数值模拟

在对跨音速压气机级的动／静双叶排中的三维定常紊流流场数值的模拟中，利用三阶高分辨率ＮＮＤ格式和ＬＵ－ＳＧＳ隐式推进迭代法，既保证了流场中激波的模拟质量和粘性流动特征的正确预估，也实现了求解过程的高效率。对某单级跨音速压气机的计算及其与实验数据的对比证实了本文方法的有效性。将计算推广至多级轴流叶轮机械的情况将会是十分简便的     

分层Reynolds应力模式

本文根据切变湍流结构的特征提出一种分层Reynolds应力模式。在切变湍流中存在湍能生成项达到狭峰的极大值层次(本文称它为活跃层),根据这一特征,在活跃层中保留微分Reynolds应力模式,但它可以简化为常微分方程。在两层活跃层之间以及活跃层和壁面之间的Reynolds应力分布,用满足衔接相容条件的多项式近似。时均湍流场将由以上得到的Reynolds应力代入Reynolds方程求出。用本方法计算了两个比较复杂的湍流例子,结果表明,本方法较常用的k-ε模式有明显的优越性,它的计算时间并不比k-ε模式多,而计算结果与实测结果的吻合则较k-ε模式好得多。另一方面本方法较之完全的微分Reynolds应力模式可大大节省计算时间。     

Bursty Bulk Flow Driven Turbulence in the Earth’s Plasma Sheet

During the last several years significant progress has been made in understanding MHD turbulence in the Earth’s plasma sheet. Due to the statistically transitory properties of fluctuations, finite size and boundary effects, however, issues of fundamental importance remain unresolved. Here we concentrate on such intrinsic features of plasma sheet turbulence as its origin and dynamical nature. In particular, we investigate bursty bulk flow driven multi-scale transfer of energy towards the dissipation scale, and provide evidence for the presence of non-linear interactions. We show that, in contrast with previous results, Alfvénic fluctuations together with 2D eddy interactions may appear as important constituents of turbulence in the plasma sheet.