非结构网格上Euler方程有限体积解法的改进

对二维非结构三角形网格上Euler方程有限体积解法的格点格式进行了一些改进，重点在于提高数值解的精度，细致处理人工粘性项的尺度因子以及对该项建立适当的边界条件；发展一种新的基于最长边剖分的三角形网格自适应加密方法。采取多步Runge—Kutta格式推进、当地时间步长、隐式残值光顺等措施加速迭代收敛。文末给出的数值结果非常接近于参考文献中结构网格上的结果，验证了所发展数值方法的精确性。     

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

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

高超声速马赫反射

通过实验和数值计算模拟了激波（马赫数７．６）绕４０度压缩拐角的流动。实验是在气动中心激波管中完成的，数值模拟采用了层流的Ｎ－Ｓ方程，对无粘项、粘性项和化学源项分别采用了迎风ＴＶＤ格式、中心差分格式和点隐式的方法进行差分离散，计算了平衡与非平衡条件下的马赫反射流动，其数值计算结果与实验结果符合甚好。     

基于混合动网格的二维非定常粘性流动数值模拟

利用基于Delaunay三角化的动网格方法,在混合网格上求解了带运动边界的二维非定常Navier-Stokes方程。在有限体积法格心格式的基础上,采用带人工耗散项的Jameson中心格式以及双时间步长推进方法对Navier-Stokes方程中的空间项和时间项分别进行了离散。湍流模型为Spalart-Allmaras方程模型。对在跨声速粘性绕流中进行俯仰振动以及后缘摆动的翼型进行了数值模拟,得到了令人满意的结果。     

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

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

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

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

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

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

一类线性发展方程的交替紧致差分格式

文章给出求解一类线性发展方程的交替方向紧致差分格式并运用能量法证明它是无条件稳定的.Richardson外推法使时间方向上有四阶精度,数值结果表明新算法的高精度和有效性.     

一类求解Hamilton-Jacobi方程的交错网格差分格式

Ｈａｍｉｌｔｏｎ－Ｊａｃｏｂｉ方程在控制论和微分对策中有广泛的应用,由于其表达形式与双曲守恒律方程极为相近,这有利于借助于求解双曲守恒律方程的差分格式来构造求解Ｈａｍｉｌｔｏｎ－Ｊａｃｏｂｉ方程的差分格式。文中将Ｈａｍｉｌｔｏｎ－Ｊａｃｏｂｉ方程变化为双曲守恒律方程,利用求解双曲守恒律方程的交错网格的Ｇａｕｓｓ型差分格式,构造了一类求解Ｈａｍｉｌｔｏｎ－Ｊａｃｏｂｉ方程的交错网格的Ｇａｕｓｓ型差分     

二维欧拉方程高分辨正性保持差分格式

针对二维欧拉方程研究了一类高分辨正性保持波尔兹曼型差分格式,首先分析了波尔兹曼方程和欧拉方程之间的关系;然后,利用一种特殊的插值技术构造了一类高分辨波尔兹曼型差分格式.最后通过数值实验验证了格式的实用性和有效性.     

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

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

非结构网格的高超声速粘性MHD流场数值模拟

应用数值方法在非结构网格上对磁场干扰下的二维高超声速钝头体粘性绕流进行了数值模拟。控制方程为N-S方程耦合Maxwell方程的粘性MHD方程组,空间离散采用有限体积方法,对流项用AUSM格式计算,粘性项用中心格式求解,时间推进用显式5步Runge-Kutta格式,引用双曲型散度清除技术加强▽.B=0的条件。计算模型为二维钝头体,在高超声速来流条件下,分别对有、无均匀磁场干扰下的流动进行了数值模拟。计算结果表明,在均匀磁场干扰下,激波脱体距离显著增加,物体表面压力急剧下降。对比表明文中发展的计算方法可以准确地进行二维粘性MHD流场的数值模拟。     

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

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

注水井地层伤害压力偏微分方程模型及求解

建立注水井地层伤害二维抛物型偏微分方程新模型,并进一步化简为带间断系数的抛物型方程初边值问题.考虑到差分方程解的收敛性,采用积分关系法建立其差分格式,并指出此模型在油田注水开发中的一些应用.     

混合网格在求解三维N—S方程中的应用

用三棱柱／四面体组成的混合网格求解三维Navier-Stokes方程，在物面附近采用三陵柱网格，而其他地方则采用四面体网格，与单一类型网格相比，混合网格易于调节粘性层内网点数目，并能减少存储量，提高解的效率，Navier-Stokes方程的求解采用有限体积法，对人工耗散项进行了改进，使得求解的结果更接近真实的物理解，本文对上述方法开展了数值实验研究。结果表明，用三棱柱／四面体组成的混合网格求解N－S方程是非常有效的。     

大迎角三角翼旋涡运动及其破碎特性的数值研究

从流体力学的基本方程出发，利用Ｈａｌｌ的涡核准柱假设，导出反映涡核运动的Ｎ－Ｓ方程。采用差分方法计算旋涡流场，进而分析三角翼上前缘分离涡的运动特点及其破碎机下。从计算结果可以看出，旋涡的轴向速度向下游逐渐下降，且涡心处于降较快，外缘下降较慢，反映了粘性作用自涡心外缘逐渐下降的特点；涡核外缘的径向速度开始为负，说明开始阶段有流体流入涡核，随着旋涡向下游运动，径向速度有所增加，到一定位置后增加迅速，说     

机翼非平行边界层稳定性研究

从Navier-Stocks(N-S)方程导出曲线坐标系下的抛物化稳定性方程(Parabolicstabilityequation,PSE),研究机翼非平行的可压缩边界层稳定性问题。发展了求解PSE的高效数值方法:引进法向变换,使得在临界层与壁面之间的扰动量变化最快的区域有更多法向网格点;采用包含边界邻域在内的完全四阶精度的法向差分格式,这对方程精确离散至关重要;以及全局法和局部法相结合的数值方法及其新的迭代公式,能大大加速收敛并得到更精确的特征值。算例分析研究了扰动增长因子和形状函数等演化曲线。     

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

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

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

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

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

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