边坡稳定分析的自适应性局部网格加密解法 被引量：1

1

作者 周桂云 李同春 《岩土力学》 EI CAS CSCD 北大核心 2007年第S1期398-402,共5页

常规的有限元强度折减法在确定边坡塑性区时,不能直观地反映出塑性区的形状和滑裂面位置,计算结果精度低。在塑性区范围内进行网格局部加密,不仅使安全系数的计算精度大幅度提高,而且可求出滑裂面的形状和位置。针对求解塑性区搜索滑裂... 常规的有限元强度折减法在确定边坡塑性区时,不能直观地反映出塑性区的形状和滑裂面位置,计算结果精度低。在塑性区范围内进行网格局部加密,不仅使安全系数的计算精度大幅度提高,而且可求出滑裂面的形状和位置。针对求解塑性区搜索滑裂面这一具体问题,提出塑性应变界限值网格加密准则,采用自适应性网格局部加密技术求解边坡的滑裂面和稳定安全系数。数值计算结果表明,所提出的加密准则能正确地确定滑裂面的形状和位置,计算所得稳定安全系数和理论稳定安全系数仅相差0.2%。 展开更多

关键词 边坡稳定 自适应性局部网格加密 有限元 强度折减 安全系数 滑裂面

在线阅读 下载PDF 职称材料

基于局部加密网格的LBM程序负载均衡方法研究

2

作者 黄晨希 李家辉 +2 位作者 颜辉 钟英 卢宇彤 《计算机科学》 北大核心 2025年第5期101-108,共8页

基于局部加密网格的格子玻尔兹曼方法(LBM),在大规模的非稳态计算流体力学问题中有着越来越广泛的使用。局部网格加密方法虽然可以有效减少计算量,但是也带来了严重的负载均衡问题。尤其是在大规模并行计算中,负载均衡算法的选择对整体... 基于局部加密网格的格子玻尔兹曼方法(LBM),在大规模的非稳态计算流体力学问题中有着越来越广泛的使用。局部网格加密方法虽然可以有效减少计算量,但是也带来了严重的负载均衡问题。尤其是在大规模并行计算中,负载均衡算法的选择对整体的计算效率有重要的影响。使用自研的LBM程序,简单对比了开源框架Palabos的实现结果,然后探讨了基于局部加密网格格子玻尔兹曼方法的静态负载均衡算法。程序中分别基于全局负载和分层负载,使用贪心算法和空间填充曲线算法实现了多种负载均衡策略,提出了计算节点负载均衡优化,并使用小球绕流算例和DrivAer算例进行测试和比较。根据计算结果可知,基于分层负载策略的性能比基于全局负载策略的性能要好;同时,贪心算法和空间填充曲线算法对比各有优劣,前者有较好的可扩展性,而后者在进程数较少时效率较高;最后,基于分层负载,结合贪心算法和空间填充曲线算法实现了混合算法,在进程数较多时获得了最佳性能。 展开更多

关键词 格子玻尔兹曼方法 负载均衡算法 局部加密网格 并行计算 节点负载均衡

在线阅读 下载PDF 职称材料

基于叉树法的自适应有限元局部网格加密研究 被引量：11

3

作者 魏红宁 周本宽 《计算力学学报》 EI CAS CSCD 1997年第3期324-331,共8页

本文以叉树法网格生成技术为基础，提出了效率高，易实现的局部网格加密步骤，同时设计了相互匹配的网格数据结构和树数据结构以及局部网格修正算法用以支持本文所提出的局部网格加密步骤。最后，给出了局部网格加密例子，用以证明本文... 本文以叉树法网格生成技术为基础，提出了效率高，易实现的局部网格加密步骤，同时设计了相互匹配的网格数据结构和树数据结构以及局部网格修正算法用以支持本文所提出的局部网格加密步骤。最后，给出了局部网格加密例子，用以证明本文所提出的方法可以很容易地生成任意密度分布的网格，这对于实施自适应ｈ－方案是十分有利的。 展开更多

关键词 有限元 自适应 局部网格加密 数据结构

在线阅读 下载PDF 职称材料

自适应网格局部加密方法在可压缩流场模拟中的应用 被引量：2

4

作者 沈建雄 黄兰洁 《空气动力学学报》 CSCD 北大核心 1991年第3期330-337,共8页

