POTENTIAL FIELD OF A STATIONARY ELECTRIC CURRENT USING FREDHOLM'S INTEGRAL EQUATIONS OF THE SECOND KIND*

An integral equation method is described for solving the potential problem of a stationary electric current in a medium that is linear, isotropic and piecewise homogeneous in terms of electrical conductivity. The integral equations are Fredholm's equations of the ‘second kind’ developed for the potential of the electric field. In this method the discontinuity-surfaces of electrical conductivity are divided into ‘sub-areas’ that are so small that the value of their potential can be regarded as constant. The equations are applied to 3-D galvanic modeling. In the numerical examples the convergence is examined. The results are also compared with solutions derived with other integral equations. Examples are given of anomalies of apparent resistivity and mise-a-la-masse methods, assuming finite conductivity contrast. We show that the numerical solutions converge more rapidly than compared to solutions published earlier for the electric field. This results from the fact that the potential (as a function of the location coordinate) behaves more regularly than the electric field. The equations are applicable to all cases where conductivity contrast is finite.     

Ground motion selection for simulation‐based seismic hazard and structural reliability assessment

This paper examines four methods by which ground motions can be selected for dynamic seismic response analyses of engineered systems when the underlying seismic hazard is quantified via ground motion simulation rather than empirical ground motion prediction equations. Even with simulation‐based seismic hazard, a ground motion selection process is still required in order to extract a small number of time series from the much larger set developed as part of the hazard calculation. Four specific methods are presented for ground motion selection from simulation‐based seismic hazard analyses, and pros and cons of each are discussed via a simple and reproducible illustrative example. One of the four methods (method 1 ‘direct analysis’) provides a ‘benchmark’ result (i.e., using all simulated ground motions), enabling the consistency of the other three more efficient selection methods to be addressed. Method 2 (‘stratified sampling’) is a relatively simple way to achieve a significant reduction in the number of ground motions required through selecting subsets of ground motions binned based on an intensity measure, IM. Method 3 (‘simple multiple stripes’) has the benefit of being consistent with conventional seismic assessment practice using as‐recorded ground motions, but both methods 2 and 3 are strongly dependent on the efficiency of the conditioning IM to predict the seismic responses of interest. Method 4 (‘generalized conditional intensity measure‐based selection’) is consistent with ‘advanced’ selection methods used for as‐recorded ground motions and selects subsets of ground motions based on multiple IMs, thus overcoming this limitation in methods 2 and 3. Copyright © 2015 John Wiley & Sons, Ltd.     

ON REVERSE-TIME MIGRATION*

Migration is a process whereby events in ‘image space’ are mapped into their correct positions in ‘object space’. The wave equations associated with this mapping may be defined and solved numerically either in image space or in object space. In the former the CMP section, which represents the initial conditions, is extrapolated toward increasing depths, and the migrated data are recovered at zero time. In the latter, the wave-field extrapolation takes place in the coordinate frame of the depth section, and the CMP data serve as boundary conditions at the surface. Computations begin at the last sample of the record section and continue ‘reverse time’ until time zero. This paper describes a reverse-time migration (RTM) method and compares its performance with that of an image-space method based on the idea of phase shift plus interpolation (PSPI). Synthetic zero-offset sections serve as examples for migration experiments with the RTM and PSPI methods. It is shown that the RTM approach to migration is rather expensive, but its robustness and accuracy are difficult to surpass.     

三维TTI介质高效射线追踪方法

TTI介质是石油地震勘探领域最常用的各向异性介质,快速计算TTI介质射线路径和走时信息有重要的研究意义.TTI介质传统运动学射线追踪方法一般基于任意弹性介质射线方程,利用Bond变换或者四阶张量变换来处理复杂的21个弹性参数,因而非常耗时.实际野外对称轴统一的TTI介质模型,一般可以看成VTI介质模型旋转一定角度获得.为此,本文推导了三维VTI介质射线追踪方程,提出先在本构坐标系中进行VTI介质射线追踪,再通过坐标旋转将射线路径旋转至观测坐标系中,获得TTI介质射线路径.数值模型计算表明该方法高效和精确,较传统方法效率提高了近4倍.在强各向异性等特殊情况下,体波波前面都与理论群速度面一致.     

Ground motions around a deep semielliptic canyon with a horizontal edge subjected to incident plane SH waves

