双曲型守恒律方程的一个大时间步长二阶TVD差分格式

本文根据Harten,A.的大时间步长差分分格式构造思想,为双曲型守恒律方程弱解计算构造了一个2K 3点大时间步长二阶显式差分格式——LTS-LF格式,得到了其在CFL限制K下为总变差不增差分格式(TVD格式)。文章按照Roe的方法推广格式到方程组情形,并就Burger’s方程和Euler方程组黎曼问题进行数值试验,结果令人满意。     

非线性双曲型守恒律的一类应用通量限制器的二步二阶格式

本文提出了计算非线性双曲型守恒律弱解的一类新的二步二阶精度有限差分格式,这类高精度非线性差分格式的构造应用了三点单调差分格式的数值通量和通量限制器。文章证明,在拟CFL条件限制下,格式为总变差衰减差分格式(TVD格式),因此不发生常用的诸如二步二阶Lax-Wendroff格式、Maccomark格式和其他古典二阶格式计算不连续解时在其周围出现的剧烈振荡现象。 为了减少Roe推广技术引起的计算复杂性,文中我们基于Davis,S.F.的思想,使用了Rusanov数值通量,把标量差分格式推广到方程组,构造了求解双曲型守恒律方程组的新的预测-修正格式,较之其他一般的预测-修正格式而言,仅增加了很少的计算复杂性,初步的数值试验表明,本文提出的格式无振荡和具有极好的激波分辨率。     

二维双曲型守恒律的一类MmB差分格式

基于通量分裂、单元平均分片线性重构及逆风特性进行空间离散,构造了二维标量非线性双曲型守恒律的一类新的二阶精度的半离散差分格式。进一步地利用二阶TVDRunge-Kutta离散方法对时间进行离散,得到了一类新的时空二阶精度的全离散差分格式,并证明了格式的MmB特性。之后,将格式按分量形式推广到二维非线性双曲型守恒方程组。该方法的一个主要优点是使用分量形式格式计算二维非线性双曲型守恒方程组,无须解黎曼问题且不用进行局部特征分解,因而形式简单、计算量小。通过计算二维可压缩流欧拉方程组的几个算例,数值结果表明,该格式具有高精度、高分辨率及计算简单的特点。     

求解一维非线性守恒律的一类MmB差分格式

本文构造了一维非线性双曲型守恒方程的一类ＭｍＢ差分格式，并推广到方程组情形。数值计算结果表明，格式具有高分辨率且不会产生数值为振荡。     

CFD问题显式有限差分方法的并行化

基于ＭＩＭＤ并行计算机模型，本文讨论了ＣＦＤ问题应用显式有限差分方法的并行化问题。利用区域分解法将计算问题分解为多个子问题，每个子问题由不同的处理器分别处理。针对ＳＣＢ计算格式的特点，尽量降低各处理器间的数据通讯，提高了并行计算效率。最后给出二维Ｅｕｌｅｒ方程组计算实例，计算结果令人满意，且算法有较好的可扩放性。     

二维复值Ginzburg-Landau方程的一个高阶紧致ADI差分格式

对二维复值金兹堡朗道(Ginzburg-Landau,GL)方程提出一个基于时间分裂的高阶紧致交替方向隐式有限差分格式。本文通过时间分裂法将GL方程分裂成一个非线性子问题及两个线性子问题,对非线性子问题以及其中一个线性子问题均通过精确积分进行计算,并对另一线性子问题构造紧致交替方向隐式差分格式进行数值计算。实际计算中,在每一时间步,利用追赶法求解一族常系数三对角线性代数方程组,从而使得算法既具有较高精度又拥有较快的计算速度。数值实验表明该算法在时间和空间方向分别具有二阶和四阶精度,并模拟了方程的一些动力学行为。     

弹尾超声速轴对称喷流的正格式数值模拟(英文)

采用二阶正格式方法对超声速轴对称喷流流动进行数值模拟。将二维守恒方程的正格式方法发展到轴对称Euler方程组的求解,并对导弹尾部超声速伴随射流进行了数值计算。数值结果与实验照片所反映的流动特征吻合较好,与三阶精度间断有限元方法计算结果吻合,相比于二阶高分辨率TVD格式的计算结果,正格式方法的计算结果在间断点处具有较大的梯度变化。这表明该方法对激波具有较强的捕捉能力,在间断处不会出现数值振荡,对弹尾声音速伴随射流的数值模拟真实、有效。     

轴对称进气道外罩跨音速非守恒型全位势流AF3算法

本文针对轴对称进气道外罩跨音速绕流问题的特点,采用了贴体坐标网格;并对这类网格拓扑构造了三种Baker的隐式近似因式分解AF3有限差分迭代算法。通过理论分析与数值实验,找到了既能计算亚临界情形,又能计算超临界情形的AF3格式;并就此格式讨论了适合于进气道绕流特性的AF3迭代格式中间变量的提法,具体提出了格式实施的相应方法,研制了可供分析任意轴对称进气道外罩跨音速特性的计算程序。由于采用了任意曲线坐标系上的轴对称非守恒型全位势方程适合于用不同方法生成的计算网格,因此,本文的工作还有可能进一步推广到计算更复杂的轴对称进气道跨音速绕流问题。计算表明,本方法收敛快,结果好。     

