一类非线性偏微分方程的有限差分格式的稳定性问题

对一类非线性偏微分方程进行完全离散,得到非线性有限差分格式.然后引入新的函数空间,并在这个空间上,采用变分近似法讨论了该格式的稳定性问题.文中在给出非线性有限差分格式的稳定性条件时,所用方法具有简单的特点.

平面弹性力学问题基于泡函数的简化的稳定化二阶混合有限元格式

该文研究平面弹性力学问题基于泡函数的稳定化二阶混合有限元格式,通过消去泡函数导出一种自由度很少的简化的稳定化二阶混合有限元格式,误差分析表明消去泡函数的简化格式与带有泡函数的格式具有相同的精度而可以节省18N_p个自由度(其中N_p是三角形剖分中的顶点数目).

黏声波方程波场模拟的通用格式自适应系数频域有限差分方法

黏声波方程常被用于描述地下介质的黏弹性及波的传播现象,频域有限差分(finite difference frequency domain,FDFD)方法是黏声波和黏弹性波波场模拟的常用工具.目前FDFD黏声波模拟常用的二阶五点方法和优化九点方法在一个波长内的网格点数小于4时误差较大.通过令FDFD系数随一个波长内的网格点数自适应从而提高FDFD方法的精度,本文针对黏声波波场模拟发展了一种适用于不同空间采样间隔之比的通用格式自适应系数FDFD方法.同时,为了验证自适应系数FDFD方法对一般黏声波模型的有效性,本文针对三个典型的黏声波模型,分别采用解析解和基于高阶FDFD的参考解验证了所提出方法的有效性.本方法的FDFD格式通过在传统的二阶FDFD格式的基础上引入相关校正项得到,其中校正项按网格点与中心点的距离进行分类选取,同时校正项对应的自适应FDFD系数不仅和空间采样间隔之比相关,还和一个波长内的采样点数相关.所需的自适应FDFD系数可通过声波方程的数值频散关系和查找表高效给出.数值频散分析表明,在空间采样间隔相等或不等的情况下,以相速度误差不超过1%为标准,通用格式自适应系数FDFD方法所需的一个波长内的采样点数均小于2.5.数值模拟实验表明,对于不同的空间采样间隔之比,相对于常用的二阶五点FDFD方法和优化九点FDFD方法,通用格式自适应系数FDFD方法均可在相似的计算量和内存需求下,有效提高黏声波模拟的精度.

一类p维扩散方程有限差分格式的稳定性问题

构造p维扩散方程一类有限差分格式,利用Taylor级数展开定理求出有限差分格式的截断误差并判断了有限差分格式是收敛性,利用Fourier积分法和Von Neumann条件得到当aλ≤1/2p时(其中λ=τ/h^2,τ和h分别为时间步长和空间步长),p维扩散方程的一类有限差分格式是稳定的。

不规则沟渠浅水流动的高阶平衡有限差分AWENO格式

在本文中,我们提出了具有不规则几何形状和非平坦底部地形的开放沟渠中浅水方程的高阶有限差分加权本质无振荡(Alternative Weighted Essentially Non-Oscillatory) AWENO格式。所提出的格式保持了静水稳态解的良好平衡性质,即在通量梯度和源项之间存在精确平衡时,在离散状态保持稳态。与具有恒定横截面的传统浅水方程相比,由于沟渠不规则几何形状引起的影响,构建良好平衡格式并不是一项简单的工作。为了保持良好平衡性质,我们首先重新构造源项,然后使用具有Lax-Friedrichs通量的静水重构方法,最后借助一种新颖的源项逼近方法离散源项。基准数值示例被应用来验证所得格式的良好性能:平衡性能,高阶精度,以及对不连续解的高分辨率。

非定常Stokes方程的稳定化CN有限体积元格式

建立二维非定常Stokes方程的时间二阶精度的稳定化Crank-Nicolson (CN)有限体积元格式, 并给出其稳定化CN有限体积元解的误差估计。数值实验说明时间二阶精度的稳定化CN有限体积元格式比时间一阶精度格式更优越, 从而表明稳定化CN有限体积元格式对于求解非定常Stokes方程的数值解是有效可行的。

