孔隙介质的时域BEM计算

根据Biot饱和孔隙介质动力方程,结合快、慢纵波解耦法得到时域Green函数U-P表达以及Somigliana表象积分,采用BEM分析了集中力作用下饱和孔隙介质时域动力响应.详细论述了孔隙介质时域边界积分方程的离散化方法与形式,它的Stokes状态解答和借用已有技术成果对计算奇异性的处理.在无量纲材料参数的数值分析计算中,以图表形式给出结果.由于孔隙介质的时域BEM计算在相关文献中较为罕见,因此文中结果会对两相饱和介质动力响应特性等相关研究提供一些新的途径.     

Effective models of stratified media containing porous biot layers

Periodic stratified media in which either two porous Biot layer, or an elastic and a porous layers, or a fluid and a porous layer alternate are considered. The effective models of these media are constructed and investigated. In the case of alternating porous layers, the effective model is a generalized transversely isotropic Biot medium. In this medium, the density of the fluid phase and the mean density acquire tensor character. It is shown that the effective model of a porous-fluid medium is, on the one hand, a generalized transversely isotropic Biot medium of special type and, on the other hand, a generalization of the effective model of a stratified elastic-fluid medium.Translated fromZapiski Nauchnykh Seminarov POMI, Vol. 239, 1997, pp. 140–163.This work was supported by the Russian Foundation for Basic Research under grant Nos. 96-01-00666 and 96-05-66207.     

On the ray method for Biot media

This paper is devoted to an investigation of wave propagation in a Biot porous medium, which consists of elastic and fluid phases. The space-time ray expansion of solutions of dynamic equations for a Biot medium is constructed (in the anisotropic inhomogeneous case). In the inhomogeneous isotropic case, a Rytov law analog is derived similarly to elasticity theory. Bibliography: 3 titles. __________ Translated from Zapiski Nauchnykh Seminarov POMI, Vol. 354, 2008, pp. 112–131.     

多孔介质平板通道传热模型的两种求解方法

基于Brinkman Darcy扩展模型和非局部热平衡模型,考虑液相和固相含有内热源的情况,建立了多孔介质平板通道传热的一般模型.分别采用直接法和间接法将液相与固相能量方程解耦,进而求得充分发展传热条件下的多孔介质温度场.与直接解耦法相比,间接解耦法可在原始边界条件下求解二阶微分方程,更加简单易行.通过对无量纲温度表达式系数以及温度分布的比较,验证了两种求解方法的等价性.在两种极限情形下,间接法所得温度分布解析解与现有文献结果相当吻合,这也在一定程度上证明了所建模型更具一般性.参数分析表明,液固两相温差随着Biot数或有效导热系数比的增大而减小,Nusselt数随着内热源比的增大而减小.     

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.     

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

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

Consistency of iterative operator‐splitting methods: Theory and applications

In this article, we describe a different operator‐splitting method for decoupling complex equations with multidimensional and multiphysical processes for applications for porous media and phase‐transitions. We introduce different operator‐splitting methods with respect to their usability and applicability in computer codes. The error‐analysis for the iterative operator‐splitting methods is discussed. Numerical examples are presented. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

Wave propagation in an isolated porous Biot layer with closed pores on the boundaries

On the boundaries of such an isolated porous Biot layer, the total stresses and normal relative displacement are equal to zero. For this layer, the symmetric and antisymmetric dispersion equations are established and investigated. The wave field consists of normal waves. In this layer, one bending wave, two plate waves, and infinitely many normal waves propagate. For all these waves, we determine dispersion curves by analytical methods. The velocities of the bending wave and the second plate wave for the infinite frequency are equal to the Rayleigh velocity. Bibliography: 7 titles. __________ Translated from Zapiski Nauchnykh Seminarov POMI, Vol. 354, 2008, pp. 173–189.     

On methods of deriving equations describing effective models of layered media

