Newton迭代法的一个修正格式

本文从等价方程组出发,导出了 Newton 迭代法的一个新修正公式。证明了它具有大范围收敛性和不低于线性收敛速度,并给出了数值例于。

延迟修正QUICK格式及其在弹状流中的应用

将延迟修正QUICK格式应用于数值传热学经典问题进行验证性计算 ,应用流体体积模型 ,将延迟修正QUICK格式及界面移动施主 受主格式与SIMPLEC算法相结合 ,计算了垂直管内弹状流泰勒气泡周围的流动特性 .通过计算 ,得到了满足气泡内部压力恒定要求的泰勒气泡形状及其周围流场分布 ,以及泰勒气泡周围液膜厚度随距离的变化规律 ,并与实验值进行了比较 ,得到的结果较为理想 .

动态中子输运方程的修正时间离散格式

传统方法很少考虑时间步长对数值格式的影响,导致关于时间微分的物理量出现振荡,数值计算结果的精度较低。针对自适应时间步长的特点,对球几何动态中子输运方程的时间离散格式进行了研究,构造了修正时间离散格式。数值算例表明:修正时间离散格式简单且精度较高,所需迭代次数较少,避免了时间步长变化带来的数值解的振荡,更加适合自适应时间步长的计算。

解热传导方程的一类新修正差分格式

利用lie乘积公式对一维热传导方程的初边值问题设计了一类新的数值方法,证明了此类方法中的修正隐方法和修正C-N方法是显示的无条件稳定的,且修正显方法的稳定域比经典显方法大.

基于ENO格式的三阶修正系数格式

不增加基点,仅摄动二阶ENO格式的系数(简记为MCENO),得到一类求解双曲型守恒律方程的三阶MCENO格式.由MCENO格式的构造过程可以看出,MCENO格式保留了ENO格式的许多性质,例如本质无振荡性、TVB性质等,且能提高一阶精度.进一步,利用MCENO格式模拟二维Rayleigh-Taylor(RT)不稳定性和Lax激波管的数值求解问题.数值结果表明,t=2.0时,MCENO格式的密度曲线处于三阶WENO格式和五阶WENO格式之间,是一个高效高精度格式.值得注意的是,三阶MCENO格式,三阶WENO格式和五阶WENO格式的CPU时间之比为0.62:1:2.19.表明相对于原始ENO格式,MCENO格式在光滑区域有较高精度,能提高格式精度.

Schrdinger方程的紧致修正交替方向格式

研究了多维Schrdinger方程的紧致修正交替方向格式.通过对J.Douglas等提出的交替方向格式进行误差分析可以发现其分裂误差远远大于时间离散的截断误差.为提高计算精度和效率,在格式中加入1个扰动项以提高分裂误差的阶数,使时间离散误差占优.数值实验验证了格式的优越性和扰动项的作用.

Rosenau-Burgers方程的修正局部Crank-Nicolson格式

本文给出了Rosenau-Burgers方程的两种修正局部Crank-Nicolson格式.首先,求解原有的偏微分方程对空间方向进行有限差分离散而得到的常微分方程.其次,利用矩阵分裂技术对这个方程的指数系数矩阵分别按行和元素进行逼近.最后,利用修正局部Crank-Nicolson方法得到了两种格式.讨论了格式的稳定性、收敛性和先验误差估计.数值实验结果表明了理论证明的正确性及格式的有效性.该格式具有结构简单、精度高的优点.

求解双曲守恒律方程的高阶TVD格式

针对精确度较低的单调格式,利用Osher-Chakravarthy线性通量修正格式在0<β≤{m2mm-1}的条件下对单调格式的数值通量进行改进,给出了一种求解双曲守恒律方程的高阶差分格式,并将其应用于Burgers方程及Buckley-Leverett问题.数值实验表明,该方法在保持高阶精确度的同时,还有效地防止了非物理数值振荡的产生,且具有较好的稳定性.

关于修正局部Crank-Nicolson法对于二维热传导方程的应用的注记

在《修正局部Crank-Nicolson法对于二维热传导方程的应用》一文中,作者把一维和二维热传导方程半离散化后,借助于泛函分析中Lie乘积公式,利用矩阵分裂得到修正局部Crank-Nicolson格式,该格式是显示差分格式,不需要直接解以大型矩阵为系数矩阵的线性方程组,从而计算简单,计算量小,在实际问题中有较大应用价值.本文作者在学习和应用该算法时,发现关于该算法的某些结论需要修正.

两种修正的Gauss级数序列及其它的收敛性

文章先构造级数计算的两种修正的G auss序列,并研究了修正格式的收敛性及应用等问题。

修正辛-准粒子法地震波数值模拟

