An analytic element approach for modeling polygonal inhomogeneities in multi-aquifer systems

An analytic element approach is presented for the modeling of steady groundwater flow through multi-aquifer systems with piecewise constant aquifer and leaky layer properties. Different properties may be specified for domains bounded by closed polygons, referred to as polygonal inhomogeneities. The boundary of these inhomogeneities is modeled with two types of high-order line elements. First, a string of single-aquifer line-doublets is used; these elements cut through all aquifers and are valid both inside and outside the inhomogeneity. Second, two strings of multi-aquifer line-sinks are used, one string that is valid inside the inhomogeneity and one string that is valid outside; the comprehensive extraction of these line-sinks is zero at any point along the string. The proposed approach results in a comprehensive flow field of which the component normal to the boundary of the inhomogeneity is continuous across the boundary at any point. Within each individual aquifer, continuity of head and the component of flow normal to the boundary are met approximately across the boundary; the accuracy increases when the order of the line elements is increased and/or when shorter line elements are used. The proposed analytic element approach produces results that are virtually identical to the exact solution for a cylindrical inhomogeneity, and a high-resolution MODFLOW2000 model of two rectangular inhomogeneities with a shared boundary. The practical application of the approach is demonstrated through the solution of a problem with an irregularly shaped inhomogeneity with rivers crossing the inhomogeneity boundary.     

On magnetic boundary conditions for non-spectral dynamo simulations

We address issues associated with non-local magnetic boundary conditions for non-spectral dynamo simulations. We introduce an integro-differential formulation for a domain bounded by an insulating outer domain. We show how to combine the flexibility of a local discretisation with a rigorous formulation of magnetic boundary conditions in arbitrary geometries. This formulation substantiates from mathematical point of view a new method for numerical solution of magnetohydrodynamic problems with non-local boundary conditions based on coupling finite volumes and boundary elements. Finally, we discuss practical efficiency of this new method.     

Analytic element modelling for strip aquifers

The Analytic Element Method (AEM) provides a convenient tool for groundwater flow analysis in unbounded continuous domains. The AEM is based on the superposition of analytic functions, known as elements, useful at both regional and local scales. In this study, analytic elements for strip aquifers are presented. Such aquifers occur in riverine or coastal deposits and in outcrop zones of confined aquifers. Local flow field is modelled indirectly, using a reference plane related to the aquifer domain through the Schwarz‐Christoffel transform. The regional flow is obtained as a solution of the one‐dimensional flow equation. The proposed methodology was tested by modelling two hypothetical situations, which were compared to exact solutions. It is shown that regional boundaries can be reproduced exactly while local fields are adequately reproduced with analytic elements. The developed elements are applied to simulate a real flow field in northeastern Brazil showing good agreement with measured water levels. Copyright © 2011 John Wiley & Sons, Ltd.     

Impedance matrices for circular foundation embedded in layered medium

A numerical scheme is developed in the paper for calculating torsional, vertical, horizontal, coupling and rocking impedances in frequency domain for axial-symmetric foundations embedded in layered media. In the scheme, the whole soil domain is divided into interior and exterior domains. For the exterior domain, the analytic solutions with unknown coefficients are obtained by solving three-dimensional (3D) wave equations in cylindrical coordinates satisfying homogeneous boundary conditions. For the interior domain, the analytical solutions are also obtained by solving the same 3D wave equations satisfying the homogeneous boundary conditions and the prescribed boundary conditions. The prescribed conditions are the interaction tractions at the interfaces between embedded foundation and surrounding soil. The interaction tractions are assumed to be piecewise linear. The piecewise linear tractions at the bottom surface of foundation will be decomposed into a series of Bessel functions which can be easily fitted into the general solutions of wave equations in cylindrical coordinates. After all the analytic solutions with unknown coefficients for both interior and exterior domains are found, the variational principle is employed using the continuity conditions (both displacements and stresses) at the interfaces between interior and exterior domains, interior domain and foundation, and exterior domain and foundation to find impedance functions.     

Analytic element modeling of embedded multiaquifer domains

An analytic element approach is presented for the modeling of multiaquifer domains embedded in a single-aquifer model. The inside of each domain may consist of an arbitrary number of aquifers separated by leaky layers. The analytic element solution is obtained through a combination of existing single-aquifer and multiaquifer analytic elements and allows for the analytic computation of head and leakage at any point in the aquifer. Along the boundary of an embedded multiaquifer domain, the normal flux is continuous everywhere; continuity of head across the boundary is met exactly at collocations points and approximately, but very accurately, in between. The analytic element solution compares well with an existing exact solution. A hypothetical example with a river intersecting two embedded domains illustrates the practical application of the proposed approach.     

