二维Euler方程组的高分辨率隐式差分格式的并行计算

本文首先对标量双曲型守恒律构造了一类新的二阶隐式ＴＶＤ差分格式，再利用Ｒｏｅ方法和隐式近似因子分解法推广格式到一维、二维守恒率方程组。为减少计算工作量，推导了二维Ｅｕｌｅｒ方程组隐格式的对角形式。最后，针对江南－Ⅱ并行机设计了适合于多处理机结构特征的并行计算程序，计算翼型的跨声速流场，结果令人满意。     

用守恒型差分方法计算底部分离流场

本文用守恒型显式格式及交替方向隐式迭代计算二维底部流场。守恒显式格式是从全守恒的Navier-Stokes方程组出发,在Allen—Cheng格式的基础上引进了“调整因子”。这对放大稳定性要求的步长,加速收敛均有益处。本文还用交替方向隐式迭代计算了同样的流场。结果表明只要边界条件逼近一致性较好,交替方向隐式迭代比显式格式收敛快,节约40％机时。     

用于二维无粘跨音速流计算的隐式TVD有限体积法(英文)

本文结合Harten的TVD(Total variation diminishing)概念和有限体公式建立一种求解守恒型Euler方程组的有效方法,采用通量分裂的迎风格式对方程组的隐式算子部分离散化,得到系数矩阵是主对角线占优的线性方程组,采用线Gauss-Seidel迭代推进法求解,提高了数值收敛速度,给出了二维跨音速内流的计算结果。     

二维Euler方程SCB格式的并行实现

本文讨论二维双曲型守恒律方程组的高分辨ＳＣＢ格式，并应用这些格式求解二维Ｅｕｌｅｒ方程组。计算过程中在江南－Ⅱ型并行计算机上进行的，计算过程运用了许多并行计算技术，计算结果令人满意。     

强迫耗散非线性发展方程与完全平方守恒格式

从描述大气和海洋运动的强迫耗散非线性发展方程出发，对强迫耗散非线性大气和海洋方程组显式差分格式的计算稳定性进行了分析，构造了一类强迫耗散非线性发展方程的显式准安全平方守恒差分格式。理论分析和数值试验证明，这类显式准完全平方守恒差分格式是计算稳定的，值得推广应用。     

求解双曲守恒律初值问题的GAUSS格式

构造了一类计算双曲守恒律弱解的 3点二阶显式格式 ,这类格式在CFL条件数为 1的限制下为TVD格式 ,并将这类格式推广到方程组的情形 ,进行了数值试验 ,结果是满意的。     

NND格式的推广及在粘流计算中的应用

本文把ＮＮＤ格式推广到了特征变量，守恒变量和原始变量形式。通过引入通量差分分裂方法，由守恒变量形式的ＮＮＤ格式得到了Ｒｏｅ的近拟Ｒｉｅｍａｎｎ解法，并讨论了限制器的作用。在隐式ＮＮＤ格式中引入了一种分解算法，提高了计算效率。本文的数值计算表明，ＮＮＤ格式对复杂波系干扰，边界层分离与再附等流动现象均有较高的分辨率。     

叶轮机转子内流动的三维跨声Navier-Stokes方程组数值计算

应用有限体积法及Baldwin-Lomax代数紊流模型求解叶轮机转子内流动的三维可压缩非定常薄层Navier-Stokes方程组,计算一个绕X轴等转速旋转的转子内流场问题,用当地时间步长方法来加速收敛。采用隐式算法对对流通量进行求解,而扩散通量采用显式算法。采用高阶TVD下的Roe格式计算对流通量,采用中心差分格式计算扩散通量。关于守恒变量Q的非线性方程组,采用牛顿迭代方法求解,线化后的方程用对称的Gauss-Seidle方法求其渐近解。应用该算法模拟了转子NASARotor67流场,与试验做了大量的对比研究,包括性能、等马赫线以及总压和总温等等。在最高效率点附近计算与试验比较接近,但在靠近失速点时计算与试验略有差别。      

