An analytical solution for the transient response of saturated linear elastic porous media

The motions of fluid and solid phases in saturated porous media are coupled by inertial, viscous and mechanical interactions as described by Biot's equations. A one-dimensional exact analytical solution of the Biot's equations for the completely general solution of the transient problem in saturated, linear, elastic, porous media is presented. The problem is solved by using the Fourier series. The transient response of porous media is shown for typical material properties of a natural granular deposit and for different degrees of viscous coupling. The analytical results show the mechanics of dispersive wave propagation in saturated porous media and they should provide a useful comparison term for the existing numerical solution methods.     

二维饱和多孔介质因点汇诱发比奥固结的解析解

给出了有限二维饱和多孔介质因点汇诱发的Biot固结的一个解析解。其中假设多孔介质为均匀各向同性和线弹性,假设孔隙压力场符合第1类边界条件,数学模型采用可压缩多孔介质模型。利用傅里叶和拉普拉斯变换及相应反演获得了双重无穷项级数和形式的精确解。然后特别探讨了定流量点汇诱发的稳态解析解,并用文献现有解析解进行了验证。所提出的解析解适合于验证数值解,并可用于深入分析有限二维多孔介质的流-固耦合行为。     

饱和土半空间中圆柱形孔洞对平面P波的散射

在 Biot饱和多孔介质动力学理论的基础上,首次建立了求解饱和土半空间中圆柱形孔洞对平面 P 波散射问题的波函数展开法。首先。分析了具有圆柱形孔洞的饱和土半空间场地在平面 P 波入射下产生散射波系,并将入射波和散射波的波函数在圆柱坐标下展开。然后,引入边界条件,求出散射波函数的待定系数,从而,得到饱和土半空间中圆柱形孔洞对平面P波的散射问题的解析解。根据所得的波函数的解,可求解区域内的位移、应力的值,同时,分析了入射波频率、入射角对柱面上的应力集中因子的影响。     

Parallel simulations of pore-scale flow through porous media

Smoothed particle hydrodynamics (SPH) is a versatile technique which can be applied to single and multiphase flow through porous media. The versatility of SPH is offset by its computational expense which limits the practicability of SPH for large problems involving low Reynolds number flow. A parallel pore-scale numerical model based on SPH is described for modeling flow phenomena in porous media. Aspects of SPH which complicate parallelization are emphasized. The speed of the method is demonstrated to be proportional to the number of processors for test cases where load balance was achieved. The parallel algorithm permits the application of SPH to more complicated porous media problems than previously considered. For such problems, best performance is achieved when several soil grains are simulated by each processor. Finally, future applications of the method and possible extensions are discussed.     

横观各向同性含液饱和多孔介质黏弹性瞬态反应

在考虑横观各向同性含液饱和多孔介质固体骨架和流体可压缩性以及固体骨架的黏弹性特征下,基于横观各向同性含液饱和多孔介质u-w形式的三维动力控制方程,以固相位移u、液相相对位移w为基本未知量,综合运用Laplace变换、双重Fourier变换等方法,在直角坐标系下通过引入中间变量,将六元2阶动力控制方程组化为两组各含4个未知变量的常微分方程组,给出了直角坐标系下横观各向同性含液饱和多孔介质三维黏弹性动力反应的积分形式一般解;作为理论推导的验证,通过引入初始条件和边界条件,对横观各向同性含液饱和多孔介质半空间黏弹性瞬态反应问题进行了求解。解答的退化验证表明,所推导的理论解是正确的。     

An analytical solution for axisymmetric tunnel problems in elasto-viscoplastic media

An analytical solution is presented for the time-dependent stresses and displacements in plane strain around a circular hole when it is loaded by an axisymmetric internal and far-field pressure. The material is assumed to be elasto-viscoplastic with dilatant plastic deformations according to a non-associated flow rule. Strain softening is considered by a modified St. Venant slider which is characterized by Mohr-Coulomb yield conditions for both the peak and the residual strengths.     

A pore-scale numerical model for flow through porous media