本文基于M.J.Berger和J.Oliger发展的自适应网格局部加密算法,它采用多个组合网格的思想,用于求解双曲型方程组。它结合显式有限差分格式,采用Richardson外推技术自动进行局部截断误差估计,对精度低的区域产生新细网格局部加密,或去除... 本文基于M.J.Berger和J.Oliger发展的自适应网格局部加密算法,它采用多个组合网格的思想,用于求解双曲型方程组。它结合显式有限差分格式,采用Richardson外推技术自动进行局部截断误差估计,对精度低的区域产生新细网格局部加密,或去除无需再加密的旧细网格,以最小的运算量达到给定的精度要求。网格可一层层加密下去,按层次覆盖,每个网格是有任意方向的矩形均匀网格。这套算法独立于求解所用的差分格式,很容易和各种格式结合。 我们从M.J.Berger和J.Oliger程序出发,参考M.J.Berger和P.Colella,对守恒型方程组,实现了粗细网格交接面上数值通量守恒,从而能计算间断解。对网格边界处理、误差估计和网格重新生成等做了一些改进,并且可以灵活地在任何感兴趣区域指定加密。 我们用不定常Euler方程组计算了带前台阶的二维管道中Mach-3流动,采用MacCormack显式附加人工粘性差分格式。共用了四层网格结构:Δ=1/10、1/20、1/80和1/320。其中基本网格是Δ=1/10;到台阶上第一次激波反射区域从Δ=1/20加密至Δ=1/80,计算到定常所需运算量仅是同等均匀网格△=1/80的28％;我们还用Δ=1/320网格对台阶拐角和切向间断处指定加密,结果表明可明显改善台阶上的激波Mach反射,切向间断清晰可辨,所需的计算量是上述均匀网格Δ=1/ 展开更多

关键词 自适应网格 局部网格 加密

在线阅读 下载PDF 职称材料

基于局部加密非结构网格的海洋可控源电磁法三维有限元正演 被引量：12

5

作者 叶益信 李予国 +2 位作者 刘颖 李刚 杨海燕 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2016年第12期4747-4758,共12页

本文基于非结构网格实现了海洋可控源电磁法三维有限元正演模拟.该算法采用完全非结构网格剖分,可以模拟任意起伏地形和复杂地电模型.为了避免场源的奇异性,采用一次场/二次场分解算法,一次场由基于Schelkunoff势函数的一维解析公式得到... 本文基于非结构网格实现了海洋可控源电磁法三维有限元正演模拟.该算法采用完全非结构网格剖分,可以模拟任意起伏地形和复杂地电模型.为了避免场源的奇异性,采用一次场/二次场分解算法,一次场由基于Schelkunoff势函数的一维解析公式得到.为了提高算法的精度和效率,采用对测点附近单元和异常体区域进行体积约束加密的方法,实现了非结构网格的局部加密.一、二维模型计算和分析表明,本文采用的局部加密方法能够明显地改善算法的精度,最大相对误差基本在1%以内.对三维模型计算及对比分析,说明了该算法对三维可控源电磁正演的实用性.复杂海底地形模型的正演模拟表明,海底地形对电磁场的影响很大,在进行海洋可控源电磁资料解释时,地形的影响有必要考虑在内. 展开更多

关键词 海洋可控源电磁 非结构网格 局部加密 有限元

在线阅读 下载PDF 职称材料

局部网格加密技术在混凝土裂缝扩展模拟中的应用 被引量：6

6

作者 韦未 李同春 +1 位作者 牛志伟 屈寒飞 《华南农业大学学报》 CAS CSCD 北大核心 2007年第4期112-116,共5页

采用有限元分析混凝土裂缝问题时,常遇到无法一步到位地确定采用什么样的剖分方案使有限元计算能到达精度要求.针对这一问题,文章提出一种局部网格加密技术,网格根据计算精度要求以及混凝土开裂位置对原有网格进行局部加密,并且利用非... 采用有限元分析混凝土裂缝问题时,常遇到无法一步到位地确定采用什么样的剖分方案使有限元计算能到达精度要求.针对这一问题,文章提出一种局部网格加密技术,网格根据计算精度要求以及混凝土开裂位置对原有网格进行局部加密,并且利用非协调网格算法进行计算,通过2个算例证明,局部网格加密在保证达到精度要求的条件下,将会大大减少计算机计算量. 展开更多

关键词 混凝土 局部网格加密 非协调网格

在线阅读 下载PDF 职称材料

基于局部加密非结构化网格的三维电阻率法有限元数值模拟 被引量：42

7

作者 任政勇 汤井田 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2009年第10期2627-2634,共8页

