一种MWENO谱体积（SVMWEN05）格式

构造一种新的谱体积（spectral volume）格式来求解双曲型守恒律，记为SVMWEN05（Spectral Volume method by Multi—Weighted Essentially Non-Oscillatory）格式。其主要思想是：第一步，将空间计算区域划分为一系列单元，称为谱体积，等分每个谱体积为一些子单元，称为控制体积（controlvolume）。第二步，在谱体积内部采用类似MWEN05（Multi-Weighted Essentially Non—Oscillatory）的格式进行重构，而在谱体积的边界处采用传统WEN05格式进行处理。第三步，利用Runge—KutmTVD离散方法对半离散格式进行时间离散，得到时空一致高精度全离散格式。最后．在文中给出几个经典数值算例用以验证本格式的计算能力。     

一种MWENO格式的构造和应用

在结构网格上提出了一种新的在混合模板集上守恒量插值的非线性权系数计算方法,构造了一种多重加权实质无波动激波捕捉格式,即所谓的MWENO(Multi-Weighted Essentially Non-Oscillatory)格式.与DWENO(Double Weighted Essentially Non-Oscillatory)格式相比非线性权系数的计算不用人为调节.随后文中给出了几个典型算例.通过MWENO格式和WENO格式的计算结果比较,表明MWENO格式具有较高的间断分辨率.     

双曲守恒律方程的一种熵相容格式

提出了一种求解双曲守恒律方程的熵相容数值通量。在熵守恒通量中添加一个二阶迎风项和一个三阶的差商项来保持熵稳定并且抵消解在跨过激波时所产生的激波强度立方倍的熵增,从而实现熵相容。新的数值通量能精确保持定常的接触间断、消除非物理的膨胀激波及负压力等现象。通过采用近年发展起来的WENO方法在单元交界面处进行高阶重构,得到高阶精度的熵相容格式。数值算例采用空间半离散格式,并结合显式三步三阶Runge-Kutta(RK3)方法进行时间推进。不同的算例结果表明,格式具有稳定性、高分辨率和无振荡性等特点。     

双曲型守恒律方程离散格式的若干特性述评

从计算数学、物理力学分析和计算经验相结合的观点分析一些对于工程计算有价值的离散格式,梳理出脉络,从改进和创新的观点评述了当代计算流体力学的若干问题.     

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

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

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

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

一种TVD方法在三维非结构网格中的应用

以三维非结构网格的显式有限体积法为基础,采用了一种TVD方法求解三维Euler方程,使用Roe通量来计算网格单元边界处的守恒量通量.为了验证方法的可行性,用该方法模拟三维爆炸问题,得出的结论是我们的方法可行,稳定且有效.     

多介质流动改进型Level-Set重新初始化方法

根据Level-Set函数在网格点处的值,预估出间断面的位置,并对间断面附近网格点上的LS值作相应调整,使其保持距离函数的性质,能够避免由于重新初始化造成间断面位置移动的问题.数值计算采用高精度WENO格式,对于多个可压缩流问题进行了数值模拟,表明了该方法的有效性.     

应用维数分裂方法推广MUSCL和WENO格式的若干问题

首先对有限差分法和有限体积法的差异进行了讨论,在已有文献观点的基础上补充了二者在边界条件处理及网格需求等方面存在差异的新论据。介绍了逐维推导的MUSCL和WENO格式计算控制体界面通量的过程,认为此类格式计算界面通量的方法直接应用于高斯积分型有限体积法不够严谨,从而得到了应用维数分裂方法构造的MUSCL格式和WENO格式不属于高斯积分型有限体积法的观点,“积分格式”这一定义更能准确反映这类格式的特点。此外,讨论了MUSCL格式和WENO格式在曲线坐标系下不能保证守恒的原因,并简单介绍了消除方法。     

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

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

求解双曲型守恒律的高精度,无波动样条逼近有限体积方法

本文从三次样条逼近出发，提出了求解流体力学双曲型守恒律的一种高精度、无波动的数值方法。针对样条函数的特性，本文提出了一种新的通量限制技术，使该方法在光滑区可以达到四阶精度，在流动参数的空间分布出现拐点或极值点时分别退化为二阶或一阶精度的格式。数值实验表明，该方法对流场中的激波和接触间断有很高的分辨率，优于二阶精度的ＴＶＤ格式。     

基于有限体积格式的自适应笛卡尔网格虚拟单元方法及其应用

基于自适应笛卡尔网格,结合虚拟单元方法处理浸入边界,发展了一种有限体积格式实现模式,数值模拟了二维无粘可压缩流动。本文发展的方法相对于传统的有限差分方法具有执行效率快,易于获得守恒性等优点。最后,利用一些标准算例和多物体流场的计算问题验证了算法的有效性。     

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

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

基于时间谱方法的飞行器动导数高效计算技术

针对动导数计算工程应用亟需的高效、高精度发展目标,提出求解非定常雷诺平均Navier-Stokes(RANS)方程的时间谱方法(TSM)的全隐格式,以改善采样点数较大时的数值稳定性,并将TSM离散推广到Menter剪切应力输运(SST)湍流模型,以提高TSM的工程实用性。将TSM应用于数值模拟NACA0015强迫振荡,所得计算结果与试验数据和双时间步(DTS)方法的计算结果均能较好地吻合,验证了TSM对周期运动的模拟能力。采用发展的TSM对高超声速HBS标模和超声速Finner标模进行计算,并分析研究攻角和马赫数对动导数的影响规律。结果表明:对于周期运动,具有与DTS方法相当的计算精度,但TSM的计算效率会随来流马赫数的增大而提高,其效率优势在高超声速范围时可达一个量级以上。