Three-dimensional finite element simulation of fracture propagation in rock specimens with pre-existing fissure(s) under compression and their strength analysis

A FORTRAN program, consistent with the commercially available finite element (FE) code ABAQUS, is developed based on a three-dimensional (3D) linear elastic brittle damage constitutive model with two damage criteria. To consider the heterogeneity of rock, the developed FORTRAN program is used to set the stiffness and strength properties of each element of the FE model following a Weibull distribution function. The reliability of the program is assessed against available experimental results for granite cylindrical specimens with a throughgoing, flat and inclined fissure. The calibration procedure of the material parameters is explained in detail, and it is shown that the compressive to tensile strength ratio can have a substantial influence on the failure response of the specimens. Numerical simulations are conducted for models with different levels of heterogeneity. The results show a smaller load bearing capacity for models with less homogeneity, representing gradual coalescence of fully damaged elements forming throughout the models during loading. The maximum load bearing capacity is studied for various combinations of inclination angles of two centrally aligned, throughgoing and flat fissures of equal length embedded in cylindrical models under uniaxial and multiaxial loading conditions. The key role of the compressive to tensile strength ratio is highlighted by repeating certain simulations with a lower compressive to tensile strength ratio. It is proven that the peak loads of the rock models with sufficiently small compressive to tensile strength ratios containing two throughgoing fissures of equal length are similar, provided that the minimum inclination angles of the models are the same. The results are presented and discussed with respect to the existing experimental findings in the literature, suggesting that the numerical model applied in this study can provide useful insight into the failure behaviour of rock-like materials.     

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.     

Dynamic modelling of retrogressive landslides with emphasis on the role of clay sensitivity

This paper presents a detailed numerical study of the retrogressive failure of landslides in sensitive clays. The dynamic modelling of the landslides is carried out using a novel continuum approach, the particle finite element method, complemented with an elastoviscoplastic constitutive model. The multiwedge failure mode in the collapse is captured successfully, and the multiple retrogressive failures that have been widely observed in landslides in sensitive clays are reproduced with the failure mechanism, the kinematics, and the deposition being discussed in detail. Special attention has been paid to the role of the clay sensitivity on each retrogressive failure as well as on the final retrogression distance and the final run‐out distance via parametric studies. Moreover, the effects of the viscosity of sensitive clays on the failure are also investigated for different clay sensitivities.     

Going underground: in search of Carboniferous coal forests

The development of coal forests during the Carboniferous is one of the best-known episodes in the history of life. Although often reconstructed as steamy tropical rainforests, these ancient ecosystems were a far cry from anything we might encounter in the Amazon today. Bizarre giant club-mosses, horsetails and tree ferns were the dominant plants, not flowering trees as in modern rainforests. At their height, coal forests stretched all the way from Kansas to Kazakhstan, spanning the entire breadth of tropical Pangaea. Most of what we know of their biodiversity and ecology has been quite literally mined out of the ground through two centuries of hard labour. Without coal mining, our knowledge would be greatly impoverished. Over the past few years, we've been exploring underground coal mines in the United States, where entire forested landscapes have been preserved intact over huge areas. Never before have geologists had the opportunity to walk out through mile upon mile of fossilized forest. In this feature article, we describe some of our recent explorations and attempt to shed new light on these old fossils.     

Finite element formulation and algorithms for unsaturated soils. Part II: Verification and application

The finite‐element formulation and integration algorithms developed in Part I are used to analyse a number of practical problems involving unsaturated and saturated soils. The formulation and algorithms perform well for all the cases analysed, with the robustness of the latter being largely insensitive to user‐defined parameters such as the number of coarse time steps and error control tolerances. The efficiency of the algorithms, as measured by the CPU time consumed, does not depend on the number of coarse time steps, but may be influenced by the error control tolerances. Based on the analyses presented here, typical values for the error control tolerances are suggested. It is also shown that the constitutive modelling framework presented in Part I can, by adjusting one constitutive equation and one or two material parameters, be used to simulate soils that expand or collapse upon wetting. Treating the suction as a strain variable instead of a stress variable proves to be an efficient and robust way of solving suction‐dependent plastic yielding. Moreover, the concept of the constitutive stress is a particularly convenient way of handling the transition between saturation and unsaturation. Copyright © 2003 John Wiley & Sons, Ltd.     

Richtmyer-Meshkov Experiments on the Omega Laser

Observations of the interstellar medium reveal a dynamic realm permeated by shocks. These shocks are generated on a large range of scales by galactic rotation, supernovae, stellar winds, and other processes. Whenever a shock encounters a density interface, Richtmyer-Meshkov instabilities may develop. Perturbations along the interface grow, leading to structure formation and material mixing. An understanding of the evolution of Richtmyer-Meshkov instabilities is essential for understanding galactic structure, molecular cloud morphology, and the early stages of star formation. An ongoing experimental campaign studies Richtmyer-Meshkov mixing in a convergent, compressible, miscible plasma at the Omega laser facility. Cylindrical targets, consisting of a low density foam core and an aluminum shell covered by an epoxy ablator, are directly driven by fifty laser beams. The aluminum shell is machined to produce different perturbation spectra. Surface types include unperturbed (smooth), single-mode sinusoids, multi-mode (rough), and multi-mode with particular modes accentuated (specified-rough). Experimental results are compared to theory and numerical simulations.     

A techno-economic comparison of power systems for autonomousunderwater vehicles

This paper compares lead-acid batteries, sodium-sulfur batteries, solid polymer fuel cells and closed-cycle diesel engines for autonomous underwater vehicle (AUV) applications. The service is described in terms of a parametric mission and life cycle. A generic AUV is used as a basis for comparison. Power systems are evaluated by two criteria: (1) submerged endurance capability and (2) life cycle cost. This study determines categories of service for which each power system is preferred. The solid polymer fuel cell can provide greater submerged endurance than other power systems examined. For extremely long duration AUV missions, the fuel cell is the required system, indicating a possible market niche for today's fuel cell technology. Considering cost projections for each power system, the results also show that the SPFC can become cost-competitive with conventional technologies, particularly for services characterized by high levels of utilization