Scattering of SH waves induced by a symmetrical V-shaped canyon: a unified analytical solution

This paper reports a series solution of wave functions for two-dimensional scattering and diffraction of plane SH waves induced by a symmetrical V-shaped canyon with different shape ratios. A half-space with a symmetrical V-shaped canyon is divided into two sub-regions by using a circular-arc auxiliary boundary. The two sub-regions are represented by global and local cylindrical coordinate systems, respectively. In each coordinate system, the wave field satisfying the Helmholtz equation is represented by the separation of variables method, in terms of the series of both Bessel functions and Hankel functions with unknown complex coefficients. Then, the two wave fields are described in the local coordinate system using the Graf addition theorem. Finally, the unknown coefficients are sought by satisfying the continuity conditions of the auxiliary boundary. To consider the phase characteristics of the wave scattering, a parametric analysis is carried out in the time domain by assuming an incident signal of the Ricker type. Surface and subsurface transient responses demonstrate the characteristics and mechanisms of wave propagating and scattering.     

New analytical solutions for groundwater flow in wedge-shaped aquifers with various topographic boundary conditions

The hydraulic head distribution in a wedge-shaped aquifer depends on the wedge angle and the topographic and hydrogeological boundary conditions. In addition, an equation in terms of the radial distance with trigonometric functions along the boundary may be suitable to describe the water level configuration for a valley flank with a gentle sloping and rolling topography. This paper develops a general mathematical model including the governing equation and a variety of boundary conditions for the groundwater flow within a wedge-shaped aquifer. Based on the model, a new closed-form solution for transient flow in the wedge-shaped aquifer is derived via the finite sine transform and Hankel transform. In addition, a numerical approach, including the roots search scheme, the Gaussian quadrature, and Shanks’ method, is proposed for efficiently evaluating the infinite series and the infinite integral presented in the solution. This solution may be used to describe the head distribution for wedges that image theory is inapplicable, and to explore the effects of the recharge from various topographic boundaries on the groundwater flow system within a wedge-shaped aquifer.     

On the well posedness of some wave formulations of the shallow water equations

According to the theory of characteristics, the number of boundary conditions required for the adequate definition of a PDE problem is equal to the number of characteristic half planes entering the domain associated with the PDE problem.This theory was applied to the primitive form of the shallow water equations two decades ago to determine the number of initial and boundary conditions required by these equations. The results of this early study are remembered here. Subsequently, the same theory is applied to some wave formulations of the shallow water equations (the wave continuity and primitive momentum equation model, and the wave continuity and wave momentum equation model).Circumstances under which the number of boundary conditions required by the mathematical model can be reduced to the number of available boundary conditions are discussed for both the primitive and wave formulations of the shallow water equations.     

Transient analysis of wave propagation in non‐homogeneous elastic unbounded domains by using the scaled boundary finite‐element method

The scaled boundary finite‐element method has been developed for the dynamic analysis of unbounded domains. In this method only the boundary is discretized resulting in a reduction of the spatial dimension by one. Like the finite‐element method no fundamental solution is required. This paper extends the scaled boundary finite‐element method to simulate the transient response of non‐homogeneous unbounded domains with the elasticity modulus and mass density varying as power functions of spatial coordinates. To reduce the number of degrees of freedom and the computational cost, the technique of reduced set of base functions is applied. The scaled boundary finite‐element equation for an unbounded domain is reformulated in generalized coordinates. The resulting acceleration unit‐impulse response matrix is obtained and assembled with the equation of motion of standard finite elements. Numerical examples of non‐homogeneous isotropic and transversely isotropic unbounded domains demonstrate the accuracy of the scaled boundary finite‐element method. Copyright © 2006 John Wiley & Sons, Ltd.     

Assessment of managed aquifer recharge potential using ensembles of local models

A simple quantitative approach for assessing the artificial recharge potential of large regions using spatial ensembles of local models is proposed. The method extends existing qualitative approaches and enables rapid assessments within a programmable environment. Spatial discretization of a water resource region into continuous local domains allows simple local models to be applied independently in each domain using lumped parameters. The ensemble results can be analyzed directly or combined with other quantitative and thematic information and visualized as regional suitability maps. A case study considers the hydraulic potential for surface infiltration across a large water resource region using a published analytic model for basin recharge. The model solution was implemented within a geographic information system and evaluated independently in >21,000 local domains using lumped parameters derived from existing regional datasets. Computer execution times to run the whole ensemble and process the results were in the order of a few minutes. Relevant aspects of the case study results and general conclusions concerning the utility and limitations of the method are discussed.     

A new analytical solution for assessing climate change impacts on subsurface temperature

