基于带波动算子非线性Schr(o)dinger方程数值分析的守恒差分算法

对一类带波动算子的非线性Schr(o)dinger方程进行了数值分析,提出了一个含参数的二阶守恒差分格式,根据参数选取的差异,该格式既可隐式计算也可显式计算.对初值条件进行了中心差分离散,使其具有二阶精度,从而与守恒格式的精度一致.利用矩阵理论证明了差分解的存在惟一性,并利用一个重要的不等式在先验估计的基础上,运用能量估计的方法证明了该格式按无穷范数以二阶精度收敛到真实解.数值实验表明该格式具有较高的计算效率.     

基于带波动算子非线性Schroedinger方程数值分析的守恒差分算法

对一类带波动算子的非线性Schrodinger方程进行了数值分析，提出了一个含参数的二阶守恒差分格式，根据参数选取的差异，该格式既可隐式计算也可显式计算。对初值条件进行了中心差分离散，使其具有二阶精度，从而与守恒格式的精度一致。利用矩阵理论证明了差分解的存在惟一性，并利用一个重要的不等式在先验估计的基础上，运用能量估计的方法证明了该格式按无穷范数以二阶精度收敛到真实解。数值实验表明该格式具有较高的计算效率。     

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

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

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

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

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

利用Hamilton-Jacobi方程与双曲型守恒律的紧密联系,借助于求解双曲型守恒律的一类无波动无自由参数的耗散差分(NND)格式构造了一类求解Hamilton-Jacobi方程的差分格式.数值实验结果表明:该格式具有计算量小且高分辨率等优点.     

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

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

一类基于非均匀剖分的二维高精度差分格式

利用dimension-by-dimension方法,将求解一维非线性双曲型守恒律的一类基于非等距单元平均值重构的高效差分格式推广到二维标量双曲型守恒律方程,得到求解二维双曲型守恒律的一类二维高精度差分格式.证明了该类格式的无振荡特性.然后,将格式推广到二维双曲型守恒方程组情形.最后,给出了几个标准数值算例,验证了格式具有高阶精度、高效捕获激波等间断的能力.     

双曲型守恒律的无振荡中心差分格式

讨论了双曲型守恒律的一类无振荡中心差分格式。Ｈ．Ｎｅｓｙａｈｕ和Ｅ．Ｔａｄｍｏｒ研究了以交错型ＬａｘＦｒｉｅｄｒｉｃｈ格式（ＬｘＦ）为模块的无振荡中心格式的构造与熵不等式。此类格式利用高阶的ＭＵＳＣＬ型插值替代一阶分片常数逼近,减少了ＬａｘＦｒｉｅｄｒｉｃｈ格式的过多数值粘性,建立了一维标量非线性双曲型守恒律的一类高分辨格式。讨论以非交错ＬａｘＦｒｉｅｄｒｉｃｈ格式为模块建立起的差分格式。证明了此格式具有二阶精度、ＴＶＤ性质并在一定条件下满足熵条件。     

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

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

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

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

一类新的二阶精度非交错形式差分格式的构造

本文对二维标量非线性双曲型守恒律，以非交错形式Ｌａｘ－Ｆｒｉｅｄｒｉｃｈｓ格式为基本模块，构造了一类新的不须解Ｒｉｅｍａｎｎ问题的差分格式，并证明其具有时空二阶精度。数值实验表明，格式具有较高的分辨激波、稀疏波及接触不连续的能力     

高精度混合型格式的构造及数值模拟

在守恒变量形式的无波动无自由参数耗散(NND)格式的基础上,引入样条插值的思想,并采用具有一致二阶精度格式的限制器函数,构造了一个在光滑区具有四阶精度,在激波附近降为二阶精度的混合型格式。在计算中采用具有总变差减少(TVD)性质的四阶R unge-K u tta方法进行时间离散,通过对多个典型算例的计算,并与NND格式作了比较,分析结果表明:混合四阶格式引起的耗散和波动较小,并且能够高精度地分辨间断。     

一类新的二阶精度非交错形式差分格工工的构造

本文对二维标量非线性双曲型守恒律，以非交错形式上Ｌａｘ－Ｆｒｉｅｄｒｉｃｈｓ格式为基本模块，构造了一类新的不须解Ｒｉｅｍａｎｎ问题的差分格式，并证明其具有时空二阶精度。数值实验表明，格式具有较高的分辨激波，稀疏波及接触不连续的能力。     

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

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

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

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

基于非确定性的翼型鲁棒反设计方法

实际的翼型是在一个非确定性的状态范围内工作的,为此应用鲁棒控制理论建立了翼型的鲁棒气动反设计方法及数值求解过程,基本的优化器为基于伴随方程的位流优化方法。对守恒型位流方程应用高效AF 2格式求解,对非守恒型的伴随方程采用了旋转差分和高效AF 3格式,提高了梯度计算的精度和效率。优化算例表明:本文方法所得到的翼型在一个较宽的工作状态范围内都具有良好的性能。     

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

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

空气动力学数值模拟的一种分量形式TVD格式

本文选取一个精度较高的ＴＶＤ限制子,把著名的二阶Ｍａｃ Ｃｏｒｍａｃｋ格式和Ｗａｒｍｉｎｇ－Ｂｅａｍ格式结合,构造了一种分量形成的二阶ＴＶＤ差分格式。通过两个典型算例的计算,并与ＥＮＯ格式的结果比较,显示了该格式精度高和计算逻辑简单的特点。     

基于混合网格的直升机旋翼/机身非定常干扰流场数值模拟方法

基于混合网格建立了一个适合于直升机旋翼/机身非定常干扰流场的数值模拟方法。其中,采用混合网格系统生成旋翼/机身干扰流场网格以降低网格生成难度;将二阶MUSCL格式和Roe格式相结合进行空间离散以提高流场的计算精度;同时引入隐式LU-SGS时间格式和OpenMP并行模式以提高干扰流场的计算效率。应用建立的计算方法,对Georgia和Robin旋翼/机身模型进行了数值模拟和分析,与试验结果良好的相关性表明了该方法的有效性。     

基于WENO-分段线性格式的旋翼流场数值模拟

为了降低二阶JST(Jameson-schmidt-turkel)格式导致的数值耗散,发展了一套适合于格心格式的基于加权本质无振荡(Weighted essentially non-oscillatory,WENO)-分段线性格式的旋翼流场数值模拟方法。采用运动嵌套网格方法生成围绕旋翼的网格系统,主控方程选择Navier-Stokes(N-S)方程。为了更加有效地对桨尖涡等流动细节进行捕捉,在笛卡尔背景网格上采用七阶Roe-WENO格式计算对流通量;在桨叶贴体非结构网格上采用二阶精度的Roe-分段线性格式计算对流通量。时间离散采用了高效的双时间隐式LU-SGS(Lower upper symmetric Gauss-Seidel)方法进行时间推进。最后,应用上述方法对悬停状态的C-T(Caradonna-tung)旋翼和Helishape 7A旋翼进行了数值模拟,将数值计算结果与实验数据进行了对比,计算值与实验值吻合较好;并将桨尖涡模拟效果与二阶JST格式的模拟效果进行了对比。对比结果表明:本文方法能有效对旋翼流场进行计算,且在相同计算条件下,WENO-分段线性格式能够更有效地捕捉旋翼涡流场的流动特性,表明在计算旋翼涡流场时WENO-分段线性格式相比传统二阶JST格式具有更低的数值耗散。