预处理方法在化学非平衡流场数值模拟中的应用

将预处理方法发展应用于化学非平衡流场的数值模拟中。采用有限差分,LU-SGS隐式方法对二维可压化学非平衡控制方程进行了耦合求解,对化学源项和无粘项隐式处理,粘性项则采用显式。通过化学反应剪切层和激波点火的算例表明,所采用的数值方法能够有效求解化学非平衡流场,证明了本算法的可行性,为下一步的工程应用奠定了基础。     

非结构同位网格SIMPLE类算法收敛性能比较

以非结构网格的SIMPLE算法为基础,将算法扩展为SIMPLEC算法.利用30°角的斜方腔流动计算成果,分析了非结构同位网格的SIMPLE/SIMPLEC算法的收敛性能;比较了因网格的非正交而引入的非正交项的取舍对该算法收敛性能的影响;并采用显式校正步法对SIMPLEC算法进行了显式校正.比较表明,在非结构同位网格SIMPLEC算法中可忽略非正交项,但有必要对压力作亚松弛.显式校正步法可显著地加速在非结构网格上求解N-S方程的收敛性能,而且在不同的松弛因子组合下,均有较好的收敛速率.     

柔性扑翼非定常流场的数值计算方法

提出一种将Delaunay图映射网格变形技术和非结构嵌套网格方法结合使用的策略,解决网格变形和嵌套网格单独用于柔性扑翼流场计算时需要网格再生的问题。该方法为嵌套网格中的每个嵌于背景网格的贴体非结构网格生成Delaunay背景图;每个时间步,根据扑翼的运动和变形规律移动背景图,再根据网格点和背景图的映射关系移动网格点,之后自动完成嵌套边界的定义和插值关系的建立。为方便嵌套关系的建立,嵌套网格进行分层管理。也研究了一种内存消耗少、效率较高的搜索算法,以及格心格式和格点格式统一的边界拓宽算法。非定常可压缩Navier-Stokes方程在非结构的动态网格上用有限体积法离散,并用预处理的双时间步推进、隐式LU-SGS迭代求解。几个扑翼算例的结果表明,该方法充分利用了Delaunay图映射网格变形方法的高效率,同时也发挥了嵌套网格处理大幅运动的优势;用于既有整体大幅扑动又有局部小变形的柔性扑翼流场计算,可取得令人满意的精度和效率。     

栅格翼导弹流场混合网格N-S方程数值计算

采用结构和非结构混合网格技术,对栅格翼导弹黏性流场进行数值模拟,预测其气动特性.计算的马赫数为0.7～2.5.用有限体积法和LU-SGS算法求解N-S方程.计算与实验数据以及栅格翼附近流场结构进行了比较.     

基于非结构网格流场计算的网格重排序

提出基于非结构网格流场计算时的网格重排序,优化了网格信息在内存中的存贮位置,使得计算过程中调用数据的处理速度加快,有效提高了计算效率。以M6机翼及DLR-F4翼身组合体的跨音速无粘流场为算例,分别通过显式四步龙格-库塔推进格式和隐式LU-SGS推进格式来验证网格重排序后对计算效率的提高作用。原始网格均用Delaunay方法生长,并分别使用了两种重排序方法进行比较。计算结果表明:重排序后的计算效率可以提高到未排序时的3到5倍,很具实用价值。     

高阶DG格式多重网格方法构造中的不同隐式策略影响研究

DG方法是一种非常具有潜力的高精度方法,但其在对复杂外形的数值模拟方面仍存在内存需求量大、计算量巨大等不足.为了进一步提高DG方法求解Euler方程的效率,在传统p型多重网格的基础上,结合LU-SGS和GMRES两种隐式迭代方法,研究其整体加速性能.p型多重网格方法通过对不同阶次多项式近似解进行递归迭代求解,来达到加速收敛的目的.高阶近似(p＞0)使用显式龙格库塔格式,最低阶近似(p=0)使用隐式格式.对NACA0012翼型和ONERA M6机翼跨音速无粘流动进行数值模拟,结果表明:与显式TVD-RKDG时间格式相比,DG(p0)层上采用LU-SGS和GMRES的p型多重网格方法收敛速度均得到明显提高,且GMRES迭代法性能最佳,LU-SGS迭代法次之.     

基于动网格的非定常粘性空气动力学计算

