首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到18条相似文献,搜索用时 109 毫秒
1.
基于非均匀网格,提出了一种求解一维定常对流扩散反应方程的高精度紧致差分格式。首先采用坐标变换方法将原方程由物理空间的非均匀网格转换为计算空间的均匀网格,然后给出一阶导数和二阶导数在均匀网格上的中心差分逼近式,并结合变换后的方程,得到了定常对流扩散反应方程具有四阶精度的紧致差分格式。最后,通过数值算例验证了该方法的精确性和高分辨率的特点。数值实验结果表明,对于所研究问题,该方法较不进行坐标变换而直接在物理域上建立的非均匀网格上的高阶紧致格式具有更高精度。  相似文献   

2.
本文主要研究一维四阶双曲方程初边值问题.首先通过引入一个中间函数将其转化为二阶方程组,然后对方程中的空间导数项采用四阶紧致差分格式离散,时间导数项采用二阶中心差分格式离散,构造出问题的隐式紧致差分格式.数值算例表明该格式具有较好的计算效果.  相似文献   

3.
针对时间分数阶扩散方程,提出了一种新的隐式差分方法,其中空间导数采用中心差分方法离散.对于时间分数阶导数,将Caputo分数阶导数转化为Riemman-Liouville分数阶导数后,写成Hadamard有限部分积分,再用分段二次多项式对该有限积分部分逼近,由此推导出Caputo分数阶导数的3-α阶离散方法,从而得到无条件稳定的和收敛的分数阶扩散方程的隐式差分格式.数值实验验证该隐式差分格式的有效性.  相似文献   

4.
研究了一类含弱非线性的改进型Boussinesq水波方程,在非交错网格下,利用有限差分法建立了混合四阶Adams-Bashforth-Moulton的预报校正格式的波浪数值模型。在数值模型中,关于空间一阶导数差分格式采用四阶精度、二阶导数差分格式采用二阶精度。针对波浪的一维、二维传播变形问题进行了数值计算,并通过与相关实验结果对比分析考察了该数值模型的适用性。  相似文献   

5.
非均匀介质地震波传播交错网格高阶有限差分法模拟   总被引:10,自引:0,他引:10  
采用常规的二阶声波方程有限差分方法对于非均匀介质进行了数值模拟时,其数值模拟精度较低。而采用一阶双曲型标量波动方程,则无须对介质的弹性常数进行空间求导。根据Taylor级数展开式,推导出了交错网格一阶空间导数的任意偶数阶精度展开式和相应差分系数计算式以及一阶双曲型标量波动方程交错网格任意偶数阶精度差分格式,并给出了该差分算法的稳定性条件。用该差分算法对均匀介质模型、非均匀介质模型和Mannousi模型进行了数值模拟试验,并与伪谱法进行了对比。结果表明,一阶双曲型标量波动方程交错网格高阶差分法的模拟精度与伪谱法的精度非常接近,计算效率高,且适合于模拟非均匀介质、复杂构造和复杂地质体的地震波场。  相似文献   

6.
本文研究了带有阻尼项的四阶梁振动方程初边值问题,基于紧致差分方法,给出了数值求解该问题的四种高阶紧致差分格式.对方程中的一阶和二阶时间导数项采用中心差分离散,对四阶空间导数项分别采用五点、七点和带紧致的五点、七点四种方法进行离散,得到四种高阶紧致差分格式,这四种格式均在时间方向达到二阶精度,在空间方向分别达到二阶、四阶、四阶和六阶精度.最后利用数值算例验证了四种格式的精度阶与理论结果一致.本文相对于之前的研究,对弹性梁的振动增加了阻尼因素,因此也更加适合对实际问题的数值计算.  相似文献   

7.
变系数分数阶反应-扩散方程的数值解法   总被引:1,自引:0,他引:1  
考虑了变系数分数阶反应一扩散方程,将一阶的时间偏导数和二阶的空间偏导数分别用Caputo分数阶导数和Riemann-Liouville分数阶导数替换,利用L1算法和G算法对方程的变系数分数阶导数进行适当的离散,给出了该方程的一种计算有效的隐式差分格式,并证明了这个差分格式是无条件稳定和无条件收敛的,且具有o(τ+h)收敛阶.最后用数值例子说明差分格式是有效的.  相似文献   

8.
时间分数阶对流-扩散方程可以用来模拟由传统的对流-扩散方程演变而来的反常扩散方程.本文针对一类时间分数阶对流-扩散方程提出了一个新的隐式差分格式,时间分数阶导数采用直接离散,空间导数采用中心差分格式离散,讨论了差分解的存在唯一性,并利用能量范数证明了该格式的无条件稳定性、收敛性,分析了收敛阶.数值试验验证了该格式的有效性.  相似文献   

9.
研究二维有限域上的空间分数阶扩散方程的数值解法,通过移位的Grunwald公式对空间分数阶导数进行离散,得到Euler隐式差分格式。利用傅里叶变换理论证明了交替差分格式的一致性。  相似文献   