Groundwater temperature is an important water quality parameter that affects species distributions in subsurface and surface environments. To investigate the response of subsurface temperature to atmospheric climate change, an analytical solution is derived for a one‐dimensional, transient conduction–advection equation and verified with numerical methods using the finite element code SUTRA. The solution can be directly applied to forward model the impact of future climate change on subsurface temperature profiles or inversely applied to produce a surface temperature history from measured borehole profiles. The initial conditions are represented using superimposed linear and exponential functions, and the boundary condition is expressed as an exponential function. This solution expands on a classic solution in which the initial and boundary conditions were restricted to linear functions. The exponential functions allow more flexibility in matching climate model projections (boundary conditions) and measured temperature–depth profiles (initial conditions). For example, measured borehole temperature data from the Sendai Plain and Tokyo, Japan, were used to demonstrate the improved accuracy of the exponential function for replicating temperature–depth profiles. Also, the improved accuracy of the exponential boundary condition was demonstrated using air temperature anomaly data from the Intergovernmental Panel on Climate Change. These air temperature anomalies were then used to forward model the effect of surficial thermal perturbations in subsurface environments with significant groundwater flow. The simulation results indicate that recharge can accelerate shallow subsurface warming, whereas upward groundwater discharge can enhance deeper subsurface warming. Additionally, the simulation results demonstrate that future groundwater temperatures obtained from the proposed analytical solution can deviate significantly from those produced with the classic solution. Copyright © 2013 John Wiley & Sons, Ltd.     

基于无单元Galerkin法探地雷达正演模拟

无单元Galerkin法采用滑动最小二乘法拟合场函数,只需节点无需单元,具有前处理简单、精度高、解高次连续等优点,被用于求解探地雷达(GPR)正问题.本文从Maxwell方程出发,推导了GPR正演需满足的波动方程;详细介绍了滑动最小二乘法形函数的构造方法.针对EFGM不满足插值条件导致强加边界条件的处理变复杂的特性,采用罚因子法对强加边界条件进行了处理;同时为了消除EFGM进行GPR正演模拟时来自截断边界处的超强反射,采用透射边界条件把GPR波在截断边界处的反射波透射出去,进而压制了来自截断边界处的反射波.然后,编制了EFGM的GPR正演模拟Matlab程序,应用该程序对典型GPR地电模型进行了正演模拟,并把该正演剖面图与基于线性插值FEM正演剖面图进行了对比,结果表明了EFGM用于GPR正演计算的正确性及有效性,并且在相同节点数条件下,EFGM比矩形剖分的FEM的精度要高,更有利于指导雷达剖面的数据解译.     

An Hermitian finite element solution of the two-dimensional saturated-unsaturated flow equation

This paper describes a Galerkin-type finite element solution of the two-dimensional saturated-unsaturated flow equation. The numerical solution uses an incomplete (reduced) set of Hermitian cubic basis functions and is formulated in terms of normal and tangential coordinates. The formulation leads to continuous pressure gradients across interelement boundaries for a number of well-defined element configurations, such as for rectangular and circular elements. Other elements generally lead to discontinuous gradients; however, the gradients remain uniquely defined at the nodes. The method avoids calculation of second-order derivatives, yet retains many of the advantages associated with Hermitian elements. A nine-point Lobatto-type integration scheme is used to evaluate all local element integrals. This alternative scheme produces about the same accuracy as the usual 9- or 16-point Gaussian quadrature schemes, but is computationally more efficient.     

A parallel domain decomposition algorithm for coastal ocean circulation models based on integer linear programming

This paper presents a new parallel domain decomposition algorithm based on integer linear programming (ILP), a mathematical optimization method. To minimize the computation time of coastal ocean circulation models, the ILP decomposition algorithm divides the global domain in local domains with balanced work load according to the number of processors and avoids computations over as many as land grid cells as possible. In addition, it maintains the use of logically rectangular local domains and achieves the exact same results as traditional domain decomposition algorithms (such as Cartesian decomposition). However, the ILP decomposition algorithm may not converge to an exact solution for relatively large domains. To overcome this problem, we developed two ILP decomposition formulations. The first one (complete formulation) has no additional restriction, although it is impractical for large global domains. The second one (feasible) imposes local domains with the same dimensions and looks for the feasibility of such decomposition, which allows much larger global domains. Parallel performance of both ILP formulations is compared to a base Cartesian decomposition by simulating two cases with the newly created parallel version of the Stevens Institute of Technology’s Estuarine and Coastal Ocean Model (sECOM). Simulations with the ILP formulations run always faster than the ones with the base decomposition, and the complete formulation is better than the feasible one when it is applicable. In addition, parallel efficiency with the ILP decomposition may be greater than one.     