带有几何源项的浅水波方程组的高精度熵稳定有限差分格式

本文建立了浅水波方程组的有限差分格式,该格式具有高阶精度和保持熵稳定性。首先,我们构造了一个二阶熵守恒格式,该格式满足给定熵对的熵等式,并能准确地保持湖泊静止稳态。主要思想是使源项的离散化与通量梯度项的离散化相匹配。其次,以具有合理熵守恒通量的二阶熵守恒格式为基础,实现了高阶熵守恒格式。第三,在已有的熵守恒格式上增加适当的耗散项,得到满足离散熵不等式的半离散熵稳定格式。特别是熵稳定格式可以避免熵守恒格式的振荡。其中,耗散项的建立采用根据熵变量构造的加权本质无振荡重构。最后,采用龙格–库塔方法进行时间离散化。文中给出了大量的数值算例,以验证该方法的高精度性及解决捕捉间断的能力。

一种新的不稳定导热有限差分格式

本文提出不稳定导热问题的一种新的有限差分格式。这种格式根据所取时间、空间步长及边界条件确定相应的温度加权因子。这样给出的差分格式与C-N格式计算量相同,但准确性却明显提高,并且是无条件稳定的。

分数阶平流-扩散方程的隐式有限差分格式的稳定性分析

文献[1]中给出了分数阶平流-扩散方程的隐式有限差分格式,并通过数值实验验证了方法的有效性。在文献[1]研究的基础上,进一步给出数值方法的稳定性分析,得到差分格式是无条件稳定的。

非定常Stokes方程一种基于POD方法的简化有限差分格式

特征正交分解(proper orthogonal decomposition,简记为POD)方法是一种可对偏微分方程的物理模型(如流体流动)做简化的技术.这种方法已经成功地用于对复杂系统模型降阶.该文推广应用POD方法,将POD方法应用于具有实际应用背景的非定常Stokes方程经典的有限差分格式,建立一种维数较低而精度足够高的简化差分格式,并给出简化差分格式解与经典差分格式解的误差估计.数值例子说明数值计算结果与理论结果相吻合.进一步表明基于POD方法的简化差分格式对求解非定常Stokes方程数值解是可行和有效的.

对流扩散方程的变步长摄动有限差分格式

摄动有限差分(PFD)方法是构造高精度差分格式的一种新方法。变步长摄动有限差分方法是等步长摄动有限差分方法的发展和推广。对需要局部加密网格的计算问题,变步长PFD格式不需要对自变量进行数学变换,且和等步长PFD格式一样,具有如下的共同特点:从变步长一阶迎风格式出发,通过把非微商项(对流系数和源项)作变步长摄动展开,展开幂级数系数通过消去摄动格式修正微分方程的截断误差项求出,由此获得高精度变步长PFD格式。该格式在一、二和三维情况下分别仅使用三、五和七个基点,且具有迎风性。文中利用变步长PFD格式对对流扩散反应模型方程,变系数方程及Burgers方程等进行了数值模拟,并与一阶迎风和二阶中心格式及其问题的精确解作了比较。数值试验表明,与一阶迎风和二阶中心格式相比,变步长PFD格式具有精度高,稳定性与收敛性好的特点。变步长PFD格式与等步长PFD格式相比,变步长PFD解在薄边界层型区域的分辨率得到了明显的提高。

基于紧致差分格式的高效时域有限差分算法

探讨一种基于紧致差分格式的高效时域有限差分算法(high-order compact-FDTD),该方法不仅提高计算精度,而且网格结点少、内存使用率和CPU时间大为降低.利用紧致格式FDTD方法实现无耗波导系统及光子晶体光纤中电磁波传播的数值模拟.通过计算实例验证算法的高效性.

双曲守恒型方程的二阶摄动有限差分格式

