有限元有限差分法二维波动逆时偏移初探

提出了采用有限元有限差分实现二维波动方程的逆时偏移算法。该方法在空间上 ,联合采用有限元法和有限差分法 ;对于地表 (水平 )方向 ,使用有限元法进行离散 ,将原方程转化为一个一维 (深度和时间 )问题的方程组 ;在深度和时间方向上 ,采用有限差分法来求解。介绍了算法的基本原理 ,给出了计算实例并与使用F K(频率波数 )域相移法、频率空间域有限差分法的结果进行了比较。与采用有限元的偏移方法相比 ,本方法可以节省大量内存 ;与采用有限差分的偏移方法相比 ,可以在一定程度上提高计算精度。本算法有可能在地震勘探数据处理中发挥一定的作用     

有限元-有限差分法二维波动逆时偏移初探

提出了采用有限元-有限差分实现二维波动方程的逆时偏移算法。该方法在空间上，联合采用有限元法和有限差分法；对于地表(水平)方向，使用有限元法进行离散，将原方程转化为一个一维(深度和时间)问题的方程组；在深度和时间方向上，采用有限差分法来求解。介绍了算法的基本原理，给出了计算实例并与使用F-K(频率-波数)域相移法、频率-空间域有限差分法的结果进行了比较。与采用有限元的偏移方法相比，本方法可以节省大量内存；与采用有限差分的偏移方法相比，可以在一定程度上提高计算精度。本算法有可能在地震勘探数据处理中发挥一定的作用。     

非达西渗流时一维固结方程的两种数值解法

Terzaghi一维固结理论对低渗饱和地基沉降的预测值与实际值之间存在较大误差,原因之一在于其假定地基的渗流符合Darcy定律.在采用考虑起始水力梯度的非Darcy渗流方程和承认Terzaghi一维固结理论其他假定的基础上,重新推导了饱和粘土一维固结方程,并分别给出了有限差分法和有限体积法的离散格式.对单面排水情况下渗流前锋和固结度的计算结果证明了两种数值方法的适用性,而有限体积法具有更好的收敛性.     

模拟自由面渗流的适体坐标变换方法

用传统的有限元方法求解复杂边界的自由面渗流是很困难的,为此提出了基于适体坐标变换的有限差分法。该方法通过求解Poisson方程自动生成计算区域的曲线网格,并将其变换至规则统一的直角网格系统下进行有限差分离散和数值求解。这种数值网格生成技术可以精确而有效率地模拟复杂几何边界,算法成熟,通用性强,可以避免有限元法的网格在迭代过程中变化的问题,对不同的几何边界可以实现统一的数值求解算法,自动化程度高。计算实例表明,基于适体坐标变换的有限差分法计算结果的精度和有限元法相当,但计算效率上较有限元方法有很大的提高,并具有简单、灵活的优点。     

煤层瓦斯流固耦合渗流的二维数值模拟

将渗流力学与弹塑性力学相结合,考虑煤层瓦斯和煤体骨架之间的相互作用,建立了煤层瓦斯运移的数学模型,并给出其数值解。对煤层瓦斯渗流方程和煤体变形场方程分别进行离散,得到其相应的泛函方程,根据有限元法原理推导出其耦合求解方法,最后进行了实例计算。     

静压气体球轴承承载力的三维计算方法

为研究环面节流静压气体球轴承的承载力特性,提出一种基于三维模型的计算方法.应用非结构同位网格技术对轴承进行网格划分,采用有限体积法对三维稳态可压缩N-S方程进行离散,将适用于可压缩流体的改进的SIMPLE算法应用于离散方程的求解.实验表明,该方法的计算结果与二维有限元法计算结果及实验数据相符,证明算法是有效的.     

泡沫油流渗流规律数学模型及其数值解

为了准确模拟泡沫油流在储层内的流动规律的从渗流力学基本方程出发,推导处泡沫油流渗流基本方程。利用有限差分方法求得其数值解。将计算结果与黑油模型以及室内实验进行了对比,结果表明,所提出的模型更接近实验结果。     

用径向基函数法计算不稳定渗流达西速度

离散地下水不稳定渗流方程的时间变量t,得到半离散差分格式,再用径向基函数配点法得到全离散格蓼,由此提出了达西速度的径向基函数计算方法,该方法较有限元法简单,节点配置较有限元和有限差分等方法灵活,在高维空间中比传统方法更具优势．模型试算表明,其相对误差到10^-5,比目前精度较高的混合有限元法缩小10^3倍．     

求解正压原始方程的一个算法

本文给出了有限元法离散正压原始方程的推导过程,并采用选点迭代法求解有限元离散后的线性方程组,以提高运算速度。数值天气预报是一种客观定量的天气预报的方法。在数值天气预报中使用最多的是有限差分法。有限元法较有限差分法来说在剖分形式上更加灵活,能够更好地刻划边界条件,在相同的情况下,精度也更高。但有限元法的计算量要比有限差分法大一些,并且对计算机的要求也比有限差分法高一些。目前,有限元法已得到了数值天气预报工作者的广泛重视,并对许多具体的局部问题展开了讨论。本文给出了一种近似求解正在原始方程组的有限元算法,主要是针对减少计算量和求介稳定性所作的工作。     

