Attenuation, transport and diffusion of scalar waves in textured random media

Most theoretical investigations of seismic wave scattering rely on the assumption that the underlying medium is statistically isotropic. However, deep seismic soundings of the crust as well as geological observations often reveal the existence of elongated or preferentially oriented scattering structures. In this paper, we develop mean field and radiative transfer theories to describe the attenuation and multiple scattering of a scalar wavefield in an anisotropic random medium. The scattering mean free path is found to depend strongly on the propagation direction. We derive a radiative transfer equation for statistically anisotropic random media from the Bethe–Salpeter formalism and propose a Monte Carlo method to solve this equation numerically. At longer times, the energy density is shown to obey a tensorial diffusion equation. The components of the diffusion tensor are obtained in closed form and excellent agreement is found between Monte Carlo simulations and analytical solutions of the diffusion equation. The theory has important potential implications for lithospheric models where scatterers are for example flat structures preferentially aligned along the surface. In this simple geometry, analytical expressions of the Coda Q parameter can be given explicitly in the diffusive regime. Our results suggest that pulse broadening and coda decay are controlled by different parameters, related to the eigenvalues of the diffusion tensor. These eigenvalues can differ by more than one order of magnitude. This theory could be applied to probe the anisotropy of length scales in the lithosphere.     

Numerical modelling of the dynamic process of oil displacement by water in sandstone reservoirs with random pore structures

In order to maintain the production rate of a reservoir and improve the displacement efficiency, it is crucial to have an in-depth understanding of the process of oil displacement by water. However, with respect to the conceptualization of porous media of a reservoir, very limited efforts have been made to the pore structures inside the reservoirs. In this paper, the pore structures of a sandstone reservoir were generated by using the method of random growth algorithm. Based on the randomly generated model, a theoretical model to describe the dynamic process of oil displacement by water in the sandstone reservoir was established, and then corresponding numerical modelling was performed. The effects of the displacement velocity, the viscosity ratio of oil-water phase and the porosity of reservoirs on the displacement performance were also analyzed. Results show that due to a great difference in the viscosity between oil and water phases, the moving interface of water phase is not uniform, and the viscous fingering occurs, tending to proceed along the direction with the minimum flow resistance. There is not a linear relationship between the displacement velocity and the displacement efficiency. Too high displacement velocities do not lead to much better displacement efficiency, while a higher pressure drop is caused.Choosing a proper displacement velocity is indispensable in practical engineering. A lower oil-water viscosity ratio is more favorable to obtain high displacement efficiency. Under the present simulation conditions, when the viscosity ratio is 1.2, the displacement efficiency reaches 96.2% at a moderate Reynolds number. The porosity is not a sole factor determining the displacement performance. Even for the same porosity, the shape and length of preferential flow paths are different and randomly distributed,causing a different displacement performance. A large tortuosity tends to result in a low hydraulic conductivity and displacement efficiency.     

Propagation of electromagnetic waves from an induction source in media with time-dependent conductivity and permittivity

Electromagnetic wave propagation from an induction source (a transmitter coil) in media with time-dependent conductivity and permittivity is simulated using the vector finite-element method. A special algorithm is suggested for source approximation.     

Contaminant transport with non‐equilibrium processes in unsaturated soils and implicit characteristic Galerkin scheme

A non‐equilibrium sorption—advection—diffusion model to simulate miscible pollutant transport in saturated–unsaturated soils is presented. The governing phenomena modelled in the present simulation are: convection, molecular diffusion, mechanical dispersion, sorption, immobile water effect and degradation, including both physical and chemical non‐equilibrium processes. A finite element procedure, based on the characteristic Galerkin method with an implicit algorithm is developed to numerically solve the model equations. The implicit algorithm is formulated by means of a combination of both the precise and the traditional numerical integration procedures. The stability analysis of the algorithm shows that the unconditional stability of the present implicit algorithm is enhanced as compared with that of the traditional implicit numerical integration procedure. The numerical results illustrate good performance of the present algorithm in stability and accuracy, and in simulating the effects of all the mentioned phenomena governing the contaminant transport and the concentration distribution. Copyright © 2000 John Wiley & Sons, Ltd.     

The intrinsic permeability of microcracks in porous solids: Analytical models and Lattice Boltzmann simulations

