高阶精度DG/FV混合方法在二维粘性流动模拟中的推广

DG/FV混合方法因其具有紧致性、易于推广至高阶及相比同阶DGM计算量、存储量小等优点,已成功应用于一维/二维标量方程和Euler方程的求解。在此基础上,将该方法推广于二维三角形/矩形混合网格上的NavierStokes方程数值模拟,将格式形式精度提高至4~5阶。物理量的空间重构及离散使用DG/FV混合重构方法;无粘通量计算采用Roe格式;粘性通量计算采用BR2格式;时间方向离散采用高阶显式R-K方法或隐式方法。利用该方法计算了有解析解的Couette流动问题以验证几种格式的数值精度阶,并计算了层流平板流动和定常、非定常圆柱绕流问题等经典算例。计算结果表明DG/FV混合方法达到了设计的精度阶,在较粗的网格上亦能得到高精度的计算结果;定性分析和数值结果表明相比同阶DG方法单步计算量减少约40%。     

二元喷管无粘流的有限体积法分析

本介绍了一种计算收敛－扩散喷管内流特性的有限体积法，它将求解域分成许多有限体积单元，将欧拉方程在小体积单元上积分，从而得到半离散方程，然后采用Mac Cormack格式进行时间推进，求得定常解。计算结果与试验值吻合良好，是一种有效的计算具有喉道小曲率半径不连续的收敛－扩散喷管特性的方法。     

矢通量分裂显式格式及二元喷管射流模拟

本文在较常用的显式两步格式的基础上,采用 Steger- Warming形式的矢通量分裂 (简称矢分或 FVS)的原理,提出了一种矢分二阶精度显式两步格式。在分析了格式的各种特性后,与Mac Cormack格式进行了实算对比。又给出了某二元喷管射流问题的模拟实例。计算结果表明,矢分与该显式格式的配合具有矢分法的基本优点和计算稳定、速度较快的优点。      

双曲守恒律方程的一种熵相容格式

提出了一种求解双曲守恒律方程的熵相容数值通量。在熵守恒通量中添加一个二阶迎风项和一个三阶的差商项来保持熵稳定并且抵消解在跨过激波时所产生的激波强度立方倍的熵增,从而实现熵相容。新的数值通量能精确保持定常的接触间断、消除非物理的膨胀激波及负压力等现象。通过采用近年发展起来的WENO方法在单元交界面处进行高阶重构,得到高阶精度的熵相容格式。数值算例采用空间半离散格式,并结合显式三步三阶Runge-Kutta(RK3)方法进行时间推进。不同的算例结果表明,格式具有稳定性、高分辨率和无振荡性等特点。     

不可压缩流动的高精度非标准格子玻尔兹曼方法

首先回顾了高精度非标准格子玻尔兹曼方法的发展历程,基于高精度通量重构格式,发展了一种通量重构格子玻尔兹曼方法（FRLBM）。采用两种求解方法：一种是将碰撞项隐式处理,直接求解离散速度玻尔兹曼方程（直接法）;另一种是先执行碰撞步,再求解纯对流方程（分步法）。通过收敛性研究,比较了这两种方法的精度和稳定性。研究结果表明,在小时间步长下两种方法误差近似,都能取得高阶精度;然而当时间步长增大,直接法误差几乎不变,分步法误差出现明显上升。由此表明,当取得近似误差时,直接法可以采用较大时间步长,计算效率更高,且直接法的稳定性略占优。接着通过模拟顶盖驱动方腔流验证了FRLBM捕捉流场细节的能力,并且比较了基于半隐格式的显式方法和一阶、二阶及三阶隐式-显式Runge-Kutta格式的时间离散,在不同雷诺数下的最大允许Courant-Friedrichs-Lewy数,数值结果表明二阶隐式-显式Runge-Kutta格式效果最优。最后数值模拟了圆柱绕流,验证了FRLBM计算复杂外形绕流的可靠性。     

不同马赫数的无粘和粘性流动高阶精度隐式计算方法

 应用隐式时间推进法对不同马赫数的无粘和粘性流动进行数值分析,给出了基于预处理方法的高阶精度隐式求解方法。方程离散采用改进的高阶精度对流迎风分裂格式。该方法通过求解曲线坐标系可压缩Euler和 Navier-Stokes方程,对低速到超音速范围内的无粘和粘性流动的典型问题进行了数值分析。计算结果与文献的数值结果或实验数据对比分析表明,该方法对低速到超音速范围内的无粘和粘性流动问题进行数值分析是可行而有效的。     

用隐式WENO格式计算悬停旋翼跨声速流场

发展了一种基于高阶迎风格式的计算悬停旋翼跨声速流场的隐式有限体积法。对流项采用Roe通量差分分裂格式,使用五阶WENO格式进行左右状态重构,并与MUSCL插值进行比较;粘性项采用中心有限体积法。为提高收敛到定常解的效率,时间推进采用LU-SGS隐式方法。数值模拟采用了一种能够有效传递网格间流场信息的重叠网格,其中使用了三层内边界和贡献边界的方法以便插值的直接进行。用该方法对一跨声速悬停旋翼粘性流场进行了数值计算,数值结果表明:与MUSCL格式相比,WENO格式对激波位置捕捉得更准确,具有更强的涡捕捉能力,同时还表明了WENO格式在很大程度上能够克服涡耗散问题。     

捕捉间断的高精度数值方法

为发展适用于捕捉超声速流场中各种间断的高精度算法,将通量限制的思想引入到紧致格式中,构造了一个传统方法与紧致格式混合组成的通量限制型差分格式.通过在时间方向上利用一阶精度格式计算的一维定常激波,以及在时间方向采用多步Runge-Kutta方法计算的一维非定常激波管问题上的数值试验与二阶精度的TVD格式所计算的结果比较,表明新方法比二阶精度方法在间断的捕捉上具有明显的优势.通过新方法的计算结果与精确解的比较,表明新方法的准度也是非常令人满意的.     

超燃发动机燃烧室三维化学反应流场的数值模拟

采用ＬＵ隐式方法强耦合求解有限体积法离散的完全N-S方程及组份方程,数值模拟超燃发动机燃烧室内的化学反应流场。计算中采用了AUSM通量分裂格式,两方程及八方程的化学反应模型。计算得到了流场的结构与组份分布,并对两种化学反应模型的结果进行了比较与分析。     

高温富油燃气超声速燃烧数值模拟

发展了二维超燃流场的数值模拟计算程序 ,用来进行高温富油燃气超燃流场的数值计算。数值解法采用了 Mac Cormack预估 -校正两步显式格式求解矢通分裂形式的 N- S方程组 ,湍流模型采用了壁面率修正的 Baldwin- Lomax代数方程模型和用于剪切混合层的普朗特混合长度模型。在超燃计算中 ,应用 C2 H4- O2 两步反应的有限速率化学反应模型 ,模拟高温富油燃气超燃试验的试验状态。计算结果表明 ,在高温富油燃气超燃流场中 ,热量释放较少