总结了目前常用的结构化网格及其局限性,分析了非结构化网格对复杂地质体边界的适应性和Delaunay三角化算法及其网格加密策略,提出了一种新的局部节点加密方法,实现了对复杂模型的完全非结构化四面体全自动剖分,给出了三维直流电阻率模... 总结了目前常用的结构化网格及其局限性,分析了非结构化网格对复杂地质体边界的适应性和Delaunay三角化算法及其网格加密策略,提出了一种新的局部节点加密方法,实现了对复杂模型的完全非结构化四面体全自动剖分,给出了三维直流电阻率模拟中四面体网格的质量评价标准和最优指标.计算和分析表明,数值解在点电流源及附近的奇异区精度最低,网格加密策略可以有效地减少其影响,极大地提高数值解的精度.本文提出的局部节点加密策略计算量最小,对精度的改善也优于经典方法.在精度要求苛刻或模型十分复杂时,局部体积加密策略和二次单元是高精度模拟的可靠保证. 展开更多

关键词 非结构化网格 DELAUNAY算法 网格局部加密 直流电阻率 二次单元

在线阅读 下载PDF 职称材料

MRT-LBM三重网格局部加密算法研究 被引量：5

8

作者 曲延鹏 张坤 陈颂英 《北京理工大学学报》 EI CAS CSCD 北大核心 2018年第5期441-448,共8页

应用MRT-LBM计算流场中运动粒子上的力和力矩时,可使用三重网格局部加密方法来提高计算结果的精度.将流场计算区域划分为互不重合的粗糙区、过渡区和加密区,相邻区域界面重合节点上的物理参数如密度、速度和应力等保持连续,不同区域的... 应用MRT-LBM计算流场中运动粒子上的力和力矩时,可使用三重网格局部加密方法来提高计算结果的精度.将流场计算区域划分为互不重合的粗糙区、过渡区和加密区,相邻区域界面重合节点上的物理参数如密度、速度和应力等保持连续,不同区域的分布函数通过过渡区边界点进行传递,运算程序采用碰撞-迁移算法,给出了数据的初始化过程和加密的网格结构.通过对含非对称放置的粒子Couette流动,分别采用标准Boltzmann、全场加密、局部二重加密和三重加密的4种网格进行计算.数值计算结果表明,三重网格加密算法减少了计算结果的波动.对方腔流顶盖下方左右奇异角落处应用三重局部网格加密,结果显示,对雷诺数1 000的方腔流动,沿腔中心线的速度分布与经典文献结果对比效果良好,压力轮廓图的噪声明显降低,应力振荡明显减少.模拟结果证实了所建立的局部加密方法的有效性. 展开更多

关键词 格子BOLTZMANN方法 多松弛时间模型 局部网格加密 COUETTE流动 方腔流动

在线阅读 下载PDF 职称材料

二维有限元网格的局部加密方法 被引量：7

9

作者 黄晓东 杜群贵 《华南理工大学学报（自然科学版）》 EI CAS CSCD 北大核心 2004年第12期56-60,共5页

设计了关键点、圆域和矩形域三种有限元网格局部加密方式 .在所构造的局部加密嵌入面域内 ,利用基于Delaunay剖分的动态节点 -单元一体化生成算法生成局部加密网格 ,并将其嵌入原网格 ,以实现分析对象关键区域有限元网格的局部加密 .算... 设计了关键点、圆域和矩形域三种有限元网格局部加密方式 .在所构造的局部加密嵌入面域内 ,利用基于Delaunay剖分的动态节点 -单元一体化生成算法生成局部加密网格 ,并将其嵌入原网格 ,以实现分析对象关键区域有限元网格的局部加密 .算例表明 ,该方法生成的网格质量良好 ,符合有限元分析要求 . 展开更多

关键词 有限元 网格生成 局部加密 DELAUNAY剖分

在线阅读 下载PDF 职称材料

基于局部加密等级网格的2.5D直流电法有限元模拟 被引量：3

10

作者 胡宏伶 肖晓 +2 位作者 潘克家 汤井田 谢维 《中南大学学报（自然科学版）》 EI CAS CSCD 北大核心 2014年第7期2259-2267,共9页

