求解非线性随机微分方程加权格式的收敛性

通过研究加权格式用于求解非线性随机微分方程的收敛性,利用随机变量服从正态分布的性质,得到了在噪声为乘性噪声时,加权格式用于求解非线性随机微分方程均值意义上的局部收敛阶为2,均方意义上的局部收敛阶为3/2,强收敛阶为1.

求解非线性随机延迟微分方程加权格式的收敛性

通过研究加权格式用于求解非线性随机微分方程的收敛性,利用随机变量服从正态分布的性质,得到了在噪声为乘性噪声时,加权格式用于求解非线性随机微分方程均值意义上的局部收敛阶为2,均方意义上的局部收敛阶为1,强收敛阶为1.

波动方程的三层加权格式

本文对波动方程的三层加权格式的稳定性和收敛性作了分析,并用数值例子进行验证并编程实现。

求解广义Rosenau-Kawahara方程的一个非线性加权守恒差分格式

对广义Rosenau-Kawahara方程的初边值问题进行数值研究。在二阶精度前提下,在空间层引入两个加权系数,构造了一个带有两个加权系数的两层非线性差分格式。该格式很好地模拟了原问题的一个守恒性质。利用离散泛函分析方法证明了该格式的二阶收敛性与无条件稳定性。数值实验表明,通过适当调整两个加权系数可使计算精度大幅度提高,证明本文提出的加权格式是有效的。

高阶加权非线性格式的拟线性频谱分析方法研究

拟线性频谱特性分析方法,能够更加准确地给出非线性空间离散格式的频谱特性,已成为非线性格式性能评估的重要手段,但时间离散方法和计算点数等因素严重影响了该方法的预估精度和使用的方便性。为了剔除这些因素的影响,首先,通过理论推导获得了与时间离散无关的频谱表达式,解释了该式中各项的物理含义,并给出了推进时间步长对频谱计算结果的影响;其次,基于时间离散无关的拟线性频谱分析方法,分析了格式频谱特性曲线在某些点发生跳跃的主要原因是选取了不恰当的计算点数,并给出了一种计算点数选取方法,当计算点数超过某整数后,该方法可以有效地消除计算点数和初始相位变化对格式频谱特性所带来的影响。在此基础上,基于两个三阶精度WCNS格式,开展了它们的频谱特性和典型算例的数值模拟研究,结果表明,所发展的方法得到的拟线性频谱特性在定性上能够正确评价非线性空间离散格式特性,但定量上仍显不足。

RLW方程的高精度加权平均差分格式

对正则长波(RLW)方程建立一种高精度加权平均差分格式,所建格式是三层隐式的,具有二阶时间精度和四阶空间精度,且能够精确地模拟原问题的质量和能量守恒.利用离散能量法证明了所建格式解的存在性、唯一性、收敛性和稳定性.数值实验验证了理论分析的有效性.

控制体有限元方法中质量加权迎风格式的迁延修正

控制体有限元方法对于处理流动问题提供了一种较好的数值模拟方法，但对于大Ｒｅｙｎｏｌｄｓ数的回流流动问题，该方法在采用各种迎风格式时仍存在着严重的假扩散．为此，对二维质量加权迎风格式，给出了一种由插值单元及其一个与插值点最靠近的相邻单元所确定的多项式插值格式迁延修正的形式．通过对Ｒｅ＝４００、１０００、５０００三种情形的二维空穴流问题的计算，表明该格式可有效地减少控制体有限元方法对回流问题数值模拟的假扩散．

双时间步长加权ENO-强紧致高分辨率格式及在叶轮机械非定常流动中的应用

在非结构网格下,针对非定常三维N-S方程组发展了一种双时间步长高精度快速迭代格式。该格式在时间上具有二阶精度,在空间上将r=3的加权ENO格式与强紧致格式相结合去处理N-S方程中的对流项以及离散方程的右端项,并用四阶精度的紧致格式去计算N-S方程中的粘性项。典型的3个算例从不同侧面对格式进行了考核。计算表明:该算法具有高效率与高分辨率的特征,所得的计算结果与相关实验数据比较吻合,初步表明了该算法可以在非结构网格下模拟非定常流动的物理过程。

非结构网格下非定常流场双时间步长的加权ENO-强紧致杂交高分辨率格式

针对三维非定常、可压缩流场的Navier-Stokes方程组,本文提出一种新的双时间步长高精度快速迭代格式。该格式在时间上具有二阶精度,在空间离散上不低于三阶。在对流项与粘性项的处理上,本格式分别采用了加权ENO-强紧致格式与紧致四阶精度格式的思想。几个典型算例的实践表明:计算结果与相关实验数据比较吻合,初步表明了该算法可以在非结构网格下具有高效率与高分辨率的特征。

扩散方程的高精度加权差分格式

基于已发展的二阶微商三次样条四阶逼近公式，提出了数值求解含源汇项扩散方程的二层三点且在空间方向上达到四阶精度的加权差分格式．通过Ｆｏｕｒｉｅｒ方法讨论了文中格式的稳定性．证明了当１／２θ１时，格式是无条件稳定的，而当０θ＜１／２时，只有０＜ｒ１／［３（１－２θ）］，格式是稳定的．其中θ是权参量，ｒ＝Ｄτ／ｈ２为Ｆｏｕｒｉｅｒ数，Ｄ为扩散系数，而τ。

