对流扩散方程定常非指数型松弛和超松弛半隐格式

考虑定常对流扩散方程的数值解，采用时间相关法，从迎风格式出发，得到半隐格式，并进一步提出松弛和超松弛改进半隐格式，以达到节省内存空间和提高定常解的收敛速度的目的。上述改进半隐格式是无条件稳定和单调的，当定常态时是相容、收敛的。上述格式为非指数型，可化为显格式，改进了过去指数型格式受计算机字长限制的弱点，文末，对一维非线性Ｂｕｒｇｅｒｓ方程做数值实验表明，本文提出的三点改进半隐格式适合非线性计算，且保持无条件稳定和单调的特性，并使收敛加快，精度提高。     

QUICK和乘方格式在顶盖驱动方腔流动数值计算中的比较

推导了非均分网络系统中QUICK格式附加源项的表达式，并用QUICK和乘方两种差分格式，分别对不同网格数下的二维、三维顶盖驱动方腔流动进行了数值计算。通过与文献中基准解的比较，考察了两种对流差分格式的数值预测性能。     

对流扩散方程定常非指数型松驰和超松驰半隐格式

考虑定常对流扩散方程的数值解，采用时间相关法，从迎风格式出发，得到半隐格式，并进一步提出松驰和超松驰改进半隐格式，以达到节省内存空间和提高定常的收敛速度的目的，上述改进半隐格是无条件稳定和单调的，当定常态时相容，收敛的，上述格式为非指数型，可化为显格式，改进了过去指数型格式受计算机字长限制的弱眯，末，对一维非线性Ｂｒｕｇｅｒｓ方程做数值实验表明，本提出的三点改进半隐格式适合非线性计算，且保持无     

有限空间年龄结构种群模型的数值分析

对一类有限空间年龄结构种群模型提出了一种数值方法．运用迎风格式和Euler向前差分格式进行求解．给出了格式的稳定性，得到了最大模误差估计，证明了数值解的非负性并给出了数值算例．     

一阶线性双曲方程组的隐-显迎风差分格式

采用显格式与隐格式交替使用的方法，针对一阶线性双曲方程组提出了一种隐－显迎风差分格式．它综合了隐格式与显格式的优点，具有稳定性好、计算简便的特性．数值计算结果表明，这种格式是实用的．     

求解双曲型方程的隐式迎风格式的构造

针对一维单个守恒律的初值问题研究了NND格式，在已有的数值方法的基础上，提出了一类隐式迎风格式，并证明了它是NND格式。数值实验结果显示，新格式具有较高的激波分辨能力，且是无波动的。     

海冰的数值模拟

在前人工工作基础上建立了海冰的数值模型，建立了动量方程的迭代差分格式及连续方程的分裂差分格式。数值试验表明所建模型及差分格式的可行性，所得结果较好地反映了冰的形成及演变过程。     

含源定常对流扩散方程的一种高精度差分格式

提出一维定常对流扩散方程的一种高精度差分格式。该格式呈现指数型，具有四阶精度，数值算例表明，该格式较其它格式具有更高精度。     

一种改进的MUSCL格式

本文主要对MUSCL格式进行了一些改进,得到了一种新的格式-Modify-MUSCL(简记为M-MUSCL),并通过对线性初边值问题、一维Burgers方程初边值问题、Sod Riemann问题的数值求解,对MUSCL格式和M-MUSCL格式进行了数值测试和定量的比较,发现M-MUSCL格式有明显的优势.     

线性传输方程的Entropy-Monotone格式

在守恒律方程组的数值计算中,线性间断的磨损问题备受关注.为了减少线性间断的磨损,针对线性传输方程,提出了Entropy-Monotone格式.该格式属于Godunov型格式,包括重构、发展和求网格平均3个步骤.与传统的Godunov格式不同,该格式同时计算数值解和数值熵,并通过它们构造分片常数的台阶函数.数值实验表明,此格式对线性间断的模拟非常有效.     

三维非定常粘性流动的数值方法及其应用

提出了适用于模拟非定常三维粘性流动问题的数值计算体系，包括高速收敛稳定性好的新的ＬＵ型隐式格式，捕捉弱间断面和滑移面的高精度、高分辨率的改良型ＭＵＳＣＬＴＶＤ格式。将基于此数值模型而发展起来的非定常粘性流动通用软件应用于叶轮机械气动热力学计算，为研究内流粘性损失、激波和非定常尾涡干涉、三维分离流动，尤其为探明分离泡内部的涡结构打下了基础。此软件亦适用于叶栅设计等方面的工程应用问题     

紧致差分格式的分辨率与精度的实例比较与讨论

通过比较先前建立的4阶最优紧致差分格式,以及传统的6阶和8阶紧致差分格式,来研究精度和分辨率之间的关系,主要比较了空间离散格式的有效波数范围、实际数值计算精度、以及对小尺度波动的模拟能力.数值试验结果表明:(1)这3种格式的计算精度都可以达到理论精度,并且此时精度越高,误差越小;(2)对于小尺度波动,最优4阶紧致格式比6阶和8阶紧致格式具有更高的分辨率;(3)对于行波问题,最优4阶紧致格式能够更加准确地模拟波动的传播行为.理论分析和数值算例的比较结果均表明,数值格式的精度和分辨率并不能互相替代,而是要根据计算问题的需要选择具有合适的精度和分辨率的数值格式.     

低耗散TVD格式的构造

