Parallel finite element algorithm with domain decomposition for three-dimensional magnetotelluric modelling     

Fabio I. Zyserman  Juan E. Santos   《Journal of Applied Geophysics》2000,44(4):351

We present a new finite element (FE) method for magnetotelluric modelling of three-dimensional conductivity structures. Maxwell's equations are treated as a system of first-order partial differential equations for the secondary fields. Absorbing boundary conditions are introduced, minimizing undesired boundary effects and allowing the use of small computational domains. The numerical algorithm presented here is an iterative, domain decomposition procedure employing a nonconforming FE space. It does not use global matrices, therefore allowing the modellization of large and complicated structures. The algorithm is naturally parallellizable, and we show results obtained in the IBM SP2 parallel supercomputer at Purdue University. The accuracy of the numerical method is verified by checking the computed solutions with the results of COMMEMI, the international project on the comparison of modelling methods for electromagnetic induction.  相似文献   

A high strain‐rate constitutive model for sand and its application in finite‐element analysis of tunnels subjected to blast     

William Higgins  Tanusree Chakraborty  Dipanjan Basu 《国际地质力学数值与分析法杂志》2013,37(15):2590-2610

The paper presents a constitutive model for simulating the high strain‐rate behavior of sands. Based on the concepts of critical‐state soil mechanics, the bounding surface plasticity theory and the overstress theory of viscoplasticity, the constitutive model simulates the high strain‐rate behavior of sands under uniaxial, triaxial and multi‐axial loading conditions. The model parameters are determined for Ottawa and Fontainebleau sands, and the performance of the model under extreme transient loading conditions is demonstrated through simulations of split Hopkinson pressure bar tests up to a strain rate of 2000/s. The constitutive model is implemented in a finite‐element analysis software Abaqus to analyze underground tunnels in sandy soil subjected to internal blast loads. Parametric studies are conducted to examine the effect of relative density and type of sand and of the depth of tunnel on the variation of stresses and deformations in the soil adjacent to the tunnels. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

A practical and efficient numerical scheme for the analysis of steady state unconfined seepage flows     

Mohammad Hossein Bazyar  Adel Graili 《国际地质力学数值与分析法杂志》2012,36(16):1793-1812

The scaled boundary finite‐element method (SBFEM), a novel semi‐analytical technique, is applied to the analysis of the confined and unconfined seepage flow. This method combines the advantages of the finite‐element method and the boundary element method. In this method, only the boundary of the domain is discretized; no fundamental solution is required, and singularity problems can be modeled rigorously. Anisotropic and nonhomogeneous materials satisfying similarity are modeled without additional efforts. In this paper, SBFE equations and solution procedures for the analysis of seepage flow are outlined. The accuracy of the proposed method in modeling singularity problems is demonstrated by analyzing seepage flow under a concrete dam with a cutoff at heel. As only the boundary is discretized, the variable mesh technique is advisable for modeling unconfined seepage analyses. The accuracy, effectiveness, and efficiency of the method are demonstrated by modeling several unconfined seepage flow problems. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

Numerical simulation of ductile shear localization: Plastic‐flow network and slip‐line field     

Sheng‐zu Wang  Jun‐xiu Sun  Jie‐yuan Ning 《国际地质力学数值与分析法杂志》2008,32(3):307-324

  相似文献   

Hydrogeochemical balance of forest umbrisol profiles (‘Sierra de Gata’, central western Spain)     

I. Menéndez  G. Moreno  J. F. Gallardo  J. Saavedra 《水文研究》2007,21(14):1949-1956

Three techniques for obtaining soil water solutions (gravitational and matrical waters extracted using both in situ tension lysimeters and in vitro pressure chambers) and their later chemical analysis were performed in order to know the evolution of the soil‐solution composition when water moves down through the soil, from the Ah soil horizon to the BwC‐ or C‐horizons of forest soils located in western Spain. Additionally, ion concentrations and water volumes of input waters to soil (canopy washout) and exported waters (drainage solutions from C‐horizons) were determined to establish the net balance of solutes in order to determine the rates of leaching or retention of ions. A generalized process of sorption or retention of most components (even Cl^?) was observed, from the soil surface to the C‐horizon, in both gravitational and matrical waters, with H₄SiO₄, Mn²⁺, Na⁺, and SO₄^2? being the net exported components from the soil through the groundwater. These results enhance the role of the recycling effect in these forest soils. The net percentages of elements retained in these forest soils, considering the inputs and the outputs balance, were 68% K⁺, 85% Ca²⁺, 58% Mg²⁺, 7% Al³⁺, 5% Fe³⁺, 34% Zn²⁺, 57% Cl^?, and 20% NO₃^?, and about 75% of dissolved organic carbon was mineralized. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