两边空间-时间分数阶扩散方程的加权有限差分格式(英文)

对于空间-时间分数阶扩散方程的初边值问题提出了一种加权差分格式.利用能量估计,得到了差分格式的稳定性.然后使用数学归纳法证明了在相同的条件下,所提出的的格式是收敛的.最后通过一个例子说明了所提出的格式是可靠的、有效的.

高精度加权紧致非线性格式的研究进展

综述了高精度加权紧致非线性格式WCNS在理论分析以及复杂流动应用方面的研究进展.首先回顾了国内外高精度格式研究的概况,然后介绍了WCNS的研究与发展历程.在对WCNS进行了Fourier分析和渐近稳定性分析后,给出了WCNS求解多维复杂流动的算例.

加权基本无振荡格式研究进展

加权型基本无振荡WENO格式是近十年发展起来的一类高阶、高精度格式,它是在ENO格式的基础上采用加权思想构造的,对流场内的间断和细致结构具有较高的分辨率,适于求解包含激波、膨胀波以及接触间断等复杂结构的流场,目前已发展成为计算流体力学中一种重要的方法。本文针对加权型基本无振荡格式近年来的进展作一简要介绍。

加权ENO格式的构造及数值模拟

将加权ENO格式推广到非结构三角形网格上 ,构造了一类加权ENO有限体积格式 ,提出的插值多项式的构造方式 ,可以减少计算时间 .对于出现的病态方程组 ,给出了解决方法 .此外还给出了插值点的选取方式及加权因子的构造方法 .结合三阶TVDRunge Kutta时间离散 。

一种通量加权型紧致格式

将通量加权的思想引入到紧致格式中,构造了一个传统方法与紧致格式混合组成的加权型差分格式.利用该格式与二阶TVD格式分别计算了一维方波、组合波问题和一些黎曼问题,如Sod问题和Shu问题,以及一维定常激波问题.计算结果的比较表明加权格式无论在捕捉各种间断,还是在分辨各种复杂波系上,都具有较大的优势,并与精确解非常吻合.

二维扩散反应方程的三次样条高精度加权隐式差分格式及其多重网格方法

利用二阶微商的三次样条四阶紧致差分逼近公式,推导出两种数值求解二维扩散反应方程的两层9点加权隐式紧致差分格式.当θ=1/2时,该格式在时间和空间方向上分别达到二阶和四阶精度.通过Fourier方法讨论知,当1/2≤θ≤1时,格式是无条件稳定的;当0≤θ<1/2时,格式是条件稳定的.为了克服传统迭代法在求解隐格式方面的困难,差分方程采用多重网格方法进行求解并将本文格式的结果与P-R格式及C-N格式下的结果进行比较.数值实验结果验证本文方法的精确性和可靠性及多重网格方法的效率.

解抛物型方程的一种高精度加权差分格式

利用二阶徽商的四阶精度紧致差分逼近公式，给出解抛物型方程精度为O[1－20）t，t2＋x4]的一种新的加权差分格式，并通过Fourier方法讨论格式的稳定性．证明了当1／2≤θ≤1时，格式是无条件稳定的；当0≤θ＜1／2时，只有r≤1／3（1－2θ），格式才是稳定的，其中θ是加权参数（因子），t，x分别为时空方向的网格长度，r=（D是二阶导数项系数）．

加权型紧致格式与加权本质无波动格式的比较

线性紧致格式和加权本质无波动格式是两种典型的高阶精度数值格式,它们各有优缺点.线性紧致格式在具有高阶精度的同时,格式的分辨率也比较高,耗散低,是计算多尺度流场结构的较好格式,但是不能计算具有强激波的流场.加权本质无波动格式是一种高阶精度捕捉激波格式,鲁棒性好,但耗散比较高,分辨率也不理想.近年来,在莱勒的线性紧致格式基础上,采用加权本质无波动格式捕捉激波思想,发展了一系列加权型紧致格式.本文较全面地比较了加权型紧致格式和加权本质无波动格式,包括构造方法、鲁棒性、分辨率、耗散特性、收敛特性以及并行计算效率.结果表明,现有的加权型紧致格式基本保持了加权本质无波动格式的性质,对于气动力等宏观量的计算,比加权本质无波动格式没有明显的优势.

双重加权实质无波动激波捕捉格式

本文构造了一种新的加权激波捕捉格式,命名为双重加权实质无波动(DoubleWeightedEssentiallyNon-Oscillatory,简称DWENO)格式。DWENO格式不需要引进最优权值,并且可以递推地达到任意高阶精度。文中给出了典型算例的计算结果,这些结果表明DWENO格式具有高于WENO格式的间断分辨率。

扩散方程高精度加权差分格式的MATLAB实现

对扩散方程混合问题,利用二阶微商三次样条四阶逼近公式构造其四阶加权差分格式.使用MATLAB软件编程,将问题的解用图像表示出来,通过数值结果验证了该方法的可行性和稳定性.