粘性流体两相多孔介质非线性动力问题的罚有限元法

采用基于混合物理论的多孔介质模型,给出粘性流体饱和两相多孔介质非线性动力问题的控制场方程以及相应边值和初值问题的提法,用Galerkin加权残值法导出罚有限元公式,并给出该非线性方程组的迭代求解方法。考虑了体积分数和渗透率与变形相关的情况。用编制的有限元程序计算分析了一维多孔柱体在脉冲载荷作用下的瞬态响应,数值结果表明文中方法正确有效。     

各向异性粘弹性多孔介质中平面波的传播

讨论了弹性多孔介质中波的传播的（或许是）最一般的模型。考虑的介质是粘弹性的、各向异性的、多孔固体骨架，其各向异性可渗透的孔隙中充满着粘性液体。考虑一般类型的各向异性，并且介质中的衰减波作为非均质波处理。对介质中4种衰减波中的每一种，将复慢矢量分解定义为相速度、均质衰减、非均质衰减和衰减角。用一个无量纲参数来度量非均质波与其均质波的区别。利用北海沙岩的数值模型，分析传播方向、非均质参数、频率范围、各向异性对称性、骨架滞弹性和孔隙流体粘度，对该类介质中波传播特性的影响。     

多孔介质瞬态分析中非分裂PML及时域有限元实现

在波场的数值模拟中,完全匹配层(perfectly matched layer,PML)已经被证明是一种十分有效的吸收技术,并得到了广泛的应用.为了解决具有无限域的多孔介质中2阶弹性波动方程数值模拟中的吸收边界问题,提出了一种非分裂格式的PML(non-splitting perfectly matched layer,NPML).首先,基于Biot多孔介质波动理论,建立了以固相和流相位移表示的2阶动力控制方程,其中考虑了固体颗粒和孔隙流体的可压缩性、惯性以及孔隙流体的粘性.其次,根据复伸展坐标变换的定义,通过Laplace变换获得了非分裂格式PML的频域表达式.然后,借助辅助函数将该方程变换到时间域内,得到了一种有效的非分裂PML.最后,基于Galerkin近似方法,给出了其时域有限元计算格式.通过数值算例分析了该非分裂格式的PML在饱和介质动力响应分析中的有效性.     

横观各向同性含液饱和多孔介质中应力波传播的特征分析

根据广义特征理论,对横观各向同性含液饱和多孔介质中应力波传播特性进行了特征分析．给出了特征曲面的微分方程以及沿次特征线的相容条件,得到了波阵面的解析表达式．详细地讨论了应力波在横观各向同性含液饱和多孔介质中传播时,其速度曲面和波阵面的形状及性质．分析结果亦表明,纯固体中应力波传播的特征方程,是含液饱和多孔介质中应力波特征方程的特例．     

孔隙介质中磁热粘弹波的传播

这一课题的内容和意义以及历史发展,请见文献[1-4]。 本文研究两相介质(饱含流体的多孔介质)在磁场作用下,当固相骨架为导电、导热的粘弹性体时地震波的传播。首先,用文[23]的思想建立数学模型,直接耦合流动、应力场和电磁场;其次,研究波的传播规律,具体     

饱和多孔介质的超粘弹性本构理论研究

为了建立能考虑固体材料、多孔固体与流体可逆和不可逆变形的饱和多孔介质超粘弹性理论,以多孔固相为参考构型,以有效应力、材料真实应力和流相真实孔压作为状态变量,结合混合物均匀化响应原理获得各项均符合热力学功共轭特征的饱和多孔介质能量平衡方程,根据非平衡热力学熵分解理论求得熵流和熵产.结果表明,超弹塑性理论是该理论的一个特例;多孔固体的总变形可分为固相间隙和材料变形两部分,间隙应变与Terzaghi有效应力构成功共轭对,材料应变与材料真实应力构成功共轭对.饱和多孔介质的自由能可分为固相和流相两部分.当固相间隙和材料变形解耦时,固相所含的自由能又可分为间隙和材料两部分.证明了Skempton有效应力不是饱和多孔介质的基本应力状态变量.     

二维各向同性多孔介质的弹性动力学通解

在二维直角坐标系下,从固体位移和流体流速满足的基本方程出发,研究了二维各向同性多孔介质的弹性动力学通解.首先引入4个物理量,对固体骨架的运动方程、流体流速运动方程、连续性方程进行整理,将方程组分解成膨胀波和扭转波两部分,并利用Lur’e算子矩阵理论,获得由3个类调和函数表示的动力学通解,该通解满足全部基本方程.最后将时间项退化获得稳态通解,并证明了稳态通解的完备性.     

