高精度IPDG湍流模拟及通量格式数值特性北大核心CSCD

为更高精度数值求解复杂湍流问题,在内罚间断伽辽金(internal penalty discontinuous Galerkin,IPDG)方法框架内发展了基于SST k-ω模型的湍流模拟方法,通过对亚/跨/超声速工况下流场的数值计算与湍流特征捕捉,验证了方法的适用性,进而系统分析了AUSM、Lax-F、HLL和Roe 4种通量格式在IPDG湍流模拟中的数值特性。结果表明:AUSM格式在超声速工况下脱体激波面“褶皱”现象明显;Lax-F格式和HLL格式数值耗散大,在激波解析方面精度较低,且Lax-F格式在激波脚后诱导产生流动分离;Roe格式具备宽速域适用性,计算精度较高且能精确解析激波结构。     

超声速横喷干扰湍流场数值模拟

本文探讨了超声速平板横向喷流干扰湍流流场的数值模拟方法。分别采用Roe、Harten-Yee及AUSM 格式对三维雷诺平均Navier-Stokes方程进行离散,湍流模拟应用Baldwin-Lomax代数模型,讨论了B-L模型中Fy最大值的不同选取方式。分析了数值格式、网格划分对计算结果的影响。将来流马赫数为5、喷流马赫数为3的实验结果与计算结果进行比较,选择了与实验符合较好的数值模拟方法,为进一步利用数值计算方法研究喷流干扰流场特性奠定基础。     

二维高亚声速空腔流激振荡的数值模拟研究

采用基于P .L .Roe的近似Riemann解的修正Osher&Chakravarthy (MOC)三阶TVD有限差分格式 ,数值求解二维雷诺平均全Navier Stokes方程 ,并用Cebeci Smith代数湍流模型 (对腔内区域作修正 )来模拟湍流效应 ,时间方向的积分采用四阶Runge Kutta方法 ,通过对GAMM超音速前台阶绕流的计算 ,验证了格式及程序的有效性 ,对高亚声速来流下的空腔流动作了数值模拟研究 ,取得了较好的效果     

斜切喷管欠膨胀超声速射流的数值模拟

采用三阶MUSCL TVD格式解三维可压缩平均雷诺纳维尔-斯托克斯方程组，湍流模型为Spalart-Allmaras代数模型，数值模拟了两种斜切角喷管在喷口压比为2.0，2．8，3．4和4．0条件下的欠膨胀超声速射流场，获得了射流场流谱和参数分布，计算结果和实验数据符合良好，并给出了喷管斜切角和喷口压强比与射流偏转角和扩张角的分布曲线。     

基于RANS/LES方法的超声速底部流场数值模拟

分别采用基于两方程k-ω剪切应力输运(SST)湍流模型的延迟DES(DDES)、更改的DDES(MDDES)和改进的DDES(IDDES)方法,并引入可压缩修正,结合三阶MUSCL-Roe和五阶WENO-Roe两种空间离散格式,针对超声速底部的复杂流动现象,开展了数值模拟研究。计算结果表明本文方法能够捕捉到超声速底部流动中丰富的湍流结构,通过分析计算结果对超声速底部的流动机理有了进一步的认识,为下一步的超声速底部流动减阻改进和雷诺平均NavierStokes/大涡模拟(RANS/LES)方法在非定常高可压缩性流动中的应用提供了参考。通过对比分析不同空间离散格式的计算结果研究了数值耗散对计算的影响,五阶WENO-Roe格式的计算结果与实验结果吻合良好;对不同RANS/LES混合方法的计算结果进行了对比分析,结果表明IDDES方法在近壁区的表现优于DDES和MDDES方法。     

求解可压缩流动的低耗散E-CUSP迎风格式性能分析

采用对比分析的方法,对求解可压缩流动的低耗散E-CUSP(LDE)迎风格式的性能进行测评,目的是进一步更好地将其用于分离涡数值模拟。对比格式采用广泛使用的Roe格式。控制方程为Euler和雷诺平均Navier-Stokes方程耦合Spalart-Allmaras一方程湍流模型,采用全隐式无分裂Gauss-Seidel线松驰时间推进方法。结果表明,LDE格式具有强的捕捉激波能力,计算结果与Roe格式一致,和实验数据吻合良好。相对于Roe格式,LDE格式具有较快的收敛速度。     