10.
为进一步发展更有效的湍流模型,提高对高负荷跨声速压气机流动的数值计算精度,以侧偏平均脉动速度和应变率张量为基础构造了新的侧偏平均脉动速度(partial average fluctuation velocity, PAFV)湍流模型。在使用Python计算机语言进行编程计算过程中,将流动控制方程在一般曲线坐标系下变换后,采用有限差分方法进行空间离散。方程对流项采用Steger-Warming通量向量分裂法,扩散项采用中心差分格式离散,时间导数项采用具有TVD(total variation diminishing)性质的二阶Runge-Kutta法离散。使用该湍流模型对德国航空航天研究试验院的L030-4叶栅亚声速和超声速流动进行了计算,在进行网格无关性研究基础上,获得了流场压力等值线云图、马赫数等值线云图,压力系数分布等结果。计算显示侧偏平均脉动速度在激波附近区域受到明显抑制,避免了激波附近区域出现过大的湍流黏性。数值计算结果与实验符合比较好,初步验证了PAFV湍流模型可以有效地模拟该类型的叶栅流动。  相似文献   

11.
为了既能精确地拟合不规则区域边界,又使有限差分的计算能在规则的网格上进行,在长江口南支河段非恒定流的平面二维数值模拟中,采用了近年来得到广泛注意的边界拟合坐标法,为了提高计算的稳定性,本文引用了二维人工粘性项和人工光滑的方法。计算结果与实测资料相当吻合。计算结果还表明,在同样的网格密度下,使用边界拟合坐标的差分计算的CPU仅比使用昔通网格的差分计算略有增加。所以,在二维水流计算中,这是一种值得推荐的方法。  相似文献   

12.
借助空间坐标变换,把移动区域模型转化为固定区域模型,通过构造显式和隐式两种有限差分格式求解热源周期振荡条件下的一维融化问题.对所构造的两种差分格式分别研究它们的数值稳定性,比较它们的计算量和计算效率;应用这两种差分格式分别数值模拟融化过程中移动边界的运动及液态介质内温度场的分布.数值实验结果表明,这两种差分格式的数值结果吻合得非常好,而隐式差分格式的计算效率要明显优于显式差分格式.  相似文献   

13.
有限长轴承非稳态油膜力自由边界问题数值解法   总被引:1,自引:0,他引:1  
讨论了不可压缩流体润滑的动载径向滑动有限长轴承油膜压力分布的自由移动边界问题的有限元计算方法。将自由边界问题转化为全球的具有不等式约束的微分方程边值问题及相应的二次规划问题,借助于牛顿非光滑算法迭代求解,首次给出了有限长轴承真实的油膜压力分布。本算法可以用紧缩式的存储方式,既节省了存储空间又减少了计算量。计算结果与通常的全周油膜理论的结果进行了比较。  相似文献   

14.
顾及2套坐标误差的三维坐标变换方法   总被引:2,自引:1,他引:1  
探讨顾及2套坐标误差改正数的三维坐标转换模型,以2套坐标改正数的加权平方和最小为准则导出了转换参数的解算公式;并根据公共点的坐标改正数以及公共点与待转换点间的协方差阵计算待转换点的坐标改正数并对其改正,从而求得近似无缝的坐标转换结果.试验表明:当公共点与转换点坐标相关性较大时,该方法能显著地提高坐标转换精度.  相似文献   

15.
采用显格式与隐格式交替使用的方法,针对一阶线性双曲方程组提出了一种隐-显迎风差分格式.它综合了隐格式与显格式的优点,具有稳定性好、计算简便的特性.数值计算结果表明,这种格式是实用的.  相似文献   

16.
提出了一种在极坐标系下快速推进带电粒子的等离子体粒子模拟算法,该算法在每个时间步长内对带电粒子位置附近的区域采用局部直角坐标系近似,并通过坐标变换,得到粒子运动方程求解的显式格式,从而避免了极坐标系下隐式格式直接迭代求解所带来的容许误差,并提高了计算效率.数值实验表明:显式粒子推进算法的计算速度比容许误差取10-3时的隐式算法提高了约1倍.该算法可用于轴向均匀等离子体器件和高功率微波器件的2 5维粒子模拟.  相似文献   

17.
采用 L U- AUSML W算法 ,对非线性双曲型 Euler方程进行求解 .该算法综合了 L U- SGS方法与 AUSMPW格式的优点 .为了验证 L U- AUSML W算法的有效性 ,在曲线坐标系中 ,对亚音速、跨音速和超音速流动进行了数值模拟 .数值试验的计算结果与文献计算结果相符很好  相似文献   

18.
采用高精度紧致有限差分-Fourie谱杂交的方法直接数值模拟了三维不可压缩的NavierStokes方程,该算法的时间离散采三阶精度混合显隐分裂格式,空间离散则结合Fourie谱方法及高精度紧天才有限差分逼近,该方法与普通的有限差分格式相比,具有很高的逼的精度及波数分辨率,针对三维平面槽道流的情况,应用该算法,直接数值模拟了三维T-S波在平面槽道流的传播问题,计算结果与流动稳定性分析结果吻合一致。  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号