弹性固体和充满两种互不相溶粘性流体的多孔固体之间界面的反射和折射及其衰减波

在充满两种互不相溶粘性流体的多孔固体中,研究弹性波的传播.用3个数性的势函数描述3个纵波的传播,用1个矢性的势函数单独描述横波的传播.根据这些势函数,在不同的组合相中,定义出质点的位移.可以看出,可能存在3个纵波和1个横波.在一个弹性固体半空间与一个充满两种互不相溶粘性流体的多孔固体半空间之间,研究其界面上入射纵波和横波所引起的反射和折射现象.由于孔隙流体中有粘性,折射到多孔介质中的波,朝垂直界面方向偏离.将入射波引起的反射波和折射波的波幅比,作为非奇异的线性代数方程组计算.进一步通过这些波幅比,计算出各个被离散波在入射波能量中所占的份额.通过一个特殊的数值模型,计算出波幅比和能量比系数随入射角的变化.超过SV波的临界入射角,反射波P将不再出现.越过界面的能量守恒原理得到了验证.绘出了图形并对不同孔隙饱和度以及频率的变化,讨论它们对能量分配的影响.     

饱和土中波传播的二相动力分析模型

将饱和土作为二相介质,考虑土骨架和孔隙流体之间的摩阻力,建立了饱和土中一维压缩波传播的二相动力分析控制微分方程,采用有限差分法求解。该模型可以计算地表面爆炸冲击荷载作用下饱和土中波的传播和波在不动障碍上的反射问题。该模型可考虑饱和土土体的非线性塑性力学性质及土性随深度变化或土介质的分层情况。     

多孔磁流体动力学

本文建立饱含可压缩电磁流体的多孔固体介质动力学。直接耦合应力场、流场和电磁场,研究了它们的相互作用和波的传播规律,纠正了某些作者在磁流体波中的误失,证明了介质中具有多种P波和S波,一般多有阻尼和频散效应。 导电流体(如海水、石油和血液等)在磁场(如地磁场和生物场等)中运动感生电流因而改变磁场;电流在磁场中运动则产生机械力,影响固体骨架的变形与流体的运动;固体骨架的变形和空隙通道中的Darcy渗流又反过来改变流体的运动。新理论的困难和可以预料的丰富的内容正来自复杂的相互作用,不妨称作(微)百幕大效应。相应的守恒律、等价原理、传输线理论相似性、电磁效应、非线性波、动应力集中、变分原理以及模型的进一步推广将在另一组论文中发表。 深入研究这一理论对地球物理、海洋科学、生命科学以及能源科学无疑是重要的。多孔磁流体介质动力学是一门年轻的学科,大量未解决的问题和实验工作期待着人们去努力!     

流体饱和多孔隙介质波动方程小波有限差分法

研究流体饱和多孔隙介质中波动方程的数值模拟．针对求解二维弹性波方程问题,提出小波有限差分法．该方法综合了小波多分辨分析计算灵活、计算效率高特性和有限差分易于实现的优点．数值模拟的结果显示,此方法对于求解流体饱和多孔隙介质方程的数值模拟是有效稳定的．     

Spectral element method for acoustic wave simulation in heterogeneous media

In this paper, we present a spectral element method for studying acoustic wave propagation in complex geological structures. Due to complexity (both lithological and stratigraphical), the use of numerical methods of higher accuracy and flexibility is needed to achieve the correct results. The spectral element method shows more accurate results compared to the low-order finite element, the conventional finite difference and the pseudospectral methods. High accuracy is reached even for rather long wave propagation times and dispersion errors are essentially eliminated; pirregular interfaces between different media can be well described so that numerical artifacts or noises are not at all introduced. The method is tested against analytical solutions both in the two-dimensional homogenous and heterogeneous media. The results of different simulations are presented.     

On the Matrix Method in the Theory of Wave Propagation in Layered Biot Media

Wave propagation in porous Biot media with homogeneous isotropic layers is investigated by the matrix method. In order to use this method, matrices describing porous layers and half-spaces are established. On the basis of these matrices, the wave fields in layered porous media are derived and examined. In this paper, matrices of inhomogeneous porous layers are also determined, and these matrices are represented by convergent series. Bibliography: 9 titles.     

