Effective block diagonal preconditioners for Biot's consolidation equations in piled‐raft foundations

The finite element (FE) simulation of large‐scale soil–structure interaction problems (e.g. piled‐raft, tunnelling, and excavation) typically involves structural and geomaterials with significant differences in stiffness and permeability. The symmetric quasi‐minimal residual solver coupled with recently developed generalized Jacobi, modified symmetric successive over‐relaxation (MSSOR), or standard incomplete LU factorization (ILU) preconditioners can be ineffective for this class of problems. Inexact block diagonal preconditioners that are inexpensive approximations of the theoretical form are systematically evaluated for mitigating the coupled adverse effects because of such heterogeneous material properties (stiffness and permeability) and because of the percentage of the structural component in the system in piled‐raft foundations. Such mitigation led the proposed preconditioners to offer a significant saving in runtime (up to more than 10 times faster) in comparison with generalized Jacobi, modified symmetric successive over‐relaxation, and ILU preconditioners in simulating piled‐raft foundations. Copyright © 2012 John Wiley & Sons, Ltd.     

A continuous/discontinuous Galerkin framework for modeling coupled subsurface and surface water flow

We consider conjunctive surface-subsurface flow modeling, where surface water flow is described by the shallow water equations and ground water flow by Richards’ equation for the vadose zone. Coupling between the models is based on the continuity of flux and water pressure. Numerical approximation of the coupled model using the framework of discontinuous Galerkin (DG) methods is formulated. In the subsurface, the local discontinuous Galerkin (LDG) method is used to approximate ground water velocity and hydraulic head; a DG method is also used to approximate surface water velocity and elevation. This approach allows for a weak coupling of the models and the use of different approximating spaces and/or meshes within each regime. A simplified LDG method based on continuous approximations to water head is also described. Numerical results that investigate physical and numerical aspects of surface–subsurface flow modeling are presented. This work was supported by National Science Foundation grant DMS-0411413.     

THE SETTLEMENT OF A RIGID CIRCULAR FOUNDATION RESTING ON A HALF-SPACE EXHIBITING A NEAR SURFACE ELASTIC NON-HOMOGENEITY

The present paper examines the elastostatic problem pertaining to the axisymmetric loading of a rigid circular foundation resting on the surface of a non-homogeneous elastic half-space. The non-homogeneity corresponds to a depth variation in the linear elastic shear modulus according to the exponential form G(z)=G₁+G₂e^-ζz. The equations of elasticity governing this type of non-homogeneity are solved by employing a Hankel transform technique. The mixed boundary value problem associated with the indentation of the half-space by the rigid circular foundation is reduced to a Fredholm integral equation which is solved via a numerical technique. The numerical results presented in the paper illustrate the influence of the near-surface elastic non-homogeneity on the settlement of the foundation.     

Tensionless–frictionless interaction of flexible annular foundation with a transversely isotropic multi‐layered half‐space

A transversely isotropic multi‐layered half‐space, with axis of material symmetry perpendicular to the free surface, supports a flexible either annular or solid circle foundation. The contact area of the foundation and the half‐space is considered to be both frictionless and tensionless. The foundation is assumed to be affected by a vertical static axisymmetric load. Detailed analysis of the interaction of these two systems with different thickness of layers is the target of this paper. With the use of ring load Green's functions for both the foundation and the continuum half‐space, an integral equation accompanied with some inequalities is introduced to model the complex BVP. With the incorporation of ring‐shape FEM, we are capable of capturing both regular and singular solution smoothly. The validity of the combination of the analytical and numerical method is proved with comparing the results of this paper with a number of benchmark cases of both linear and nonlinear interaction of circular and annular foundation with half‐space. Some new illustrations are presented to portray the aspect of the anisotropy and layering of the half‐space. Copyright © 2014 John Wiley & Sons, Ltd.     

Na-K-Mg三角图修正与Na-K温标选取

Na-K-Mg三角图作为水-岩平衡状态初步判定的方法,未考虑具体温标公式和矿物成分对平衡状态的影响,因此通过基于不同Na-K温标公式修正后的Na-K-Mg三角图判断水岩平衡状态,并利用适合的Na-K温标公式计算地下热储温度,对划分地热系统成因类型和开发利用地热资源具有重要意义。研究发现,修正后的Na-K-Mg三角图中任意一点,对某种温标公式确定的曲线处于平衡状态,对另一温标公式确定的曲线则有可能处于部分平衡状态;且温度越低,这种现象越显著。最贫Mg上边界和最富Mg下边界是由不同的温标公式确定的,随温度具有明显的分段性。利用各Na-K温标公式作出T-lg（Na/K）曲线图,绘制曲线交点对应的基于各\     