A pore-scale numerical model based on Smoothed Particle Hydrodynamics (SPH) is described for modelling fluid flow phenomena in porous media. Originally developed for astrophysics applications, SPH is extended to model incompressible flows of low Reynolds number as encountered in groundwater flow systems. In this paper, an overview of SPH is provided and the required modifications for modelling flow through porous media are described, including treatment of viscosity, equation of state, and no-slip boundary conditions. The performance of the model is demonstrated for two-dimensional flow through idealized porous media composed of spatially periodic square and hexagonal arrays of cylinders. The results are in close agreement with solutions obtained using the finite element method and published solutions in the literature. Copyright © 1999 John Wiley & Sons, Ltd.     

Effective two-phase flow through highly heterogeneous porous media: Capillary nonequilibrium effects

We consider the two-phase flow through a dual-porosity medium, characterized by a period of heterogeneity ω, a ratio of global permeabilities ∈_K, and a ratio of the order of capillary forces ∈_c. The limit when ω tends to zero at different values of ∈_K and ∈_c gives four classes of global behavior, differing by the type of elementary flows at the one-cell level. We propose a diagram of their predominance. A macro-scale model is constructed by formal homogenization techniques for one of these classes; it shows a nonlinear kinetic relationship for the averaged capillary pressure functions, and leads to a decomposition for the effective phase permeability tensors. A capillary relaxation time is explicitly determined. This revised version was published online in August 2006 with corrections to the Cover Date.     

Flow through porous media: A procedure for locating the free surface

A finite element procedure is developed to accurately locate the free surface of unconfined seepage flow through porous media. The free surface is taken as the boundary between wet and dry soils, with flow in the saturated region characterized by Darcy's law. The method involves equations and meshing which are fully consistent with a general formulation for geotechnical engineering problems involving simultaneous solution of pore fluid pressures and soil skeleton displacements. Accuracy and versatility of the proposed procedure are demonstrated by solving various unconfined seepage flow problems through earth structures. Free surfaces and flownets are presented for the calculated flow fields.     

A frequency-domain analytical solution of Maxwell's equations for layered anisotropic media

An analytical solution of Maxwell's equations for layered anisotropic media is presented in a form which allows estimating the sought parameters by layer stripping without round-off accumulation. The solution in each layer is reduced to the standard procedures of solving a fourth-order algebraic equation, multiplication, addition, and inversion of second-order non-singular matrices. The algorithm has no limitations on layer thickness and is applicable to both very thick and very thin layers. The new numerical code is straightforward and can be easily parallelized.     

An analytical solution for the transient response of saturated porous elastic solids

The general forms for the field equations governing the transient response of poroelastic media given by Biot and by Zienkiewicz are compared and relations between the material constants are obtained. A one-dimensional analytical solution is presented for the situation where the solid and fluid materials satisfy Biot'S dynamic compatibility relation. The transient response of porous media is illustrated for varying degrees of solid and fluid compressibility when subjected to step, cyclic and short duration spike surface tractions. The results obtained (for the special situation where the materials are dynamically compatible) exhibit the overall characteristics of wave propagation in porous media and will provide representative test problems which allow a quantitative evaluation of the accuracy of various numerical solution methods (e.g. finite element models).     

A mixed solution for two-dimensional unsteady flow in fractured porous media

A mixed finite element–boundary element solution for the analysis of two-dimensional flow in porous media composed of rock blocks and discrete fractures is described. The rock blocks are modelled implicitly by using boundary elements whereas finite elements are adopted to model the discrete fractures. The computational procedure has been implemented in a hybrid code which has been validated first by comparing the numerical results with the closed-form solution for flow in a porous aquifer intercepted by a vertical fracture only. Then, a more complex problem has been solved where a pervious, homogeneous and isotropic matrix containing a net of fractures is considered. The results obtained are shown to describe satisfactorily the main features of the flow problem under study. © 1997 by John Wiley & Sons, Ltd.     

随机介质热弹性力学模型球坐标问题的解析解

针对随机介质热弹性力学模型,在随机韦泊分布和指数分布下,提出了一种解决随机概率分布用连续函数表示的方法,并推导出该模型球坐标问题的解析解.通过算例,对随机模型与均质模型热应力的区别做了详细研究,结果表明随机概率分布参数m对随机模型的热应力有重要影响.     

Experimental study on the velocity-dependent dispersion of the solute transport in different porous media

