Fully coupled elastoplastic hydro-mechanical analysis of unsaturated porous media using a meshfree method

This paper presents the first application of an advanced meshfree method, ie, the edge-based smoothed point interpolation method (ESPIM), in simulation of the coupled hydro-mechanical behaviour of unsaturated porous media. In the proposed technique, the problem domain is spatially discretised using a triangular background mesh, and the polynomial point interpolation method combined with a simple node selection scheme is adopted for creating nodal shape functions. Smoothing domains are formed on top of the background mesh, and a constant smoothed strain, created by applying the smoothing operation over the smoothing domains, is assigned to each smoothing domain. The deformation and flow models are developed based on the equilibrium equation of the mixture, and linear momentum and mass balance equations of the fluid phases, respectively. The effective stress approach is followed to account for the coupling between the flow and deformation models. Further coupling among the phases is captured through a hysteretic soil water retention model that evolves with changes in void ratio. An advanced elastoplastic constitutive model within the context of the bounding surface plasticity theory is employed for predicting the nonlinear behaviour of soil skeleton. Time discretisation is performed by adopting a three-point discretisation method with growing time steps to avoid temporal instabilities. A modified Newton-Raphson framework is designed for dealing with nonlinearities of the discretised system of equations. The performance of the numerical model is examined through a number of numerical examples. The state-of-the-art computational scheme developed is useful for simulation of geotechnical engineering problems involving unsaturated soils.     

Assessment of alternative adsorption models and global sensitivity analysis to characterize hexavalent chromium loss from soil to surface runoff

We investigate our ability to assess transfer of hexavalent chromium, Cr(VI), from the soil to surface runoff by considering the effect of coupling diverse adsorption models with a two‐layer solute transfer model. Our analyses are grounded on a set of two experiments associated with soils characterized by diverse particle size distributions. Our study is motivated by the observation that Cr(VI) is receiving much attention for the assessment of environmental risks due to its high solubility, mobility, and toxicological significance. Adsorption of Cr(VI) is considered to be at equilibrium in the mixing layer under our experimental conditions. Four adsorption models, that is, the Langmuir, Freundlich, Temkin, and linear models, constitute our set of alternative (competing) mathematical formulations. Experimental results reveal that the soil samples characterized by the finest grain sizes are associated with the highest release of Cr(VI) to runoff. We compare the relative abilities of the four models to interpret experimental results through maximum likelihood model calibration and four model identification criteria (i.e., the Akaike information criteria [AIC and AIC_C] and the Bayesian and Kashyap information criteria). Our study results enable us to rank the tested models on the basis of a set of posterior weights assigned to each of them. A classical variance‐based global sensitivity analysis is then performed to assess the relative importance of the uncertain parameters associated with each of the models considered, within subregions of the parameter space. In this context, the modelling strategy resulting from coupling the Langmuir isotherm with a two‐layer solute transfer model is then evaluated as the most skilful for the overall interpretation of both sets of experiments. Our results document that (a) the depth of the mixing layer is the most influential factor for all models tested, with the exception of the Freundlich isotherm, and (b) the total sensitivity of the adsorption parameters varies in time, with a trend to increase as time progresses for all of the models. These results suggest that adsorption has a significant effect on the uncertainty associated with the release of Cr(VI) from the soil to the surface runoff component.     

AN EVALUATION OF SAS/GRAPH SOFTWARE FOR COMPUTER CARTOGRAPHY1

The SAS® computer software system, widely used and respected for its capabilities in statistical analysis and data base management, now includes a new set of graphic and cartographic procedures called SAS GRAPH?. We have used these cartographic procedures in research on mapping ethno-cultural census data from metropolitan areas in Ontario and in undergraduate and graduate courses in computer cartography. On the basis of that experience, we describe and evaluate SAS/GRAPH'S cartographic capabilities and illustrate with maps drawn by various devices.     

SOME IMPORTANT CHARACTERISTICS OFSURFACE ALBEDO OF CHINA IN RECENT YEARS

Sdsealbedoisdefinedasaratioofreflectivesolarradiationtototalacradiation,whichcanindicatetheaborptionandreflectivecapabilityofthe~hforsolarradiation.Itisan~tfactortoinfluencethes~radiationbalance,espeiallyfortheformationandvariationoflocalndcroclirnate.Generally,theactualdatafromo~tionalstationSforrebationareedintheanalysisofsolacealal.HOwever,therepresentativenessOfthedateisIratedduetothescarcityofobservationalstationandunevendistributionofs~condition.Afterthe1960s,theanalysisfors~ealbedoise…     

Surface Photometry of Comet Hale-Bopp Carbon Images

A very long series of photographic observations of the comet Hale-Bopp has been made during January–April 1997 at the double astrograph (400/2000) of the Main Astronomical Observatory (Kyiv, Ukraine). Some of the cometary photos were obtained with two wide-band filter combinations. One of these combinations isolates C₂ emission, another — the nearby dust continuum. The images were digitized by means of AMDPH-XY machine and then calibrated following the standard procedure. After subtraction of the dust continuum the distribution of surface brightness in the C₂ emission coma of comet Hale-Bopp was studied. We found an asymmetric brightness distribution both pre- and post-perihelion. On 21.77 April 1997 a secondary brightness peak is found at the distance of 1.03 × 10⁵ km from the nucleus. It is possible that this peak is related to the extended source of the C₂ molecules. This revised version was published online in July 2006 with corrections to the Cover Date.     

The second order characteristic condition for plastic flow

In this paper the second order characteristic (discontinuous bifurcation) condition is derived for the granular flow (fully plastic) equations. This second order bifurcation equation is shown to be formally identical to the first order localization requirement during steady elastoplastic deformation provided the elastic compliance tensor is substituted for the product of the plastic multiplier with the flow Hessian. For isotropic yield and flow functions the invariant form of the characteristic condition is given in detail, as well as an alternative expression in adapted co‐ordinates. The characteristic condition can be regarded as defining a hardening function which is maximized to identify the critical angles. When the method is applied to 3D Coulomb flow, Mohr's 3D fracture plane conditions are obtained uniquely. Copyright © 2003 John Wiley & Sons, Ltd.     

Long-chain alkenones and U37′ variability along a south–north transect in the Western Pacific Ocean

In this paper, we present 50 surface water samples collected during the IMAGES III cruise (June–July 1997) along a transect from New Zealand to the China Sea (42°S–178°E, 21°N–120°E) covering a temperature range from 13.3 to 30.4 °C. A very worthwhile aspect of this study is a coupling of both biomarker (alkenone) and coccolithophorid counting. We show that the U₃₇^k′–temperature relationship is very similar to the Prahl et al. [Geochim. Cosmochim. Acta 52 (1988) 2203] culture calibration and to the global core top calibration of Müller et al. [Geochim. Cosmochim. Acta 62 (1998) 1757]. However, in the warmest surface waters of the Western Pacific ocean (>26.4 °C) where Gephyrocapsa oceanica is likely the most widespread species, the associated U₃₇^k′ has a constant value of 1.0. The consequence is that above this temperature threshold, U₃₇^k′ cannot be used as an accurate paleothermometer.