Accuracy of numerical simulations of contaminant transport in heterogeneous aquifers: A comparative study     

F.A. Radu  N. Suciu  J. HoffmannA. Vogel  O. Kolditz  C.-H. ParkS. Attinger 《Advances in water resources》2011,34(1):47-61

This work deals with a comparison of different numerical schemes for the simulation of contaminant transport in heterogeneous porous media. The numerical methods under consideration are Galerkin finite element (GFE), finite volume (FV), and mixed hybrid finite element (MHFE). Concerning the GFE we use linear and quadratic finite elements with and without upwind stabilization. Besides the classical MHFE a new and an upwind scheme are tested. We consider higher order finite volume schemes as well as two time discretization methods: backward Euler (BE) and the second order backward differentiation formula BDF (2). It is well known that numerical (or artificial) diffusion may cause large errors. Moreover, when the Péclet number is large, a numerical code without some stabilising techniques produces oscillating solutions. Upwind schemes increase the stability but show more numerical diffusion. In this paper we quantify the numerical diffusion for the different discretization schemes and its dependency on the Péclet number. We consider an academic example and a realistic simulation of solute transport in heterogeneous aquifer. In the latter case, the stochastic estimates used as reference were obtained with global random walk (GRW) simulations, free of numerical diffusion. The results presented can be used by researchers to test their numerical schemes and stabilization techniques for simulation of contaminant transport in groundwater.  相似文献   

Perturbations in the Anomalistic Period of Artificial Satellites Caused by the Direct Solar Radiation Pressure Centre for Astronomy and Space Sciences,Satellite Tracking Station No. 1132     

V. Mioc  E. Radu 《Astronomische Nachrichten》1979,300(6):313-315

Approximate formulae for the evaluation of the influence of direct solar radiation pressure (neglecting the POYNTING -ROBERTSON effect) on the anomalistic period of artificial satellites with quasi-circular and completely illuminated orbits are derived.  相似文献   

Enzymological classification of salt lakes in Romania     

Vasile Muntean  Radu Crisan  Daniela Paşca  Stefan Kiss  Mihail Drăgan-Bularda 《International Journal of Salt Lake Research》1996,5(1):35-44

Sediments of 56 salt lakes from Romania were studied enzymologically. The following 7 enzymatic and nonenzymatic catalytic activities have been measured: phosphatase, H₂O₂-splitting in nonautoclaved (catalase) and autoclaved samples, TTC reduction in nonautoclaved (dehydrogenase) and autoclaved samples, without or with glucose addition. A formula is proposed for the evaluation of the enzymatic potential of salt lake sediments or other habitats. A hierarchy of the studied lakes has been established on the basis of the enzymatic potential of their sediments expressed as enzymatic indicator.  相似文献   

Modelling chemo‐hydro‐mechanical behaviour of unsaturated clays: a feasibility study     

Z. Liu  N. Boukpeti  X. Li  F. Collin  J.‐P. Radu  T. Hueckel  R. Charlier 《国际地质力学数值与分析法杂志》2005,29(9):919-940

Effective capabilities of combined chemo‐elasto‐plastic and unsaturated soil models to simulate chemo‐hydro‐mechanical (CHM) behaviour of clays are examined in numerical simulations through selected boundary value problems. The objective is to investigate the feasibility of approaching such complex material behaviour numerically by combining two existing models. The chemo‐mechanical effects are described using the concept of chemical softening consisting of reduction of the pre‐consolidation pressure proposed originally by Hueckel (Can. Geotech. J. 1992; 29 :1071–1086; Int. J. Numer. Anal. Methods Geomech. 1997; 21 :43–72). An additional chemical softening mechanism is considered, consisting in a decrease of cohesion with an increase in contaminant concentration. The influence of partial saturation on the constitutive behaviour is modelled following Barcelona basic model (BBM) formulation (Géotech. 1990; 40 (3):405–430; Can. Geotech. J. 1992; 29 :1013–1032). The equilibrium equations combined with the CHM constitutive relations, and the governing equations for flow of fluids and contaminant transport, are solved numerically using finite element. The emphasis is laid on understanding the role that the individual chemical effects such as chemo‐elastic swelling, or chemo‐plastic consolidation, or finally, chemical loss of cohesion have in the overall response of the soil mass. The numerical problems analysed concern the chemical effects in response to wetting of a clay specimen with an organic liquid in rigid wall consolidometer, during biaxial loading up to failure, and in response to fresh water influx during tunnel excavation in swelling clay. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

Comparison between VELACS numerical ‘class A’ predictions and centrifuge experimental soil test results     

