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

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

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

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

高超声速边界层感受性问题高精度数值模拟

通过五阶WENO格式和六阶对称紧致格式以及三阶TVD R-K结合的方法,对存在强激波和小扰动相互干扰的高超声速边界层感受性问题进行研究.结果表明:此方法能够模拟边界层内不稳定波的产生和发展以及小扰动和激波相互干扰等现象,因此能广泛应用于含激波的感受性等问题的可压缩湍流直接数值模拟以及具有激波和边界层干扰等复杂流场的计算.     

用CSCM格式求解激波反射问题

本文对Euler方程采用CSCM格式,数值求解了激波反射问题,与目前流行的交替方向隐式(ADI)格式以及总变差递减(TVD)格式的计算结果进行了比较,并作了简要分析讨论。     

高分辨率格式在非定常无粘可压缩流动中的应用研究

本文发展了一种计算非定常无粘可压缩流动的高分辨率格式。通量计算采用AUSMPW+Rie-mann求解器,其中使用了基于单调性的压力加权函数f,并与Roe Riemann求解器进行比较。为提高精度,使用五阶WENO格式进行左右状态插值,且时间推进采用三阶TVD Runge-Kutta方法。用该方法对移动接触间断问题、Sod激波管问题、二维激波反射问题和双马赫反射问题进行了数值计算,数值结果表明基于五阶WENO插值的AUSMPW+格式有很高的激波及接触间断分辨能力,并比基于五阶WENO插值的Roe格式有更高的计算效率。     

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

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

高精度ENO格式在射流数值模拟中的应用

本文分析了ENO格式的特点，并应用于全NS方程的迁移项和压力顶，模拟了轴对称和三维射流。首先对轴对称自由射流为算例比较了高阶TVD格式和ENO格式，证实了ENO格式捕获复杂激波波系的能力；其次计算了四喷管湍射流干扰流场，给出了相应的计算结果；最后计算了四股射流流入圆管内的流场，给出了物理量沿壁面和截面的分布。     

一种通量加权型紧致格式

将通量加权的思想引入到紧致格式中,构造了一个传统方法与紧致格式混合组成的加权型差分格式.利用该格式与二阶TVD格式分别计算了一维方波、组合波问题和一些黎曼问题,如Sod问题和Shu问题,以及一维定常激波问题.计算结果的比较表明加权格式无论在捕捉各种间断,还是在分辨各种复杂波系上,都具有较大的优势,并与精确解非常吻合.      

管道中复杂马赫反射流场的高精度数值模拟

本文分析了一类高精度迎风TVD欧拉方程差分格式的精度,并用240×80个网格点进行马赫反射的计算,获得了目前实验和计算不易得到的两个激波三叉点的复杂波系结构。本文还研究了激波管上壁面对马赫反射的影响,并计算了轴对称扩张管道中激波的运动及激波与半球形障碍物的相互作用。计算结果分析指出,本文所采用的高精度迎风TVD格式与计算方法可以用来模拟具有复杂波系的工程问题研究。     

TVD格式计算一维Euler方程间断解存在的几个理论问题

本文讨论了TVD格式按本文第二节A段意义推广到拟线性方程组的几个理论问题。指出:1.对于Riemann间断分解,Euler方程的解和Riemann不变量都没有TVD性质,因此不能将TVD格式自然推广;2.TVD格式是提高分辨率的一种途径,它不能保证得到物理解,能得到物理解的TVD格式一定是耗散守恒格式,因此将TVD格式推广时,必须保证格式是耗散的,否则仍会得到非物理解。并给出了一阶迎风格式,TVD2,UNO2得到的若干非物理解算例;3.非物理解现象只是发生在具有剧烈膨胀的单波区的计算上,而某些TVD格式若会得到非物理解,则一定发生在特征值变号的单波区的计算上。     

无波动、无自由参数的耗散差分格式

本文通过对“模型问题”的研究发现,在激波附近差分解的虚假波动与差分方程的修正方程式的三阶色散项相关。如果适当地调节色散项的系数使之过激波时该系数按要求变号,则可完全抑制激波附近虚假的波动。 在以上发现的基础上,提出了一个有效的无自由参数、无波动的差分格式。文中证明,这个格式就是“TVD格式”,亦可视为格式的推广。 文中对许多二维、三维流动作了计算,结果表明,本文的格式应用方便且可得到满意的结果,对激波有较高的分辨能力。     

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

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

不同半径圆柱诱导CH4/空气预混燃烧的计算研究