超声速二维湍流分离流动的数值解法

本文采用层流中引导简化N-S方程的方法,建立了简化的湍流运动方程组,分别用零方程、一方程模型计算了绕二维压缩拐角的湍流分离流动。数值求解用显隐式的反扩散格式。     

叶轮机转子内流动的三维跨声速Navier-Stokes方程组数值计算及试验验证

应用有限体积法及Baldwin-Lomax代数湍流模型求解叶轮机转子内流动的三维可压缩非定常薄层Navier-Stokes方程组,计算一个绕X轴等转速旋转的转子内流场问题,用当地时间步长方法来加速收敛.采用隐式算法对对流通量进行求解,而耗散通量采用显式算法.采用高阶TVD下的Roe格式计算对流通量,采用中心差分格式计算耗散通量.关于守恒变量Q的非线性方程组,采用牛顿迭代方法求解,线化后的方程用对称的Gauss-Seidle方法求其渐近解.应用该算法模拟了转子Rotor67流场,与试验做了大量的对比研究,包括性能、等马赫线以及总压和总温等.在最高效率点附近计算与试验比较接近,但在靠近失速点时计算与试验有较大的差别.     

用虚拟黏性法构造Navier-Stokes方程黏性项的隐式格式

提出的虚拟黏性法是一种关于Navier-Stokes方程黏性项(黏性应力项和传热项)的隐式新方法。通过引入虚拟时间和虚拟黏性项,将隐式格式的构造大大简化,从而避免了大型复杂隐式差分方程组的常规求解。在虚拟时间推进过程中所需求解的方程组的系数矩阵是一个三对角矩阵,它具有计算简单且计算量小的优点。用模型方程和Navier-Stokes方程进行了数值仿真,研究各种参数对计算的影响,并在精度和效率上与显式方法进行比较,证实了算法的正确性和优势。最后分析了本文格式的适用范围。     

三维Euler方程组隐式LU分解的高性能并行计算

许多非定常无粘流体力学问题的数值模拟都需要利用Euler方程组来进行计算,而由于在隐格式下,所选取的时间步长可以比在显格式下时大得多,所以隐格式越来越受到重视,其中隐式LU分解是最常用的方法之一。对三维Euler方程组,采用隐式LU分解进行计算时,网格点所在的各个对角阵面之间存在数据依赖关系,本文分析了采用区域分解且边界上用显格式代替隐格式进行计算的高效性,在长方体建筑物内的爆炸模拟表明,在有112个CPU的某MPP巨型机上,并行计算效率超过60%。本文还分析了计算结果与串行计算时的差异,以及利用区域重叠减小这种差异的方法,同时考虑了对处理器进行合理的逻辑组织,将计算网格映射到处理器网格,以最大限度减少通信开销的方法。文中最后以一个爆炸毁伤的例子实际说明了所述方法的可行性与高效性。     

气动计算中一个简单有效的差分格式

本文采用算子附加修正方法,提出了一个气动计算中简单而有效的差分格式。此差分格式的特点是保持原显式格式的精度,计算过程有显式格式的简单,又有隐式格式的特点,即对于任意步长此格式都是稳定的。本文利用此格式求解了N-S方程,给出了二维激波与边界层干扰的分离流计算结果,并与MacCormack的隐式格式进行了比较。计算结果表明,本文提出的格式与MacCormack隐式格式相比,大大简化了计算过程,节省了计算机时。     

扩散抛物化(DP)NS方程组的意义及其在计算流体力学中的应用