We study scattering of antiplane shear waves induced by a deep semielliptic canyon with a horizontal edge. We employ the region-point-matching technique to cope with the problem considered. Through an auxiliary boundary, a part of the circumference of a semiellipse, the whole analyzed region is divided into two subregions. We express the displacement fields in terms of Mathieu functions. We unify two distinct elliptic coordinates via a simple coordinate transformation relation. Integration of the coordinate transformation relation into the region-point-matching technique simplifies the procedure for constructing simultaneous equations. Imposing the continuity conditions and traction-free ones, we obtain the expansion coefficients. Frequency-domain results demonstrate ground motion variability based on several key factors. Ground surface responses under seismic shaking are also simulated in the time domain.     

结构动力学方程的显式积分格式

本文从空间解耦有限元常微分方程组出发，探讨了结构动力学方程的高精度显式积分格式。通过被积函数的拉格朗日多项式内插和分部积分导出了波动数值模拟的一组显式时步积分公式。这组公式是时间和空间解耦的，即波场内任一离散节点在任一时刻的波动数据可以用这组公式依据该节点及其邻近节点在该时刻之前的n＋1个时刻的波动数据显式地算出（n为非负整数），阐明了这组公式的如下特点：第一，其截断误差的量级不超过0（Δt^n＋3），Δt为时间步距。第二，它不仅可用于线性波动的数值模拟，而且可用于本构方程具有强非线性情形。第三，这组公式也可推广应用于一系列数学物理暂态问题的数值求解。针对一个简单的时不变系统初步分析了此组积分格式的稳定性。但是，对其稳定性尚需作进一步研究。     

Exact solution for two-dimensional flow to a well in an anisotropic domain

Although most current applications of the analytic element method are formulated for isotropic hydraulic conductivity, anisotropic domains can be modeled with analytic elements using the well-known coordinate transformation where one coordinate axis is scaled by the square root of the anisotropy ratio. If the standard analytic solution for steady radial flow to a well is used with this coordinate transformation, the resulting solution correctly models the far field but it does not meet the constant head boundary condition at the well radius. This could be a significant shortcoming if you are interested in the flow field close to the well or want to estimate the head at the pumping well. A new solution for two-dimensional steady flow to a well in an anisotropic domain is presented. This solution satisfies the governing equations exactly and meets the constant head boundary condition at the well radius exactly. It was derived using a conformal mapping.     

Seismic response of pile groups under oblique-shear and Rayleigh waves

A simple analytical solution is presented to calculate the pile-soil-pile interaction and eventually the dynamic response of pile groups when excited by the passage of Rayleigh waves and obliquely incident SH-waves. A dynamic Winkler model, with realistic frequency-dependent ‘springs’ and ‘dashpots’ in conjunction with physically motivated approximations is utilized to compute the wave field radiating from an oscillating pile and the effect of this field on an adjacent pile. The coupled rocking motion of the pile group resulting from Rayleigh waves and the torsional motion of the pile group resulting from SH-waves is accurately predicted by a simple mathematical expression. The results of the presented method can be obtained with ‘hand calculations’ and are in excellent agreement with results from ‘rigorous’ solutions based on integral equation formulations. It is found that the group response is primarily affected from the phase difference of the input seismic motion at the location of each pile (wave-passage effect). Pile-soil-pile interaction has insignificant effect and can be neglected.     

A perspective on numerical analysis of the diffusion equation

The equation of transient groundwater motion is founded on the principle of mass conservation and can be mathematically described by the diffusion equation. Recently, powerful integral formulations have been developed for numerically solving the diffusion equation under complex conditions. In the literature, it is customary to formulate the integral equations by integrating point differential equations. Instead, in this paper, we shall employ a direct method of formulation, starting from the concepts of set and measure, the notion of partitions and the definition of set-averages.When the direct approach is applied to formulate the well-known finite element (FEM) equations, it is seen that the ‘Galerkin’ weighting function, which is mathematically treated as an artifice for weighting residuals, is but an appropriate spatial partition function. The logical framework of the direct approach is then applied to study the properties of ‘lumped’ and ‘consistent’ matrices arising in the use of the FEM. The lumped matrix, stemming naturally from the direct approach, seeks to conserve mass locally as well as globally, while the consistent matrix, which results only when the differential equation is integrated in a specific fashion, attempts only to preserve global mass balance.It is concluded that the direct approach is simple and complete and, in so far as the integral formulation is concerned, there is little to be gained in starting with the differential equation. Further, in formulating integral equations, it is common practice to evaluate only the time-dependent changes in the mass content of the system and ignore the evaluation of the mass content of the system at any given instant of time. In order to be complete in itself, a true integral approach should evaluate both the time-dependent changes in the mass content of the system as well as the instantaneous mass content at any given time.     

