一类演化方程的高稳定性的差分格式

考虑一类演化方程u_t=au∂_2k+1(其中a是常数,u∂_2k+1=∂^2k+1u/∂x^2k+1,k=1,2……)的有限差分解法。构造了两类具有高稳定性的显式差分格式。并用引入耗散项的方法建立了两类半显式差分格式,它们是无条件稳定的且可显式地进行计算。     

求解奇异摄动问题的一个耦合差分格式

本文考察奇异摄动问题(1.1).在一特殊的非均匀网格上,将不稳定、二阶精度的中心差格式和稳定、一阶精度的Abrahamsson-Keller-Kreiss箱子格式相耦合,得到了一个二阶一致收敛的差分格式.最后给出了数值结果.     

拟线性抛物型方程奇异摄动问题的数值解法

本文讨论拟线性抛物型方程奇异摄动问题的差分解法,在非均匀网格上建立了线性三层差分格式,并证明了在离散的L₂范数意义下格式的一致收敛性,最后给出了一些数值例子.     

三维预处理技术初探

本文应用高阶近似LU分解法对三维椭圆型方程的正规七点差分矩阵给出了预处理阵,并对典型问题给出了预处理共轭梯度法的计算结果。文中所采用的分阶方论以及阶矩阵等概念,适用于一般的非奇异对角优势稀疏矩阵。对不同网格点数的计算结果表明,三维预处理共轭梯度法仍具有超线性的收敛速率,在高阶情形下超线性的特征尤为突出,由于非零对角线数增长过快,高阶的三维预处理方法不宜采用,零阶方法和一阶方法是值得推荐的。     

磁瑞利-泰勒不稳定性非线性演化数值模拟

基于理想磁流体力学模型,采用5阶WENO有限差分格式结合磁场伪散度清除投影算法,发展二维数值模拟程序,利用线性理论对程序进行检验,在此基础上研究磁瑞利-泰勒(MRT)不稳定性非线性演化物理过程和发展规律.结果表明:垂直于加速度方向磁场在线性和非线性阶段对MRT扰动发展都有较强的制稳作用,而平行于加速度方向磁场在线性阶段...     

三维对流扩散方程非等距网格上的四阶紧致格式及其多重网格方法

提出了数值求解三维变系数对流扩散方程非等距网格上的四阶精度19点紧致差分格式,为了提高求解效率,采用多重网格方法求解高精度格式所形成的大型代数方程组。数值实验结果表明本文方法对于不同的网格雷诺数问题,在精确性、稳定性和减少计算工作量方面均明显优于7点中心差分格式。     

高分辨率TVD格式的改进及应用

详细分析对比了已有的Harten TVD格式在三维任意曲线坐标系下的形式,尤其是其校正流量的构造。通过几何意义及物理意义的分析,指出了其缺陷,即校正流量中包含不必要的几何参数。由此,对校正流量的构造进行了几何改进,使其满足几何与物理要求。这一改进简单、意义明确、而且不增加计算量与存储量,在网格变化剧烈处可以明显改善计算的精度。同时,也考虑了预处理下的Harten TVD格式的几何修正,并通过计算证明了在预处理修正下几何改进仍然有效。     

三点强紧致加权高分辨率高阶格式及其应用

本文分别给出了采用有限体积计算时采用高阶格式遇到的网格单元体界面上三阶、四阶与五阶格式精度下的数值通量表达式,并且给出了确定权函数过程中所遇到的光滑因子表达式．文中首先对模型方程进行了格式分辨率方面的检验,然后将格式用于某涡轮级的实际流场计算．数值结果表明:所给出的高分辨率、高阶格式具有较高的激波分辨率并且具有较高的数值精度,在计算流场时,它可以用较少网格点去取代普通低阶精度格式下所采用的较密网格．     

非线性湿气迁移方程的超收敛分析

本文用混合有限元方法研究一般的非线性湿气迁移方程.利用双线性元Q₁₁和零阶Raviart-Thomas元(Q₁₀×Q₀₁)证明方程的超收敛性.利用这两个单元插值算子的性质和平均值技巧,得到了方程半离散格式的O(h²)阶超收敛结果.对于方程线性化的Crank-Nicolson(C-N)全离散格式,得到了具有O(h²+τ²)阶的超收敛结果,这里h是空间剖分参数,τ是时间步长.该方法说明如果线性化问题有超收敛性,那么对应的非线性问题有同样的超收敛性.最后,给出数值算例,证实了理论分析的正确性和方法的有效性.     

三阶非线性KdV方程的交替分段显-隐差分格式

对三阶非线性KdV方程给出了一组非对称的差分公式,用这些差分公式与显、隐差分公式组合,构造了一类具有本性并行的交替分段显-隐格式A·D2证明了格式的线性绝对稳定性．对1个孤立波解、2个孤立波解的情况分别进行了数值试验．数值结果显示,交替分段显-隐格式稳定,有较高的精确度．