共查询到20条相似文献,搜索用时 21 毫秒
1.
为高效求解扩展型Boussinesq水波方程,建立了基于有限差分和有限体积方法的混合数值格式。将一维控制方程写为守恒形式,方程中通量部分采用有限体积方法求解,剩余部分采用有限差分方法求解。其中,有限体积方法采用Godunov类高分辨率格式,并结合HLL(Harten-Lax and van Leer)式黎曼问题近似解求界面数值通量,黎曼问题界面左右变量通过高精度状态插值方法(MUSCL)构筑。有限差分方法则采用具有二阶精度的中心差分公式进行。采用具有TVD(Total Variation Diminishing)性质的三阶龙格-库塔多步积分法进行时间积分。对数值模式进行了验证,数值结果同解析解或实验数据吻合良好。 相似文献
2.
非饱和土热湿耦合传输问题的数值解——理论及一维问题的解 总被引:1,自引:0,他引:1
本文介绍了求解非饱和土壤中热量和水分耦合传输问题的一种数值方法——积分有限差分方法(IFDM)。基于菲利普-迪弗瑞斯(Philip-De Vries)多孔介质热湿耦合流动模型,采用积分有限差分方法,编制了求解热湿传输问题的计算程序 HM1,该程序可用于求解各类工程中遇到的多孔介质一维传热传湿问题。文中还给出了计算结果与实测结果相比较的实例。 相似文献
3.
4.
基于球坐标系下有限差分的地磁测深三维正演 总被引:2,自引:0,他引:2
为了计算全球尺度电磁感应的响应,本文介绍地磁测深频率域三维正演。正演算法采用球坐标系下的交错网格有限差分方法,从Maxwell方程的积分形式出发,采用PARDISO对离散后的方程组求解,避免了迭代求解的散度校正。为了验证本文结果的正确性和精度,与前人的有限元和有限差分方法进行了对比,一维层状模型的三维交错网格有限差分数值结果和解析解相对误差小于5%,双半球模型的计算结果与前人的计算结果完全吻合。三维"棋盘模型"计算表明磁场分量对异常体的大小和位置具有很好的分辨能力。 相似文献
5.
常用的求解地下水渗流模型有限差分法,三角形单元Δikj内水文地质参数相同,以相同参数的三角形单元进行参数分区,这种方法,对某些情形可以达到解决地下水渗流场模拟和预报问题,但对非均质各向异性有一定的局限性。本文阐述了以三角形单元的棱边控制面积作为参数分区的最小单元的参数分区方法,建立了单元非等参有限差分方程,给出了实际应用例子。该方法可更准确刻画非均质各向异性问题,同时兼容以往的差分方程,可退化成一般有限差分格式。 相似文献
6.
波动方程数值模拟是研究地震波传播机理的重要工具,有限差分求解波动方程是当前地震波数值模拟的主要方法之一。当地下介质中的地震波速度较低或地震波高频成分丰富时,常规有限差分技术常常产生严重的数值频散误差,这种误差会降低数值模拟的精度,影响对地震波传播机理的分析。为压制地震波数值模拟时产生的数值频散误差,提高波场模拟精度,提出了基于NAD算子的时间四阶精度波动方程差分格式。根据对应的差分格式,分析了该差分格式的数值频散关系。与常规四阶精度差分算法的频散曲线相比,基于NAD时间四阶精度差分方法不但能够实现时间频散的有效压制,同时其基于更多网格点的位移分量和位移梯度分量空间微分求解方法还能够实现空间频散的有效压制。另外在相同模型条件下,基于NAD算法的声波方程时间四阶差分解法可采用大网格对模拟空间进行差分离散,减少网格数,提高计算效率。 相似文献
7.
介绍了有限差分方法的原理,叙述了砂岩型铀矿体上氡气场的分布,并建立了氡气场分布的理论方程,通过有限差分算法求解氡气场分布方程,为砂岩型铀矿体上氡气场的正演研究提供了数值计算方法。针对砂岩型铀矿体上氡气场低值区对应矿体而高值异常场对应矿体两端的特点,以矿体端点高值异常峰值点为中心进行二维等步长剖分,得出差分方程,为问题的进一步研究提供一种算法思路。 相似文献
8.
研究了电法测井不同测量方式(井—地,地—井,井—井)下点源场井中电法的三维有限差分数值模拟。采用六面体网格剖分方式来对模型进行剖分,运用一维非零元素行压缩存储模式来存储系数矩阵,减少了内存需求和计算量;采用不完全 Cholesky 共轭梯度(ICCG)方法来求解线性方程组,提高了求解效率;编制了相应的程序实现了井—地、地—井、井—井、倾斜井条件下的电法测井三维有限差分数值模拟。设计的算例结果验证了该算法的正确性和效率性,并且分析了各种情形下的异常特征,为进一步的反演工作打下了基础。 相似文献
9.
建立了土质陡坡降雨侵蚀的数学模型,模型包括陡坡细沟间坡面流控制方程、变沟宽陡坡细沟流控制方程、陡坡细沟间坡面侵蚀泥沙连续方程及陡坡变沟宽细沟侵蚀泥沙连续方程。通过紊流冲击分布概率确定土壤剥蚀方程,确定了泥沙源/汇项的表达形式及细沟高程与宽度的表达式。通过有限元和有限差分相结合的方法求解该数学模型,利用有限单元法对数学模型的水动力方程、运动方程及泥沙连续方程进行空间上的离散,利用有限差分法解决时间域的问题,得出了顺序求解的数值计算公式及模型数值求解的具体方法和步骤。 相似文献
10.
11.
12.
时域有限差分法是麦氏方程的差分形式的数值解法,用于求解任意复杂目标与电磁场的相互作用问题。本文把这一理论应用于地质雷达的数值模拟,得到满意的结果。 相似文献
13.
根据汉口新业大厦水文地质与工程地质条件,建立开采地下水非稳定流定解问题三维数学模型,利用现有承压含水层非稳定流定解问题数值方案求解,重点讨论了源汇项的具体数学表达式。又简介根据达西定律和水均衡原理,建立了三维有限差分数值方程及其解法。根据土的有效应力原理和压密原理推导出了抽水沉降非线性压密方程,并研制了相应的三维有限元和有限差分计算程序。最后,对该环境岩土工程问题进行了分析和评价,计算结果与实测结果符合较好。 相似文献
14.
有限差分法求解分析抗滑桩(或预应力锚索抗滑桩)常采用Winkle弹性地基梁法。在对有限差分法分析的基础上,文章改进了有限差分法求解抗滑桩全桩内力和变形的计算方法。对全桩整体采用有限差分方法求解其内力,避免了将桩身分为受荷段和锚固段及利用滑面处连续条件求解的复杂过程。该方法对滑坡推力与地基系数取值,可依据实际情况假定为广义形式来模拟土体对抗滑桩的作用工况。能很好地求解预应力锚索作用于桩身或桩身设多道锚索工况下的锚索拉力值。并利用Excel对该方法编程,能便捷地计算抗滑桩内力。分析计算结果,当差分节点单元长度适当时,计算结果能够较好地满足工程应用。对比有限元计算分析结果,相同节点单元长度下两者计算相差较小,且节点长度越小,越接近有限元结果。 相似文献
15.
16.
殷文 《吉林大学学报(地球科学版)》2008,38(1):144-151
在弹性波频率空间域有限差分数值模拟方面,差分网格及边界条件是影响弹性波模拟成功与否的关键,为了压制数值模拟中的网格频散,采用25点有限差分算子,建立了有限差分矩阵方程,且借鉴匹配层衰减边界条件思想,设计了弹性波频率空间域有限差分数值模拟算法。由于采用高阶有限差分法来提高差分格式的精度,将会导致计算量显著增加,为此,对频率空间域有限差分弹性波数值模拟方法,采用流水线技术与分治策略进行了并行算法研究,提高了计算效率,使得在合理的计算时间内更精确地模拟弹性波在弹性介质中的传播过程。 相似文献
17.
孙讷正 《水文地质工程地质》1982,(4)
在解决实际问题时主要依靠数值方法求解。带有初边值条件的水动力弥散方程的数值解法与地下水流动方程的数值解法基本相同,但也有一些不同的地方。目前这方面的研究文献较多。我们将选择一些最基本的方法进行介绍,有关解水流方程数值方法方面的基本知识详见文献[1]。 有限差分方法与特征法 (1)有限差分方法的基本思想是:按时间步长△t和空间格距△x、△y、△z将时间和空间区域剖分成若干网格,用未知函数在网格点上的值所构成的差分 相似文献
18.
交错网格有限差分算法以其高效、精确、实用等优点在地震波数值模拟中得到广泛应用。目前交错网格有限差分的精度已达到时间4阶、空间2N阶;然而在求空间三次导数时,差分格式实际上并未达到所谓的2N阶精度,而是采用了低阶的差分格式及差分系数,这样有利于提高大尺度空间正演时的计算效率;但从计算精度的角度考虑,有必要推导出准确的满足2N阶精度的交错网格有限差分格式及差分系数,以得到更高精度的正演结果。笔者利用Taylor公式展开首次推导出了可导函数任意次导数的任意偶数阶精度的差分近似式及相应的差分系数,从而完善了常规高精度交错网格有限差分算法。采用新推导的交错网格有限差分格式得到的正演波形与解析解进行了对比,证明了新推导的差分格式的正确性,并与常规差分格式的正演波形进行了比较,结果显示,新推导出的交错网格有限差分格式模拟结果稳定性好,精度更高。 相似文献
19.