A global method for coupling transport with chemistry in heterogeneous porous media     

Laila Amir  Michel Kern 《Computational Geosciences》2010,14(3):465-481

Modeling reactive transport in porous media, using a local chemical equilibrium assumption, leads to a system of advection–diffusion PDEs coupled with algebraic equations. When solving this coupled system, the algebraic equations have to be solved at each grid point for each chemical species and at each time step. This leads to a coupled non-linear system. In this paper, a global solution approach that enables to keep the software codes for transport and chemistry distinct is proposed. The method applies the Newton–Krylov framework to the formulation for reactive transport used in operator splitting. The method is formulated in terms of total mobile and total fixed concentrations and uses the chemical solver as a black box, as it only requires that one be able to solve chemical equilibrium problems (and compute derivatives) without having to know the solution method. An additional advantage of the Newton–Krylov method is that the Jacobian is only needed as an operator in a Jacobian matrix times vector product. The proposed method is tested on the MoMaS reactive transport benchmark.  相似文献   

Estimating wave crest distributions using the method of L-moments     

Amir H. Izadparast  John M. Niedzwecki   《Applied Ocean Research》2009,31(1):37-43

The design of fixed or floating offshore structures requires accurate information of the met-ocean data at the intended offshore site. In the design process it is recognized that this environmental data is modified in the near-field by the interaction with the particular geometrical configuration of the offshore structure. This transformation of the incident wave field around and beneath an offshore structure presents a challenge for ocean engineers when specifying the wave gap elevation to avoid impact loads on the underside of the deck and inundation of the topsides. Thus, the accurate estimation of the wave crest distributions from measurements at various locations near and under the offshore structure during model test studies is essential. A semi-empirical approach is presented herein that builds upon the findings of previous studies and introduces the Method of L-moments. A three parameter model for a wave crest probability distribution function is presented and explicit relationships between the parameters of the distribution and its’ first three L-Moments are established. Furthermore, three narrow-band models from earlier research studies are reviewed and compared with the new model. Wave measurements from a mini-TLP model test program are used as the basis for comparison of the four distributions. The root-mean-square error is used as a metric to quantify the overall fit of the data and its accuracy in the high end tail of the data. The L-Moment model is shown to be more robust in representing the data in both the far-field and beneath the deck of the mini-TLP where the wave field demonstrates increased non-linear behavior.  相似文献   

Elasticity of a beryllium silicate (phenacite: Be2SiO4)     

Amir Yeganeh-Haeri  Donald J. Weidner 《Physics and Chemistry of Minerals》1989,16(4):360-364