Methods of deriving equations describing effective models of layered periodic media are presented. Elastic and fluid media, as well as porous Biot media, may be among these media. First, effective models are derived by a rigorous method, and then some operations in the derivation are replaced by simpler ones providing correct results. As a consequence, a comparatively simple and justified method of deriving equations of an effective model is established. In particular, this method allows us to simplify to a degree and justify the derivation of an effective model for media containing Biot layers; this method also produces equations of an effective model of a porous layered medium intersected by fractures with slipping contacts. Bibliography: 15 titles. Translated fromZapiski Nauchnykh Seminarov POMI, Vol. 250, 1998 pp. 219–243. Translated by L. A. Molotkov.     

Théorie spectrale de la propagation des ondes acoustiques dans un milieu stratifié perturbé

In this article the spectral analysis of the self adjoint operator governing the propagation acoustic waves in a perturbed stratified medium is given. Both a medium whose sound speed is a short range perturbation of that of a stratified medium and a stratified medium exterior to a compact set are considered. The basic tool is a division theorem for the self adjoint operator governing the propagation of acoustic waves in a pure stratified medium.     

An analytic method for the initial value problem of the electric field system in vertical inhomogeneous anisotropic media

The time-dependent system of partial differential equations of the second order describing the electric wave propagation in vertically inhomogeneous electrically and magnetically biaxial anisotropic media is considered. A new analytical method for solving an initial value problem for this system is the main object of the paper. This method consists in the following: the initial value problem is written in terms of Fourier images with respect to lateral space variables, then the resulting problem is reduced to an operator integral equation. After that the operator integral equation is solved by the method of successive approximations. Finally, a solution of the original initial value problem is found by the inverse Fourier transform.     

Propagation of nonlinear thermoelastic waves in porous media within the theory of heat conduction with memory: physical derivation and exact solutions

In this paper, we study a mathematical model of nonlinear thermoelastic wave propagation in fluid‐saturated porous media, considering memory effect in the heat propagation. In particular, we derive the governing equations in one dimension by using the Gurtin–Pipkin theory of heat flux history model and specializing the relaxation function in such a way to obtain a fractional Erdélyi–Kober integral. In this way, we obtain a nonlinear model in the framework of time‐fractional thermoelasticity, and we find an explicit analytical solution by means of the invariant subspace method. A second memory effect that can play a significant role in this class of models is parametrized by a generalized time‐fractional Darcy law. We study the equations obtained also in this case and find an explicit traveling wave type solution. Copyright © 2016 John Wiley & Sons, Ltd.     

A comparison of the String Gradient Weighted Moving Finite Element method and a Parabolic Moving Mesh Partial Differential Equation method for solutions of partial differential equations

We compare numerical experiments from the String Gradient Weighted Moving Finite Element method and a Parabolic Moving Mesh Partial Differential Equation method, applied to three benchmark problems based on two different partial differential equations. Both methods are described in detail and we highlight some strengths and weaknesses of each method via the numerical comparisons. The two equations used in the benchmark problems are the viscous Burgers’ equation and the porous medium equation, both in one dimension. Simulations are made for the two methods for: a) a travelling wave solution for the viscous Burgers’ equation, b) the Barenblatt selfsimilar analytical solution of the porous medium equation, and c) a waiting-time solution for the porous medium equation. Simulations are carried out for varying mesh sizes, and the numerical solutions are compared by computing errors in two ways. In the case of an analytic solution being available, the errors in the numerical solutions are computed directly from the analytic solution. In the case of no availability of an analytic solution, an approximation to the error is computed using a very fine mesh numerical solution as the reference solution.     

A Runge-Kutta/WENO method for solving equations for small-amplitude wave propagation in a saturated elastic porous medium

A high-accuracy Runge-Kutta/WENO method of up to fourth order with respect to time and fifth order with respect to space is developed for the numerical modeling of small-amplitude wave propagation in a steady fluid-saturated elastic porous medium. A system of governing equations is derived from a general thermodynamically consistent model of a compressible fluid flow through a saturated elastic porous medium, which is described by a hyperbolic system of conservation laws with allowance for finite deformations of the medium. The results of numerical solution of one- and two-dimensional wave fields demonstrate the efficiency of the method.     

Numerical analysis of wave propagation in fluid-filled deformable tubes