The paper presents a synthesis of analytical modeling and computational simulations of the intrinsic permeability of microcracks, embedded in porous materials taking into account the interaction of the fluid flow in the microcrack with the surrounding porous material. In the first part of the paper, using the DARCY , STOKES , BRINKMAN , and the BEAVERS–JOSEPH approximations, we derive the intrinsic permeability of a plain non‐rough microcrack in terms of the microcrack geometry and the permeability of the porous material surrounding the microcrack. In the second part of the paper, the intrinsic permeability of a microcrack is determined by means of computational simulations using the framework of the lattice Boltzmann method with partial bounceback conditions. The comparison of predictions from the analytical model and the numerical simulations show an excellent agreement. Copyright © 2017 John Wiley & Sons, Ltd.     

Solution of quasi‐static large‐strain problems by the material point method

A procedure for solving quasi‐static large‐strain problems by the material point method is presented. Owing to the Lagrangian–Eulerian features of the method, problems associated with excessive mesh distortions that develop in the Lagrangian formulations of the finite element method are avoided. Three‐dimensional problems are solved utilizing 15‐noded prismatic and 10‐noded tetrahedral elements with quadratic interpolation functions as well as an implicit integration scheme. An algorithm for exploiting the numerical integration procedure on the computational mesh is proposed. Several numerical examples are shown. Copyright © 2010 John Wiley & Sons, Ltd.     

列车荷载下冻土路基非平稳随机振动分析的时域显式法北大核心CSCD

随着我国冻土地区铁路运营里程的不断提升,冻土路基在随机列车荷载作用下的动力响应分析成为了急需解决的工程问题。本文以青藏铁路某路基横断面为例,采用时域显式蒙特卡罗模拟法,计算其在随机列车荷载作用下响应的统计特征。首先,提出列车驶过施加在路基顶面荷载的快速计算方法,并引入轨道不平顺及列车行驶速度两类随机参数,以生成随机分析所需的大量荷载样本。然后,采用数值积分方法将运动方程在时域上进行离散,建立任意离散时刻冻土路基动力响应关于列车荷载的显式表达式。基于该显式表达式,可以高效地实施蒙特卡罗模拟,得到在随机列车荷载作用下路基关键响应的均值、标准差和峰值等统计量。采用该方法,分析了路基在夏冬两个季节不同深度处的随机动力响应,发现在夏季路基浅层中的位移和速度响应最为剧烈。数值算例表明,时域显式蒙特卡罗模拟法在冻土路基随机振动分析中具有理想的计算精度和计算效率。     

基于薄层法和移动坐标系法的地铁竖向振动在成层土层中传播

针对地铁竖向振动在成层地基中的传播,提出了移动荷载作用下层状地基的分析模型。基于该模型,利用狄拉克函数及三重傅里叶变换将时空域内单个移动简谐荷载转换为频率-波数域内的荷载,结合薄层法和移动坐标系法推导了单个移动简谐荷载和移动简谐线荷载作用下三维层状地基动力响应的解析解,并给出了移动简谐线荷载动力响应解析解中参数n的经验取值范围;分析了荷载的移动速度对层状地基动力响应的影响以及弹性模量、泊松比、阻尼比和荷载频率对土体临界速度的影响。结果表明：荷载的移动速度对不同频带的动力响应的影响范围不同;移动速度对低频响应的影响程度大于对高频响应的影响程度;相比于泊松比和阻尼比,弹性模量对土体临界速度的影响最大;频率越接近荷载振动频率的振动响应,其幅值越大,临界速度越小。     

Vertical harmonic vibration of piled raft foundations in layered soils

This paper develops a method to analyze the piled raft foundation under vertical harmonic load. This method takes into account the interactions among the piles, soil, and raft. The responses of the piles and raft are formulated as a series of equations in a suitable way and that of layered soils is simulated with the use of the analytical layer‐element method. Then, according to the equilibrium and continuity conditions at the piles–soil–raft interface, solutions for the piled raft systems are obtained and further demonstrated to be correct through comparing with the existing results. Finally, some examples are given to investigate the influence of the raft, pile length‐diameter ratio, and layering on the response of the piled raft foundations. Copyright © 2017 John Wiley & Sons, Ltd.     

Effect of property variation and modelling on convection in a fluid overlying a porous layer

An accurate computational analysis is presented for the onset of thermal convection in a two‐layer system which is comprised of a saturated layer of porous material described by Darcy's law, over which lies a layer of the same saturating fluid. The two‐layer system is heated from below and the upper (fluid) surface is allowed to be fixed or stress free. The onset of convection may have a bi‐modal nature in which convection may be dominated by the porous medium or by the fluid depending on the depths of the relative layers, but this is strongly controlled by material parameters. The effect of variation of relevant fluid and porous material properties is investigated in detail, as is the effect of the interface boundary condition between the fluid and the porous medium. Copyright © 2001 John Wiley & Sons, Ltd.