一种可压缩流动的高阶加权紧致非线性格式(WCNS)的加速收敛方法

以高超声速粘性绕流的数值模拟为例,研究LU-SGS、高斯-赛德尔点松弛、线松弛以及GMRES等隐式求解方法在空间项采用高阶精度格式WCNS离散时的收敛性,并对GMRES(generalized minimal residual)方法中的子迭代影响作了对比计算.结果表明,采用准确的解析雅克比矩阵的点、线松弛的收敛速度优于LU-SGS,以线松弛为预处理的GMRES算法具有良好的收敛特性.     

三维非结构网格Euler方程的LU-SGS算法及其改进

将LU-SGS隐式时间推进格式运用到非结构网格Euler方程的求解中,并对传统LU-SGS格式进行改进,结合网格重排序,发展了一套效率更高的三维Euler方程求解器.以M6机翼及超临界LANN机翼的跨音速无粘流场为算例,将改进LU-SGS格式与四步龙格-库塔显式格式及传统LU-SGS格式进行了比较.计算结果表明:所有格式的计算结果与实验结果都符合很好;传统LU-SGS格式计算效率为显式格式的3倍多,而改进LU-SGS格式计算效率为显式格式的7倍多.     

PCCSAP-3D程序压力场算法改进

介绍大型先进压水堆安全壳专用分析程序PCCSAP-3D计算采用的方法,引入GMRES(Generalized Minimal RESidual)方法改进该程序的压力场算法.使用GMRES算法的实用变形,并采用合适的预处理技术,比较GMRES算法和ML-ADI算法在求解压力方程时的收敛速度.结果表明,利用压力矩阵结构化和稀疏性的特点,采用预处理GMRES算法能够更快速地求解压力方程.当压力矩阵规模变大时,GMRES算法相对于ML-ADI方法能够节省更多的时间.     

不可压缩流基于块预处理的并行有限元计算

发展了一个模拟非定常不可压缩粘性流的并行有限元求解器,时间离散使用具有二阶精度的隐式中点格式,基于三维非结构四面体网格剖分,使用高阶混合有限元离散速度场（P2）和压力场（P1）.全离散格式产生的代数方程组是大型、稀疏、非对称和病态的,基于修正的压力对流扩散预处理（PCD）和精心设计的子问题迭代执行策略,采用预处理的GMRES迭代法来高效求解线性方程组.利用相同的子问题迭代策略,同时给出基于最小二乘交换子（LSC）预处理的并行效率对比.大量数值算例验证了算法的精度、可扩展性和可靠性.三维驱动方腔流模拟结果（Re=3200.0）清晰地显示了方腔流中主涡（PE）、下游二次涡（DSE）、上游二次涡（USE）、侧壁涡（EWV）和TGL涡的存在.     

一种隐式预处理方法及其在定常和非定常流动数值模拟中的应用

将Choi-Merkle矩阵预处理方法与LU-SGS隐式方法、双时间法以及多重网格方法结合,发展适用于绕飞行器定常和非定常粘性流动的高效隐式预处理计算方法和程序.介绍一种针对定常和非定常流动的LU-SGS隐式预处理方法的统一表述方法.在不改变流动解的前提下,对Navier-Stokes方程的伪时间导数项实施Choi-Merkle矩阵预处理,从而改善可压缩控制方程在低速情况下的系统刚性,使基于LU-SGS时间推进格式的数值模拟方法同时适用于从极低马赫数到可压缩范围内的数值模拟.对Jameson中心格式的人工粘性进行相应的修改,以提高低速流动的计算精度.翼型、机翼以及翼身组合体绕流的数值模拟研究表明,隐式预处理方法获得了很高的计算效率,可使马赫数0.1左右的低速流动计算时间减少50%以上;通过对现有可压缩计算程序进行小量改动,便可使其均匀覆盖整个低速流动范围,提高CFD程序在飞行器绕流数值模拟中的实用性.     

三维笛卡儿坐标系中Lagrange流体力学的显式相容有限元方法(英文)

将Caramana等人提出的相容算法思想和有限元方法相结合,提出三维笛卡儿坐标系中Lagrange流体力学的显式相容有限元方法.采用三线性六面体单元和交错网格进行空间离散,利用质量集中进行显式求解,无需求解线性代数方程组.时间离散可采用两步显式Runge-Kutta格式.用边人工粘性消除激波振荡,用子网格扰动压力抑制网格的非物理变形.给出若干标准算例.数值算例表明,该方法具有较高的计算精度和计算效率,同时具有很好的对称性和总能量守恒性,总能量计算误差为计算机浮点计算截断误差.     

一种可扩展的三维半导体器件并行数值模拟算法

在基于漂移-扩散模型的三维半导体器件数值模拟中,通过有限体积法进行数值离散,采用完全耦合的牛顿迭代求解非线性代数方程组,并使用基于代数多重网格预条件子的GMRES方法求解牛顿迭代中的线性方程组,构造一种稳健且高度可扩展的非结构四面体网格上求解半导体方程的并行算法.基于PHG平台实现该算法的并行计算程序,并对PN结和MOS场效应晶体管等问题进行了最大网格规模达到5亿单元、最大并行规模达到1 024进程的大规模数值模拟实验,结果表明,该算法计算效率高,可扩展性好.     

Arnoldi方法简述及其在流动稳定性中的应用

流动稳定性问题常常归结于巨型非对称矩阵特征值问题。多数求解巨型非对称矩阵特征问题的算法均是经基本的Arnoldi算法演化而来。首先简述基本的Arnoldi算法; 其次简述基于Arnoldi算法的几类变体, 如显式重启Arnoldi算法,隐式重启Arnoldi算法与多重隐式重启Arnoldi算法; 最后基于Arnoldi算法及其变体结合谱位移技术求解计算流动稳定性问题, 并通过数值实验比较可知结合谱位移技术的多重隐式重启Arnoldi算法的求解效率最高。      

低速流动的预处理多重网格求解

本文采用预处理的时间推进方法,求解可压形式的Euler/Navier-Stokes方程,数值模拟低马赫数范围内的无粘、粘性流动.为了进一步加速收敛,并使用了多重网格技术.计算表明预处理能够保持可压缩求解方法在高速时的收敛特征,对于低速流动计算可以加速收敛并提高精度.     

基于LU-SGS的非结构弹簧网格迭代算法

针对运动间断拟合中需频繁更新网格点位置的特点,提出一种基于LU-SGS(lower-upper symmetricGauss-Seidel)迭代方法的非结构弹簧网格运动算法.根据弹簧网格原理构建与网格拓扑关系相对应的稀疏系数矩阵,将LU-SGS思想成功引入动网格迭代算法,并辅以合理的网格运动管理策略,实现动网格的快速迭代.研究表明,在非结构网格下,LU-SGS算法可以满足运动间断拟合的需求,在流场隐式时间推进时,仍能保证获得稳定解;与传统的SOR方法相比,计算时耗减少20%以上.