期刊文献+
共找到572篇文章
< 1 2 29 >
每页显示 20 50 100
The θ-Methods in Numerical Solution of Systems of Differential Equations with Two Delay Terms 被引量:2
1
作者 Tian Hongjiong & Kuang Jiaoxun (Department of Mathematics, Shanghai Normal University, Shanghai 200234, China) 《Journal of Systems Engineering and Electronics》 SCIE EI CSCD 1994年第3期32-40,共9页
This paper deals with the numerical solution of initial value problems for systems of differential equations with two delay terms. We investigate the stability of adaptations of the θ-methods in the numerical solutio... This paper deals with the numerical solution of initial value problems for systems of differential equations with two delay terms. We investigate the stability of adaptations of the θ-methods in the numerical solution of test equations u'(t) = a 11 u(t) + a12v(t) + b11 u(t - τ1) + b12v(t-τ2,v'(t) = a21 u(t) + a22 v(t) + b21 u(t -τ1,) + b22 v(t -τ2), t>0,with initial conditionsu(t)=u0(t),v(t) =v0(t), t≤0.where aij, bij∈C, τj >0, i,j = 1,2,, and u0(t), v0(t)are continuous and complex valued. Sufficient conditions for the asymptotic stability of test equation are derived. Furthermore, with respect to an appropriate definition of stability for the numerical method, it is proved that the linear θ-method is stable if and only if 1/2≤θ≤1 and the one-leg θ-method is stable if and only if θ= 1. 展开更多
关键词 Delay differential equations numerical solution Θ-methodS Asymptotic stability Schur polynomial.
在线阅读 下载PDF
The Numerical Simulation of Flow Field of Jet System
2
作者 Guo Jinji, Zhang Sheng and Chen TongMechanics Dept. of Zhongshan University,Guangzhou 510275. China 《Journal of Systems Engineering and Electronics》 SCIE EI CSCD 1993年第3期62-69,共8页
In this paper, the internal fluid motion of a jet system is described by the Navier Stokes mechanics equations. For the simulation of the motion, the penalty function finite element method is used, and the velocity ve... In this paper, the internal fluid motion of a jet system is described by the Navier Stokes mechanics equations. For the simulation of the motion, the penalty function finite element method is used, and the velocity vectors and stream function curves are obtained. Using the Prandtl theory, this paper derives the free jet velocity and the jet bunch width in a half-space, the latter of which is amended by experiment. The results obtained in this paper are applied to micro-type high pressure water jet cleaner and the ejector of rocket engine. 展开更多
关键词 Fluid mechanics equations Penalty finite element method numerical simulation Free jet Jet bunch width.
在线阅读 下载PDF
独立覆盖流形法的通用计算公式和通用程序设计——(三)算例验证 被引量:1
3
作者 苏海东 宋文硕 +2 位作者 龚亚琦 韩陆超 韦玉霞 《长江科学院院报》 北大核心 2025年第4期211-218,共8页
在前两篇提出的通用计算公式、通用程序设计方法的基础上,通过求解弹性力学运动微分方程、传导方程、波动方程(含稳态和瞬态分析),给出了一至三维包括位移场、温度场、渗流场、声场、静电场和势流场等算例进行全面验证,涵盖了任意形状... 在前两篇提出的通用计算公式、通用程序设计方法的基础上,通过求解弹性力学运动微分方程、传导方程、波动方程(含稳态和瞬态分析),给出了一至三维包括位移场、温度场、渗流场、声场、静电场和势流场等算例进行全面验证,涵盖了任意形状和任意连接的网格、精确几何边界的模拟及本质边界条件的准确施加、高阶级数逼近、裂纹尖端附近解析级数的应用等独立覆盖流形法的特色功能。最后对全文进行总结,并提出“级数流形元(级数元)”的新名称。 展开更多
关键词 数值流形方法 独立覆盖 偏微分方程 级数解 网格剖分 精确几何 级数流形元
在线阅读 下载PDF
独立覆盖流形法的通用计算公式和通用程序设计——(一)通用计算公式 被引量:1
4
作者 苏海东 《长江科学院院报》 北大核心 2025年第4期193-201,210,共10页
独立覆盖流形法是偏微分方程数值计算的新方法,在近似函数构造的基础层面上形成了偏微分方程的“分区级数解”,实现了有限元法及一些新方法的主要分析功能,在网格划分的灵活性、计算稳定性等方面还具有一些独特优势,但也意味着其计算公... 独立覆盖流形法是偏微分方程数值计算的新方法,在近似函数构造的基础层面上形成了偏微分方程的“分区级数解”,实现了有限元法及一些新方法的主要分析功能,在网格划分的灵活性、计算稳定性等方面还具有一些独特优势,但也意味着其计算公式和程序设计不同于现有方法。总结了近年来该方法在固体计算中的主要研究成果,归纳了一套简洁的通用计算公式,局部近似函数中的形函数表达为单位分解函数、坐标转换矩阵和级数矩阵的乘积,具体讨论了各种情况下的形函数及其求导方法,给出各种矩阵的表达式以及时间积分方法,用于求解弹性力学运动微分方程、传导方程、波动方程,包括一维至三维的稳态和瞬态分析和三类边界条件,涵盖了高阶级数、任意形状网格、精确几何边界模拟及本质边界条件的准确施加、裂纹尖端附近采用解析级数等特色功能。利用这些公式即可开展独立覆盖流形法通用程序的开发工作。 展开更多
关键词 偏微分方程 级数解 网格剖分 精确几何 独立覆盖 数值流形方法
在线阅读 下载PDF
独立覆盖流形法的通用计算公式和通用程序设计——(二)通用程序设计
5
作者 苏海东 杨震 +2 位作者 颉志强 祁勇峰 龚亚琦 《长江科学院院报》 北大核心 2025年第4期202-210,共9页
在独立覆盖流形法通用计算公式的基础上,给出了计算程序的整体流程。对一维至三维各种几何形体(包括分区、条带和边界面)的积分方式进行总结,基于点、线、面、体的单纯形几何元素开发积分程序,实现网格形状的通用性。提出将积分模块与... 在独立覆盖流形法通用计算公式的基础上,给出了计算程序的整体流程。对一维至三维各种几何形体(包括分区、条带和边界面)的积分方式进行总结,基于点、线、面、体的单纯形几何元素开发积分程序,实现网格形状的通用性。提出将积分模块与被积函数模块分开考虑的编程思路,然后再将两者任意组合,使程序具备了扩展性,有望实现偏微分方程求解的通用性。通过级数公式和相应的各种坐标以及坐标转换矩阵、级数矩阵的确定,实现了级数的通用性。所有计算参数都可以通过用户子程序输入公式,实现输入参数的通用性。最终可用较少的程序代码,实现弹性力学运动微分方程、传导方程、波动方程的一维至三维稳态和瞬态分析(含一类至三类边界条件)。 展开更多
关键词 偏微分方程 级数解 网格剖分 精确几何 独立覆盖 数值流形方法
在线阅读 下载PDF
页岩气产能评价研究进展:内涵、方法和方向
6
作者 朱苏阳 彭真 +3 位作者 邸云婷 彭小龙 刘东晨 官文洁 《油气藏评价与开发》 北大核心 2025年第3期488-499,共12页
页岩气的生产呈现早期高产而后快速递减的动态,快速递减表明初期配产可能过高,不合理配产导致气藏产能衰竭过快,从而影响气井的经济可采储量(EUR),因此,合理评价页岩气井产能对保障气藏稳产具有重要意义。为明确目前页岩气产能评价方法... 页岩气的生产呈现早期高产而后快速递减的动态,快速递减表明初期配产可能过高,不合理配产导致气藏产能衰竭过快,从而影响气井的经济可采储量(EUR),因此,合理评价页岩气井产能对保障气藏稳产具有重要意义。为明确目前页岩气产能评价方法存在的问题和相应的解决思路,研究分析了页岩气产能的特殊内涵,综述了产能流动方程解析、流动方程数值模拟、人工智能3种方法在页岩气产能评价研究中的进展。研究发现页岩气产能明显呈现出分阶段特征,不同生产阶段产能主控因素、流动机制、流动状态都不同,早期产能和后期产能的主控因素差异较大,不同评价方法存在差异化认识。其中,流动方程解析法依赖于对流动机理的认识程度,流动方程数值模拟法结果的验证需要大量可靠数据支撑和气藏工程经验判断,人工智能方法则存在高度不透明、不可解释和泛化能力差的问题。基于此,未来发展应从深化微观和宏观2个层面结合点的页岩气渗流机理着手研究,深化地质建模—应力物性演化—裂缝扩展—多相流动数值模拟—递减分析多维度产能描述融合,并加强机器学习算法的无关机理约束、因果推断的透明度、可解释程度,以规避现有产能评价方法存在的制约性,从而构建合理的页岩气产能评价模型和方法,为准确评价气井产能、实现稳产和高效开发提供理论支撑。 展开更多
关键词 页岩气 产能评价 流动方程解析法 流动方程数值模拟法 人工智能方法 动态评价
在线阅读 下载PDF
基于格子Boltzmann法的柱群布置对污染物分布影响研究
7
作者 夏鹏飞 杜海川 +1 位作者 杜璇 曹德生 《水资源与水工程学报》 北大核心 2025年第2期118-124,133,共8页
在河流、渠道和湖泊等浅水流环境中,污染物的输运现象广泛存在。为了保护水生态环境,研究污染物在这些环境中的分布规律具有重要意义。采用基于D2Q16模型的格子Boltzmann方法(LBM)模拟浅水流动,并利用D2Q5模型对污染物扩散过程进行模拟... 在河流、渠道和湖泊等浅水流环境中,污染物的输运现象广泛存在。为了保护水生态环境,研究污染物在这些环境中的分布规律具有重要意义。采用基于D2Q16模型的格子Boltzmann方法(LBM)模拟浅水流动,并利用D2Q5模型对污染物扩散过程进行模拟,探讨了圆柱群不同间距比、不同布置形式以及水流流速对柱后污染物分布的影响。结果表明:圆柱间距比对于柱后污染物浓度分布的影响较大,当进口流速一定时,随着圆柱间距比的增大,柱后各横断面的污染物分布更加集中,当圆柱间距比L/D为1.0时,各横断面的污染物浓度极值均大于其他圆柱间距比下的污染物浓度极值;当圆柱间距比一定时,随着进口流速的增大,柱后各断面污染物浓度的极值总体上呈现不断下降的趋势;柱群布置形式对于柱后污染物的分布影响较大,正方形布置形式相比于梅花形、梯形和菱形柱后同位置横断面的污染物浓度更小。 展开更多
关键词 污染物分布 污染物浓度 柱群 格子Boltzmann法 浅水方程 数值模拟
在线阅读 下载PDF
非Lipschitz条件下Lévy噪声扰动的随机比例型微分方程的数值解
8
作者 梁飞 张丽洁 《河南科技大学学报(自然科学版)》 北大核心 2025年第2期95-104,M0008,共11页
针对满足非Lipschitz条件的Lévy噪声扰动的随机比例型微分方程,首先证明了方程的精确解在非Lipschitz条件下以大概率存在于紧集中;其次运用Euler方法构造出方程的数值解,并证明了数值解在均方意义下依概率收敛于精确解;最后通过一... 针对满足非Lipschitz条件的Lévy噪声扰动的随机比例型微分方程,首先证明了方程的精确解在非Lipschitz条件下以大概率存在于紧集中;其次运用Euler方法构造出方程的数值解,并证明了数值解在均方意义下依概率收敛于精确解;最后通过一个例子验证了结论的有效性。 展开更多
关键词 随机比例型微分方程 Lévy噪声 非LIPSCHITZ条件 EULER方法 数值解
在线阅读 下载PDF
非线性时间分布阶双曲波动方程的广义BDF2⁃θ有限元方法
9
作者 侯雅馨 刘洋 李宏 《应用数学和力学》 北大核心 2025年第1期114-128,共15页
构造了一种基于带有位移参数θ的广义向后差分公式(广义BDF2-θ)的有限元(FE)方法,用于求解非线性时间分布阶双曲波动方程.时间方向由广义BDF2-θ近似进一步得到FE全离散格式.将具有高阶时间导数的模型转化为包括两个低阶方程的耦合系统... 构造了一种基于带有位移参数θ的广义向后差分公式(广义BDF2-θ)的有限元(FE)方法,用于求解非线性时间分布阶双曲波动方程.时间方向由广义BDF2-θ近似进一步得到FE全离散格式.将具有高阶时间导数的模型转化为包括两个低阶方程的耦合系统.证明了格式的稳定性以及两个函数u和p的最优误差估计结果.最后,通过数值算例验证了格式的可行性和有效性. 展开更多
关键词 非线性时间分布阶双曲波动方程 有限元方法 广义BDF2⁃θ 稳定性 误差估计 数值模拟
在线阅读 下载PDF
Some Equivalence Transform Methods for Simulating Models in Control Systems
10
作者 Degui, Liu 《Journal of Systems Engineering and Electronics》 SCIE EI CSCD 1993年第4期32-42,共11页
In this paper two classes of equivalence transform methods for solving ordinary differential equations are proposed. One class of method is the equivalence integral transform method for special differential algebraic ... In this paper two classes of equivalence transform methods for solving ordinary differential equations are proposed. One class of method is the equivalence integral transform method for special differential algebraic problems. The advantage of this class of method is such that the amount of work calculating one integration with parameters becomes that of two interpolations, when the system of nonlinear equations is solved on the right hand side function. The other class of method is the equivalence substitution method for avoiding calculating derivative on the right hand side function. In order to avoid calculation derivatives, two equivalence substitution methods are proposed here. The application instances of some special effect of the equivalence substitution methods are given. 展开更多
关键词 ALGEBRA Calculations Computer simulation Control systems Differential equations Integration Interpolation Mathematical models Mathematical transformations Nonlinear equations numerical methods
在线阅读 下载PDF
综合能源系统建模与仿真:综述、思考与展望 被引量:3
11
作者 张苏涵 顾伟 +3 位作者 俞睿智 陆帅 姚帅 庄文楠 《电力系统自动化》 EI CSCD 北大核心 2024年第17期1-21,共21页
建模和仿真将真实系统及其物理性质抽象化为适当的数学形式,并采用数学方法求解以模拟真实系统行为或性能。为综合能源系统在各时间尺度的物理过程建立恰当的数学模型并进行仿真验证,是分析系统运行规律、优化运行策略以及评估运行性能... 建模和仿真将真实系统及其物理性质抽象化为适当的数学形式,并采用数学方法求解以模拟真实系统行为或性能。为综合能源系统在各时间尺度的物理过程建立恰当的数学模型并进行仿真验证,是分析系统运行规律、优化运行策略以及评估运行性能的基础。文中围绕电气热综合能源系统建模和仿真问题,首先介绍了电、气、热子系统和耦合设备的机理模型,总结了建模的研究进展与技术难点。进一步,从仿真框架、算法以及改进策略3个方面介绍了综合能源系统仿真技术所取得的研究进展。特别地,总结了仿真算法中数值法、半解析法、解析法的原理、应用现状及其难点。最后,围绕现存技术难点,从建模和仿真算法2个方面对未来研究进行了展望。 展开更多
关键词 综合能源系统 建模 能流计算 动态仿真 微分方程 数值方法
在线阅读 下载PDF
综合能源系统分析——从数值到解析(一):气网解析法 被引量:1
12
作者 张苏涵 顾伟 +2 位作者 陆帅 俞睿智 庄文楠 《中国电机工程学报》 EI CSCD 北大核心 2024年第16期6432-6443,I0012,共13页
为解决传统数值方法冗余变量多、近似误差大、稳定性难以保证的困难,该文提出面向综合能源系统分析的解析法,以满足运行分析的高精度、低复杂度、可行性需求。作为系列论文的第一部分,该文以天然气网络为研究对象,推导其动态模型的通用... 为解决传统数值方法冗余变量多、近似误差大、稳定性难以保证的困难,该文提出面向综合能源系统分析的解析法,以满足运行分析的高精度、低复杂度、可行性需求。作为系列论文的第一部分,该文以天然气网络为研究对象,推导其动态模型的通用解析式和实用表征形式,实现了偏微分方程向代数方程的无误差变换。并根据天然气网络解析解的性质,建立天然气管道和网络的等值模型,进一步削减模型复杂度。算例仿真分析表明,天然气网络分析的解析法具有显著的精度和计算效率优势。在相近的计算开销下,解析法相较于数值法能实现数十倍的精度提升。 展开更多
关键词 综合能源系统 天然气系统 数值法 解析法 偏微分方程
在线阅读 下载PDF
基于CFD-PBE耦合的连续结晶过程数值模拟
13
作者 潘立杰 刘伯潭 +1 位作者 赵文立 盖晓龙 《化学工程》 CSCD 北大核心 2024年第1期41-47,共7页
通过CFD(计算流体力学)与PBE(粒数衡算方程)耦合对连续结晶过程进行数值模拟。采用QMOM(积分矩方法)求解PBE。在考虑成核、生长、聚并和破碎的情况下,通过网格独立性验证确定计算网格数量。在此基础上研究挡板设置、停留时间、搅拌转速... 通过CFD(计算流体力学)与PBE(粒数衡算方程)耦合对连续结晶过程进行数值模拟。采用QMOM(积分矩方法)求解PBE。在考虑成核、生长、聚并和破碎的情况下,通过网格独立性验证确定计算网格数量。在此基础上研究挡板设置、停留时间、搅拌转速以及搅拌桨类型等因素对CSD(晶体粒径分布)演变的影响,结果表明:挡板的设置增强了结晶器内湍流程度;高转速和低停留时间会获得平均粒径较小的颗粒。对比3种搅拌桨对CSD演变的影响,轴流式的螺旋桨功耗小,适用于大颗粒晶体结晶过程,而径流式的圆盘涡轮桨适用于微粒结晶过程。 展开更多
关键词 连续结晶 数值模拟 粒数衡算方程 积分矩方法 晶体粒径分布
在线阅读 下载PDF
求解广义复Sylvester张量方程的子空间迭代法
14
作者 马昌凤 李清雅 《应用数学》 北大核心 2024年第3期601-609,共9页
本文研究求解广义复Sylvester张量方程的修正共轭梯度法(MCG)和共轭梯度最小二乘法(CGLS).它们都基于求解线性方程组的经典Krylov子空间方法.在不计舍入误差的情况下,理论分析表明对于任意初始张量,所提出的算法收敛到Sylvester张量方... 本文研究求解广义复Sylvester张量方程的修正共轭梯度法(MCG)和共轭梯度最小二乘法(CGLS).它们都基于求解线性方程组的经典Krylov子空间方法.在不计舍入误差的情况下,理论分析表明对于任意初始张量,所提出的算法收敛到Sylvester张量方程的解.数值结果验证了MCG算法和CGLS算法的有效性. 展开更多
关键词 Sylvester张量方程 MCG方法 CGLS方法 收敛性分析 数值实验
在线阅读 下载PDF
一类含时Poisson-Nernst-Planck方程的时间自适应计算 被引量:1
15
作者 韦鹏 沈瑞刚 《吉林大学学报(理学版)》 CAS 北大核心 2024年第6期1352-1358,共7页
基于相邻时间步的时间自适应有限元方法求解一类含时Poisson-Nernst-Planc k(PNP)方程,以加速求解含有多种离子的二维PNP方程的长时间数值模拟效率.数值实验结果表明,该方法在两种数值格式下都可以有效加速计算,对于长时间PNP方程的数... 基于相邻时间步的时间自适应有限元方法求解一类含时Poisson-Nernst-Planc k(PNP)方程,以加速求解含有多种离子的二维PNP方程的长时间数值模拟效率.数值实验结果表明,该方法在两种数值格式下都可以有效加速计算,对于长时间PNP方程的数值计算有效. 展开更多
关键词 Poisson-Nernst-Planck方程 有限元方法 时间自适应 数值计算
在线阅读 下载PDF
一种数值阻尼耗散可控的结构动力方程积分方法
16
作者 刘伟 童小龙 金蓉 《应用数学和力学》 CSCD 北大核心 2024年第7期922-935,共14页
数值耗散是数值积分方法的重要特性,直接影响数值仿真结果的准确性.对于含有虚假高频振动的动力系统,数值耗散能够改善数值仿真结果,但是对于具有真实高频振动的动力系统,数值耗散则会导致计算结果失真.该研究针对结构动力系统的求解,... 数值耗散是数值积分方法的重要特性,直接影响数值仿真结果的准确性.对于含有虚假高频振动的动力系统,数值耗散能够改善数值仿真结果,但是对于具有真实高频振动的动力系统,数值耗散则会导致计算结果失真.该研究针对结构动力系统的求解,提出了一种数值耗散可控的两子步隐式数值积分方法.基于理论推导,详细介绍了新积分方法的谱半径、稳定性、振幅衰减和周期延长等数值特性.新隐式积分方法通过算法参数α能够对高频虚假振动数值耗散完全可控,相应的耗散比例为1-α,其中-1≤α≤1.通过单自由度动力系统、高频虚假振动系统和多自由度非线性弹簧质量系统三个典型算例,分别证明了新隐式积分方法在计算精确性、高频数值耗散和非线性求解能力方面的优越性. 展开更多
关键词 结构动力方程 隐式积分方法 数值阻尼 放大矩阵 谱半径
在线阅读 下载PDF
应用保辛Stereo-modeling方法的振幅补偿逆时偏移成像
17
作者 李静爽 张向佳 +1 位作者 贺茜君 黄雪源 《石油地球物理勘探》 EI CSCD 北大核心 2024年第5期1037-1047,共11页
地震波数值模拟的精度很大程度上决定了逆时偏移结果的分辨率,文中采用保辛近似解析离散化方法,在Birkhoffian系统下求解黏滞声波方程,并进行振幅补偿的逆时偏移(RTM)。首先,将黏滞声波方程构造成一个Birkhoffian系统,对于该系统采用二... 地震波数值模拟的精度很大程度上决定了逆时偏移结果的分辨率,文中采用保辛近似解析离散化方法,在Birkhoffian系统下求解黏滞声波方程,并进行振幅补偿的逆时偏移(RTM)。首先,将黏滞声波方程构造成一个Birkhoffian系统,对于该系统采用二阶保辛Runge-Kutta法进行时间推进,分别采用Stereo-modeling法(STEM)和传统有限差分法进行空间算子离散,相应得到SSM(Symplectic Stereo-modeling Method)和CSM(Conventional Symplectic Method)两种方法。之后,对SSM和CSM进行一系列数值性质分析,包括数值频散分析、与解析解比较、效率对比、稳定性测试等。结果显示,SSM的最大数值频散误差约为9%,而CSM约为26%;SSM的计算效率高且长时计算稳定,计算效率约为CSM的两倍。最后,应用三个模型进行了正演数值模拟和逆时偏移试验。对于衰减数据,与声波RTM相比,黏滞声波RTM成像振幅更高,能获得与无衰减数据声波RTM相近的成像振幅。对于油气藏模型,SSM比CSM的成像精度更高、数值频散更小。 展开更多
关键词 黏滞声波方程 振幅补偿 保辛方法 逆时偏移 数值频散
在线阅读 下载PDF
求解Stein张量方程的张量格式BCGSTAB算法
18
作者 马昌凤 谢亚君 卜凡 《数学物理学报(A辑)》 CSCD 北大核心 2024年第6期1652-1664,共13页
双共轭梯度稳定化方法(BCGSTAB)是双共轭梯度方法的快速和光滑收敛的变形.该文将BCGSTAB算法推广到求解Stein张量方程,给出了张量格式的算法以及解的存在性的具体证明过程,并得到了该算法的收敛性结果.数值实验证实了该算法用于求解Stei... 双共轭梯度稳定化方法(BCGSTAB)是双共轭梯度方法的快速和光滑收敛的变形.该文将BCGSTAB算法推广到求解Stein张量方程,给出了张量格式的算法以及解的存在性的具体证明过程,并得到了该算法的收敛性结果.数值实验证实了该算法用于求解Stein张量方程是有效且可行的. 展开更多
关键词 STEIN 张量方程 BCGSTAB 算法 收敛性分析 数值实验
在线阅读 下载PDF
连续Sylvester矩阵方程的参数化单步HSS迭代法
19
作者 马昌凤 《工程数学学报》 CSCD 北大核心 2024年第3期481-493,共13页
对连续Sylvester矩阵方程的数值算法进行了深入研究,并创新性地提出了一种参数化单步HSS迭代方法。该方法具有独特的求解思路,并证明了其收敛性。为提升性能,通过最小化迭代矩阵谱半径上界寻找拟最优参数。数值实验验证了新方法的有效... 对连续Sylvester矩阵方程的数值算法进行了深入研究,并创新性地提出了一种参数化单步HSS迭代方法。该方法具有独特的求解思路,并证明了其收敛性。为提升性能,通过最小化迭代矩阵谱半径上界寻找拟最优参数。数值实验验证了新方法的有效性和稳健性,展示了其在求解连续Sylvester矩阵方程时的高效和稳定,为相关数值计算提供新工具。 展开更多
关键词 连续Sylvester矩阵方程 PSHSS迭代法 收敛性分析 数值实验
在线阅读 下载PDF
基于顶点中心有限元算法的重力场矢量和重力梯度张量高精度模拟
20
作者 童孝忠 孙娅 +1 位作者 黄基文 柳建新 《Journal of Central South University》 SCIE EI CAS CSCD 2024年第5期1659-1670,共12页
密度非均质性引起的重力异常由三维泊松方程控制,而目前大多数正演模拟方法都依赖于其积分解和以单元为中心的数值方法。当利用重力位计算重力场时,这些数值策略将不可避免地失去准确性。为了缓解这一问题,本文提出了一种高效、准确的... 密度非均质性引起的重力异常由三维泊松方程控制,而目前大多数正演模拟方法都依赖于其积分解和以单元为中心的数值方法。当利用重力位计算重力场时,这些数值策略将不可避免地失去准确性。为了缓解这一问题,本文提出了一种高效、准确的高阶顶点中心有限元方法来模拟三维重力异常。首先,通过具有六面体网格的顶点中心有限元来建立正演算法,并选用ILU-BICGSTAB迭代方法求解大型对称稀疏线性方程组。其次,为了获得重力位的一阶导数和二阶导数,采用了高阶拉格朗日插值技术。最后,采用三维立方体密度模型测试了顶点中心有限元算法的准确性,并利用薄垂直矩形棱镜模型和实测模型测试了算法的灵活性。数值结果表明,高阶顶点中心有限元算法能获得高精度的重力场矢量和重力梯度张量。与精确积分解和顶点中心算法相比,高阶顶点中心有限元格式在模拟三维重力异常方面表现出更高的效率和准确性。同时,相较于单元中心数值解,高阶顶点中心有限元算法在模拟三维重力异常表现出更高的效率和准确性。 展开更多
关键词 重力异常 三维泊松方程 顶点中心有限元算法 数值模拟 ILU-BICGSTAB迭代法
在线阅读 下载PDF
上一页 1 2 29 下一页 到第
使用帮助 返回顶部