Preliminary spectral characterization of snow in a high altitude tropical glacier and potential effects of impurities in snow on albedo of tropical glaciers     

Jean Emmanuel Sicart  Yves Arnaud 《水文研究》2007,21(26):3642-3644

This paper presents some preliminary measurements of snow spectral reflectance on the tropical Bolivian Zongo glacier. Measurements show a correct agreement with theoretical spectral albedo of pure snow in the near infrared region, but lower values in the visible region (by 10–20%) probably due to aerosols contained in snow. Impurity contents ranged from 10 to 100 ppmw in one‐week‐old snow collected from the Zongo glacier, but measurements are scarce. Large amounts of snowfall partly compensate the proximity of dust sources in mid‐latitude glaciers, whereas on outer‐tropical glaciers precipitations are not abundant and are very seasonal, and sources of aerosols are proximate. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

Numerical study of soil heterogeneity effects on contaminant transport in unsaturated soil (model development and validation)     

M. Mousavi Nezhad  A. A. Javadi  A. Al‐Tabbaa  F. Abbasi 《国际地质力学数值与分析法杂志》2013,37(3):278-298

The movement of chemicals through soil to groundwater is a major cause of degradation of water resources. In many cases, serious human and stock health implications are associated with this form of pollution. The study of the effects of different factors involved in transport phenomena can provide valuable information to find the best remediation approaches. Numerical models are increasingly being used for predicting or analyzing solute transport processes in soils and groundwater. This article presents the development of a stochastic finite element model for the simulation of contaminant transport through soils with the main focus being on the incorporation of the effects of soil heterogeneity in the model. The governing equations of contaminant transport are presented. The mathematical framework and the numerical implementation of the model are described. The comparison of the results obtained from the developed stochastic model with those obtained from a deterministic method and some experimental results shows that the stochastic model is capable of predicting the transport of solutes in unsaturated soil with higher accuracy than deterministic one. The importance of the consideration of the effects of soil heterogeneity on contaminant fate is highlighted through a sensitivity analysis regarding the variance of saturated hydraulic conductivity as an index of soil heterogeneity. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

A virtual experiment technique on the elementary behaviour of granular materials with discrete element method     

Xia Li  Hai‐Sui Yu  Xiang‐Song Li 《国际地质力学数值与分析法杂志》2013,37(1):75-96

Multi‐scale investigations aided by the discrete element method (DEM) play a vital role for current state‐of‐the‐art research on the elementary behaviour of granular materials. Similar to laboratory tests, there are three important aspects to be considered carefully, which are the proper stress/strain definition and measurement, the application of target loading paths and the designed experiment setup, to be addressed in the present paper. Considering the volume sensitive characteristics of granular materials, in the proposed technique, the deformation of the tested specimen is controlled and measured by deformation gradient tensor involving both the undeformed configuration and the current configuration. Definitions of Biot strain and Cauchy stress are adopted. The expressions of them in terms of contact forces and particle displacements, respectively, are derived. The boundary of the tested specimen consists of rigid massless planar units. It is suggested that the representative element uses a convex polyhedral (polygonal) shape to minimize possible boundary arching effects. General loading paths are described by directly specifying the changes in the stress/strain invariants or directions. Loading can be applied in the strain‐controlled mode by specifying the translations and rotations of the boundary units, or in the stress‐controlled mode by using a servo‐control mechanism, or in the combination of the two methods to realize mixed boundary conditions. Taking the simulation results as the natural consequences originated from a complex system, virtual experiments provide particle‐scale information database to conduct multi‐scale investigations for better understanding in granular material behaviours and possible development of the constitutive theories provided the qualitative similarity between the simulation results from virtual experiments and observations on real material behaviour. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

Flat boundaries and their effect on sand testing     

G. Marketos  Malcolm D. Bolton 《国际地质力学数值与分析法杂志》2010,34(8):821-837

  相似文献   

Pore scale study of permeability and tortuosity for flow through particulate media using Lattice Boltzmann method     

Ali Ghassemi  Ali Pak 《国际地质力学数值与分析法杂志》2011,35(8):886-901

In this paper, Lattice Boltzmann method (LBM) has been used to study the effects of permeability and tortuosity on flow through saturated particulate media and identify the relationships between permeability and tortuosity with other parameters such as particles diameter, grain specific surface, and porosity. LBM is a simple kinematic model that can incorporate the essential physics of microscopic and mesoscopic processes involved in flow through granular soils. The obtained results indicate that the 2D LB model, due to its inherent theoretical advantages, is capable of demonstrating that the porosity and specific surface are the most influential parameters in determining the intrinsic permeability of granular media. The obtained results show that particles' size diameter has a two‐fold effect on the coefficient of permeability: one is through specific surface and the other is by tortuosity factor. Numerical study also reveals that tortuosity of granular soils decreases almost linearly with increasing the porosity. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献