A coupling of mixed and discontinuous Galerkin finite-element methods for poroelasticity     

Phillip Joseph Phillips  Mary F. Wheeler 《Computational Geosciences》2008,12(4):417-435

In this paper, we formulate a finite-element procedure for approximating the coupled fluid and mechanics in Biot’s consolidation model of poroelasticity. We approximate the flow variables by a mixed finite-element space and the displacement by a family of discontinuous Galerkin methods. Theoretical convergence error estimates are derived and, in particular, are shown to be independent of the constrained specific storage coefficient, c _o . This suggests that our proposed algorithm is a potentially effective way to combat locking, or the nonphysical pressure oscillations, which sometimes arise in numerical algorithms for poroelasticity.  相似文献   

西藏札达盆地控盆断裂有限元数值模拟   总被引：1，自引：0，他引：1  

段岩  孟宪刚  邵兆刚  朱大岗 《地质通报》2008,27(10)

札达盆地是喜马拉雅构造带中的一个山间断陷盆地,其演化过程与盆地两侧的控盆构造密切相关。对控盆断裂的构造应力场进行模拟计算,将有助于进一步深化对本区构造控盆的认识。因此,在对盆地构造地质进行详细调查的基础上,结合本区的深部地质与地球物理资料,对札达盆地控盆断裂的构造应力场进行了模拟。计算结果表明,札达盆地的演化明显受盆地两侧边界断裂的控制,札达盆地是在整体南北向挤压应力的作用下,不同块体差异隆升作用的结果。其南侧的控盆断裂为北倾的正断层,北侧的控盆断裂为南倾的逆断层,二者共同形成了南降北升的翘板式断陷盆地运动过程,是喜马拉雅地块在陆内汇聚挤压构造环境中构造应力场调整的一种方式。  相似文献   

基于二次场二维起伏地形MT有限元数值模拟   总被引：4，自引：1，他引：3  

赵广茂  李桐林  王大勇  李建平 《吉林大学学报(地球科学版)》2008,38(6):1055-1059

通过计算二次场来进行二维大地电磁数值模拟;导出了二维大地电磁二次场的微分方程,利用有限单元法来解微分方程;对矩形网格进行对角线的二次剖分,更容易且真实地模拟起伏地形。对几个典型模型进行了试算,与前人总场法的计算结果做了比较,两者视电阻率曲线一致,证明本文算法是正确的;通过2个简单的算例说明复杂地表下2种极化模式的MT观测资料都有明显的异常,视电阻率在TM模式下比TE模式更易受地形影响,TE模式下视电阻率曲线形态与地形呈"正相关",TM模式下反之。  相似文献   

弹性波的三维有限元模拟   总被引：5，自引：0，他引：5       下载免费PDF全文

邓玉琼  张之立 《地球物理学报》1990,33(1):44-53

本文讨论了弹性波瞬态传播问题的三维有限元计算方法及当前存在的实际困难.针对要求计算机内存大和计算时间长的问题,采取了改进措施.由于采用了集中质量矩阵和修正的中心差分时间积分显格式相结合的方法,可以使计算机内存和计算时间大为减少;由于采用结点定位法,最适合用于目前发展的并行计算机系统,可使计算速度大大增快;还采用了有效激发震源法,有效激发区是随时间步进的增加而逐步增大,这不仅能节省计算时间,而且使波场的传播过程一目了然,本文计算了由两种介质组成的三维楔形问题,得到若干典型剖面的瞬时波场图及随时间变化的合成地震图.  相似文献   

Semi-analytical solutions of a contaminant transport equation with nonlinear sorption in 1D     

Peter Frolkovič  Jozef Kačur 《Computational Geosciences》2006,10(3):279-290

A new method to determine semi-analytical solutions of one-dimensional contaminant transport problem with nonlinear sorption is described. This method is based on operator splitting approach where the convective transport is solved exactly and the diffusive transport by finite volume method. The exact solutions for all sorption isotherms of Freundlich and Langmuir type are presented for the case of piecewise constant initial profile and zero diffusion. Very precise numerical results for transport with small diffusion can be obtained even for larger time steps (e.g., when the Courant-Friedrichs-Lewy (CFL) condition failed).  相似文献   

Analysis of a rigid footing lying on three-layered soil using the finite difference method     

Priti Maheshwari  M. R. Madhav 《Geotechnical and Geological Engineering》2006,24(4):851-869