马赫数1.95超声速欠膨胀喷流声辐射数值分析

针对马赫数为1.95的欠膨胀超声速喷流声辐射特性,采用高精度计算格式的大涡模拟（LES）方法进行数值研究。通过对喷流平均流及湍流脉动统计结果的对比,确定数值方法的精确性。细致分析了超声速喷流流动特征,特别关注了其中激波胞格和湍流相互作用现象。基于流动与声场辐射间的关联性,分析了超声速喷流马赫波、宽频激波噪声辐射特性及形成机理。数值结果表明超声速喷流的剪切层演化及激波胞格与湍流间的相互作用构成欠膨胀超声速喷流的主要噪声源。     

新型高低速流统一算法及其叶轮机应用

 对于可压缩/不可压缩混合流动或高低速混合流动问题，目前没有较好的解决方法。基于Roe格式及预处理方法，提出了全速度Roe格式。全速度Roe格式是将Roe格式与低速Roe格式通过一个与本地马赫数相关的函数联合而得到的。与标准Roe格式相比，对亚音速流动，全速度Roe格式中声速对数值耗散的影响随着马赫数变化而变化。与预处理Roe格式相比，全速度Roe格式简单易于实现，而且没有计算稳定性的限制，不需要采用全局截断处理，从而在理论上对高低速流动具有一致精度。事实上，全速度Roe格式意味着可以通过直接改变Roe格式的非线性特征值，解决可压不可压流动统一计算精度的问题。这是对统一算法的新观点。数值模拟表明,全速度Roe格式具有与预处理Roe格式一样的激波捕获能力、极低速流动准确模拟能力。     

高阶气体动理学格式在湍流数值模拟中的应用

本文回顾了高阶气体动理学格式在湍流数值模拟中的应用。与传统的Riemann求解器相比,气体动理学格式可以提供时空耦合的演化过程,这对发展高精度格式十分重要。因此,基于两步四阶时间离散和高精度WENO重构,发展了具有四阶时间精度的气体动理学格式。该格式有更高的数值精度和稳定性,并且具有更好的处理复杂流动问题的能力。目前,两步四阶格式已经成功地应用到低雷诺数湍流直接数值模拟和高雷诺数工程湍流RANS模拟中,包括低速槽道湍流、超声速均匀各向同性衰减湍流、二维亚声速翼型湍流和三维跨声速翼身湍流等。数值结果表明该格式对湍流直接数值模拟和湍流RANS模拟具有高数值精度和高数值稳定性。下一步将利用高阶气体动理学格式研究更具有挑战性的可压缩湍流问题,例如超声速湍流边界层和激波边界相互作用等。     

基于非结构/混合网格的脱体涡模拟算法

为了提高二阶精度有限体积算法的湍流数值模拟能力,在原始Roe格式基础上建立了与脱体涡模拟(DES)方法相匹配的二阶混合耗散自适应格式,根据流场信息自动调节格式耗散,并分别基于Spalart-Allmaras一方程湍流模型和k-ω剪切应力输运(SST)两方程湍流模型,发展了基于非结构/混合网格的DES方法。采用该方法计算了雷诺数为3 900的圆柱绕流和NACA 0021翼型60°大迎角分离流两个典型算例,通过与试验数据以及其他数值结果的对比验证了该方法的可行性。同时开展了不同数值格式、湍流模型的对比分析,研究结果表明:采用混合格式的DES算法能够解析更小尺度的湍流涡结构、计算数据更接近试验值;本文的DES类算法受其基准湍流模型影响较小。     

LU-AUSMPW混合格式在可压缩无粘和粘性流动数值分析中的应用

本文首先将 AUSMPW格式与三阶 MUSCL格式融合, 给出了其在任意曲线坐标下的三维形式, 并将其与 LU-SGS格式结合, 应用于可压缩 Euler和 Navier-Stokes方程的求解。其次, 构造了一种新的通量限制器。最后, 为了验证 LU -AUSMPW混合格式的性能, 将平面叶栅跨音速无粘流动以及喷管超音速粘性流动作为算例。本文计算结果与文献计算结果和实验数据相符很好, 表明采用 LU-AUSMPW混合格式数值模拟可压缩流场, 具有较高的计算精度、较快的收敛速度和良好的稳定性。      

