基于二维泊松方程六阶紧致格式的多重网格方法

利用六阶紧致差分格式、结合多重网格V循环算法求解了二维泊松方程的Dirichlet边值问题，并用不同的松驰算子与四阶精度格式的多重网格方法进行了比较，计算结果表明，该方法在不明显增加计算量的前提下较四阶精度格式的多重网格方法具有更好的精确度和收敛阶，且ZLGS迭代不论对四阶精度还是对六阶精度格式的多重网格算法，都是一种较其他松弛算子更加有效的“光滑剂”。     

求解一维线性双曲型方程的高精度紧致差分格式

【目的】双曲型方程是一类重要的偏微分方程,由于寻求问题本身的精确解比较困难,数值方法来求解此类方程有极具深远的意义和实际应用价值。【方法】首先对于一维的线性双曲型方程,在空间上采用Kreiss提出的四阶紧致差分公式进行逼近,时间上采用Taylor级数展开及截断误差修正的方法,推导出一个隐式的紧致差分格式。【结果】该格式在时间和空间上都有四阶精度,截断误差为O(τ4+h4)。【结论】采用Fourier方法分析了该格式的稳定性。数值实验证明提出的格式具有较好的稳定性和精确性。
     

三维热传导方程恒稳定的高精度半显式差分方法

提出了数值求解三维热传导方程的一种无条件稳定的高精度半显式差分方法,该方法可以显式计算且计算量小,截断误差为O(τ2+h4).数值算例验证了方法的精确性和可靠性.     

饰纹姬蛙Dmrt基因研究过程中异常结果的分析

对DMRT基因的研究需要运用到一系列的分子生物学实验技术,如DNA提取、简并PCR、DNA切胶纯化、克隆、SSCP分析及测序等.在实验的各个环节都可能会出现非预期的异常结果.对这些异常结果的分析研究也是实验的一个重要方面,不仅有助于促进实验的顺利进行,提高整个实验的效率,并且还有可能从中发现一些新的问题.     

边界层对流扩散方程在自适应网格上的高精度紧致格式

对流扩散方程在非均匀网格上的高精度紧致格式具有精度高、模板小等优势,然而现有方法往往需要事先指定边界层或大梯度的位置,利用网格分布函数生成非均匀网格并调整网格分布参数,这严重影响方法的适用性.本文提出了一种在二维计算区域上生成正交网格的自适应方法,根据物理解的特征对网格的分布进行自适应调整,并结合非均匀网格上的高精度紧致格式对一维及二维边界层对流扩散问题进行求解.数值结果表明,自适应网格方法结合高阶紧致格式可以有效求解边界层问题,提高数值解的精度,减少计算所需的网格和计算量.     

求解二维扩散方程的一种高精度紧致差分格式

扩散方程通常用来描述扩散现象中的物质密度的变化或者与扩散相类似的现象,针对二维扩散方程提出了一种高精度紧致差分格式,该格式基于四次样条函数对空间变量进行离散,对时间导数采用(2,2)Padé逼近,从而得到了时间和空间均为四阶精度的紧致差分格式.然后证明了该格式是无条件稳定的.最后通过数值实验,验证方法的精确性和稳定性.     

数值求解一维波动方程的四阶紧致差分方法

针对一维波动方程提出了一种有限差分方法.首先,采用泰勒级数展开公式和原方程代入的方法推导出了第一个时间层未知函数值的四阶紧致差分格式.然后,用四阶紧致差分公式近似空间导数项,采用中心差分格式截断误差余项修正的方法处理时间导数项,推导出了第二个时间层以后未知函数的四阶紧致差分格式.该方法时间和空间具有整体四阶精度.利用Fourier方法分析了所提格式的稳定性.由于本文格式在未知时间层仅涉及3个网格点,因此可采用追赶法求解离散化后所得到的线性方程组.最后,用数值算例验证了本文格式的精确性和稳定性.     

求解扩散方程的一种高精度隐式差分方法

利用一阶微商和二阶微商的四阶紧致差分逼近公式,推导出了数值求解一维扩散方程的两种新的高精度隐式紧致差分格式,其截断误差分别为O(τ^2 h^4)和O(τ^4 h^4)．通过Fourier分析方法证明了格式O(τ^2 h^4)是无条件稳定的,而格式O(r^4 h^4)是无条件不稳定的．并且由于每一时间层上只用到了3个网格点,所以差分方程可采用追赶法直接进行求解．