提出了一种将阵面推进法和半扭转弹簧法相结合的动网格生成方法。运用阵面推进法决定附面层动网格的运动规律,半扭转弹簧法生成附面层外的非结构动网格。使得算法能够同时适用于非定常粘性与非粘性空气动力学计算。文中对于非定常可压缩N-S方程,在非结构动态网格上用有限体积法离散,并用预处理的双时间步推进、隐式LU-SGS迭代,计算了二维NACA0012翼型俯仰简谐运动的低雷诺数非定常空气动力,计算结果与参考文献提供的结果吻合良好。此外,本文还对三维1/4椭圆薄板的低雷诺数非定常空气动力进行了计算。     

动能BGK格式模拟高雷诺数定常流的收敛性研究

为克服动能BGK格式计算耗时多、收敛慢的缺点,在格式中引入当地时间步长、隐式LU-SGS方法和多重网格技术,基于RAE2822翼型粘性跨声速绕流的数值模拟,对其收敛性的改善进行了研究。结果表明:上述3种加速手段对求解定常问题时的收敛性有明显的改善,达到定常解所需的迭代步数随着CFL数的提高而减少;因而求解某类问题时,隐式格式比显式格式更具优势;多重网格既适用于显式格式,也适用于隐式格式,加速性更加突出。     

旋转坐标系下分区计算的LU隐式方法

将隐式时间迭代方法应用于并行计算是CFD研究中的热点问题,研究发现LU-SGS格式在旋转坐标系下对分区边界比较敏感,对旋翼悬停状态进行分区并行计算时碰到了计算发散的问题。针对此问题,对基于LU分解的3种隐式时间迭代格式（LU-SGS、DP-LUR和HLU-SGS）进行了对比研究,设计了静止流场下的旋转网格算例对计算方法进行测试。结果表明,LU-SGS格式在网格边界采用简化处理方法,当边界的逆变速度增大时会引起数值误差放大,误差的不断积累导致计算发散。DP-LUR和HLU-SGS格式通过在边界单元采用雅克比迭代算法,能有效消除分区边界影响,使计算格式在大CFL数条件下保持稳定,其中HLU-SGS继承了LU-SGS迭代效率高的特点。在此基础上,采用并行化的LU隐式方法对Caradonna-Tung旋翼进行了并行计算,针对两套稀密程度不同的背景网格,在产生复杂分区边界的计算网格条件下,均获得了气动力和旋翼尾迹都充分收敛的流场。通过计算,分析了背景网格密度对预测气动力性能和捕捉尾迹流场的影响,计算结果与试验值吻合良好,验证了当前隐式计算方法适用于旋转坐标系并行计算,适合于推广至其他大规模并行分区的流场计算。     

运动边界非定常全速流场的数值模拟

将预处理引入到ALE格式的非定常可压缩N-S方程中,发展了一种能有效模拟含运动边界的三维非定常全速流场的数值方法.控制方程在非结构动态网格上用双时间步推进并用隐式、LU-SGS迭代求解;提出了一种特征变量的边界条件以适合预处理后的特征系统.几个典型算例的结果表明,预处理方法将可压缩流方法的计算能力拓展到了低马赫数范围,有效地提高了收敛速度和计算精度.     

Revisiting the space-time gradient method: A time-clocking perspective,high order difference time discretization and comparison with the harmonic balance method

This paper revisits the Space-Time Gradient (STG) method which was developed for efficient analysis of unsteady flows due to rotor–stator interaction and presents the method from an alternative time-clocking perspective. The STG method requires reordering of blade passages according to their relative clocking positions with respect to blades of an adjacent blade row. As the space-clocking is linked to an equivalent time-clocking, the passage reordering can be performed according to the alternative time-clocking. With the time-clocking perspective, unsteady flow solutions from different passages of the same blade row are mapped to flow solutions of the same passage at different time instants or phase angles. Accordingly, the time derivative of the unsteady flow equation is discretized in time directly, which is more natural than transforming the time derivative to a spatial one as with the original STG method. To improve the solution accuracy, a ninth order difference scheme has been investigated for discretizing the time derivative. To achieve a stable solution for the high order scheme, the implicit solution method of Lower-Upper Symmetric Gauss-Seidel/Gauss-Seidel (LU-SGS/GS) has been employed. The NASA Stage 35 and its blade-count-reduced variant are used to demonstrate the validity of the time-clocking based passage reordering and the advantages of the high order difference scheme for the STG method. Results from an existing harmonic balance flow solver are also provided to contrast the two methods in terms of solution stability and computational cost.     