The adiabatic single-crystal elastic moduli of a beryllium silicate (phenacite: Be₂SiO₄, trigonal, have been determined at atmospheric pressure and 22° C by Brillouin spectroscopy. The elastic stiffness moduli in gigapascals are: C ₁₁=341.9 C ₃₃=391.0 C ₄₄=91.4 C ₆₆= 96.9 C ₁₂=148.0 C ₁₃=136.0 C ₁₄= 0.1 C ₁₅= 3.5Overall, the elastic stiffness moduli for phenacite parallel and perpendicular to the c axis are comparable (i.e., it is almost cubic in its elastic signature). The elastic moduli can be rationalized in terms of division of the structure into two types of coordination polyhedra (1Si+2Be) with slightly different stiffnesses, which are linked to form a three dimensional framework. Values of the isothermal bulk modulus and the linear compressibilities, as determined from hydrostatic compression experiments of Hazen and Au (1986), are in good agreement with those obtained here. Combining the two studies indicates a low pressure derivative of the bulk modulus for phenacite.  相似文献   

An Assessment of Changes in Winter Cold and Warm Spells over Canada   总被引：1，自引：1，他引：1  

Shabbar  Amir  Bonsal  Barrie 《Natural Hazards》2003,29(2):173-188

The recent Third Assessment Report (TAR) of the Intergovernmental Panel onClimate Change (IPCC) indicated that observed 20th century changes in severalclimatic extremes are qualitatively consistent with those expected due to increasedgreenhouse gases. However, a lack of adequate data and analyses make conclusiveevidence of changing extremes somewhat difficult, particularly, in a global sense.In Canada, extreme temperature events, especially those during winter, can havemany adverse environmental and economic impacts. In light of the aforementionedIPCC report, the main focus of this analysis is to examine observed trends andvariability in the frequency, duration, and intensity of winter (Jan–Feb–Mar) cold and warm spells over Canada during the second half of the 20th century.Cold spell trends display substantial spatial variability across the country. From1950–1998, western Canada has experienced decreases in the frequency, duration, and intensity of cold spells, while in the east, distinct increases in the frequency and duration have occurred. These increases are likely associated with morefrequent occurrences of the positive phase of the North Atlantic Oscillation (NAO)during the last several decades. With regard to winter warm spells, significantincreases in both the frequency and duration of these episodes were observedacross most of Canada. One exception was found in the extreme northeasternregions, where warm spells are becoming shorter and less frequent. The resultsof this study are discussed within the context of climate warming expectations.  相似文献   

Spatial Modeling of Geotechnical Information Using GML   总被引：1，自引：0，他引：1  

Jean-Pierre Bardet  Amir Zand 《Transactions in GIS》2009,13(1):125-165

Geotechnical data characterizes the geometry and physical properties of soil deposits, which are determined from field observations, in-situ and laboratory tests, and engineering analyses. Until now, geotechnical information was mainly described and exchanged as relational data, which occasionally included spatial characteristics. Hereafter, the spatial features of geotechnical data are investigated by invoking the geotechnical data model of the Association of Geotechnical and Geoenvironmental Specialists (AGS), which is widely used in the United Kingdom to exchange geotechnical data in engineering practice. The basic geometric objects in AGS geotechnical data are first identified. These geometric objects are mapped to basic geometric features of the Geography Markup Language (GML). AGS data can therefore be fully rendered using GML-conformant schemas, which make geotechnical data readily importable into GML-aware applications. The data can be also imported to mainstream GIS applications using Extensible Markup Language (XML) transformations. This article demonstrates the rendition of AGS data format to a GML-conformant schema and illustrates the implementation of the new format through a few geotechnical examples.  相似文献   

Prediction and validation of traffic vibration reduction due to cement column stabilization     

Andrew T. Peplow  Amir M. Kaynia   《Soil Dynamics and Earthquake Engineering》2007,27(8):793-802

The transmission and reduction of vibrations in the far-field of the surface of the ground due to a surface load is investigated theoretically and validated with given field measurement data. The performance of a given stabilization column, located directly underneath the load, at a number of receiver positions is studied and measured in terms of insertion loss. A numerical model is presented, which enables the wave-field in the region of the column to be determined, based on an integral equation formulation of the problem which is solved using a boundary element approach. It is shown that the column has a beneficial effect at low frequencies especially in certain frequency bandwidths and is validated with field data. However, when the Rayleigh wavelength becomes short compared with the depth and width of the column adverse effects occur at some frequencies which are also observed in the far-field. Various depths of columns and material properties of the surrounding soil medium are studied and results presented so that some preliminary physical conclusions may be derived.  相似文献   

Estimating ground water recharge from topography, hydrogeology, and land cover   总被引：2，自引：0，他引：2  

Cherkauer DS  Ansari SA 《Ground water》2005,43(1):102-112

Proper management of ground water resources requires knowledge of the rates and spatial distribution of recharge to aquifers. This information is needed at scales ranging from that of individual communities to regional. This paper presents a methodology to calculate recharge from readily available ground surface information without long-term monitoring. The method is viewed as providing a reasonable, but conservative, first approximation of recharge, which can then be fine-tuned with other methods as time permits. Stream baseflow was measured as a surrogate for recharge in small watersheds in southeastern Wisconsin. It is equated to recharge (R) and then normalized to observed annual precipitation (P). Regression analysis was constrained by requiring that the independent and dependent variables be dimensionally consistent. It shows that R/P is controlled by three dimensionless ratios: (1) infiltrating to overland water flux, (2) vertical to lateral distance water must travel, and (3) percentage of land cover in the natural state. The individual watershed properties that comprise these ratios are now commonly available in GIS data bases. The empirical relationship for predicting R/P developed for the study watersheds is shown to be statistically viable and is then tested outside the study area and against other methods of calculating recharge. The method produces values that agree with baseflow separation from streamflow hydrographs (to within 15% to 20%), ground water budget analysis (4%), well hydrograph analysis (12%), and a distributed-parameter watershed model calibrated to total streamflow (18%). It has also reproduced the temporal variation over 5 yr observed at a well site with an average error < 12%.  相似文献   

Vertical eddy diffusivity estimations in Swan river estuary   总被引：1，自引：0，他引：1  

Amir Etemad-Shahidi  Jorg Imberger 《Dynamics of Atmospheres and Oceans》2005,39(3-4):175-187

A field experiment is described in which diapycnal diffusivity is estimated by direct and indirect methods in Swan river estuary, Perth, Western Australia. The microstructure profiles were collected in a narrow and straight part of this estuary using a portable flux profiler (PFP). The profiles were segmented into stationary parts and the rate of vertical eddy diffusivity was estimated for the segments within the pycnocline. The direct measurement showed that instantaneous flux could be positive or negative with a low net rate of vertical eddy diffusivity for mass of about 6.5 × 10⁻⁸ m² s⁻¹. All the indirect methods overestimated this rate. However, within the indirect methods, the method of Osborn yielded the poorest estimate while the method of Osborn and Cox gave the closest estimate.  相似文献   

Non-stationary return levels of CMIP5 multi-model temperature extremes     

Linyin Cheng  Thomas J. Phillips  Amir AghaKouchak 《Climate Dynamics》2015,44(11-12):2947-2963

  相似文献   

Axisymmetric interaction of a rigid disc with a transversely isotropic half‐space     

Amir Aabbas Katebi  Ali Khojasteh  Mohammad Rahimian  Ronald Y. S. Pak 《国际地质力学数值与分析法杂志》2010,34(12):1211-1236

A theoretical formulation is presented for the determination of the interaction of a vertically loaded disc embedded in a transversely isotropic half‐space. By means of a complete representation using a displacement potential, it is shown that the governing equations of motion for this class of problems can be uncoupled into a fourth‐order partial differential equation. With the aid of Hankel transforms, a relaxed treatment of the mixed‐boundary value problem is formulated as dual integral equations, which, in turn, are reduced to a Fredholm equation of the second kind. In addition to furnishing a unified view of existing solutions for zero and infinite embedments, the present treatment reveals a severe boundary‐layer phenomenon, which is apt to be of interest to this class of problems in general. The present solutions are analytically in exact agreement with the existing solutions for a half‐space with isotropic material properties. To confirm the accuracy of the numerical evaluation of the integrals involved, numerical results are included for cases of different degrees of the material anisotropy and compared with existing solutions. Further numerical examples are also presented to elucidate the influence of the degree of the material anisotropy on the response. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献