The hydrodynamic dispersion is an important factor influencing the reactive solute transport in the porous media, and many previous studies assumed that it linearly varied with the average velocity of the groundwater flow. Actually, such linear relationship has been challenged by more and more experimental observations, even in homogeneous media. In this study, we aim to investigate the relationship between hydrodynamics dispersion and the flow velocity in different types of porous media through a laboratory-controlled experiment. The results indicate that(1) the dispersion coefficient should not be a linear function of the flow velocity when the relationship between the flow velocity and the hydraulic gradient can be described by Darcy's law satisfactorily;(2) Power function works well in describing the dispersion coefficient changing with the flow velocity for different types of porous media, and the power value is between 1.0-2.0 for different particle sizes.     

Experimental study on the velocity-dependent dispersion of the solute transport in different porous media

The hydrodynamic dispersion is an important factor influencing the reactive solute transport in the porous media, and many previous studies assumed that it linearly varied with the average velocity of the groundwater flow. Actually, such linear relationship has been challenged by more and more experimental observations, even in homogeneous media. In this study, we aim to investigate the relationship between hydrodynamics dispersion and the flow velocity in different types of porous media through a laboratory-controlled experiment. The results indicate that (1) the dispersion coefficient should not be a linear function of the flow velocity when the relationship between the flow velocity and the hydraulic gradient can be described by Darcy’s law satisfactorily; (2) Power function works well in describing the dispersion coefficient changing with the flow velocity for different types of porous media, and the power value is between 1.0-2.0 for different particle sizes.     

Finite element evaluation of diffusion and dispersion tensors in periodic porous media with advection

This paper presents three-dimensional finite element simulations to evaluate diffusion and dispersion tensors in periodic porous media in the presence of an advective velocity field. These tensors are evaluated in the framework of the double-scale expansion technique. Two problems, a Newtonian flow and a vector-valued advection–diffusion equation, have to be sequentially solved at the pore scale. Finite element techniques to approximate these problems are proposed and analyzed. Numerical results in three-dimensional networks of spheres are presented to quantitatively assess the impact of the pore morphology and of the advection velocity on the diffusion and dispersion tensors.     

单层非饱和多孔介质一维瞬态响应半解析解

基于Zienkiewicz提出的非饱和多孔介质波动理论,考虑两相流体和固体颗粒的压缩性以及惯性、黏滞和机械耦合作用,采用半解析的方法获得了一类典型边界条件下单层非饱和多孔介质一维瞬态响应解。首先推导出无量纲化后以位移表示的控制方程,并将其写成矩阵形式;然后,将边界条件齐次化,求解控制方程所对应的特征值问题,得到了满足齐次边界条件的特征值和相对应的特征函数。根据变异系数法并利用特征函数的正交性,得到了一系列仅黏滞耦合的关于时间的二阶常微分方程及相应的初始条件。在此基础上,运用精细时程积分法给出了常微分方程组的数值解。最后,通过若干算例验证了结果的正确性并探讨了单层非饱和多孔介质一维瞬态动力响应的特点。该方法可推广应用于其他典型的边界条件。     

Explicit meshfree solution for large deformation dynamic problems in saturated porous media

In this paper a new methodology to simulate saturated soils subjected to dynamic loadings under large deformation regime (locally up to 40% in equivalent plastic strain) is presented. The coupling between solid and fluid phases is solved through the complete formulation of the Biot’s equations. The additional novelty lies in the employment of an explicit time integration scheme of the \(u-w\) (solid displacement–relative fluid displacement) formulation which enables us to take advantage of such explicit schemes. Shape functions based on the principle of maximum entropy implemented in the framework of Optimal Transportation Meshfree schemes are utilized to solve both elastic and plastic problems.     

基于梯度塑性模型的多孔介质应变局部化分析的基本理论

讨论了引入梯度塑性模型情况下的含液多孔介质双重内尺度律的特征以及相互作用问题。在引入梯度塑性本构模型的情况下,一方面材料应变局部化分析数值结果的正则性得到保证;而另一方面,含液多孔介质会出现双重内尺度律参数的共同作用问题。本文给出了此时内尺度律预测的一个基本方法,并对稳定性问题进行了分析,讨论了不同情况下实波速存在的条件,给出了对于给定的渗透系数情况下实波速存在的波数区间,并对相关现象进行了解释。