本文首先讨论扩散抛物化（DP）NS方程组的早期研究工作：它的提出、数学性质、意义和在CFD中的应用，然后讨论扩散抛物化理论的一些新发展。这些新发展是对NS方程组数值计算进行物理分析的基础上得到的，其中包括NS方程组差分计算时，粘性剪切流对网格间距和格式精度的要求；粘性项只保留剪切项的广义扩散抛物化（GDP）NS方程组，它的性质和应用。由于高Re数流动在NS方程组的差分计算中，网格Re数彼此相差悬殊的特点，产生了计算离散单元守恒方程组的新的算法思路，即离散流体力学（DFD）算法。在DFD算法中需要同时计算三种不同的守恒方程组（Euler,DPNS和NS方程组）。本文讨论了DFD算法中需要同时计算三种不同的守恒方程组（Euler,DPNS和NS方程组）。本文讨论了DFD格式的构造、它的优点和应用。并以超声速绕前后台阶流动为算例，来说明GDPNS方程组的用处和DFD算法的优点。DPNS方程组、GDPNS方程组、DFD算法是高智提出的，对这些问题，他的合作者从理论、算法、数值验证和某些应用又作了系统的研究。     

超燃冲压发动机三维喷流干扰流场的数值研究

用隐式有限体积法求解三维N－S方程数值研究超燃冲压发动机三维平行、垂直喷流干扰流场。计算中采用了具有当地极值递减（LED）性质的JST和SLIP中心差分格式。为了减缓计算网格数很多时对计算机内存的压力以及提高计算效率，运用了通量守恒的分区计算方法。     

一个基于通量分裂的高精度MmB差分格式

本文研究双曲型守恒律的高精度差分方法.在新的算法中,首先将计算区间划分为互不相交的小区间,再根据精度要求等分小区间;其次,根据流动方向进行通量分裂,重构小区间交界面上的正、负数值通量,并进行校正;然后,采用高阶Runge-Kutta TVD方法进行时间离散,构造了一维非线性双曲型守恒律方程的一个高精度、高分辨率的守恒型差分格式.推广到二维双曲型守恒律方程,证明了格式的MmB特性.进而推广到二维守恒型方程组情形.最后对二维Burgers方程及Euler方程进行了数值试验,数值结果令人满意.     

欧拉方程隐式迭代矢通量分裂跨音流场计算

 本文用矢通量分裂格式的有限体积法求解时间相关欧拉方程组。文中证明,矢通量分裂的有限体积法不能再采用物面反射边界条件。二维问题用隐式迭代法解代数方程组,第三维方向上则用强隐格式往返扫描。此法CFL数可达100的数量级,收敛快,占内存少。跨音流场计算结果与实验符合良好,激波只占一两个网格。     

平面叶栅跨音流动的数值计算及变工况性能的研究

提出一种隐式矢通量分裂差分格式并用来直接求解Ｒｅｙｎｏｌｄｓ平均Ｎ－Ｓ方程组。该方法避开了近似因子分解及矩阵运算，具有精度高、稳定性好、计算量少等优点。在平面叶栅跨音流场的计算中，较好地捕获了激波，与实验比较，结果令人满意。     

求解双曲型守恒律的五阶松弛格式

给出了一种求解一维双曲型守恒律的五阶松弛格式.该格式以五阶WENO重构和显隐式Runge-Kutta方法为基础.本文格式保持了松弛格式简单的优点,即不用Riemann解算器和计算非线性通量函数的雅可比矩阵.用该格式对一维Euler方程进行了数值试验,并与三阶和四阶松弛格式的计算结果进行了比较,结果表明本文的格式具有更低的数值耗散和更高的分辨率.     

基于HDCS-E8T7格式的气动噪声数值模拟方法研究进展

发展了气动噪声高精度数值模拟方法,空间导数离散采用了HDCS-E8T7格式及精度与之匹配的边界格式,时间离散方法采用了高精度多步龙格库塔方法和高精度隐式双时间步方法,发展了高精度对接边界算法,将方法推广至多块对接网格以满足解决复杂几何构型问题的需要,采用隐式大涡的概念处理可能出现的湍流问题。在此基础上,研究了几何守恒律对计算结果的影响,展示了复杂网格中高精度计算满足几何守恒律的重要性,完成了等熵涡、双圆柱散射、串列柱翼构型和喷嘴射流等典型噪声问题的求解,所得计算结果展示了所发展的模拟方法具有良好的预测精度和解决复杂构型气动噪声问题的潜力。