The theory of Biot describing wave propagation in fluid saturated porous media is a good effective approximation of a wave induced in a fluid-filled deformable tube. Nonetheless, it has been found that Biot's theory has shortcomings in predicting the fast P-wave velocities and the amount of intrinsic attenuation. These problems arises when complex mechanical interactions of the solid phase and the fluid phase in the micro-scale are not taken into account. In contrast, the approach proposed by Bernabe does take into account micro-scopic interaction between phases and therefore poses an interesting alternative to Biot's theory. A Wave propagating in a deformable tube saturated with a viscous fluid is a simplified model of a porous material, and therefore the study of this geometry is of great interest. By using this geometry, the results of analytical and numerical results have an easier interpretation and therefore can be compared straightforward. Using a Finite Difference viscoelastic wave propagation code, the transient response was simulated. The wave source was modified with different characteristic frequencies in order to gain information of the dispersion relation. It was found that the P-wave velocities of the simulations at sub-critical frequencies closely match those of Bernabe's solution, but at over-critical frequencies they come closer to Biot's solution. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Self-similar Solutions of Nonlinear Elasticity and Poroelasticity Equations

Recent advances in nonlinear wave propagation in elastic and porous elastic (poro-elastic) material have presented new nonlinear evolutionary equations. The derivation of these equations in three-dimensional space is based on the semilinear Biot theory. The nonlinear elastodynamic equations are derived form the more general model of poro-elastodynamic using consistency arguments. For simplicity, we discuss and carry out the analysis for the nonlinear elastic model. It is found in this article that the methods of symmetry groups and self-similar solutions can furnish solutions to the nonlinear elastodynamic wave equation. It is also found that these models lead to shock wave development in finite time. Necessary conditions for the existence of the solution are given and well-posedness of the Cauchy problem is discussed.     

On the propagation of elastic waves in two-phase media

At the present time a number of papers has been already devoted to the dynamics of two-phase media. One may mention the papers by Frenkel' [1], Rakhmatulin [2], Biot [3,4], Zwikker and Kosten [5], and others. However, the basic problem of the setting up of the equations of motion in two-phase media still cannot be considered solved and requires additional study and experimental verification.

This paper is concerned with the study of the simplest case of motion, which is the propagation of elastic waves in a homogeneous isotropic medium consisting of a solid and a fluid phase. The problems of the reflection of plane waves and surface waves at the free boundary of the half-space are solved. It is shown that the stress-strain relations established by Frenkel' are equivalent to the analogous relations proposed by Biot and that the equations of motion of the latter are more general.     

Reconstructing parameters of a medium from incomplete spectral information

The problem of the synthesis of a stratified medium with specified amplitude and phase properties is investigated. The wave propagation in the medium is described by a system of differential equations. The synthesis problem considered in the paper relates to inverse problems of spectral analysis with incomplete spectral information. Using the contour integral method we study properties of spectral characteristics and obtain algorithms for the solution of the synthesis problem for differential equations with singularities.     

Numerical modeling of 1D transient poroelastic waves in the low-frequency range

Propagation of transient mechanical waves in porous media is numerically investigated in 1D. The framework is the linear Biot model with frequency-independent coefficients. The coexistence of a propagating fast wave and a diffusive slow wave makes numerical modeling tricky. A method combining three numerical tools is proposed: a fourth-order ADER scheme with time-splitting to deal with the time-marching, a space-time mesh refinement to account for the small-scale evolution of the slow wave, and an interface method to enforce the jump conditions at interfaces. Comparisons with analytical solutions confirm the validity of this approach.     

Superconvergence results for Galerkin methods for wave propagation in various porous media

Finite-element methods are considered for numerically solving the equations describing wave propagation in various porous media such as inhomogeneous elastic media, fluid saturated media, composite isotropic inhomogeneous elastic media, composite anisotropic media, etc. Quasi-projection analyses based on an asymptotic expansion to high order of finite-element solutions are given to obtain error estimates in Sobolev spaces of nonpositive index for the approximate solution. Superconvergence phenomena for the finite-element methods under consideration are also investigated. © 1996 John Wiley & Sons, Inc.