求解二维Hamilton-Jacobi方程的一类无波动的耗散差分格式

将二维Hamilton-Jacobi方程转换成双曲型守恒律方程组,然后利用求解一维H-J方程的一类无波动无自由参数的耗散差分格式的思想,构造了一类求解二维H-J方程的差分格式.数值实验结果表明:该格式易于计算且分辨率较高.     

二维非定常化学非平衡流动高精度数值模拟及其与实验的比较

本文成功地将二阶迎风 TVD 数值格式运用到弱电离、化学非平衡斜激波反射流场的数值计算中。气体力学方程与组元质量守恒方程组全耦合求解是目前非平衡流场数值计算中的新方向,文中详细给出了气体力学方程和组元质量守恒方程组全耦合求解过程,最后给出了三种不同情况斜激波反射流场数值结果并进行了详细讨论,本文所得结果与文献〔4〕提供的实验结果和文献〔1〕提供的数值结果相符合。     

高精度ENO格式的有效实现

基本无振荡（ＥＮＯ）有限差分格式是近几年出现的一类高精度差分格式，重构多项式的形成是ＥＮＯ格式构造的核心部分，一般均以网格平均值作为插值点。文中由通量值作为插值点，构造了一个简单的高精度格式，在求解高维欧拉方程时，较节省计算时间。引入三 阶ＴＶＤＲｎｇｅＫｕｔｔａ时间离散，就二维欧拉方程的激波反射与双马赫反射问题进行了数值试验，试验结果表明格式所具有的高精度和基本无振荡性质，能够满足实际计算的需要     

WEIGHTED COMPACT SCHEME FOR SHOCK CAPTURING

INTRODUCTIONRecently compact schemes have been widelyused in the simulation of complex flows,espe-cially in the directnumerical simulation of turbu-lent flows[1～ 2 ] . Standard finite differenceschemes have explicit forms and need to be atleastone pointwider than the desired approxima-tion order. It is also difficult to find suitable andstable boundary closure for high order schemes.Compared to the standard finite difference ap-proximations,the compact schemes can achievehigher order acc…     

激波运动与反射流场的ENO算法

在ＴＶＤ格式的基础上，以修正限制函数的方法，将通常的ＴＶＤ格式改造为性能更为优越的ＥＮＯ格式，并利用这种有限差分格式求解欧拉方程和薄层Ｎ－Ｓ方程。分别数值模拟了激波管中运动激波与三种楔形物体相互作用流场；给出了激波在不同物体上的反射及波涡干扰过程中非定常流场的变化情况；算出的波系结构与在ＣＡＲＤＣ小型激波管中做的相关实验结果一致。     

Euler方程的动力学通矢量分裂格式

文中基于气体分子动力学理论，研究和发展了一类新型逆风差分格式——动力学通矢量分裂格式（ＫＦＶＳ），并论证了格式的逆风性质和Ｌ１稳定性。最后给出了高精度格式的构造，并就一维激波管问题作了数值试验，还与精确解作了比较，结果令人满意。     

一类基于非等距单元平均重构的高效差分格式

构造了一维非线性双曲型守恒律的一类基于非等距单元平均值的点值重构的高精度高分辨率守恒型差分格式。其构造思想是:首先,将计算区间划分为若干个互不重叠的小区间,再根据格式精度的要求利用Gauss-Lobatto点和Gauss-Chebyshev点划分小区间,通过各非等距细小区间上的单元平均值,重构各细小区间交界面上的点值,并加以校正;其次,利用近似Riemann解计算细小区间交界面上的数值通量,并结合高阶Runge-Kutta TVD时间离散方法,得到了高精度的全离散方法。证明了该格式的无振荡特性。然后,将格式推广到一维双曲型守恒方程组情形。最后,给出了几个标准数值算例,验证了格式的高效性。     

一类分数阶薛定谔方程孤立解的对称性研究

在有界环形区域上,研究了一类分数阶薛定谔方程孤立解的对称性问题。首先将分数阶薛定谔方程转化为包含Bessel位势和Riesz位势的积分方程组,然后利用移动平面法和Hardy-Littlewood-Sobolev不等式,证明了当方程边值为常数时,环形区域必为同心球,方程正解是径向对称的,且随着到对称点的距离增大而单调递减。     

非光滑动力系统周期响应的数值解法

工程中有重要意义的非线性动力系统周期响应分析,一般因系统维数高、特性复杂而需要用数值方法进行。但目前有影响的打靶法、增量谐波平衡法均不能处理含弹塑性、间隙、干摩擦等非光滑因素的工程系统,计算光滑动力系统的强非线性周期响应时也常不收敛。本文分析了上述两种方法的欠缺,指出其原因在于方法中直接或隐含的Newton迭代格式。文中提出了用拟Newton迭代格式和无约束优化方法改进打靶法和增量谐波平衡法的8种方案,并对这些方案就收敛性、精度和效率进行了考核比较,给出了不同方案的适用对象,精选出了几种方案推荐工程界应用。