合理截取半圆形计算区域,采取局部加密的?-等级网格,结合对称行索引存贮格式(CSR)及并行稀疏直接求解器PARDISO,提出一种高效、高精度的2.5D直流电阻率法有限元正演方案,并编制相应的Fortran程序,对具有解析解的3个典型地电模型进行计... 合理截取半圆形计算区域,采取局部加密的?-等级网格,结合对称行索引存贮格式(CSR)及并行稀疏直接求解器PARDISO,提出一种高效、高精度的2.5D直流电阻率法有限元正演方案,并编制相应的Fortran程序,对具有解析解的3个典型地电模型进行计算与分析。研究结果表明:圆形截断边界不仅便于在径向方向上采取?-等级网格剖分,而且能大大简化有限元模拟中单元刚度矩阵的计算;结构化的等级网格避开了通常非结构化网格有限元计算时繁琐的网格剖分及总体刚度阵的集成过程,且能在不增大问题规模的前提下,显著提高2.5D直流电法正演源点附近的模拟精度;Intel MKL的PARDISO求解器能在普通PC机上5 s内求解电法正演有限元离散得到的100万阶稀疏线性方程组,可广泛用于各种地球物理正演问题。 展开更多

关键词 直流电阻率 等级网格 有限元 局部加密 并行稀疏直接求解器 行索引存贮格式

在线阅读 下载PDF 职称材料

基于加密物理片的数值流形法中局部网格加密 被引量：2

11

作者 刘治军 郑宏 +2 位作者 董苇 葛修润 孙冠华 《岩土力学》 EI CAS CSCD 北大核心 2017年第4期1211-1217,1225,共8页

数值流形法的数学网格不需要适应求解域内各种不连续界面和边界,因此,总可以用规则的结构化网格建立数学覆盖。但对于大多数问题,在整个求解域上布置统一密度的网格显得浪费。因此,需要研究在结构化网格上实施局部加密,建立了加密物理... 数值流形法的数学网格不需要适应求解域内各种不连续界面和边界,因此,总可以用规则的结构化网格建立数学覆盖。但对于大多数问题,在整个求解域上布置统一密度的网格显得浪费。因此,需要研究在结构化网格上实施局部加密,建立了加密物理片的方法用以解决这一问题。具体的实施过程中,确定了需要加密的区域后,先在这些区域布置规则的精细网格,然后找到这些区域中包含的原始网格中的物理片,用精细网格上建立的插值代替被加密物理片上的局部近似,从而提高了局部近似的阶次。数值算例结果表明,该方法收敛性良好。另外,如果所有物理片上的局部近似都采用0阶多项式(常数),那么将会得到正定的刚度矩阵。 展开更多

关键词 数值流形法 加密物理片 结构化网格 局部加密

在线阅读 下载PDF 职称材料

改进的栅格法三维六面体网格局部加密算法 被引量：2

12

作者 孙劲光 周勃 《计算机应用与软件》 CSCD 2016年第1期49-52,共4页

有限元网格的相容性及质量是网格划分中必须考虑的问题。针对现有加密模板生成的新增单元数目较多,网格质量不佳,容易产生大范围加密外延,且在使用时加密情况考虑较为简单的问题,在栅格法生成六面体网格的基础上提出了一套基于8分法的... 有限元网格的相容性及质量是网格划分中必须考虑的问题。针对现有加密模板生成的新增单元数目较多,网格质量不佳,容易产生大范围加密外延,且在使用时加密情况考虑较为简单的问题,在栅格法生成六面体网格的基础上提出了一套基于8分法的加密模板,改进了面模板中新增小网格质量较差的问题,增加了拐角模板。为了使得所有加密情况均能直接加密,给出加密模板与网格单元中加密节点所有可能的排列情况之间的对应法则,并给出加密模板的集体应用方法。通过实验结果及对比分析,该套模板能够减少新增单元数量,改善新增单元的质量,减少加密外延现象。 展开更多

关键词 局部加密 网格质量优化 三维六面体网格 加密模板 栅格法 自适应

在线阅读 下载PDF 职称材料

气泡堆积法生成局部加密非结构化网格 被引量：2

13

作者 武利龙 陈斌 《西安交通大学学报》 EI CAS CSCD 北大核心 2009年第9期19-22,共4页

对生成非结构化网格的气泡堆积法(BPM)进行改进,提出了一种生成局部加密非结构化网格的算法.该算法通过向多边形顶点添加任意大小气泡以及在内部设置点源,并利用反距离加权平均法传递节点信息,简单地实现了局部加密/稀疏网格的生成.同时... 对生成非结构化网格的气泡堆积法(BPM)进行改进,提出了一种生成局部加密非结构化网格的算法.该算法通过向多边形顶点添加任意大小气泡以及在内部设置点源,并利用反距离加权平均法传递节点信息,简单地实现了局部加密/稀疏网格的生成.同时,还提出了气泡数量的控制方法,可对过于稠密或者疏松的气泡群进行气泡的删除和添加,使所有的气泡更好地与计算区域相匹配,消除了不规则三角形网格的生成,提高了生成网格的质量.用改进的BPM算法对矩形区域生成了边界加密、中心稀疏的非结构化网格,在雷诺数为1 000的工况下对方腔顶盖驱动流进行了数值模拟,与基准解对比非常吻合,证明生成的网格质量较好. 展开更多