Linearization and first‐order expansion of the rocking motion of rigid blocks stepping on viscoelastic foundation

In structural mechanics there are several occasions where a linearized formulation of the original non‐linear problem reduces considerably the computational effort for the response analysis. In a broader sense, a linearized formulation can be viewed as a first‐order expansion of the dynamic equilibrium of the system about a ‘static’ configuration; yet caution should be exercised when identifying the ‘correct’ static configuration. This paper uses as a case study the rocking response of a rigid block stepping on viscoelastic supports, whose non‐linear dynamics is the subject of the companion paper, and elaborates on the challenge of identifying the most appropriate static configuration around which a first‐order expansion will produce the most dependable results in each regime of motion. For the regime when the heel of the block separates, a revised set of linearized equations is presented, which is an improvement to the unconservative equations published previously in the literature. The associated eigenvalues demonstrate that the characteristics of the foundation do not affect the rocking motion of the block once the heel separates. Copyright © 2008 John Wiley & Sons, Ltd.     

SH波入射时柔性基础上等腰三角形坝体结构的出平面反应

本文利用复变函数和坐标移动方法研究了SH波入射对柔性基础上等腰三角形坝体结构的影响。首先建立问题的数学模型并根据分区和辅助函数法将模型分割为2部分,其1为等腰三角形和半圆形组成的区域Ⅰ,其余为区域Ⅱ;其2在区域Ⅰ内构造1个满足等腰三角形两边应力自由的驻波解,在区域Ⅱ内构造满足水平边界应力自由的散射波;通过移动坐标在区域Ⅰ、Ⅱ的公共边界实现位移和应力的连续,建立起求解该问题的无穷代数方程组;最后,本文给出了例题和数值结果并对其进行了讨论,并通过算例强调了与文献[9]的本质区别。     

Simulation of synthetic ground motions for specified earthquake and site characteristics

A method for generating a suite of synthetic ground motion time‐histories for specified earthquake and site characteristics defining a design scenario is presented. The method employs a parameterized stochastic model that is based on a modulated, filtered white‐noise process. The model parameters characterize the evolving intensity, predominant frequency, and bandwidth of the acceleration time‐history, and can be identified by matching the statistics of the model to the statistics of a target‐recorded accelerogram. Sample ‘observations’ of the parameters are obtained by fitting the model to a subset of the NGA database for far‐field strong ground motion records on firm ground. Using this sample, predictive equations are developed for the model parameters in terms of the faulting mechanism, earthquake magnitude, source‐to‐site distance, and the site shear‐wave velocity. For any specified set of these earthquake and site characteristics, sets of the model parameters are generated, which are in turn used in the stochastic model to generate the ensemble of synthetic ground motions. The resulting synthetic acceleration as well as corresponding velocity and displacement time‐histories capture the main features of real earthquake ground motions, including the intensity, duration, spectral content, and peak values. Furthermore, the statistics of their resulting elastic response spectra closely agree with both the median and the variability of response spectra of recorded ground motions, as reflected in the existing prediction equations based on the NGA database. The proposed method can be used in seismic design and analysis in conjunction with or instead of recorded ground motions. Copyright © 2010 John Wiley & Sons, Ltd.     

Discrete crack modelling for dynamically loaded,unreinforced concrete structures

The dynamic response of unreinforced concrete structures is studied taking account of initiation, extension, closing and reopening of so-called discrete cracks. The computational procedure is based on the finite-element method and is at present restricted to two-dimensional situations. The discrete cracks are simulated by separation of originally adjacent finite elements. An equivalent tensile-strength criterion is used for the initiation and extension of the cracks which are assumed to propagate perpendicularly to the principal tensile stress. If this direction does not coincide with the interelement boundaries of the finite-element mesh, the latter is automatically altered. Between elements being separated by a crack special ‘crack elements’ are introduced, which take account of the stress transfer by aggregate interlock. The equations of motion are integrated numerically using an explicit formulation. The procedures outlined are demonstrated on a simplified cross-section of a concrete gravity dam subjected to horizontal earthquake excitation.     