对双曲守恒型方程,将其一阶迎风格式空间差商的常系数摄动展开为时间步长和空间步长的幂级数,通过确定幂级数系数而获得二阶精度的摄动有限差分(PFD)格式。进而从双曲守恒型方程的通量分裂型一阶迎风格式出发,通过类似的摄动展开方法,获得空间精度为二阶的通量分裂形式的摄动有限差分(FPFD)格式。这两类格式保留了一阶守恒迎风格式的简洁结构形式,使用三节点即可达到二阶精度,又避免了三点二阶格式的非物理数值振荡。并将这两类格式推广应用到双曲守恒型方程组,最后通过模型方程和一维激波管流动的数值算例验证了格式的高精度、高分辨率性质。

有限差分-有限元混合方法的隐式格式研究

利用文［１］提出的有限差分—有限元混合方法，给出一种新的有限元隐式格式，它避免了在传统有限元方法中采用隐式格式时需要求解大型带状稀疏矩阵以及存储量大等问题，同时利用差分法中近似因式分解方法和ＰｕｌｉａｍＴＨ等人为避免块对角矩阵的求解提出的对角化技术，以期提高计算效率。

高精度有限差分法模拟Kelvin-Helmholtz不稳定性

WENO有限差分格式有较高的分辨精度,适合复杂流场的计算,在国际上被广泛采用。本文利用WENO有限差分格式求解2维守恒型欧拉方程,实现了对无粘流体中Kelvin-Helmholtz不稳定性的数值模拟。速度剪切方向采用周期边界条件;扰动增长方向采用嵌边出流边界条件,一个不稳定波长分布64个网格。数值模拟给出的扰动幅值线性增长率与线性稳定性分析给出的结果很好符合,显示了该格式的有效性和精度。数值模拟给出了清晰的密度等值线,表明该方法还具有较好的界面变形捕捉能力。

基于优化方法的时间-空间域隐格式有限差分算子确定方法

声波方程数值模拟构成了地震逆时偏移成像技术和全波形反演的基础。对于有限差分法而言,在满足一定的稳定性条件时,普遍存在着因网格化而形成的数值频散效应。如何有效地压制数值频散是有限差分方法研究的关键所在。为了进一步抑制数值频散,利用隐式有限差分比显式有限差分更能压制数值频散的特点,采用前人提出的新的有限差分模板(在保持相同精度的情况下增大了时间步长),应用信赖域优化方法在时间-空间域确定隐格式有限差分系数。频散分析和数值模拟试算的结果表明,这种新模板隐格式有限差分优化方法既提高了声波数值模拟精度又提高了计算效率。

地下水溶质运移定解问题差分格式的稳定性分析

溶质在地下水中的运移问题用对流扩散方程来描述，因而应用该方程求解定解问题时必须对数学模型进行稳定性条件分析，本文采用矩阵理论证明了地下水污染物运移问题迎风差分格式稳定性条件，其显式格式的稳定性条件为，隐式格式无条件稳定，并以承压含水层为例，结合目前计算能力和数据可视化极强的MATLAB软件，对所给出的稳定性条件进行了分析和验证，结果表明，在满足稳定性条件的前提下，其计算结果与实际相符合。

三维双曲方程非齐次边值问题的推广型LOD有限差分格式

针对三维非齐次双曲方程第一边值问题提出了一种新型的LOD有限差分格式,此格式能够将高维问题完全分解为一系列一维问题进行求解,克服了LOD格式源项难以分解、过渡层条件不易确定的缺陷.证明了该LOD有限差分格式按照离散L2模具有二阶收敛精度,与抛物型方程相比,源项的扰动达到了Δt4,从而使Δt的取法有更大的灵活性.

RLW方程基于混合有限元法的差分格式及其数值模拟

给出了RLW方程ut+auux-δuxxt=0的基于混合有限元法的一种最低阶的差分格式 ,并给出数值解的例子 ,与以往的处理RLW方程的有限差分法不同之处是该方法能同时求出流体的波峰及其流通量的近似解 ,而且得到的数值解具有很好的稳定性 .

二维双曲方程非齐次边值问题的推广型LOD有限差分格式

针对二维非齐次双曲方程第一边值问题提出了一种新型的LOD有限差分格式,此格式能够将高维问题完全分解为一系列一维问题进行求解,克服了LOD格式源项难以分解、过渡层条件不易确定的缺陷.证明了此种LOD有限差分格式按照离散L2模具有二阶收敛精度.数值算例表明计算效果良好.