Radu Popescu  Jean H. Prevost 《Soil Dynamics and Earthquake Engineering》1995,14(2)

The numerical ‘class A’ predictions performed within the framework of the VELACS Project are compared to the experimental results recorded in the centrifuge experiments. The comparisons are made in terms of: (1) the root mean square error of the predictions with respect to the mean of the experimental results; and (2) the size of a confidence interval centered at the predicted value which contains the estimated true value of the experimental results with a 75% probability. An assessment of the capability of various groups of constitutive soil models to predict excess pore pressures induced by dynamic loading is also presented.  相似文献   

Isoseismal map of the Vrancea earthquake of March 4, 1977     

G. Radu  V. Krník  G. Polonic  D. Prochzkov  Z. Schenkov 《Tectonophysics》1979,53(3-4)

The iaoseismal map of the intermediate Vrancea earthquake of March 4, 1977 was compiled using the contributions from many European countries within the shaken area. Characteristic features of the macroseiamic field are discussed. Mean radii of isoseismals and intensity attenuation are determined.  相似文献   

Dynamics of nonlinear porous media with applications to soil liquefaction   总被引：2，自引：0，他引：2  

Radu Popescu  Jean H. Prevost  George Deodatis  Pradipta Chakrabortty 《Soil Dynamics and Earthquake Engineering》2006,26(6-7):648-665

  相似文献   

Orbits in an anisotropic radiation field     

V. Mioc  E. Radu 《Astronomische Nachrichten》1992,313(6):353-357

The elliptic-type motion around a source whose luminosity is anisotropic is being perturbatively treated. Using a Fourier expansion for the perturbing force, exact expressions describing the evolution of each orbital element are determined from Newton-Euler equations. For near-unit frequency the resonant solutions are pointed out. An approximate expression for the variation of the nodal period is determined, too. The behaviour of the radius vector of the unstable orbit (resonance case) is pointed out, showing periodic variations of constantly increasing amplitude. The (both angular and physical) time scale for escape is estimated. Concrete astronomical situations modellable in this way are mentioned.  相似文献   

An earthquake-source-based metric for seismic fragility analysis     

Alin Radu  Mircea Grigoriu 《Bulletin of Earthquake Engineering》2018,16(9):3771-3789

The seismic fragility of a system is the probability that the system enters a damage state under seismic ground motions with specified characteristics. Plots of the seismic fragilities with respect to scalar ground motion intensity measures are called fragility curves. Recent studies show that fragility curves may not be satisfactory measures for structural seismic performance, since scalar intensity measures cannot comprehensively characterize site seismicity. The limitations of traditional seismic intensity measures, e.g., peak ground acceleration or pseudo-spectral acceleration, are shown and discussed in detail. A bivariate vector with coordinates moment magnitude m and source-to-site distance r is proposed as an alternative seismic intensity measure. Implicitly, fragility surfaces in the (m, r)-space could be used as graphical representations of seismic fragility. Unlike fragility curves, which are functions of scalar intensity measures, fragility surfaces are characterized by two earthquake-hazard parameters, (m, r). The calculation of fragility surfaces may be computationally expensive for complex systems. Thus, as solutions to this issue, a bi-variate log-normal parametric model and an efficient calculation method, based on stochastic-reduced-order models, for fragility surfaces are proposed.  相似文献   

Investigation on regional attenuation of Vrancea (Romania) intermediate-depth earthquakes     

Florin Pavel  Radu Vacareanu 《地震工程与工程振动(英文版)》2018,17(3):501-509

This paper aims at investigating possible regional attenuation patterns in the case of Vrancea(Romania) intermediate-depth earthquakes.Almost 500 pairs of horizontal components recorded during 13 intermediate-depth Vrancea earthquakes are employed in order to evaluate the regional attenuation patterns.The recordings are grouped according to the azimuth with regard to the Vrancea seismic source and subsequently,Q models are computed for each azimuthal zone assuming similar geometrical spreading.Moreover,the local soil amplification which was disregarded in a previous analysis performed for Vrancea intermediate-depth earthquakes is now clearly evaluated.The results show minor differences between the four regions situated in front of the Carpathian Mountains and considerable differences in attenuation of seismic waves between the forearc and backarc regions(with regard to the Carpathian Mountains).Consequently,an average Q model of the type Q(f) = 115×f~(1.25) is obtained for the four forearc regions,while a separate Q model of the type Q(f) = 70×f~(0.90) is computed for the backarc region.These results highlight the need to evaluate the seismic hazard of Romania by using ground motion models which take into account the different attenuation between the forearc/backarc regions.  相似文献