A numerical procedure is described for the analysis of the vertical deformation and the stress distribution of the strip footings on layered soil media. Three layers of soil with different stiffness are considered with the middle soil layer the thinnest and most stiff layer. The soil media is discretized and using the theory of elasticity, the governing differential equations are obtained in terms of vertical and horizontal displacements. These equations along with appropriate boundary and continuity conditions are solved by using the finite difference method. The vertical and horizontal displacements, strains and stresses are found at various nodes in the soil media. Parametric studies are carried out to study the effect of the placement depth of the middle soil layer, the relative ratios of the moduli of deformation of the soil layers on the vertical displacement of the footing and the vertical stress distribution. These studies reveal that the middle thin but very stiff layer acts like a plate and redistributes the stresses on the lower soft soil layer uniformly. The displacement on the top and bottom of the middle soil layer is almost the same showing that the compression of the middle layer is negligible as it is very stiff.  相似文献   

Performance of Mortar and Chemical Grout Injection into Surrounding Soil When Slurry Pipe-jacking Method is Used   总被引：1，自引：0，他引：1  

Hideki Shimada  Takashi Sasaoka  Saeid Khazaei  Yasuhiro Yoshida  Kikuo Matsui 《Geotechnical and Geological Engineering》2006,24(1):57-77

Small-diameter shallow tunnels are often being built by using the slurry pipe-jacking method. This system involves the pushing or thrusting of a drivage machine and concrete pipes into the ground. Chemical grout injection into the surrounding soil around the tunnel is carried out after the drivage and pushing processes are finished. The purpose of the chemical grout injection is to maintain permanent stability of the surrounding soil. However, the behavior of the chemical grouting material in the surrounding soil around the tunnel and the amount of optimum injection is not clearly understood. From these points of view, this paper discusses the performance of the chemical grouting material, when it is injected into the surrounding soil around the tunnel, by means of 2-D Eulerian–Lagrangian seepage analysis. Moreover, the effectiveness of the chemical grout injection was evaluated by using the non-linear finite element method. This investigation show when the range of the grouted zone is designed; it is necessary that the relationship between Young’s modulus of the soil/grouted zone and the confining stress be taken into consideration in order to establish effective, economical and safe chemical grout injection system. Understanding the performance of the seepage/dispersion behavior of the chemical grout and the characteristics of soil/ grouted zone is also important.  相似文献   

Time domain Gauss–Newton seismic waveform inversion in elastic media     

Dong-Hoon Sheen  Kagan Tuncay  Chang-Eob Baag  Peter J. Ortoleva 《Geophysical Journal International》2006,167(3):1373-1384

  相似文献   

Dual simulations of fluid flow and seismic wave propagation in a fractured network: effects of pore pressure on seismic signature   总被引：1，自引：0，他引：1  

S. Vlastos  E. Liu  I. G. Main  M. Schoenberg  C. Narteau  X. Y. Li  B. Maillot 《Geophysical Journal International》2006,166(2):825-838

  相似文献   

Comparison of methods to model the gravitational gradients from topographic data bases   总被引：7，自引：0，他引：7  

Christopher Jekeli  Lizhi Zhu 《Geophysical Journal International》2006,166(3):999-1014

A number of methods have been developed over the last few decades to model the gravitational gradients using digital elevation data. All methods are based on second-order derivatives of the Newtonian mass integral for the gravitational potential. Foremost are algorithms that divide the topographic masses into prisms or more general polyhedra and sum the corresponding gradient contributions. Other methods are designed for computational speed and make use of the fast Fourier transform (FFT), require a regular rectangular grid of data, and yield gradients on the entire grid, but only at constant altitude. We add to these the ordinary numerical integration (in horizontal coordinates) of the gradient integrals. In total we compare two prism, two FFT and two ordinary numerical integration methods using 1" elevation data in two topographic regimes (rough and moderate terrain). Prism methods depend on the type of finite elements that are generated with the elevation data; in particular, alternative triangulations can yield significant differences in the gradients (up to tens of Eötvös). The FFT methods depend on a series development of the topographic heights, requiring terms up to 14th order in rough terrain; and, one popular method has significant bias errors (e.g. 13 Eötvös in the vertical–vertical gradient) embedded in its practical realization. The straightforward numerical integrations, whether on a rectangular or triangulated grid, yield sub-Eötvös differences in the gradients when compared to the other methods (except near the edges of the integration area) and they are as efficient computationally as the finite element methods.  相似文献