等效水力隙宽和水力梯度对岩体裂隙网络非线性渗流特性的影响

等效水力隙宽和水力梯度是影响岩体裂隙网络渗流特性的重要因素。制作裂隙网络试验模型,建立高精度渗流试验系统;求解纳维-斯托克斯方程,模拟流体在裂隙网络内的流动状态,研究等效水力隙宽和水力梯度对非线性渗流特性的影响。结果表明,当水力梯度较小时,等效渗透系数保持恒定的常数,流体流动属于达西流动区域,流量与压力具有线性关系,可采用立方定律计算流体流动;当水力梯度较大时,等效渗透系数随着水力梯度的增加而急剧减少,流体流动进入强惯性效应流动区域,流量与压力具有强烈的非线性关系,可采用Forchheimer方程计算流体流动。随着等效水力隙宽的增加,区别线性和非线性流动区域的临界水力梯度呈幂函数关系递减。当水力梯度小于临界水力梯度时,控制方程可选立方定律;当水力梯度大于临界水力梯度时,控制方程可选Forchheimer方程,其参数A和B可根据经验公式计算得到。其研究结果可为临界水力梯度的确定及流体流动控制方程的选取提供借鉴意义。     

Application and verification of fuzzy algebraic operators to landslide susceptibility mapping

The aim of this study was to apply and to verify the use of fuzzy logic to landslide susceptibility mapping in the Gangneung area, Korea, using a geographic information system (GIS). For this aim, in the study, a data-derived model (frequency ratio) and a knowledge-derived model (fuzzy operator) were combined. Landslide locations were identified by changing the detection technique of KOMPSAT-1 images and checked by field studies. For landslide susceptibility mapping, maps of the topography, lineaments, soil, forest, and land cover were extracted from the spatial data sets, and the eight factors influencing landslide occurrence were obtained from the database. Using the factors and the identified landslide, the fuzzy membership values were calculated. Then fuzzy algebraic operators were applied to the fuzzy membership values for landslide susceptibility mapping. Finally, the produced map was verified by comparing with existing landslide locations for calculating prediction accuracy. Among the fuzzy operators, in the case in which the gamma operator (λ = 0.975) showed the best accuracy (84.68%) while the case in which the fuzzy or operator was applied showed the worst accuracy (66.50%).     

二维圆柱坐标下FDTD法对多频电磁波测井的数值模拟

在介质围绕Z轴旋转对称的情况下,三维直角坐标下的麦克斯韦方程可变成两组分别对应TM波和TE波的二维柱坐标下的偏微分方程.采用线圈作为电磁波测井的发射器和接收器的时候,可采用TE波对应的方程,这时相当于磁极子的情况.而采用电偶极子作为发射和接收器的时候,则对应TM波的方程.采用时间和空间均为二阶精确的有限差分方法,将偏微分方程进行差分化.这样,空间的电磁场可由时间域有限差分法(FDTD)来求解.     

Modeling of electromagnetic waves in axisymmetric geometry using a Laguerre spectral method

A spectral method for modeling high-frequency electromagnetic waves in axisymmetric geometry is proposed.The method is based on the expansion of the solutions of Maxwell’s equations in Laguerre functions in the time region.The spectral method is used to solve Maxwell’s equations for both 2D media and stratified media. In the case of stratified media, a Fourier–Bessel expansion in the radial variable is used. The effectiveness of the spectral and finite-difference methods is compared. Harmonic solutions and solitary solutions by the Laguerre method are considered, and the dynamics of monochromatic and broadband electromagnetic pulses are examined.     

关于频率电磁测深几个问题的探讨(八)——频率电磁测深与瞬变电磁测深的关系

从探测原理出发,说明人工源频率电磁测深与瞬变电磁测深的关系密切,这是一个问题的两个方面,可以互相转换。两种探测方法实质是一样的,具有共性;但毕竟从频率域和时间域两个侧面反映问题,又具有各自特性。两种方法的原理所依据的麦克斯韦方程,其正演公式可通过傅立叶变换相互转换;从实际探测所发收的信号特征,场区划分与探测深度等主要方面阐明两种探测方法异同,以加深对探测理论与方法的认识,结合实际择优应用。