Kinematic transformations for planar multi‐directional pseudodynamic testing

The pseudodynamic (PSD) test method imposes command displacements to a test structure for a given time step. The measured restoring forces and displaced position achieved in the test structure are then used to integrate the equations of motion to determine the command displacements for the next time step. Multi‐directional displacements of the test structure can introduce error in the measured restoring forces and displaced position. The subsequently determined command displacements will not be correct unless the effects of the multi‐directional displacements are considered. This paper presents two approaches for correcting kinematic errors in planar multi‐directional PSD testing, where the test structure is loaded through a rigid loading block. The first approach, referred to as the incremental kinematic transformation method, employs linear displacement transformations within each time step. The second method, referred to as the total kinematic transformation method, is based on accurate nonlinear displacement transformations. Using three displacement sensors and the trigonometric law of cosines, this second method enables the simultaneous nonlinear equations that express the motion of the loading block to be solved without using iteration. The formulation and example applications for each method are given. Results from numerical simulations and laboratory experiments show that the total transformation method maintains accuracy, while the incremental transformation method may accumulate error if the incremental rotation of the loading block is not small over the time step. A procedure for estimating the incremental error in the incremental kinematic transformation method is presented as a means to predict and possibly control the error. Copyright © 2009 John Wiley & Sons, Ltd.     

曲线坐标系程函方程的求解方法研究

笛卡尔坐标系中经典的程函方程在静校正、叠前偏移、走时反演、地震定位、层析成像等许多地球物理工作都有应用,然而用其计算起伏地表的地震波走时时却比较困难.我们通过把曲线坐标系中的矩形网格映射到笛卡尔坐标系的贴体网格推导出了曲线坐标中的程函方程,此时,曲线坐标系的程函方程呈现为各向异性的程函方程(尽管在笛卡尔坐标系中介质是各向同同性的).然后尝试用求解各向同性程函方程的快速推进法和Lax-Friedrichs快速扫描算法来分别求解该方程.数值试验表明未加考虑各向异性程函方程与各向同性程函方程的差别而把求解各向同性程函方程的快速推进法直接拓展到曲线坐标中的程函方程的做法是错误的,而Lax-Friedrichs快速扫描算法总能稳定地求解曲线坐标系的程函方程,进而有效地处理了地表起伏的情况,得到稳定准确的计算结果.     

Solution of the transport equations using a moving coordinate system

A convection-diffusion equation arises from the conservation equations in miscible and immiscible flooding, thermal recovery, and water movement through desiccated soil. When the convection term dominates the diffusion term, the equations are very difficult to solve numerically. Owing to the hyperbolic character assumed for dominating convection, inaccurate, oscillating solutions result. A new solution technique minimizes the oscillations. The differential equation is transformed into a moving coordinate system which eliminates the convection term but makes the boundary location change in time. We illustrate the new method on two one-dimensional problems: the linear convection-diffusion equation and a non-linear diffusion type equation governing water movement through desiccated soil. Transforming the linear convection diffusion equation into a moving coordinate system gives a diffusion equation with time dependent boundary conditions. We apply orthogonal collocation on finite elements with a Crank-Nicholson time discretization. Comparisons are made to schemes using fixed coordinate systems. The equation describing movement of water in dry soil is a highly non-linear diffusion-type equation with coefficients varying over six orders of magnitude. We solve the equation in a coordinate system moving with a time-dependent velocity, which is determined by the location of the largest gradient of the solution. The finite difference technique with a variable grid size is applied, and a modified Crank-Nicholson technique is used for the temporal discretization. Comparisons are made to an exact solution obtained by similarity transformation, and with an ordinary finite difference scheme on a fixed coordinate system.     

Modelling of fluid–solid interactions using an adaptive mesh fluid model coupled with a combined finite–discrete element model