Interaction of an axisymmetric body with obliquely incident seismic waves by global local finite elements

Asymmetric steady-state structure-media interaction due to obliquely incident body waves is investigated via a version of the global local finite element method. In the present version, a local region that houses an axisymmetric structure is modelled by conventional finite elements, while the behaviour in the remaining portion of the homogeneous semi-infinite medium is presented by the spherical harmonics that are the eigensolutions of the entire space problem. The solution scheme involves (1) full displacement and traction continuity along the boundary between the local and the exterior regions and (2) satisfaction of the traction-free requirement on the surface of the half-space beyond the discretized region by virtue of a sequence of integral constraints of the non-zero weighted surface tractions of the spherical harmonics. The numerical results presented are for a perfectly bonded rigid circular foundation resting on the surface of the half-space and subjected to obliquely incident body waves. Dependence of the displacement response of the footing upon incident angles and dimensionless wave numbers is thoroughly studied.     

卫星磁异常的理论模型

本文介绍了计算卫星磁异常理论模型的数学方法 ,即球谐分析方法、冠谐分析方法、矩谐分析方法和等效源方法 .根据相同的 MAGSAT资料 ,计算的卫星磁异常冠谐模型、矩谐模型和等效源模型都能很好地表示卫星磁异常的分布 .由于在整个研究区域 ( 1 0°N～60°N,70°E～ 1 40°E)都有卫星资料 ,所以这些理论模型没有所谓的“边界效应”.这一结论对计算地磁场的区域模型是很有意义的 .     

Heat diffusion and size reduction of a spherical magma chamber

An analytical solution is derived for the size reduction of a spherical magma chamber cooling by conduction. The use of moving boundary conditions and the constraint of a spherical symmetry allow one to ignore the details of the heat redistribution processes which take place within the magma chamber. The dependence of the solution on the initial conditions is investigated. A simple solution is found for short time, which is shown to be valid for times long enough to make it useful in the volcanological context. Moreover, the general solution confirms that the hydrothermal contribution to heat transfer in Phlegraean Fields cannot be extremely important.     

A hybrid modelling of soil–structure interaction problems for deeply embedded structures in A multilayered medium

Dynamic response of deeply embedded structures, such as underground tunnels and deep foundations, in a multilayered elastic half-space are analysed when the structure is excited by a plane P or SV wave propagating at some angle. The scattered field is represented by the sum of three Green's functions, corresponding to two oscillating forces and one oscillating moment at the centroid position of the buried structure. The amplitudes of these two forces and one moment are a priori unknown and are obtained by satisfying displacement and stress continuity conditions across the near-field/far-field boundary. The distinguishing feature of this technique from direct or indirect boundary integral techniques is that in these techniques a distribution of sources of unknown amplitude are considered at the near-field/far-field boundary, and a large number of sources are needed for different combinations of source-receiver arrangements. But in this technique the sources of unknown amplitude are placed at the location of the structure, not at the near-field/far-field boundary and, using the Saint Venant's principle, the scattered field is modelled. Thus, the number of sources required is reduced to only three. Two example problems are solved. The first one is for a deeply embedded footing in a three-layer soil mass and the second one is for a rectangular tunnel in a two-layer soil mass.     

考虑刚度及边界条件的短索索力求解与试验研究

建立了一个非齐次边界拉索力学模型,基于拉索线性振动理论,推导了考虑拉索抗弯刚度及边界条件的拉索频率方程;根据双曲函数的数理特性对频率方程进行了合理简化,得到了实用的频率简化方程,提出了逆过程求解索力的方法和步骤,并探讨了不同支承刚度对索力的影响。设计了室内拉索试验模型并进行了27种工况索力试验,比较了不同工况下的计算索力与实际索力,验证了本文提出的频率简化方程和逆过程求解方法的实用可靠性。     

Hydroelastic vibration and buckling of rectangular Mindlin plates on Pasternak foundations under linearly varying in-plane loads

Hydroelastic vibration and buckling analysis of horizontal rectangular plates resting on Pasternak foundation and subjected to linearly varying in-plane loads is investigated for different boundary conditions. Fluid is assumed to be inviscid and incompressible. To analyze the interaction of the Mindlin plate with the elastic foundation and fluid system, three displacement components of the plate are expressed in the Ritz method by adopting a set of static Timoshenko beam functions satisfying geometric boundary conditions. Convergence study is first carried out to confirm the stability of the present method. Then, a few comparisons are made for limited cases between present results and the available data in the literature. Finally, parametric study is conducted to highlight the effects of loading factors, fluid levels, foundation stiffness parameters, aspect ratios, thickness to width ratios and different boundary conditions on the critical buckling loads and wet natural frequencies of rectangular plates.