Surface waves in a general anisotropic poroelastic solid half-space   总被引：1，自引：0，他引：1  

M. D. Sharma 《Geophysical Journal International》2004,159(2):703-710

  相似文献   

Scattering of plane SV waves by cylindrical canyons in saturated porous medium   总被引：4，自引：0，他引：4  

Li Weihua  Zhao Chenggang 《Soil Dynamics and Earthquake Engineering》2005,25(12):981-995

On the basis of Biot dynamic theory, an analytic solution of two-dimensional scattering and diffraction of plane SV waves by circular cylindrical canyons in a half space of saturated porous media is presented in this paper for the first time. The solution is obtained by employing the Fourier–Bessel series expansion technique. Parametric studies had been carried out, which includes: the angle of incidence, the frequency of the incident SV wave, the porosity of saturated porous medium and the stiffness and Poisson's ratio of the solid-skeleton. All the outcomes are useful for the seismic analysis of the surface topography conditions.  相似文献   

Numerical modeling of toxic nonaqueous phase liquid removal from contaminated groundwater systems: mesh effect and discretization error estimation          下载免费PDF全文

Chongbin Zhao  Thomas Poulet  Klaus Regenauer‐Lieb 《国际地质力学数值与分析法杂志》2015,39(6):571-593

Numerical modeling has now become an indispensable tool for investigating the fundamental mechanisms of toxic nonaqueous phase liquid (NAPL) removal from contaminated groundwater systems. Because the domain of a contaminated groundwater system may involve irregular shapes in geometry, it is necessary to use general quadrilateral elements, in which two neighbor sides are no longer perpendicular to each other. This can cause numerical errors on the computational simulation results due to mesh discretization effect. After the dimensionless governing equations of NAPL dissolution problems are briefly described, the propagation theory of the mesh discretization error associated with a NAPL dissolution system is first presented for a rectangular domain and then extended to a trapezoidal domain. This leads to the establishment of the finger‐amplitude growing theory that is associated with both the corner effect that takes place just at the entrance of the flow in a trapezoidal domain and the mesh discretization effect that occurs in the whole NAPL dissolution system of the trapezoidal domain. This theory can be used to make the approximate error estimation of the corresponding computational simulation results. The related theoretical analysis and numerical results have demonstrated the following: (1) both the corner effect and the mesh discretization effect can be quantitatively viewed as a kind of small perturbation, which can grow in unstable NAPL dissolution systems, so that they can have some considerable effects on the computational results of such systems; (2) the proposed finger‐amplitude growing theory associated with the corner effect at the entrance of a trapezoidal domain is useful for correctly explaining why the finger at either the top or bottom boundary grows much faster than that within the interior of the trapezoidal domain; (3) the proposed finger‐amplitude growing theory associated with the mesh discretization error in the NAPL dissolution system of a trapezoidal domain can be used for quantitatively assessing the correctness of computational simulations of NAPL dissolution front instability problems in trapezoidal domains, so that we can ensure that the computational simulation results are controlled by the physics of the NAPL dissolution system, rather than by the numerical artifacts. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

A model for moisture and heat flow in fractured unsaturated porous media          下载免费PDF全文

A.M. Varnosfaderani  B. Gatmiri  E.M. Haghighi 《国际地质力学数值与分析法杂志》2017,41(6):828-858

The present study investigates propagation of a cohesive crack in non‐isothermal unsaturated porous medium under mode I conditions. Basic points of skeleton deformation, moisture, and heat transfer for unsaturated porous medium are presented. Boundary conditions on the crack surface that consist of mechanical interaction of the crack and the porous medium, water, and heat flows through the crack are taken into consideration. For spatial discretization, the extended finite element method is used. This method uses enriched shape functions in addition to ordinary shape functions for approximation of displacement, pressure, and temperature fields. The Heaviside step function and the distance function are exploited as enrichment functions for representing the crack surfaces displacement and the discontinuous vertical gradients of the pressure and temperature fields along the crack, respectively. For temporal discretization, backward finite difference scheme is applied. Problems solved from the literature show the validity of the model as well as the dependency of structural response on the material properties and loading. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

The effective stress concept and the evaluation of changes in pore air pressure under jacketed isotropic compression tests for dry sand based on the 2‐phase porous theory          下载免费PDF全文

Sayuri Kimoto  Fusao Oka  Yasuhiro Morimoto 《国际地质力学数值与分析法杂志》2017,41(18):1894-1907

The effective stress concept for solid‐fluid 2‐phase media was revisited in this work. In particular, the effects of the compressibility of both the pore fluid and the soil particles were studied under 3 different conditions, i.e., undrained, drained, and unjacketed conditions based on a Biot‐type theory for 2‐phase porous media. It was confirmed that Terzaghi effective stress holds at the moment when soil grains are assumed to be incompressible and when the compressibility of the pore fluid is small enough compared to that of the soil skeleton. Then, isotropic compression tests for dry sand under undrained conditions were conducted within the triaxial apparatus in which the changes in the pore air pressure could be measured. The ratio of the increment in the cell pressure to the increment in the pore air pressure, m, corresponds to the inverse of the B value by Bishop and was obtained during the step loading of the cell pressure. In addition, the m values were evaluated by comparing them with theoretically obtained values based on the solid‐fluid 2‐phase mixture theory. The experimental m values were close to the theoretical values, as they were in the range of approximately 40 to 185, depending on the cell pressure. Finally, it was found that the soil material with a highly compressible pore fluid, such as air, must be analyzed with the multi‐phase porous mixture theory. However, Terzaghi effective stress is practically applicable when the compressibilities of both the soil particles and the pore fluid are small enough compared to that of the soil skeleton.  相似文献   