关键词 气泡堆积法 非结构化网格 局部加密

在线阅读 下载PDF 职称材料

局部网格粗化与加密技术在大底水油藏数值模拟中的应用 被引量：5

14

作者 刘广天 李保振 《科学技术与工程》 北大核心 2012年第13期3207-3210,共4页

大型底水油藏数值模拟模型的网格数量庞大,其中底水区域的网格数目占绝大部分,占用很大数值计算量。同时,对于储层较厚的底水油藏,常规的网格系统通常无法准确描述油水分布及开发动态。针对大型底水油藏数值模拟中存在的问题,提出了底... 大型底水油藏数值模拟模型的网格数量庞大,其中底水区域的网格数目占绝大部分,占用很大数值计算量。同时,对于储层较厚的底水油藏,常规的网格系统通常无法准确描述油水分布及开发动态。针对大型底水油藏数值模拟中存在的问题,提出了底水区域实施适当网格粗化,生产井附近区域适当加密的处理方法。采用此方法既能保证计算精度,同时减少了模型网格数,大幅提高模拟计算效率,同时为该类油藏生产井的射孔位置、厚度等指标的设计优化提供了技术保障。 展开更多

关键词 大底水油藏 数值模拟 局部网格加密 网格粗化

在线阅读 下载PDF 职称材料

多孔介质中混溶驱动问题的时空局部网格加密有限差分格式 被引量：2

15

作者 刘伟 《工程数学学报》 CSCD 北大核心 2006年第1期139-146,共8页

给出了多孔介质中一维混溶驱动问题在时间和空间上进行局部网格加密的有限差分格式,压力方程采用中心差分格式近似,饱和度方程采用修正迎风格式,且在交界面上采用线性插值,并利用极大值原理给出了误差估计。最后给出了数值算例。

关键词 耦合方程组 局部网格加密 混溶驱动 误差估计

在线阅读 下载PDF 职称材料

三维 CAD 曲面四边形网格剖分及局部加密技术 被引量：1

16

作者 李华 《大连理工大学学报》 CAS CSCD 北大核心 1998年第6期738-741,共4页

利用ＣＡＤ及ＣＡＧＤ系统中三维自由型曲面造型的数学方法，根据曲面上网格剖分的密度要求及曲率变化，先在定义曲面的参数平面域上生成密度可变的全四边形网格，然后利用曲面方程将参数域上的网格映射到曲面上．实现了任意自由型曲面... 利用ＣＡＤ及ＣＡＧＤ系统中三维自由型曲面造型的数学方法，根据曲面上网格剖分的密度要求及曲率变化，先在定义曲面的参数平面域上生成密度可变的全四边形网格，然后利用曲面方程将参数域上的网格映射到曲面上．实现了任意自由型曲面密度可变的全四边形网格剖分．为有限元分析系统与ＣＡＤ几何造型系统的集成提供了一种高效、可靠的前处理手段．输入数据是ＣＡＤ或ＣＡＧＤ系统中曲面的控制顶点或插值边界．该方法运算速度快，效率高，生成单元的质量好． 展开更多

关键词 有限元 四边形 曲面 局部加密 CAD 网格划分

在线阅读 下载PDF 职称材料

多孔介质中两相不可压缩流体的时间局部网格加密有限差分格式

17

作者 刘伟 《高校应用数学学报（A辑）》 CSCD 北大核心 2005年第2期170-176,共7页

引入Charent压力变量,对于多孔介质中两相不可压缩流体的非混溶驱动问题,其模型表现为耦合的非线性偏微分方程组,一个是压力方程,另一个为饱和度方程.文中考虑一维问题且假定达西速度u已知,建立了在时间上进行局部加密的有限差分格式,... 引入Charent压力变量,对于多孔介质中两相不可压缩流体的非混溶驱动问题,其模型表现为耦合的非线性偏微分方程组,一个是压力方程,另一个为饱和度方程.文中考虑一维问题且假定达西速度u已知,建立了在时间上进行局部加密的有限差分格式,给出了饱和度的最大模误差估计.最后给出了数值算例. 展开更多