几种典型迎风格式的分析与比较

重点分析了AUSM 格式的构造方法,同时将其与Van Leer, Roe格式等构造方法进行比较和分析,通过求解一维激波管问题,以及二维和三维跨音速流场的计算比较了AUSM 格式,VanLeer格式,Roe格式和中心格式的数值模拟结果和网格适应性,以及在经典四步Runge-Kutta显式时间推进和LU-SGS隐式时间推进下上述四种格式的计算效率.获取了对于基础研究和工程应用很有意义的结论.     

非结构网格上新型的NND有限元格式

张涵信发展的ＮＮＤ差分格式是由中心差分格式,二阶迎风格式和一阶迎风格式混合组成的杂交型格式。众所周知,和中心差分格式相对应的是Ｇａｌｅｒｋｉｎ有限元格式。通过对中心型有限元格式加修正项的方法本文成功地构造出二阶迎风型有限元格式和一阶迎风型有限元格式。     

烧蚀外形转捩流动的数值模拟

本文采用雷诺平均Navier-Stokes方程和修正Baldwin-Lomax代数湍流模型,以Roe二阶流通量差分分裂格式进行离散,数值模拟了典型烧蚀外形的高超声速流场,重点讨论了流动转捩的判定及其与复杂流态的关联.结果表明:数值模拟出的球钝锥外形转捩点位置与工程估算结果符合较好;对工程估算中存在困难的奶头状凹陷烧蚀外形,数值模拟也可给出符合流场结构的转捩点位置.     

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

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

冲压加速器燃烧流场的数值模拟

本文利用ＥＮＮ格式,通过求解包含非平衡效应的完全Ｎ－Ｓ方程,数值模拟了超声速燃烧流动中很有兴趣的Ｈ２／Ａｉｒ斜爆轰波冲压加速器的绕流。     

MLP高阶重构格式应用

在有限体积框架下,利用MLP(multi-dimensional limiting process)系列重构格式结合HLLHLLC(Harten-Lax-van Leer with contact)近似黎曼求解器,同时引入激波探测函数进一步降低MLP在光滑流动区域的数值耗散,数值模拟了超声速前台阶流动、结冰翼型的非定常流动、高超声速双楔流动和DLR F6-WB跨声速流场,研究了MLP系列格式在可压缩复杂流场中的表现.结果表明:在多维空间中,MLP格式能够在如强斜激波与网格线不重合等复杂流场数值模拟中保持严格的流场单调性;具有和传统MUSCL(monotone upstream-centered schemes for conservation laws)格式类似的计算效率,可以实现5阶,甚至更高阶重构;数值耗散更低,捕获更准确的激波位置,对航空工程数值模拟具有重要意义.     

典型标模音爆的数值预测与分析

精确预测音爆对超声速民机的研制具有重要意义。主流的音爆强度预测方法分为两步,首先通过风洞试验或CFD方法得到近场音爆过压（Over-pressure）分布,再运用修正线化理论或非线性声学理论将近场过压传播至地面,最终获得地面音爆的声压信号。本文运用典型标模对当前音爆数值预测方法的精度进行了验证和确认。在近场音爆过压分布的数值预测方面,分别考察了超声速尖点构型前缘修形尺度、不同空间离散格式和无黏/有黏流动控制方程求解对近场过压计算结果的影响。远场音爆预测方面,以LM1021全机构型近场过压分布为输入,使用基于波形参数法的远场传播工具分别考察了不同离散格式和有/无黏性计算的近场过压分布差异对地面音爆结果的影响。算例结果表明,尖点构型近场音爆预测中进行几何修形是十分必要的,使用相对合理的过渡球半径可以保证近场音爆预测精度,过大的修形尺度会对激波形状、激波和膨胀波的峰值均产生较明显的影响;就近场波形而言,熵相容格式计算得到的结果与试验测量值吻合最好,但不同离散格式导致的近场预测波形差异对传播到远场的波形关键指标（主要是最大过压和上升时间）的影响很小;是否计入黏性对近场波形结果尽管仅有小幅的影响,但将近场信号传播到远场得到地面波形时,这些细微差异会在远场波形的音爆评价关键指标上表现出明显的区别。