油藏热流固耦合渗流问题的有限元方法

提出了一种有效实用的求解油藏热流固耦合渗流问题的数值计算方法。该方法以有限元法为主,结合有限差分法,用有限元法求解耦合温度场方程和岩石耦合变形方程,用有限差分法求解流体耦合渗流方程,发挥有限元法网格技术和单元划分灵活的特点及处理复杂的油藏边界优势,兼顾了有限差分法在流场分析方面的成熟应用,使复杂的热流固耦合数学模型得以完整求解,取得了单由有限元法或有限差分法难以取得的效果,是一种新型的油藏数值模拟方法。     

两相渗流驱动问题有限元—特征线修正差分格式及理论分析

研究多孔介质中可压缩可混溶两相渗流驱动问题的计算方法：压力方程用有限元方法求解，饱合度方程用特征线修正有限差分方法求解，构造了全离散数值计算格式，证明了最佳收敛阶     

有限体积法和有限元方法之间的比较

有限体积法现在已经成为和有限元方法并驾齐驱的一种求解偏微分方程的数值方法。与有限元方法相比,有限体积法保持物理量的局部守恒性质,并且计算更加简单。本文主要介绍有限体积法和有限元法之间的一些相同点和不同点。     

裂缝性气水两相地层井筒组合数值模拟

气水井试井对天然气开发具有重大意义。借助数值模拟方法了解地下气流和水流的变化,对于提高采收率是极为重要的。通常的数值模拟中仅考虑了储层的渗流动态,而把井筒视为线汇处理,这不能准确的反映出井筒和储层之间的相互影响。建立了气水两相组合模型,借助有限体积法导出了组合模型的离散方程,并讨论了求解的边界条件。实算表明了这种方法能更有效地反映出试井过程中的各种现象。     

两相渗流驱动问题有限元——特征线修正差分格式及理论分析

研究多孔介质中可压缩可混溶两相渗流驱动问题的计算方法，压力程用有限元方法求解，饱合度方程用特征线修正有限差分方法求解，构造了全离散数值计算格式，证明了最佳收敛阶。     

四边形网上双曲型方程的有限体积元法

有限体积元法已引起国内外学者和专家的广泛关注,该方法与有限元方法有着相同的收敛阶,具有计算简单,保持物理量守恒性等优点.讨论一类双曲型方程在四边形网上的有限体积元法,在四边形网格单元满足h2拟平行四边形条件下,给出了双曲型方程半离散有限体积元格式下最优的H1模和L2模误差估计以及两个全离散有限体积元格式下的误差估计.     

一种用于结构动态过程数值模拟的广义有限元法

提出了一种能同时处理连续体和离散体动态过程的有限元法. 该法通过构造接触力模型, 将传统的离散元当做标准的一节点有限单元来看待. 新构造的一节点单元与常规有限单元具有同样物理特征, 包括应力和应变. 从而, 将离散元算法与有限元算法紧密结合在一起, 形成了统一的、与现有有限元系统兼容的广义有限元计算模型. 这种模型不仅简化了计算过程, 还可最大限度地发挥现有有限元代码的功能. 数值算例表明, 新方法可以更加有效地模拟大量离散体与连续体相互作用的复杂动力学过程.     

二维方腔流动的有限元模拟

采用了Ｔａｙｌｏｒ－Ｇａｌｅｒｋｉｎ有限元法离散Ｎ－Ｓ方程，该法保留了高阶的空间离散精度，隐含了流线迎风的耗散作用，采用了压力校正法求解流场中各原始变量，推导了压力修正Ｐｏｉｓｓｏｎ方程的有限元离散格式，最后给出了二维不可压缩方腔流动计算结果。     

基于非结构化网格有限体积的LBM

提出了一种基于非结构化网格有限体积的LBM.采用Cell-vertex有限体积法离散控制方程.该方法在时间上采用伪、实二时间步,其中伪时间步采用向前差分,实时间步采用二阶向后差分方法;空间上采用edge-based通量计算方法,采用高阶TVD格式计算控制体边界通量.离散后的控制方程采用隐式迭代,控制变量采用五层二阶Runge-Kutta方法求解.二维同心圆环内圆柱间Couette流与顶盖驱动方腔流的数值结果显示该方法为一种有效求解不规则边界流体动力特性的实用工具.     

薄壁铝型材挤压成形的一种有效模拟方法

基于大变形弹塑性有限元理论和有限体积法基本原理，建立了金属塑性成形的弹塑性UL有限元列式以及塑性流动中的有限体积控制方程．提出了有限元模拟系统到有限体积模拟系统的数据传递和信息继承方法。建立了铝型材挤压成形有限元/有限体积法复合模拟系统，对铝型材挤压过程进行了数值模拟，预示金属在成形中的塑性变形行为，从而为模具设计及工艺参数选取提供理论依据．     

轴对称问题中的无网格Galerkin法

采用无网格Galerkin法分析轴对称问题，得到弹性力学中的对称问题的无网格离散方程．将这一方法与有限元耦合，即在边界处布置有限单元，这样就可以用传统有限元方法方便地处理力学边界条件．算例考察表明：本文方法通过了分片检验，计算结果达到了较高的精度，最大误差不超过5％．