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

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

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

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

关于NND格式的两点注记

对差分格式在激波上计算准确性的标准提出了新的看法,并证明ＮＮＤ格式完全符合这个标准。     

双曲守恒型方程的二阶摄动有限差分格式

对双曲守恒型方程，将其一阶迎风格式空间差商的常系数摄动展开为时间步长和空间步长的幂级数，通过确定幂级数系数而获得二阶精度的摄动有限差分(PFD)格式。进而从双曲守恒型方程的通量分裂型一阶迎风格式出发，通过类似的摄动展开方法，获得空间精度为二阶的通量分裂形式的摄动有限差分(FPFD)格式。这两类格式保留了一阶守恒迎风格式的简洁结构形式，使用三节点即可达到二阶精度，又避免了三点二阶格式的非物理数值振荡。并将这两类格式推广应用到双曲守恒型方程组，最后通过模型方程和一维激波管流动的数值算例验证了格式的高精度、高分辨率性质。     

预估校正的高分辨率格式及其应用

本文从基于通量分裂的一阶迎风格式出发,构造了一类校正型具有高分辨瘁的二阶全变差递减(TVD)格式,并用此格式进行了一维激波管问题及二维问题的计算。结果表明格式的解为二阶精度,激波处无明显的振荡,激波过渡区为1～2个点。     

差分计算中激波上、下游解出现波动的探讨

本文探讨了差分计算中激波上、下游解出现波动的原因。研究指出:在NS激波方程中,附加二阶耗散项,可平滑激波;但附加三阶弥散项或四阶耗散项,在一定条件下,可引起解在激波上、下游波动。因此在求解NS激波方程时,若采用二阶或三阶差分格式,由于格式弥散项和耗散项的存在,将引起解的波动。     

声速再入体表面热流数值模拟研究

本文采用计算效率较高的标量对角化隐式NND格式，通过求解Navier-Stokes方程对影响再入体表面热流计算准确的诸因素进行了综合分析。研究了Steger-Warming通矢量分裂、VanLeer通矢量分裂和通量差分分裂方法及相应熵修正方法对热流的分辨能力，并阐明了在物面边界上采用二阶中心格式、二阶中心和二阶迎风混合格式、以及一阶迎风格式等不同边界格式对热流计算的影响。在此基础上，采用通量差分形式NND格式对钝锥和钝双锥高超声速粘性绕流进行了数值模拟，计算给出了与试验结果相吻合的热流分布。     

NS方程激波计算的摄动有限差分方法

摄动有限差分(PFD)方法从一阶迎风差分格式出发,将差分系数展开为网格步长的幂级数,通过提高修正微分方程的逼近精度来获得更高精度的差分格式。由于格式基于一阶迎风格式,因此具有迎风效应、网格节点少等特点。本文首先通过对Burgers方程的摄动差分格式的推导,将摄动有限差分格式引入时间相关法的计算,并构造了守恒形式的摄动有限差分格式,然后推广到一维Navier-Stokes方程组的计算。数值比较研究表明:本文构造的NS方程摄动有限差分格式具有比一阶迎风较高的精度和分辨率,而且保持了一阶迎风格式的无振荡性质。     

基于5阶精度格式WCNS-E-5的p-multigrid方法研究

p-multigrid方法的基本思想是:在保证收敛结果为高阶精度的同时,利用低阶精度格式耗散大的特点,来改善高精度有限差分格式在迭代计算时收敛速度慢的弱点.本文基于5阶精度WCNS-E-5差分格式,引入1阶精度迎风格式和3阶精度加权格式,构造了p-multigrid方法,在迭代过程中采用了V循环、W循环、S循环、PreV和FMG循环等不同方式来应用这三种格式,并通过典型算例考察了这些循环方式对收敛速度的影响,初步数值试验表明,采用恰当的循环方式,本文所设计的p-multigrid方法能够加快收敛速度,并保证了最终收敛结果与5阶精度WCNS-E-5差分格式的一致性.     

高精度两步TVD格式的构造及数值检验