Fluid–structure interactions are modelled by coupling the finite element fluid/ocean model ‘Fluidity-ICOM’ with a combined finite–discrete element solid model ‘Y3D’. Because separate meshes are used for the fluids and solids, the present method is flexible in terms of discretisation schemes used for each material. Also, it can tackle multiple solids impacting on one another, without having ill-posed problems in the resolution of the fluid’s equations. Importantly, the proposed approach ensures that Newton’s third law is satisfied at the discrete level. This is done by first computing the action–reaction force on a supermesh, i.e. a function superspace of the fluid and solid meshes, and then projecting it to both meshes to use it as a source term in the fluid and solid equations. This paper demonstrates the properties of spatial conservation and accuracy of the method for a sphere immersed in a fluid, with prescribed fluid and solid velocities. While spatial conservation is shown to be independent of the mesh resolutions, accuracy requires fine resolutions in both fluid and solid meshes. It is further highlighted that unstructured meshes adapted to the solid concentration field reduce the numerical errors, in comparison with uniformly structured meshes with the same number of elements. The method is verified on flow past a falling sphere. Its potential for ocean applications is further shown through the simulation of vortex-induced vibrations of two cylinders and the flow past two flexible fibres.     

Contamination of a well in a uniform background flow

In this paper we describe the transport of pollution in groundwater in the neighbourhood of a well in a uniform background flow. We compute the rate at which contaminated particles reach the well as a function of the place of the source of pollution. The motion of a particle in a dispersive flow is seen as a random walk process. The Fokker-Planck equation for the random motion of a particle is transformed using the complex potential for the advective flow field. The resulting equation is solved asymptotically after a stretching transformation. Finally, the analytical solution is compared with results from Monte Carlo simulations with the random walk model. The method can be extended to arbitrary flow fields. Then by a numerical coordinate transformation the analytical results can still be employed.     

Increasing numerical efficiency of iterative solution for total least-squares in datum transformations

Cartesian coordinate transformation between two erroneous coordinate systems is considered within the Errors-In-Variables (EIV) model. The adjustment of this model is usually called the total Least-Squares (LS). There are many iterative algorithms given in geodetic literature for this adjustment. They give equivalent results for the same example and for the same user-defined convergence error tolerance. However, their convergence speed and stability are affected adversely if the coefficient matrix of the normal equations in the iterative solution is ill-conditioned. The well-known numerical techniques, such as regularization, shifting-scaling of the variables in the model, etc., for fixing this problem are not applied easily to the complicated equations of these algorithms. The EIV model for coordinate transformations can be considered as the nonlinear Gauss-Helmert (GH) model. The (weighted) standard LS adjustment of the iteratively linearized GH model yields the (weighted) total LS solution. It is uncomplicated to use the above-mentioned numerical techniques in this LS adjustment procedure. In this contribution, it is shown how properly diminished coordinate systems can be used in the iterative solution of this adjustment. Although its equations are mainly studied herein for 3D similarity transformation with differential rotations, they can be derived for other kinds of coordinate transformations as shown in the study. The convergence properties of the algorithms established based on the LS adjustment of the GH model are studied considering numerical examples. These examples show that using the diminished coordinates for both systems increases the numerical efficiency of the iterative solution for total LS in geodetic datum transformation: the corresponding algorithm working with the diminished coordinates converges much faster with an error of at least 10^-5 times smaller than the one working with the original coordinates.     

Modelling of fluid–solid interactions using an adaptive mesh fluid model coupled with a combined finite–discrete element model

Fluid–structure interactions are modelled by coupling the finite element fluid/ocean model ‘Fluidity-ICOM’ with a combined finite–discrete element solid model ‘Y3D’. Because separate meshes are used for the fluids and solids, the present method is flexible in terms of discretisation schemes used for each material. Also, it can tackle multiple solids impacting on one another, without having ill-posed problems in the resolution of the fluid’s equations. Importantly, the proposed approach ensures that Newton’s third law is satisfied at the discrete level. This is done by first computing the action–reaction force on a supermesh, i.e. a function superspace of the fluid and solid meshes, and then projecting it to both meshes to use it as a source term in the fluid and solid equations. This paper demonstrates the properties of spatial conservation and accuracy of the method for a sphere immersed in a fluid, with prescribed fluid and solid velocities. While spatial conservation is shown to be independent of the mesh resolutions, accuracy requires fine resolutions in both fluid and solid meshes. It is further highlighted that unstructured meshes adapted to the solid concentration field reduce the numerical errors, in comparison with uniformly structured meshes with the same number of elements. The method is verified on flow past a falling sphere. Its potential for ocean applications is further shown through the simulation of vortex-induced vibrations of two cylinders and the flow past two flexible fibres.