用半拉格朗ElENO和WEN0格式解Vlasov方程

借鉴ENO和WENO格式的活动模板选择的思想，构造了一种新的semi-Lagrange类格式，并且用它计算了Vlasov方程．该格式在空间插值函数的构造上，并不使用固定的模板来构造插值函数，而是利用Newton差商来选取不同的模板，以达到对振荡的抑制作用．算例表明，该格式对振荡具有很好的抑制作用，同时又具有semi—Lagrange类格式的可以用大CFL条件的特点．该格式对Vlasov方程的计算结果也是令人满意的．     

重心插值配点法求解Volterra积分方程

文章提出了求解Volterra积分方程的一种高精度数值方法：重心插值配点法（包括重心Lagrange插值配点法和重心有理插值配点法）。该方法分为两步：首先对Volterra积分方程采用两种重心插值配点法进行离散，构造出Volterra积分方程的数值求解格式；然后，依次选取第二类Chebyshev节点和等距节点进行数值计算。文章主要研究积分项中含有未知函数的一阶导函数的Volterra积分方程的离散格式构造及数值实现。数值实验结果表明：在使用第二类Chebyshev节点时，用重心Lagrange插值配点法较好；在使用等距节点时，使用重心有理插值配点法较好。     

对流占优扩散问题的一种特征差分方法

用基于一般的 L agrange插值的特征差分方法求解对流占优扩散问题 ,会出现较大的数值扩散或者数值振荡等困难 ,高阶单调插值又计算复杂。该文采用 A.A .Sam arskii构造差分格式的方法 ,建立了一种新的特征差分方法。先对对流扩散方程的扩散项进行修改 ,然后再进行特征差分。此方法具有较高精度 ,并消除了非物理振荡。证明了方法的无条件稳定性。数值结果表明 ,该方法可成功求解对流占优扩散问题。     

强间断多介质流的高精度伪弧长方法

考虑到双曲守恒律方程随着时间的发展,产生的解会包含强间断. 伪弧长算法可以削弱方程的奇异性,因此结合伪弧长算法和高精度权基本无振荡格式发展了高阶伪弧长算法,有效提高计算求解的精度和分辨率.由于在变形的网格中直接构造高精度格式较为复杂,因此通过坐标变换将控制方程映射至正交均匀的弧长空间,然后在弧长空间中完成计算. 结合Level Set技术和虚拟流体法界面处理,将算法拓展到多介质流的计算中,针对网格移动以后Level Set距离函数的插值,提出了三阶非守恒插值格式.计算结果表明,高精度伪弧长算法具有较高的收敛阶,可以有效降低间断处的数值震荡,提高间断分辨率.      

插值样条δ-序列求解非线性对流扩散方程

提出了一种用广义函数δ序列求解偏微分方程的数值方法.首先对一阶B样条函数N1(x)进行卷积得到四阶B样条函数N4(x),用N4(x)的线性组合构造出三次样条插值基函数;然后用样条插值基序列逼近δ函数,利用δ函数的性质构造插值样条δ序列,该δ序列具有对称、Riesz基和插值性质.以非线性对流扩散方程(伯格方程)为例,用插值样条δ序列离散该方程的空间形式,用四阶龙格库塔方法描述发展过程,取得了较好的精度.为减少计算量,加快插值函数的收敛速度,进一步提高求解精度,对δ序列进行了改进,对同一算例进行数值实验,结果表明,改进后的算法求解过程稳定发展,能够有效描述局部快速变化的情况.     

一类求解浅水波方程的基本无振荡熵稳定格式

针对浅水波方程,提出了一类低耗散基本无振荡熵稳定格式.在Roe型熵稳定通量中添加熵守恒格式的熵数值黏性绝对值的量来抵消解在跨过激波时所产生的熵增,从而抑制伪振荡;并且,利用通量限制器函数构造出相应的高分辨率熵稳定格式.利用新格式模拟一维和二维经典问题,数值结果表明,该格式具有低耗散、高分辨率、基本无振荡性等特点,是求解浅水波方程较为理想的方法.     

一种Lagrange格式及重映算法

从积分形式的二维Lagrange流体力学方程组出发，用有限体积格式进行计算，考虑压力梯度分布对速度和能量改变的影响，构造了在两个控制体上的动量方程的计算格式。在重映算法上，采用积分重映的方法，针对不规则的四边形网格，根据非结构网格的ENO插值的思想，构造线性插值多项式，由于利用其插值点自适应选取的特性，在物理量变化剧烈区域选取最光滑区域的点来插值，虽然不能保持单调性，但只允许出现非常小的振荡。数值结果表明了该方法的可行性。     

一种基于弦截法的预估校正格式

利用弦截法计算预估值,并在预估值附近用抛物插值方程代替非线性方程,再在预估值处对抛物方程使用牛顿公式计算校正值,构造了一种无需计算导数值的新型预估校正格式,从理论上分析了格式的近三阶收敛性.该方法每步迭代需要计算三次函数值,因此其效率指数为1.719.最后通过数值算例验证所提格式的可行性和有效性.     

Chebyshev-Gauss-Lobatto节点的一个注记

利用Lagrange插值基函数和Chebyshev多项式的性质,推导以Chebyshev-Gauss-Lobatto点为插值点构造的插值基函数的一阶、二阶微分矩阵的显示格式,并由插值点的性质得出两者之间的关系.通过对具有解析解的一维对流扩散方程进行数值求解,验证了一阶、二阶微分矩阵显式格式的正确性.数值结果表明:由微分矩阵显式格式可以方便地构造配置点谱方法中的拟谱算子,利用其求解微分方程,在较少的网格点时,即可得到快速收敛的高精度的数值结果.研究工作对配置点谱方法的应用具有一定的理论指导意义.     

求解非线性伪抛物方程的重心Lagrange插值配点法

提出了重心Lagrange插值配点法求解一类非线性伪抛物方程。首先,介绍了重心Lagrange插值并给出了微分矩阵表达式。其次,构造了求解非线性伪抛物方程的直接线性化迭代格式、部分线性化迭代格式、Newton线性化迭代格式。再次,未知函数和初边值条件利用重心Lagrange插值函数来近似,利用配点法得到离散方程,获得了方程的矩阵表达式。最后,数值算例表明,重心Lagrange插值配点法具有高精度和高效率的优点。