基于Hamilton力学,从分子动力学角度建立了针对地震波传播的空间离散的准粒子体系,在该体系中各个粒子仅与上、下、左、右和4个对角粒子发生相互作用,且认为其作用力和相对位移呈近似线性关系,给出了粒子间的相互作用系数,并将该体系变换至Hamilton系统。在时间离散上,通过在二阶辛格式的基础上构建了一个新的修正辛格式,该格式具有三阶时间精度,并在理论上分析了其数值频散和稳定性。修正辛格式所有的系数均为正系数,符合时间演进方向,且具有长时间计算能力。为了测试修正辛格式和空间准粒子体系的准确性,利用拉梅问题验证对弹性波数值模拟的计算精度和效率,为了测试稳定性,选取Sigsbee 2B速度模型进行验证。结果表明,修正辛格式在数值频散压制和数值稳定性提升等方面具有明显的优势。

弱条件下列修正Broyden法的收敛性

在弱条件下,利用优序列技巧,给出了求解非线性方程组的列修正Broyden法的存在收敛性定理。

一类非线性抛物方程的直接间断有限元方法

文章针对一类非线性抛物方程构造全离散数值格式,在空间方向采用直接间断有限元方法,在时间方向利用修正的Crank-Nicolson格式。基于双线性算子强制性与连续性的详细讨论,得到全离散格式数值解的存在唯一性。随后的误差分析表明在能量范数下的最优结果为O(h^(2)+hΔt+Δt^(2)),并用数值算例验证该方法的有效性和理论结果。

非线性对流扩散问题的1种新数值方法

建立了非线性对流扩散问题的1种新数值方法——再生核函数法,利用再生核函数,直接给出了每个离散时间层上近似解的显式表达式,验证了算法的稳定性及收敛性,并给出了一些数值结果。

荷载属于H^(-1)(Ω)时不完全双二次非协调板元的新的误差估计式

在真解u∈H3(Ω),荷载f∈H-1(Ω)时,通过修正的离散格式,给出了不完全双二次非协调板元的能量模及L2-模新的误差估计式,它改善了以往相应的结果且减少了计算工作量.

US-FE-LSPIM四边形单元及其在几何非线性问题中的应用

为了提高在网格畸变时的数值计算精度,基于非对称有限单元的概念,提出US-FE-LSPIM四边形单元。该单元是利用传统的四节点等参元形函数集和FE-LSPIM四边形单元形函数集分别作为检验函数和试函数而构成。前者用于满足单元间和单元内的位移连续性要求,后者用于满足位移完备性要求。该单元结合了有限单元法和无网格法的优点,能方便地施加整段长度的位移边界条件。在分析几何非线性问题时,使用修正拉格朗日格式建立有限元方程,采用牛顿迭代法求解,编制了FORTRAN程序。数值算例表明,在规则网格和畸变网格下,US-FE-LSPIM四边形单元都具有很高的计算精度对网格畸变不敏感,性能优于传统的四节点等参元和QM6单元。

积分微分方程各向异性有限元的收敛性分析

本文研究具有各向异性特征的双二次元对具有积分型边界条件的积分微分方程的逼近问题,通过采用积分恒等式和插值后处理技术,在不需要Ritz-Volterra投影及任何修正格式情况下,利用该单元的特殊性质,在各向异性网格下得到了相应的超逼近和超收敛结果。

半导体瞬态问题的修正迎风有限体积格式

半导体器件的瞬时状态由包含三个拟线性偏微分方程所组成的方程组的初边值问题来描述.其中电子位势方程是椭圆型的,电子和空穴浓度方程是对流扩散型的.对电子位势方程采用一次元有限体积法来逼近,对电子浓度和空穴浓度方程采用修正的迎风有限体积方法来逼近,并进行详细的理论分析,关于位势得到O(h+Δt)阶的H^1模误差估计结果,关于浓度得到O(h^2+Δt)阶的L^2模误差估计结果.最后,给出数值例子.

曲率障碍下四阶变分不等式问题的Morley元逼近一个新的误差估计式

对简支情形的曲率障碍问题,采用修正的离散格式给出了Morley非协调元逼近的一个新的误差估计,改善了已有的结果.

大变形空间曲梁有限元分析(英文)

为了进行空间曲梁大变形有限元分析,将直梁单元位移分量插值的思想改进为位移矢量插值用以建立曲梁单元的位移场,分别建立了适用于任意曲线形式的全拉格朗日和修正拉格朗日增量格式空间曲梁有限元列式.通过采用等参数曲线代替实际曲线的策略,使得修正拉格朗日增量格式曲梁单元可以应用于更广泛的场合.算例对比结果表明,曲梁单元的建立过程正确,曲梁单元的精度要明显高于直梁单元.一般情况下,仅用直梁单元数量1/5的曲梁单元就可以达到相同的计算精度.