关键词 局部网格加密 非混溶驱动 误差估计

在线阅读 下载PDF 职称材料

局部网格加密数模技术研究剩余油分布 被引量：7

18

作者 丁增勇 李微 +3 位作者 刘建华 田礼妍 廖家彬 张德海 《油气地质与采收率》 CAS CSCD 2001年第2期43-46,共4页

针对稠油油藏热采开发的特点 ,采用在吞吐井周围进行局部加密的方法 ,既兼顾稠油油藏热采特点 ,保证了拟合计算的合理性 ,又减少网格节点数 ,节省计算时间 ,从而实现了锦 90块兴I组油层全块的模拟。从宏观上和整体上查明了锦 90块兴I组... 针对稠油油藏热采开发的特点 ,采用在吞吐井周围进行局部加密的方法 ,既兼顾稠油油藏热采特点 ,保证了拟合计算的合理性 ,又减少网格节点数 ,节省计算时间 ,从而实现了锦 90块兴I组油层全块的模拟。从宏观上和整体上查明了锦 90块兴I组油层的目前油层动用状况和剩余油分布规律 。 展开更多

关键词 稠油油藏 蒸汽吞吐 数据模型 局部网格加密 剩余油分布

在线阅读 下载PDF 职称材料

数值流形法中基于适合分析T样条的局部网格加密算法 被引量：2

19

作者 刘登学 张友良 +2 位作者 丁秀丽 黄书岭 裴启涛 《岩土力学》 EI CAS CSCD 北大核心 2019年第4期1584-1595,共12页

数值流形方法中一般采用有限元网格或规则网格作为其数学覆盖系统,而规则的网格突出的优点是不需要适应求解域的边界和各种不连续面。采用规则的矩形网格作为数值流形方法中的数学网格,并借助适合分析的T样条实现了数值流形方法中的局... 数值流形方法中一般采用有限元网格或规则网格作为其数学覆盖系统,而规则的网格突出的优点是不需要适应求解域的边界和各种不连续面。采用规则的矩形网格作为数值流形方法中的数学网格,并借助适合分析的T样条实现了数值流形方法中的局部加密。适合分析的T样条定义在一个限制的T网格上,其基函数具有线性无关、单位分解、局部加密等许多重要性质,使得其非常适合用于工程设计及分析。当对一个适合分析的T网格加密后,所产生的新的网格往往不再是适合分析的T网格。基于此,提出了一种简单的数学网格加密算法,该算法能保证局部加密后的数学网格仍然是适合分析的。算例结果表明:在应力集中区和裂纹尖端等应力梯度较大区域,该算法均具有较强的适用性。 展开更多

关键词 数值流形方法 规则网格 局部加密算法 适合分析T样条

在线阅读 下载PDF 职称材料

局部网格加密方法计算不可压N-S方程

20

作者 沈建雄 《空气动力学学报》 CSCD 北大核心 1991年第1期67-74,共8页

本文发展了二维非定常不可压N-S方程的局部网格加密方法(LMR):需加密区域预先指定;仅对空间方向加密,分粗细二层网格,细网格覆盖在粗网格上;粗和细网格上分别用显式和隐式差分格式。结合压力修正法类的余量型差分格式,恰当地处理了粗细... 本文发展了二维非定常不可压N-S方程的局部网格加密方法(LMR):需加密区域预先指定;仅对空间方向加密,分粗细二层网格,细网格覆盖在粗网格上;粗和细网格上分别用显式和隐式差分格式。结合压力修正法类的余量型差分格式,恰当地处理了粗细网格之间的信息传递,使得粗细网格交接面上质量守恒,即满足连续性条件。粗网格通过插值给细网格边界值影响细网格;细网格通过粗细网格压力Poisson方程耦合求解影响粗网格,并且由于压力Poisson方程从动量方程形成,在细网格覆盖下的粗网格上,压力Poisson方程是细网格方程的组合,从而耦合求解时粗细网格压力Poisson方程不需迭代。 本文中计算了二维方腔内的自然对流模型问题,Rayleigh数是10~6。粗网格用显式ULWC格式,细网格用隐式余量型近似因式分解格式。Poisson方程用快速直接算法PO1STG(在FISHPACK中)。粗网格距△=1/16,细网格距△=1/64,加密1/4区域,计算到定常局部网格加密方法所需CPU时间比同等均匀网格(△=1/64)节约一倍以上。 展开更多

关键词 局部网格加密 不可压缩流 N-S方程

在线阅读 下载PDF 职称材料