首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到10条相似文献,搜索用时 15 毫秒
1.
针对结构自适应加密网格(SAMR)上扩散方程的求解,分析几种有限体格式的逼近性,同时设计和分析一种两层网格算法.首先,讨论一种常见的守恒型有限体格式,并给出网格加密区域和细化/粗化插值算子的条件;接着,通过在粗细界面附近引入辅助三角形单元,消除粗细界面处的非协调单元,设计了一种保对称有限体元(SFVE)格式,分析表明,该格式具有更好的逼近性,且对网格加密区域和插值算子的限制更弱;最后,为SFVE格式构造一种两层网格(TL)算法,理论分析和数值实验表明该算法的一致收敛性.  相似文献   

2.
针对交替方向显式离散格式,提出一个基于结构网格局部加密技术(SAMR)的求解流体力学方程组的自适应时间积分算法;基于该算法,在JAsMIN框架上研制多介质流体力学并行自适应数值模拟程序;在512个处理器上模拟惯性约束聚变中的二维内爆模型.数值模拟结果和并行性能分析显示了算法的正确性和并行实现的高效率.  相似文献   

3.
研究自适应Runge-Kutta间断Calerkin(RKDG)方法求解双曲守恒律方程组,并提出两种生成相容三角形网格的自适应算法.第一种算法适用于规则网格,实现简单、计算速度快.第二种算法基于非结构网格.设计一类基于间断界面的自适应网格加密策略,方法灵活高效.两种方法都具有令人满意的计算效果.而且降低了RKDG的计算量.  相似文献   

4.
张荣培  蔚喜军  崔霞  冯涛 《计算物理》2012,29(5):647-653
提出一种求解二维非平衡辐射扩散方程的数值方法.空间离散上采用加权间断Galerkin有限元方法,其中数值流量的构造采用一种新的加权平均;时间离散上采用隐-显积分因子方法,将扩散系数线性化,然后用积分因子方法求解间断Galerkin方法离散后的非线性常微分方程组.数值试验中在非结构网格上求解了多介质的辐射扩散方程.结果表明:对于强非线性和强耦合的非线性扩散方程组,该方法是一种非常有效的数值算法.  相似文献   

5.
为提高油藏数值模拟算法的计算效率,在求解单向稳态渗流的有限分析算法基础上,构建二维非均匀多孔介质中不可压两相渗流的有限分析算法.算法中,网格界面上的平均渗透率不是简单地取为相邻网格渗透率的调和平均值,而是通过奇点邻域解析解积分求得.相比于传统的数值算法,有限分析算法随着网格的加密,能够很快地收敛(仅需将原始网格细分至2×2或3×3),并且其计算精度和收敛性不依赖于介质的非均匀强度,从而计算效率得到提高.  相似文献   

6.
结构网格自适应程序需要使用矩形区域求差集算法计算网格层间数据依赖关系和网格层嵌套关系.原有的矩形区域求差集算法时间复杂度较高,成为该类应用大规模并行计算可扩展性能瓶颈.本文利用分而治之的方法,构造近似线性时间复杂度的矩形区域求差集快速算法,并利用区域分解实现该算法的并行计算.分别针对规则矩形区域和多层自适应网格的非规则矩形区域求差集问题,验证该算法的效率.结果表明,该算法具有近似线性计算复杂度,对于大规模计算问题,加速效果显著.  相似文献   

7.
二维保单调保守恒插值算子   总被引:1,自引:0,他引:1  
基于一个一维保单调保守恒插值算子,利用不完全双二次插值提出一个二维保单调保守恒插值算子.从插值逼近角度,通过几个数值实验验证该插值算子有效.用得到的二维插值算子作为结构网格自适应加密(structured adaptive mesh refinement,SAMR)算法中的细化插值算子,求解几个二维Euler方程数值例子,结果表明,提出的二维插值算子有效.  相似文献   

8.
面向结构网格并行应用的一类快速通信算法   总被引:1,自引:0,他引:1  
通信算法需要在相邻子区域间传输数据.通过求解子区域间的相交问题可以寻找相邻区域.针对子区域的求交问题,基于区间树,结合结构网格应用的特点,构造近似线性时间复杂度的算法.数值实验表明该算法具有较高的计算效率和可扩展性,能够支持百万量级矩形子区域的并行计算.  相似文献   

9.
建立求解稳态对流-扩散方程的一种稳定、高效的无单元Galerkin方法.该方法计算积分时采用基于局部Taylor展开的节点积分,并根据对流占优的程度对积分节点进行自适应迎风偏移.与传统的使用稳定化的无单元Galerkin方法相比,该方法是一种不依赖于背景网格积分的纯无网格方法,具有更好的稳定性和较高的计算效率,其程序实施更为简便.  相似文献   

10.
分数阶混沌系统的Adomian分解法求解及其复杂性分析   总被引:3,自引:0,他引:3       下载免费PDF全文
贺少波  孙克辉  王会海 《物理学报》2014,63(3):30502-030502
根据分数阶微分定义,采用Adomian分解算法,研究了分数阶简化Lorenz系统的数值解.研究发现,该算法与预估-校正算法相比,求解结果更准确,所耗计算资源和内存资源更少,求解整数阶系统时较Runge-Kutta算法更准确;利用Adomian算法得到的分数阶简化Lorenz系统出现混沌的最小阶数为1.35,比利用预估-校正算法得到的最小阶2.79更小.采用相图、分岔图分析了该系统的动力学特性,基于谱熵算法(SE)和C0算法分析了该系统的复杂度.结果表明,复杂度结果和分岔图一致,说明系统的复杂度同样能反映出系统动力学特性;复杂度随阶数q的增加呈总体减小的趋势,而混沌态时系统参数c变化对系统复杂度影响不大.为分数阶混沌系统应用于信息加密、保密通信领域提供了理论与实验依据.  相似文献   

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

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