Why asymptotic limit of the acid dissolution capacity can lead to a sharp dissolution front in chemical dissolution of porous rocks?          下载免费PDF全文

Chongbin Zhao  B.E. Hobbs  A. Ord 《国际地质力学数值与分析法杂志》2017,41(15):1590-1602

The use of the asymptotic limit can greatly simplify the theoretical analysis of chemical dissolution front instabilities in fluid‐saturated rocks and therefore make it possible to obtain mathematical solutions, which often play a crucial role in understanding the propagation behavior of chemical dissolution fronts in chemical dissolution systems. However, there has been a debate in recent years that the asymptotic limit of the acid dissolution capacity (i.e., the acid dissolution capacity number approaching zero) alone cannot lead to a sharp dissolution front of the Stefan type in the acidization dissolution system, in which the dissolvable minerals of carbonate rocks are chemically dissolved by the injected acid flow. The acid dissolution capacity number is commonly defined as the ratio of the volume of the carbonate rock dissolved by an acid to that of the acid. In this paper, we use four different proof methods, including (i) direct use of the fundamental concepts; (ii) use of the mathematical governing equations of an acidization dissolution system; (iii) use of the different time scaling approach; and (iv) use of a moving coordinate system approach, to demonstrate that the asymptotic limit of the acid dissolution capacity can indeed lead to sharp dissolution fronts of the Stefan type in acidization dissolution systems on a much larger time scale (than the dissolution time scale). Our new finding is that on the reaction time scale, the condition of the conventional time derivative of porosity approaching zero alone can ensure that the acidization dissolution front has a sharp shape of the Stefan type. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献   

Analytical solutions for advective–dispersive solute transport in double‐layered finite porous media     

Peter John Cleall 《国际地质力学数值与分析法杂志》2011,35(4):438-460

Analytical solutions for advection and dispersion of a conservative solute in a one‐dimensional double‐layered finite porous media are presented. The solutions are applicable to five scenarios that have various combinations of fixed concentration, fixed flux and zero concentration gradient conditions at the inlet and outlet boundaries that provide a wide number of options. Arbitrary initial solute concentration distributions throughout the media can be considered via explicit formulations or numerical integration. The analytical solutions presented have been verified against numerical solutions from a finite‐element‐based approach and an existing closed‐form solution for double‐layered media with an excellent correlation being found in both cases. A practical application pertaining to advective transport induced by consolidation of underlying sediment layers on contaminant movement within a capped contaminated sediment system is presented. Comparison of the calculated concentrations and fluxes with alternative approaches clearly illustrates the need to consider advection processes. Consideration of the different features of contaminant transport due to varying pore‐water velocity fields in primary consolidation and secondary consolidation stages is achieved via the use of non‐uniform initial concentration distributions within the proposed analytical solutions. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

Interaction of cnoidal waves with an array of vertical concentric porous cylinders     

《Applied Ocean Research》2016

An array of large concentric porous cylinder arrays is mounted in shallow water exposed to cnoidal waves. The interactions between waves and cylinders are studied theoretically using an eigenfunction expansion approach. Semi-analytical solutions of hydrodynamic loads and wave run-up on each cylinder are obtained using first approximation to cnoidal waves. The square array configuration of four-legged identical concentric porous cylinder is investigated in present study. Numerical results reveal the variation of dimensionless wave force and wave run-up on individual cylinder with angle of incidence, porosity parameter, spacing between outer and inner cylinders, spacing between concentric porous cylinders and wave parameter. Different mechanism of wave force is found under different range of scattering parameter.  相似文献   

九寨沟“8.8”地震背景下火花海钙华坝体无损探测分析     

罗尧东  代群威  董发勤  王富东  赵学钦  廖传康  崔杰  李琼芳  王安平 《中国岩溶》2021,40(1):147-156

为研究震后九寨沟钙华景区的溶洞、裂隙、暗流等地质特性,丰富钙华景区岩溶洞穴无损探测技术,探讨了多物探技术在多孔疏松钙华地区作业的可靠性。首次运用地质雷达法及高密度电法相结合的无损探测技术,在九寨沟景区火花海特定区域进行了现场勘查和无损探测试验。结合对比两种无损探测成果与测区内出露地质剖面3组勘探数据,结果发现:高密度电法与地质雷达法相结合,两者能优劣互补,相互佐证,在解决钙华区域浅部溶洞、暗流、裂隙等工作中,对探测目标体具有极高的识别度。两种方法探测数据与震后出露剖面地质特征保持极高的一致性,验证了两种无损检测方法在多孔疏松介质条件下数据的真实可靠性,解决了多物探技术在多孔疏松介质地区施工的疑点和难点,丰富了工程物探领域技术。   相似文献   

基于粒径分维值检验的松散介质渗透系数空间估值——以河北平原为例     

罗晓娟  徐恒力  叶含春 《地质科技情报》2007,26(1):87-90

为了克服松散含水介质非均质性各向异性对介质渗透系数的求取所带来的困难,以水动力条件为线索,利用互相关法寻找同一时代地层及同一时代地表水的主、次流方向.在主流方向上,建立了石家庄-束鹿-衡水段介质颗粒粒径随距离递变的数学模型,并以此为基础推求无钻孔地段的粒径.同时,运用分形理论,计算了中、上更新统地层主流方向平均粒径空间分布的分维值,并将其作为检验粒径估值正确与否的判据.通过不同岩性的渗透系数-深度模型,实现了用平均粒径预测渗透系数的全过程.  相似文献