基于通量分裂的思想,利用Taylor级数理论和第二步限制器函数的控制,构造了在光滑区域空间为三阶、时间为一阶精度的高分辨率全变差递减(TVD)格式.该格式在激波过渡点降为一阶迎风格式.并从单个线形方程推广到非线性方程及方程组情形.通过几个典型算例的计算,并与二阶TVD格式作了比较,表明该方法对流场中的激波有较高的分辨率,且是无波动的.     

关于建立高阶差分格式的问题

为了能在不太密的网格上捕捉到流场的细致结构。通常采用高阶精度的差分 式进行数值模拟。为能抑制计算中和在激波附近产生的虚假波动,本文从物理构思出发,提出了建立高阶格式的两个基本原则,作者称之为抑制波动的原则和稳定性原则。     

高超声速再入体表面热流数值模拟研究

本文采用计算效率较高的标量对角化隐式NND格式，通过求解Navier-Stokes方程对影响再入体表面热流计算准确的诸因素进行了综合分析，研究了Steger-Warming通矢量分裂、Van Leer通矢量分裂和通量差分分裂方法及相应熵修正方法对热流的分辨能力，并阐明了在物面边界上采用二阶中心格式、二阶中心和二阶迎风混合格式、以及一阶迎风格式等不同边界格式对热流计算的影响。在此基础上，采用通量差分形式NND格式对钝锥和钝双锥高超声速粘性绕流进行了数值模拟，计算给出了与试验结果相吻合的热流分布。     

非定常数值模拟方法的发展及其在动态绕流中的应用

基于混合通量分裂的思想,通过应用Gauss-Seidel迭代求解差分方程,构造了一种时空二阶精度、无条件稳定的隐式迭代NND算法,并讨论了时间精度与稳定性,亚迭代收敛判则,几何守恒律的应用以及动壁边界条件等相关问题.通过引入一种简便易行的加权函数来综合刚性动网格生成技术和超限插值生成动网格这两种方法的优点,发展了一种计算量小、比较实用的加权动网格生成技术.作为应用实例,本文给出了多个动态物体绕流的数值模拟算例,计算结果表明了本文数值方法的成功.     

NND格式在多维理想磁流体方程组中的应用

采用修正的四步Runge-Kutta方法求解三维一般曲线坐标系下的理想磁流体方程组,为克服数值振荡,加特征型NND格式进行后处理.特征型NND格式推广到求解三维磁流体(MHD)问题需要知道雅可比通量的左右特征矩阵,在具体计算时需要克服矩阵的奇性.本文用三维程序采用推广的特征NND格式计算了一维MHD激波管和二维(MHD)喷管流动,计算结果表明,特征NND格式保持了TVD格式高精度的优点,又具有计算简单的特点,在包括强弱间断等复杂波系的定常和非定常MHD流场数值模拟中是成功的.     

跨音速压气机级三维紊流流场数值模拟

在对跨音速压气机级的动／静双叶排中的三维定常紊流流场数值的模拟中，利用三阶高分辨率ＮＮＤ格式和ＬＵ－ＳＧＳ隐式推进迭代法，既保证了流场中激波的模拟质量和粘性流动特征的正确预估，也实现了求解过程的高效率。对某单级跨音速压气机的计算及其与实验数据的对比证实了本文方法的有效性。将计算推广至多级轴流叶轮机械的情况将会是十分简便的     

高速粘性内流的高分辨率高精度迎风型杂交格式

本文在有限体积离散和ＬＵ分解的基础上，构造出一个新的隐式迎风型杂交格式，并用于求解定常流动的稳态解。     

高超声速滚转阻尼导数数值模拟

采用非定常Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程描述物体简谐振动流场,并在Ｅｔｋｉｎ理论下给出绕定轴转动时滚转阻尼导数的计算公式,定常流场的计算采用空间二阶精度的交替方向隐式分解的ＮＮＤ格式,非定常流场的计算采用时、空二阶精度的Ｒｕｎｇｅ－Ｋｕｔｔａ多步格式,采用代数和方法生成物体静、动网格。最后给出高超声速下钝体外形滚转阻尼导数的计算结果,以及滚转力矩系数随瞬时振幅的变化曲线。     

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

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