本文研究不同半径圆柱诱导CH4/空气预混燃烧流场。采用保自由流5阶WENO格式求解贴体坐标变换后的多组分Euler方程,用基元反应模型描述CH4/空气燃烧。不同于标准WENO格式通量构造方法,该WENO格式数值通量由方程的解进行WENO插值得到,在曲线坐标系下具有保自由流性质。首先给出了保自由流WENO格式精度和保自由流的数值验证,然后计算圆柱诱导CH4/空气预混气燃烧流场,并考察不同半径圆柱的影响,给出燃烧流场压力与温度等值线和云图、压力和温度沿过驻点线分布。结果表明:在考核计算结果网格无关性基础上,该WENO格式可准确地捕捉激波和火焰阵面形状、激波和火焰面驻点距离,得到的计算结果和文献结果相符。增大圆柱半径,激波和火焰面被推向来流方向,激波和火焰面之间距离也减小。和TVD格式相比,5阶WENO格式采用四分之一的网格数可得到近似相同的计算结果。     

超声速主流中横向喷流干扰流场数值模拟

本文用二阶迎风 TVD格式求解层流 N- S方程 ,数值模拟了超声速主流与横向喷流相互作用复杂干扰流场 ,得到了合理的涡系结构 ,计算波系结构及沿物面压力分布也与实验结果一致。计算还表明 ,TVD格式不但具有较强的捕捉激波的能力 ,而且用来模拟复杂涡系结构也是有效的     

收敛-扩张喷管中运用激波矢量控制法的计算研究

采用二阶TVD格式的有限体积法耦合RNG k-ε湍流模型和非平衡壁面函数求解二维守恒型雷诺平均N-S方程，数值模拟了运用激波矢量控制法的二维收敛一扩张喷管中的流动现象，很好的模拟出了由于次流喷射所产生的激波及上下壁面的压力分布，并将计算结果与国外的实验结果进行了对比。     

全隐式TVD格式求解二维Euler方程的对角化算法

本文采用有限体积法,Chakravarthy提出的基于近似Riemann解的全隐式TVD差分格式求解二维非定常Euler方程。采用牛顿法对半离散化后的全隐式差分方程进行线性化,然后对线性化后的方程采用近似因式分解法处理,得到了两组三对角块矩阵方程组。在求解该方程组时,采用了矩阵对角化思想使原来的4×4三对角块矩阵方程组转化为4组分离的三对角代数方程组,因而大大节省了求解方程组所耗费的计算时间,提高了计算效率。通过对翼型的跨声速绕流问题及圆柱超声速流动问题的计算,证实了本文算法不但简便可靠,而且还具有较强的通用性。采用TVD格式捕捉到的激波分辨率高,上、下游无任何波动。     

激波捕捉法中信号传输的监控和检测

依据本文引入的模拟流场中信息传输的监测因子θ,能十分容易地判定流场的激波。为了提高捕捉激波的质量,文中提出了通过选择合适的限控函数R(θ)构造高精度TVD差分格式的原理及利用由R(θ)提供的信息,用最小二乘法确定激波位置的设想。     

高精度差分求解气动方程的几个问题

本文探讨了发展高精度格式的必要性，研究了高精度格式与网格以及熵增条件的关系，并发展一种半离散化的空间为三阶精度的格式。从此半离散化的格式出发，可建立多步显式格式及隐式格式。模型问题的计算表明，该三阶精度的格式具有好的精度，且激波附近基本上没有波动。     

双重加权实质无波动激波捕捉格式

本文构造了一种新的加权激波捕捉格式，命名为双重加权实质无波动(Double Weighted Essentially Non—Oscillatory，简称DWENO)格式。DWENO格式不需要引进最优权值，并且可以递推地达到任意高阶精度。文中给出了典型算例的计算结果，这些结果表明DWENO格式具有高于WENO格式的间断分辨率。     

非结构网格上的TV-van Leer混合迎风格式

为了改进迎风型对流通量格式的激波捕捉计算稳定性,对迎风格式产生数值激波不稳定现象的原因进行分析,并且在分析结果基础上构造新型迎风型格式。为此,首先设计了数值算例,分析非结构网格上迎风格式出现数值激波不稳定现象的主要影响因素。结果表明,在一定条件下,网格和流动状态共同作用导致产生不合理的数值结果,并影响计算的稳定。其中,网格间的动量通量对数值激波不稳定现象影响显著。因此,基于TV格式构造了一种自适应调节动量通量耗散项的混合格式,从而在激波附近引入耗散,抑制扰动的产生和传播。数值计算结果显示新型混合格式在非结构网格离散条件下能够稳定地捕捉激波,同时混合格式保持了TV格式原有的黏性流动计算分辨率。