期刊文献+
共找到143篇文章
< 1 2 8 >
每页显示 20 50 100
耦合非线性薛定谔方程组孤立子解的局部间断Petrov-Galerkin方法数值模拟
1
作者 赵国忠 蔚喜军 《工程数学学报》 CSCD 北大核心 2024年第6期1109-1132,共24页
耦合非线性薛定谔方程组在量子物理、非线性光学、晶体物理、波色–爱因斯坦凝聚和水波动力学等很多物理领域有着重要的应用价值。提出了一种局部间断PetrovGalerkin方法。首先,将耦合非线性薛定谔方程组改写为一阶微分方程组。空间离... 耦合非线性薛定谔方程组在量子物理、非线性光学、晶体物理、波色–爱因斯坦凝聚和水波动力学等很多物理领域有着重要的应用价值。提出了一种局部间断PetrovGalerkin方法。首先,将耦合非线性薛定谔方程组改写为一阶微分方程组。空间离散采用间断Petrov-Galerkin方法,时间离散采用三阶总变差不增Runge-Kutta方法。数值实验表明,该算法对线性元和二次元都能达到最优收敛阶。通过数值算例计算了质量、动量和能量守恒量,该算法可以很好地模拟单孤立子传输、双孤立子碰撞和三孤立子碰撞现象。此外,该算法可以在较长的时间间隔内模拟复杂波型的相互作用或传播,还可以模拟孤子传输和孤子产生现象。 展开更多
关键词 局部间断Petrov-galerkin方法 耦合非线性薛定谔方程 孤立子碰撞 守恒量
在线阅读 下载PDF
间断Galerkin有限元隐式算法GPU并行化研究
2
作者 高缓钦 陈红全 +1 位作者 贾雪松 徐圣冠 《空气动力学学报》 CSCD 北大核心 2024年第2期21-33,I0001,共14页
为了提高间断伽辽金(discontinuous Galerkin,DG)有限元方法的计算效率,围绕求解Euler方程,构建了基于图形处理器(graphics processing unit,GPU)并行加速的隐式DG算法。算法结合Roe格式进行空间离散,采用人工黏性法处理激波等间断问题... 为了提高间断伽辽金(discontinuous Galerkin,DG)有限元方法的计算效率,围绕求解Euler方程,构建了基于图形处理器(graphics processing unit,GPU)并行加速的隐式DG算法。算法结合Roe格式进行空间离散,采用人工黏性法处理激波等间断问题,时间推进选用下上对称高斯-赛德尔(lower-upper symmetric Gauss-Seidel,LU-SGS)隐式格式。为了克服传统隐式格式固有的数据关联依赖问题,借助于本文提出的面向任意网格的单元着色分组技术,先给出了LUSGS隐式格式的并行化改造,使得隐式时间推进能按颜色组别依次并行,由于同一颜色组内算法已不存在数据关联,可以据此实现并行化。在此基础上,再结合DG算法局部紧致等特点,基于统一计算设备架构(compute unified device architecture,CUDA)编程模型,设计了依据单元的核函数,并构建了对应的线程与数据结构,给出了DG有限元隐式GPU并行算法。最后,发展的算法通过了多个二维和三维典型流动算例考核与性能测试,展示出隐式算法GPU加速的效果,且获得的计算结果能与现有的文献或实验数据接近。 展开更多
关键词 间断伽辽金方法 LU-SGS隐式格式 GPU并行化 单元着色分组 EULER方程
在线阅读 下载PDF
Caputo型时间分数阶变系数扩散方程的局部间断Galerkin方法
3
作者 代巧巧 李东霞 《上海大学学报(自然科学版)》 CAS CSCD 北大核心 2024年第1期174-190,共17页
提出一种带有Caputo导数的时间分数阶变系数扩散方程的数值解法.方程的解在初始时刻附近通常具有弱正则性,采用非一致网格上的L1公式离散时间分数阶导数,并使用局部间断Galerkin(local discontinuous Galerkin,LDG)方法离散空间导数,给... 提出一种带有Caputo导数的时间分数阶变系数扩散方程的数值解法.方程的解在初始时刻附近通常具有弱正则性,采用非一致网格上的L1公式离散时间分数阶导数,并使用局部间断Galerkin(local discontinuous Galerkin,LDG)方法离散空间导数,给出方程的全离散格式.基于离散的分数阶Gronwall不等式,证明了格式的数值稳定性和收敛性,且所得结果关于α是鲁棒的,即当α→1^(-)时不会发生爆破.最后,通过数值算例验证理论分析的结果. 展开更多
关键词 局部间断galerkin方法 非一致时间网格 α-鲁棒 弱正则性 变系数
在线阅读 下载PDF
任意单元间断Galerkin有限元计算方法研究 被引量:16
4
作者 贺立新 张来平 张涵信 《空气动力学学报》 EI CSCD 北大核心 2007年第2期157-162,共6页
基于龙格库塔间断Galerkin(RKDG)有限元法的构造思想,通过局部坐标变换,发展了非正交单元DG有限元计算方法;借鉴非结构网格有限体积隐式计算方法,发展了适应于DG有限元方法的隐式计算方法;借鉴一维和二维(三角形单元)DG有限元限制器构... 基于龙格库塔间断Galerkin(RKDG)有限元法的构造思想,通过局部坐标变换,发展了非正交单元DG有限元计算方法;借鉴非结构网格有限体积隐式计算方法,发展了适应于DG有限元方法的隐式计算方法;借鉴一维和二维(三角形单元)DG有限元限制器构造方法,提出了非正交三棱柱单元限制器方法;利用上述方法数值模拟了球头和双椭球的高超声速粘性绕流,得到了清晰的流场结构,并得到了较好的压力和热流分布,表明该方法在复杂高超声速流动的数值模拟方面具有广阔的应用前景。 展开更多
关键词 非正交单元 间断Glerkin有限元 隐式计算方法
在线阅读 下载PDF
起伏地表弹性波传播的间断Galerkin有限元数值模拟方法 被引量:17
5
作者 薛昭 董良国 +1 位作者 李晓波 刘玉柱 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2014年第4期1209-1223,共15页
间断Galerkin有限元法(DG-FEM)作为一种有效的高阶有限元法受到了国内外学者的广泛关注.本文基于任意高阶间断Galerkin有限元法对弹性波方程进行空间离散,并将离散后所得的非齐次线性常微分方程系统齐次化,最后结合针对齐次问题的强稳... 间断Galerkin有限元法(DG-FEM)作为一种有效的高阶有限元法受到了国内外学者的广泛关注.本文基于任意高阶间断Galerkin有限元法对弹性波方程进行空间离散,并将离散后所得的非齐次线性常微分方程系统齐次化,最后结合针对齐次问题的强稳定性保持龙格库塔(SSP Runge-Kutta)算法,将DG-FEM推广至时间任意高阶精度.另外,借鉴近最佳匹配层(NPML)的思想,基于复频移(CFS)拉伸坐标变换推导了一种新的PML吸收边界条件(简称为CFS-NPML),该CFS-NPML能够与DG-FEM算法很好地结合,形成有效的起伏地表地震波传播数值模拟技术.数值试验结果表明,DG-FEM具有高阶精度,可以适应任意复杂起伏地表和复杂构造情况下的弹性波传播数值模拟.同时,CFS-NPML对包括面波等震相的人为边界反射都具有良好的吸收效果. 展开更多
关键词 间断galerkin有限元法 起伏地表 弹性波传播 任意高阶Runge-Kutta时间格式 CFS-NPML
在线阅读 下载PDF
基于雅可比矩阵精确计算的GMRES隐式方法在间断Galerkin有限元中的应用 被引量:5
6
作者 龚小权 贾洪印 +2 位作者 陈江涛 赵辉 周桂宇 《空气动力学学报》 CSCD 北大核心 2019年第1期121-132,共12页
为改善高阶间断Galerkin有限元方法(DG)时间推进效率,在三维非结构网格下针对该方法建立了并行广义最小残差(Generalized Minimal Residual,GMRES)隐式时间迭代方法,GMRES方法基于科学计算工具包PETSc中的Krylov子空间求解器实现。为进... 为改善高阶间断Galerkin有限元方法(DG)时间推进效率,在三维非结构网格下针对该方法建立了并行广义最小残差(Generalized Minimal Residual,GMRES)隐式时间迭代方法,GMRES方法基于科学计算工具包PETSc中的Krylov子空间求解器实现。为进一步提高GMRES的计算效率,发展了方程组右端项残值雅可比精确计算方法,针对无黏通量Roe格式和黏性通量BR2(Bassi Rebay 2)黏性计算方法,分别解析给出其对守恒变量多项式自由度的雅可比矩阵。基于建立的方法首先采用NACA0012翼型研究了GMRES的重启次数及收敛参数对方法收敛性影响,然后采用无黏及黏性算例对比研究了基于雅可比矩阵不同计算方法的GMRES计算效率,同时对比研究了雅可比矩阵完全近似求解下GMRES和LU-SGS(Lower Upper-Symmetric Gauss-Seidel)的计算效率。结果表明,建立的基于右端项残值雅可比矩阵精确求解的GMRES方法能够大幅提高不同精度DG方法的CFL(CourantFriedrichs-Lewy)数,相比前面提到的其它方法具有更高的计算效率,其收敛速度实现量级以上的提高。 展开更多
关键词 隐式迭代 间断galerkin有限元方法 GMRES 精确雅可比矩阵 计算效率
在线阅读 下载PDF
一种适用于非结构网格的间断Galerkin有限元LU-SGS隐式方法 被引量:5
7
作者 马明生 龚小权 +1 位作者 邓有奇 赵辉 《西北工业大学学报》 EI CAS CSCD 北大核心 2016年第5期754-760,共7页
具有TVD性质的显式Runge-Kutta间断Galerkin(RKDG)格式在CFD领域得到广泛应用,但是显式计算稳定性差、计算效率低。为改善时间推进效率,基于高阶间断Galerkin有限元方法,采用欧拉一阶后差(BDF1),发展了一套高效的隐式LU-SGS(lower upper... 具有TVD性质的显式Runge-Kutta间断Galerkin(RKDG)格式在CFD领域得到广泛应用,但是显式计算稳定性差、计算效率低。为改善时间推进效率,基于高阶间断Galerkin有限元方法,采用欧拉一阶后差(BDF1),发展了一套高效的隐式LU-SGS(lower upper-symmetric Gauss-Seidel)求解方法,方法基于MPI并行实现,适合于不同计算精度。针对非线性系统左端项矩阵,对比了简化前后LU-SGS的计算效率。建立的间断Galerkin有限元方法基于非结构网格,采用Taylor基函数,计算精度最高达到四阶精度。通过NACA0012翼型以及M6机翼算例对发展的LU-SGS方法进行了考察,与显式算法相比,隐式格式的迭代步数和CPU时间均较大程度减小,效率能够提高1个量级以上。最后将隐式算法用于复杂外形翼身组合体F4的流场计算,结果表明所发展的隐式方法具有较好的鲁棒性,能够用于复杂外形计算。 展开更多
关键词 间断galerkin有限元 欧拉方程 Taylor基函数 LU-SGS 计算效率 非结构网格
在线阅读 下载PDF
Navier-Stokes方程间断Galerkin有限元方法研究 被引量:24
8
作者 于剑 阎超 《力学学报》 EI CSCD 北大核心 2010年第5期962-970,共9页
通过引入全局提升算子和局部提升算子,发展了求解Navier-Stokes方程的间断Galerkin(discontinuousGalerkin,DG)有限元方法的一般框架,并在此框架下给出了几种典型黏性离散格式的具体表达形式.对局部提升算子的求解给出了详细的计算步骤... 通过引入全局提升算子和局部提升算子,发展了求解Navier-Stokes方程的间断Galerkin(discontinuousGalerkin,DG)有限元方法的一般框架,并在此框架下给出了几种典型黏性离散格式的具体表达形式.对局部提升算子的求解给出了详细的计算步骤.同时还给出了一种简单有效的计算方法来对物面边界进行高阶近似.为了能够对NS方程进行精度测试,采用对原始系统添加源项的方法构造精确解.二维Euler和NS系统的精度测试表明该方法达到了DG方法的理论精度.二维圆柱无黏绕流的计算结果表明关于物面边界的高阶近似方法能够保持DG方法原有的精度.卡门涡街数值模拟则进一步验证了该方法的正确性并且显示出DG方法较高的计算精度和分辨率. 展开更多
关键词 间断galerkin有限元方法 NAVIER-STOKES方程 黏性项 计算流体力学
在线阅读 下载PDF
不连续温度场问题的间断Galerkin方法 被引量:5
9
作者 刘冬欢 郑小平 刘应华 《力学学报》 EI CSCD 北大核心 2010年第1期74-82,共9页
针对不连续温度场问题建立了一种间断Galerkin有限元方法,该方法的主要特点是允许插值函数在单元边界上存在跳变.在建立有限元方程时,通过在单元边界上引入数值通量项和稳定性项来处理间断效应,并且数值通量可以直接由接触热阻的定义式... 针对不连续温度场问题建立了一种间断Galerkin有限元方法,该方法的主要特点是允许插值函数在单元边界上存在跳变.在建立有限元方程时,通过在单元边界上引入数值通量项和稳定性项来处理间断效应,并且数值通量可以直接由接触热阻的定义式导出.数值算例表明该方法可以很方便且准确地捕捉到结构内部由于接触热阻而引起的温度跳变,同时在局部高梯度温度场的模拟方面也比常规连续Galerkin有限元方法效率明显要高.该方法也为研究由接触热阻引起的温度场与应力场之间的耦合问题提供了一种新的数值模拟手段. 展开更多
关键词 间断galerkin 接触热阻 高梯度 数值通量 耦合非线性
在线阅读 下载PDF
高阶间断Galerkin方法求解三维欧拉方程的研究 被引量:2
10
作者 郝海兵 杨永 左岁寒 《西北工业大学学报》 EI CAS CSCD 北大核心 2011年第1期128-132,共5页
在三维非结构网格上,对高阶间断Galerkin方法求解定常三维欧拉方程进行研究。文中使用Roe格式通量函数计算网格单元边界上的数值通量;时间方向采用显式Runge-Kutta方法推进;并引入激波探测器和斜率限制器技术,成功地抑制流场解在间断处... 在三维非结构网格上,对高阶间断Galerkin方法求解定常三维欧拉方程进行研究。文中使用Roe格式通量函数计算网格单元边界上的数值通量;时间方向采用显式Runge-Kutta方法推进;并引入激波探测器和斜率限制器技术,成功地抑制流场解在间断处的数值振荡。对M6机翼跨音速无粘流场进行数值模拟,结果表明:计算结果和实验值吻合较好,和同等精度的有限体积法相比,具有更低的数值耗散和更强的激波捕捉能力。 展开更多
关键词 间断galerkin有限元 激波探测器 斜率限制器 欧拉方程
在线阅读 下载PDF
基于简单WENO-间断Galerkin的Euler方程自适应计算 被引量:3
11
作者 吴泽艳 王立峰 武哲 《北京航空航天大学学报》 EI CAS CSCD 北大核心 2016年第4期806-814,共9页
为了得到Euler方程的高精度、高分辨率数值解,介绍了间断Galerkin方法、三角形单元上简单WENO限制器的基本原理以及基于自适应网格加密的激波捕捉方法。将简单WENO限制器-间断Galerkin方法应用到曲边四边形单元上,通过单元边界上高斯积... 为了得到Euler方程的高精度、高分辨率数值解,介绍了间断Galerkin方法、三角形单元上简单WENO限制器的基本原理以及基于自适应网格加密的激波捕捉方法。将简单WENO限制器-间断Galerkin方法应用到曲边四边形单元上,通过单元边界上高斯积分点的坐标来搜索相邻单元从而得到相邻单元的单元编号,实现了基于"问题单元"的局部网格加密自适应计算。对若干典型问题进行编程计算,结果表明,简单WENO限制器可以应用到曲边四边形单元上,且可适用于局部网格加密时具有"悬挂节点"的非结构网格上的激波捕捉。 展开更多
关键词 间断galerkin方法 简单WENO限制器 EULER方程 自适应计算 曲边四边形单元
在线阅读 下载PDF
一维溃坝问题的间断Galerkin方法 被引量:3
12
作者 张飞 汪继文 《安徽大学学报(自然科学版)》 CAS 北大核心 2005年第2期5-8,共4页
将间断Galerkin方法应用于一维溃坝问题中,采用了两种不同数值流函数,对其进行比较。使用TVDM,TVBM型的限制器,同时利用一种改进的限制器来消除振荡。与一种高分辨率TVD性质的差分方法做了比较。结果表明间断Galerkin方法得到了更加逼... 将间断Galerkin方法应用于一维溃坝问题中,采用了两种不同数值流函数,对其进行比较。使用TVDM,TVBM型的限制器,同时利用一种改进的限制器来消除振荡。与一种高分辨率TVD性质的差分方法做了比较。结果表明间断Galerkin方法得到了更加逼真的图像。 展开更多
关键词 间断galerkin方法 一维 溃坝 差分方法 限制器 振荡 性质 问题 TVD 消除
在线阅读 下载PDF
固体非傅立叶温度场的时域间断Galerkin有限元法 被引量:3
13
作者 武文华 李锡夔 《计算力学学报》 EI CAS CSCD 北大核心 2007年第2期219-223,共5页
运用时域间断Galerkin有限元法[1],对高频非傅立叶热波动问题[2-3]进行分析。其主要特点是:取温度及温度的时间导数为基本未知量,对其分别进行3次Hermite插值和线性插值。在保证节点温度自动保持连续的基础上,温度的时间导数在离散时域... 运用时域间断Galerkin有限元法[1],对高频非傅立叶热波动问题[2-3]进行分析。其主要特点是:取温度及温度的时间导数为基本未知量,对其分别进行3次Hermite插值和线性插值。在保证节点温度自动保持连续的基础上,温度的时间导数在离散时域存在间断。数值结果表明所提出的方法能够滤掉虚假的数值震荡,能够良好地模拟固体中的非傅立叶热波动行为。 展开更多
关键词 非傅立叶 热波动 间断galerkin有限元
在线阅读 下载PDF
对流占优的扩散问题的局部间断Galerkin方法 被引量:3
14
作者 王阿霞 马逸尘 《西安交通大学学报》 EI CAS CSCD 北大核心 2008年第2期234-237,共4页
针对具有周期性边界条件对流占优的扩散问题中的二阶导数,引入辅助变量,构造了局部间断Galerkin(LDG)方法,并给出了方法的稳定性结果和误差估计式.局部间断Galerkin方法是Runge-Kutta间断Galerkin方法的推广,具有高阶精度,能够灵活处理... 针对具有周期性边界条件对流占优的扩散问题中的二阶导数,引入辅助变量,构造了局部间断Galerkin(LDG)方法,并给出了方法的稳定性结果和误差估计式.局部间断Galerkin方法是Runge-Kutta间断Galerkin方法的推广,具有高阶精度,能够灵活处理复杂区域,易于处理复杂边界的边值问题,能够有效去除近似解在间断、大梯度处产生的虚假振荡.数值实验说明,当有限元空间取为一次多项式空间时,LDG方法具有二阶收敛,误差满足理论估计式.该方法可以推广到更高阶的方程,如Korteweg-de Vries方程、重调和方程等. 展开更多
关键词 局部间断galerkin方法 对流占优的扩散问题 高阶精度 误差估计
在线阅读 下载PDF
基于局部间断Galerkin方法的p型有限元 被引量:2
15
作者 李子然 吴长春 《中国科学技术大学学报》 CAS CSCD 北大核心 2003年第3期324-329,共6页
将基于三变量能量原理的局部间断Galerkin方法 (localdiscontinuousGalerkin ,LDG)应用于 p型单元的构造 .该方法采用间断的单元试解 ,不需要满足普通有限元所必须的协调条件 ,就能使构造高阶的插值函数变得更加灵活和容易 .在此基础上 ... 将基于三变量能量原理的局部间断Galerkin方法 (localdiscontinuousGalerkin ,LDG)应用于 p型单元的构造 .该方法采用间断的单元试解 ,不需要满足普通有限元所必须的协调条件 ,就能使构造高阶的插值函数变得更加灵活和容易 .在此基础上 ,对应力和应变场采用Legendre正交多项式进行插值 ,避免了柔度矩阵的求逆过程 .数值算例表明这种方法构造出的 p型单元不仅升阶过程简单 ,而且具有较高的精度 . 展开更多
关键词 p型有限元 局部间断galerkin方法 Legendre正交多项式
在线阅读 下载PDF
一种新型的节点间断Galerkin格式滤波器设计 被引量:2
16
作者 郭永恒 杨永 +1 位作者 刘逸鸥 张强 《西北工业大学学报》 EI CAS CSCD 北大核心 2013年第1期94-97,共4页
使用高阶数值格式求解双曲守恒律时,在间断解附近产生的Gibbs振荡很容易导致非物理解的产生,从而使计算过程中断。文章以节点间断Galerkin格式为背景,使用广义Fourier级数这一分析工具,提出了相对振荡强度的概念,并以此作为滤波器的筛... 使用高阶数值格式求解双曲守恒律时,在间断解附近产生的Gibbs振荡很容易导致非物理解的产生,从而使计算过程中断。文章以节点间断Galerkin格式为背景,使用广义Fourier级数这一分析工具,提出了相对振荡强度的概念,并以此作为滤波器的筛查对象;构造了新的指数型罚函数,它能够根据振荡的强弱,自动识别非物理解的广义Fourier分量并加以过滤。在这项算法的支持下,利用节点间断Galerkin格式顺利求解了Euler方程,数值模拟了超临界翼型RAE2822翼型的跨声速外流场和超燃冲压发动机的超声速内流场。结果表明,在这种新型的滤波器基础上,不仅很好地捕捉到了激波的正确位置,而且有效地抑制了激波附近的非物理振荡。 展开更多
关键词 节点间断galerkin 相对振荡强度 广义Fourier级数 滤波器
在线阅读 下载PDF
谐振管内非线性驻波的二维间断Galerkin计算方法 被引量:1
17
作者 郭琪磊 马尧 宁方立 《机械设计与制造》 北大核心 2017年第3期162-165,169,共5页
基于柱坐标下的二维瞬态可压缩Navier-Stokes方程、理想气体状态方程以及保强稳定的Runge-Kutta方法,提出了求解谐振管内非线性驻波的二维间断Galerkin计算方法。通过处理边界条件以及采取局部时间步长,利用所提出的二维间断Galerkin计... 基于柱坐标下的二维瞬态可压缩Navier-Stokes方程、理想气体状态方程以及保强稳定的Runge-Kutta方法,提出了求解谐振管内非线性驻波的二维间断Galerkin计算方法。通过处理边界条件以及采取局部时间步长,利用所提出的二维间断Galerkin计算方法对圆柱形谐振管内非线性驻波进行求解,得到了谐振管内压力空间分布以及压力波形、速度波形、密度波形等物理特性,证明时间域、空间域内均产生激波,从而使压力幅值和密度幅值无法进一步提升;同时研究了激励加速度幅值对谐振管内非线性驻波压力波形及非稳态进程的影响。 展开更多
关键词 谐振管 非线性驻波 二维 间断galerkin 强声密封
在线阅读 下载PDF
时域间断Galerkin有限元法在激光热加工过程中应用 被引量:1
18
作者 吴志刚 郭攀 武文华 《大连理工大学学报》 EI CAS CSCD 北大核心 2012年第1期1-5,共5页
针对半无限体和薄膜结构受多种激光热源作用下的非傅里叶热传导过程,采用时域间断Galerkin有限元法进行数值仿真.其主要特点是在时域内对温度及其时间导数分别进行三次Hermite插值和线性插值.对于半无限激光热源热传导问题,其计算结果... 针对半无限体和薄膜结构受多种激光热源作用下的非傅里叶热传导过程,采用时域间断Galerkin有限元法进行数值仿真.其主要特点是在时域内对温度及其时间导数分别进行三次Hermite插值和线性插值.对于半无限激光热源热传导问题,其计算结果与解析解吻合良好.算例表明,时域间断Galerkin有限元法在高频激光脉冲问题中,没有虚假的数值振荡,具有广泛的工程应用性. 展开更多
关键词 非傅里叶热传导 时域间断galerkin有限元法 激光热源 数值仿真 薄膜
在线阅读 下载PDF
基于模拟有限差分法的两相渗流流线模拟新方法 被引量:1
19
作者 饶翔 何旭鹏 +1 位作者 赵辉 刘怡娜 《计算力学学报》 北大核心 2025年第1期26-36,共11页
首次提出了基于模拟有限差分法的流线模拟新方法,该方法利用适用于全张量渗透率及非结构网格的模拟有限差分法获取高精度的压力及渗流速度分布,采用简易的流线追踪方法在二维三角网格或三维四面体网格情况下快速获取流线分布并计算流线... 首次提出了基于模拟有限差分法的流线模拟新方法,该方法利用适用于全张量渗透率及非结构网格的模拟有限差分法获取高精度的压力及渗流速度分布,采用简易的流线追踪方法在二维三角网格或三维四面体网格情况下快速获取流线分布并计算流线上的飞行时间分布,使用基于简化龙格库塔格式的间断伽辽金法在各一维流线上获取低数值耗散的含水饱和度分布。理论上来说,相较于传统各类流线模拟方法,该方法能够同时适用于复杂油藏几何和各向异性储层,相较于仅计算网格平均饱和度的有限体积、模拟有限差分等方法,能够显著降低饱和度计算的数值耗散,相较于计算网格顶点饱和度的间断伽辽金法等高阶方法,该方法能避免强稳定性条件对时间步的限制而提高计算效率,且能进一步降低数值耗散。数值算例也表明该方法能够在各向异性渗透率、非结构网格、复杂几何情况下具有计算精度和计算效率的双重优势,且能够提供直观的流动诊断结果。 展开更多
关键词 模拟有限差分 流线模拟 间断伽辽金 油水两相渗流
在线阅读 下载PDF
复杂地质条件的间断有限元地震波数值模拟及GPU加速
20
作者 韩德超 刘卫华 +2 位作者 张春丽 袁媛 白鹏 《石油物探》 北大核心 2025年第4期639-652,共14页
间断Galerkin有限元方法(DGFEM)是一种具有较高模拟精度的有限元方法,但其算法编程难度大,其针对各类复杂介质的波动方程的算法目前未见统一的计算格式。为此,基于三角形非结构化网格以及局部Lax-Friedrichs数值流,构建了针对复杂介质... 间断Galerkin有限元方法(DGFEM)是一种具有较高模拟精度的有限元方法,但其算法编程难度大,其针对各类复杂介质的波动方程的算法目前未见统一的计算格式。为此,基于三角形非结构化网格以及局部Lax-Friedrichs数值流,构建了针对复杂介质波动方程模拟的DGFEM编程计算矩阵,并进一步得出了适用于各类复杂介质模拟的单一波场分量的通用计算格式。该通用计算格式能够有效提升DGFEM算法编程的可拓展性。基于该格式给出了DGFEM的通用CUDA核函数的构建方法,并形成CPU+GPU的二维DGFEM并行计算程序框架。通用CUDA核函数可以将DGFEM算法进一步延伸到更加复杂的介质以及三维情况。理论模型和复杂山地模型的数值实验结果表明,构建的通用计算格式和CUDA核函数可以准确模拟声波、弹性波、粘弹性波、孔隙弹性波方程描述的纵波、横波以及慢纵波等波现象。相比单核CPU模拟,二维DGFEM弹性波GPU计算加速比平均在100倍左右。同时,弹性波、粘弹性波、孔隙弹性波模拟耗时约为声波模拟的1.7,2.3,3.0倍。此结果可以指导复杂介质耦合条件模拟时多进程的负载平衡。 展开更多
关键词 间断galerkin有限元方法 弹性波 粘弹性波 孔隙弹性波 数值模拟 GPU并行计算.
在线阅读 下载PDF
上一页 1 2 8 下一页 到第
使用帮助 返回顶部