目的　为进一步提高双曲型方程数值解的精度,减少数值耗散,同时保证稳定性,处理有强间断的流动问题; 方法　从总变差减小（ＴＶＤ）定义出发,分析保证数值解ＴＶＤ的条件,找出满足研究目的的可能性和途径; 结果　首次提出和证明双曲型方程数值解ＴＶＤ的充分必要条件,并据此简化导出适用于时间显示格式的ＴＶＤ充分条件,并据此构造了一种低耗散ＴＶＤ 格式; 结论　进行的数值实验证明,高分辨格式可以进一步减小数值耗散,提高数值解精度;     

双曲守恒律方程WENO格式的优化方法

Weighted Essentially Non-Oscillatroy(WENO)是求解双曲守恒律方程的高精度高分辨率数值格式．论文讨论了双曲守恒律方程WENO格式的一些优化策略，减少了非线性权的计算次数和特征分解的次数，通过数值算例证明了这些策略的可行性，并比较了优缺点．     

浓度梯度对瓦斯爆炸影响的数值模拟

针对煤矿巷道中存在瓦斯浓度梯度的问题,基于时间上的TVD Runge-Kutta格式,空间上的5阶加权本质无振荡(WENO)格式离散控制方程组,自主研发了高精度大规模的并行计算程序. 利用该程序模拟了煤矿巷道中的爆轰波的传播过程. 研究结果得出了纵向和横向浓度梯度对瓦斯爆炸的影响规律,模拟结果与理论分析结果基本吻合.      

三对角四阶紧致差分格式的优化和初步应用

 提高数值解的精度和分辨率,有助于更精确地求解日趋复杂的工程问题。本文依据差分格式的伪波数应该在尽可能大的波数范围内接近物理波数的思想,构造了满足四阶精度的具有高分辨率的三对角紧致差分格式。一方面,它可以与近些年发展的求解(循环)三对角方程组的高效算法相结合,以更高的分辨率、更小的计算量来计算一阶导数;另一方面,与传统格式相比,该格式的最大精确求解波数可以达到2.5761,大于传统格式的1.13097。因此,优化格式更适合模拟小尺度波动。数值计算结果表明：(1) 虽然优化格式仍然是四阶精度,但要比传统四阶紧致差分格式的计算误差小,尤其对于小尺度波动,优化格式的计算误差会更小;(2) 对于行波问题,优化格式能够更加准确地模拟波动的传播行为,其优势也更加明显。理论分析和数值算例的比较结果均表明,优化的紧致差分格式更适合求解小尺度波动问题。     

Riemann问题基于双激波近似的激波捕捉格式

Riem ann问题迭代求解的计算量很大 ,为了提高效率 ,减小计算量 ,提出一种 Riemann问题的近似解法。把Riem ann问题中的膨胀波看成是“膨胀激波”,认为左右状态之间存在由两个激波围成的相似解区域。该方法的优点是保留了 Riem ann问题的非线性特征 ,且保证左右波之间的自相似解区域的熵不减少。该方法在一维及多维无粘可压缩流动求解中取得了较好的结果。     

无结构三角形网格下运动界面追踪的WENO有限体积法

在处理运动界面追踪问题的流体体积函数(VOF)法的基础上,给出了一种无结构三角形网格下的高分辨率的运动界面捕捉方法.该方法采用高精度的加权本质无振荡(WENO)有限体积格式离散VOF函数的空间导数,采用三阶TVD Runge-kutta方法离散时间导数,采用Lax-Friedrichs通量作为数值流通量.数值试验结果表明,用该方法来进行旋转速度场和剪切速度场的运动界面追踪,可以得到与理论解非常一致的追踪结果.     

Low Dispersion and Boundedness Practical Discretization Scheme

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｄｅｓｉｒｅｆｏｒａｎｏｐｔｉｍａｌｄｉｓｃｒｅｔｉｚａｔｉｏｎｓｃｈｅｍｅｒｅｍａｉｎｓｏｎｅｏｆｔｈｅｃｒｉｔｉｃａｌｉｓｓｕｅｓｉｎｃｏｍｐｕｔａｔｉｏｎａｌｆｌｕｉｄｄｙｎａｍｉｃｓａｎｄｎｕｍｅｒｉｃａｌｈｅａｔｔｒａｎｓｆｅｒ[1,2].Ａｎｏｐｔｉｍａｌｓｃｈｅｍｅｐｏｓｓｅｓｓｅｓａｃｃｕｒａｃｙ,ｓｔａｂｉｌｉｔｙ:ｂｏｕｎｄｅｄｎｅｓｓ(ｎｏｕｎｐｈｙｓｉｃａｌｏｓｃｉｌｌａｔｉｏｎｓ)ａｎｄａｌｇｏｒｉｔｈｍｉｃｓｉｍｐｌｉｃｉｔｙ.Ｔｈｅｍａｉｎｃｏｎ…     

利用小波求解二维非线性双曲型守恒律

在C.W.Shu和S.Osher提出的基于点值的高阶精确基本无振荡格式的基础上，使用插值小波变换对函数进行多尺度分解，直接利用点值插值误差来探测具有奇异性的区域范围（比如包括激波的区域）。在这些奇性区域，利用代价昂贵的高阶基本无振荡格式计算单元边界的数值通量；而在解光滑的区域，则用廉价的多项式插值根据先前低分辨尺度上得到的值来计算出其数值散度，从而减少计算工作量。一些其它的技巧也被用来提高其计算效率。最后这种方法被应用到求解二维非线性双曲型守恒律，同时还可以很容易地通过在一维基础上计算数值散度的方式扩展到高维情形。