非结构混合网格高雷诺数流动计算中一种改进的LU-SGS隐式格式(英文)

The lower-upper symmetric Gauss-Seidel (LU-SGS) implicit relaxation has been widely used because it has the merits of less dependency on grid topology,low numerical complexity and modest memory require...     

预处理法求解定常/非定常混合网格的全速流场

运用时间导数预处理法在混合网格上求解三维定常/非定常N-S方程。预处理后的N-S方程,用二阶迎风格式的有限体积法离散,隐式、LU-SGS迭代求解,双时间步推进求解非定常流。本文以一维预处理欧拉方程为基础进行分析,提出了一种特征变量的边界条件以适应预处理后的特征系统。从极低速到跨声速范围的众多算例表明,本文的边界条件处理稳定有效,预处理法在各种速度下都具备同样快速的收敛性,在较宽的使用范围内均能得到理想的计算结果。     

一种高效的隐式间断Galerkin方法研究

基于线性化处理,在时间方向上对间断Galerkin方程进行了隐式离散,从整体上对迭代过程进行了合理的优化,并以此求解了计算流体力学中的二维Euler方程。其中,LU-SGS方法得到了进一步的推广,被用来高效求解隐式格式对应的大型稀疏线性系统。数值实验表明,无论对于亚声速问题还是跨声速问题,该格式都是无条件稳定的;与显式的Runge-Kutta间断Galerkin格式相比,当残值下降到相同量级时,隐式格式所需的迭代步数和CPU时间均在很大程度上得到了减少。     

高阶加权紧致非线性格式（WCNS）在二维流动计算中的加速收敛研究

以二维圆柱超声速无粘和粘性绕流的数值模拟为例,空间项采用高阶精度格式WCNS离散,对比研究了LU-SGS、高斯-赛德尔点松弛、线松弛以及GMRES等隐式求解方法的收敛性,并对右端项中的无粘通矢量、GMRES方法中的预处理和子迭代等影响作了对比计算.结果表明,右端项采用Steger-Warming无粘通矢量的收敛性优于其他通失量方法,采用了准确的解析雅克比矩阵的点、线松弛的收敛速度优于LU-SGS,以线松弛为预处理的GMRES算法具有良好的收敛特性.     

高超声速压缩拐角湍流计算的收敛性分析

应用基于shear stress transport(SST)湍流模型的隐式紧耦合方法对高超声速压缩拐角流动进行了数值计算.考察了差分格式、限制器以及湍流时间步长等7个因素对湍流计算过程稳定性和收敛性的影响.结果表明:差分格式和限制器对收敛性有很大影响,总体而言耗散性越大的差分格式和限制器收敛性越好;增加湍流时间步长对收敛性影响很小;适当降低湍流变量最低限制值有利于湍流方程的收敛;增大Courant-Friedrichs-Lewy(CFL)数和隐式内迭代次数可显著加速收敛,但均存在最优值,超过最优值则加速效果不明显.      

用隐式多重网格法计算实用外形无粘绕流

网格展弦比的好坏严重地影响多重网格的效率。文中细致地处理了多重网格算法的细节,了网格展弦比的影响,发展了一种计算实用外形绕流的隐式多重网格解法,计算效率提高了５倍。     

一种适用于三维混合网格的GMRES加速收敛新方法

为提高流场计算收敛效率,发展了一套适用于三维混合网格Naiver-Stokes方程求解的并行广义最小残差（GMRES）隐式时间推进方法。该方法由科学计算可移植扩展工具包（PETSc）中的Krylov子空间求解器实现,线性方程系统中的系数矩阵直接以显式给出以提高算法的稳定性。为进一步提高GMRES方法的收敛速度,对非结构网格的序号进行了重排序,使得系数矩阵的非零元素尽量向主对角线靠近。利用所发展的GMRES方法,完成了对ONERA-M6机翼、AIAA阻力预测会议通用研究模型（CRM）等算例的计算,计算结果与试验结果吻合良好。通过与其他隐式推进方法进行比较,对算法的收敛特性进行了研究。结果表明,所发展的GMRES方法计算更加稳定,残差下降速度相对LU-SGS（Lower-Upper Symmetric Gauss-Seidel）方法更快,尤其是气动力系数向着收敛解逼近的速度更加明显,提高了计算效率。