含液体非弹性多孔介质中波传播过程的失稳与逸散性

在基于Biot理论的饱和-非饱和多孔介质的动力-渗流模型中计及流固惯性耦合效应。对单轴应变的一维情况讨论了饱和和非饱和多孔介质中波传播过程的驻值失稳和逸散性,分析了流固粘性耦合,流固惯性耦合,流固混合体的压缩性,孔隙水饱和度,及固体骨架在高应变速率下材料粘弹塑性本构行为等因素的影响。该工作将对克服饱和与非饱和多孔介质在强动荷载下波传播过程的数值求解困难提供理论上的依据和启示。     

Infinite elements in a poroelastodynamic FEM

Wave propagation phenomena in unbounded domains occur in many engineering applications, e.g., soil structure interactions. When simulating unbounded domains, infinite elements are a possible choice to describe the far field behavior, whereas the near field is described through conventional finite elements. Finite element formulations for porous materials in terms of Biot's theory [1] have been published, e.g., by Zienkiewicz [2]. For infinite elements, several approaches are described in [3,4]. Infinite elements are based on special shape functions to approximate the semi-infinite geometry as well as the Sommerfeld radiation condition, i.e., the waves decay with distance and are not reflected at infinity. If there is only one wave traveling in the media, a formulation in time-domain can be established. But in poroelastodynamics, there are three body waves and eventually also a Rayleigh wave. Unfortunately, the extension to more than one wave is not straight forward. Here, an infinite element is presented which can handle all wave types, as it is needed in poroelasticity. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Three-dimensional finite/infinite elements analysis of fluid flow in porous media

Fluid flow in naturally fractured porous media can always be regarded as an unbounded domain problem and be better solved by finite/infinite elements. In this paper, a three-dimensional two-direction mapped infinite element is generated and combined with conventional finite elements and one direction infinite element to simulate poroelasticity. Therefore, the entire semi-infinite domain can be included in the numerical analysis. Both single- and dual-porosity porous media are considered. For the purpose of validation, we compare the results of finite/infinite elements with those of finite elements under two extreme boundary conditions. The comparison indicated that mapped infinite element is an appropriate approach to model fluid flow in porous media and provides an intermediate solution.     

固体中短波传播的单位分解有限元法

提出了固体中短波传播数值模拟的单位分解有限元法．有限元空间由形成单位分解的标准等参有限元形函数乘以定义为局部子空间基函数的特殊形函数构成．特殊形函数使试空间中包含了关于波动方程的已有知识,因而在单个单元内能近似地再现高度振荡性质．数值例题显示了所提出单位分解有限元在计算精度和效率上的良好性能．     

A fully conservative block-centered finite difference method for simulating Darcy-Forchheimer compressible wormhole propagation

Numerical Algorithms - In this paper, a block-centered finite difference method is derived for the Darcy-Forchheimer compressible wormhole propagation in porous media by introducing an auxiliary...     

Extended finite element modeling of deformable porous media with arbitrary interfaces

In this paper, an enriched finite element method is presented for numerical simulation of saturated porous media. The arbitrary discontinuities, such as material interfaces, are encountered via the extended finite element method (X-FEM) by enhancing the standard FEM displacements. The X-FEM technique is applied to the governing equations of porous media for the spatial discretization, followed by a generalized Newmark scheme used for the time domain discretization. In X-FEM, the material interfaces are represented independently of element boundaries and the process is accomplished by partitioning the domain with some triangular sub-elements whose Gauss points are used for integration of the domain of elements. Finally, several numerical examples are analyzed, including the dynamic analysis of the failure of lower San Fernando dam, to demonstrate the efficiency of the X-FEM technique in saturated porous soils.     

Energetic BEM-FEM coupling for the numerical solution of the damped wave equation

Time-dependent problems modeled by hyperbolic partial differential equations can be reformulated in terms of boundary integral equations and solved via the boundary element method. In this context, the analysis of damping phenomena that occur in many physics and engineering problems is a novelty. Starting from a recently developed energetic space-time weak formulation for the coupling of boundary integral equations and hyperbolic partial differential equations related to wave propagation problems, we consider here an extension for the damped wave equation in layered media. A coupling algorithm is presented, which allows a flexible use of finite element method and boundary element method as local discretization techniques. Stability and convergence, proved by energy arguments, are crucial in guaranteeing accurate solutions for simulations on large time intervals. Several numerical benchmarks, whose numerical results confirm theoretical ones, are illustrated and discussed.