Cavity expansion analyses of crushable granular materials with state‐dependent dilatancy

Crushability is one of the important behaviors of granular materials particularly under high stress states, and affects both the deformability and strength of the materials that are in essence associated with state‐dependent dilatancy. In this presentation, first, a new critical state model is proposed to take into account the three different modes of compressive deformation of crushable granular materials, i.e. particle rearrangement, particle crushing and pseudo‐elastic deformation. Second, the governing equations for cavity expansion in crushable granulates are introduced, in which the state‐dependent dilatancy as well as the bounding surface plasticity model are used. Then, the procedure to obtain semi‐analytical solutions to cavity expansion in the material is described in detail, in which a commercial differential equation solver is employed. Finally, cavity expansion analyses are carried out on Toyoura sand, a well‐documented granular material, to demonstrate the effects of crushability and state‐dependent dilatancy. The study shows that particle crushing does occur at both high stress and critical states and affects the stress fields and the deformation behavior of the material surrounding the cavity in association with state‐dependent dilatancy. This leads to conclusion that particle crushing and state‐dependent dilatancy have to be taken into account when cavity expansion theory is used to interpret cone penetration tests and pressuremeter tests. Copyright © 2011 John Wiley & Sons, Ltd.     

Consolidation in spatially random unsaturated soils based on coupled flow‐deformation simulation

This paper integrates random field simulation of soil spatial variability with numerical modeling of coupled flow and deformation to investigate consolidation in spatially random unsaturated soil. The spatial variability of soil properties is simulated using the covariance matrix decomposition method. The random soil properties are imported into an interactive multiphysics software COMSOL to solve the governing partial differential equations. The effects of the spatial variability of Young's modulus and saturated permeability together with unsaturated hydraulic parameters on the dissipation of excess pore water pressure and settlement are investigated using an example of consolidation in a saturated‐unsaturated soil column because of loading. It is found that the surface settlement and the pore water pressure profile during the process of consolidation are significantly affected by the spatially varying Young's modulus. The mean value of the settlement of the spatially random soil is more than 100% greater than that of the deterministic case, and the surface settlement is subject to large uncertainty, which implies that consolidation settlement is difficult to predict accurately based on the conventional deterministic approach. The uncertainty of the settlement increases with the scale of fluctuation because of the averaging effect of spatial variability. The effects of spatial variability of saturated permeability k_sat and air entry parameters are much less significant than that of elastic modulus. The spatial variability of air entry value parameters affects the uncertainties of settlement and excess pore pressure mostly in the unsaturated zone. Copyright © 2016 John Wiley & Sons, Ltd.     

A coupled discrete element model for the simulation of soil and water flow through an orifice

Soil erosion around defective underground pipes can cause ground collapses and sinkholes in urban areas. Most of these soil erosion events are caused by fluidization of the surrounding soil with subsequent washing into defective sewer pipes. In this study, this soil erosion process is simplified as the gradual washout of sand particles mixed with water through an orifice. The discrete element method is used to simulate the large deformation behavior of the sand particles, and the Darcy fluid model is coupled with this approach to simulate fluid flow through porous sand media. A coupled 3D discrete element model is developed and implemented based on this scheme. To simulate previous experiments using this coupled model considering the current computing capacity, we incorporated a ‘supply layer’ to study the continuous erosion process. The coupled model can predict the erosion flow rates of sand and water and the shape of erosion void. Thus, the model can be used as an effective and efficient tool to investigate the soil erosion process around defective pipes. Copyright © 2017 John Wiley & Sons, Ltd.     

Elastic fields in two joined transversely isotropic media of infinite extent as a result of rectangular loading

This paper presents the closed‐form solutions for the elastic fields in two bonded rocks induced by rectangular loadings. Each of the two bonded rocks behaves as a transversely isotropic linear elastic solid of semi‐infinite extent. They are completely bonded together at a horizontal surface. The rectangular loadings are body forces along either vertical or horizontal directions and are uniformly applied on a rectangular area. The rectangular area is embedded in the two bonded rocks and is parallel to the horizontal interface. The classical integral transforms are used in the solution formulation, and the elastic solutions are expressed in the forms of elementary harmonic functions for the rectangular loadings. The stresses and displacements in the rocks induced by both the horizontal and vertical body forces are also presented. The numerical results illustrate the important effect of the anisotropic bimaterial properties on the stress and displacement fields. The solutions can be easily implemented for numerical calculations and applied to problems encountered in rock mechanics and engineering. Copyright © 2011 John Wiley & Sons, Ltd.     

Modeling three‐dimensional hydraulic fracture propagation using virtual multidimensional internal bonds

Propagation of fractures, especially those emanating from wellbores and closed natural fractures, often involves Mode I and Mode II, and at times Mode III, posing significant challenges to its numerical simulation. When an embedded inclined fracture is subjected to compression, the fracture edge is constrained by the surrounding materials so that its true propagation pattern cannot be simulated by 2D models. In this article, a virtual multidimensional internal bond (VMIB) model is presented to simulate three‐dimensional (3D) fracture propagation. The VMIB model bridges the processes of macro fracture and micro bond rupture. The macro 3D constitutive relation in VMIB is derived from the 1D bond in the micro scale and is implemented in a 3D finite element method. To represent the contact and friction between fracture surfaces, a 3D element partition method is employed. The model is applied to simulate fracture propagation and coalescence in typical laboratory experiments and is used to analyze the propagation of an embedded fracture. Simulation results for single and multiple fractures illustrate 3D features of the tensile and compressive fracture propagation, especially the propagation of a Mode III fracture. The results match well with the experimental observation, suggesting that the presented method can capture the main features of 3D fracture propagation and coalescence. Moreover, by developing an algorithm for applying pressure on the fracture surfaces, propagation of a natural fracture is also simulated. The result illustrates an interesting and important phenomenon of Mode III fracture propagation, namely the fracture front segmentation. Copyright © 2012 John Wiley & Sons, Ltd.     

The environmental impacts of pollutants generated by routine shipping operations on ports

The environmental impacts generated by shipping operations have increasingly become an important research topic, where its pollutants often pose negative externalities to natural habitats and economic losses to coastal areas. While the environmental impact costs generated by shipping disasters, notably large scale accidental oil spills, have been widely studied, hitherto, works dedicated to the assessment of environmental impact costs of pollutants generated by routine shipping operations remain scarce due to their relatively implicit nature and thus delays of consequential risks. Hence, by proposing an economic model and calibrating it to Port of Rotterdam, this paper assesses the environmental impact costs generated by routine shipping operations on ports. By shedding light on this important, but under-researched, issue critical to the well-being of the global shipping industry, this paper provides a decent framework for further research in sustainable shipping and port management for future generations.     

Modelling and observation of hedgerow transpiration effect on water balance components at the hillslope scale in Brittany

Vegetation has a major influence on the water and energy balance of the earth's surface. In the last century, human activities have modified land use, inducing a consequent change in albedo and potential evapotranspiration. Linear vegetation structures (hedgerows, shelterbelts, open woodland, etc) were particularly abundant but have declined considerably over the past several decades. In this context, it is important to quantify their effect on water and energy balance both on a global scale (climate change and weather prediction) and on a local scale (soil column, hillslope and watershed). The main objective of this study was to quantify the effect of hedgerows on the water cycle by evaluating spatial and temporal variations of water balance components of a hillslope crossed by a hedgerow. Water flow simulation was performed using Hydrus‐2D to emphasize the importance of transpiration in the water balance and to evaluate water extraction from groundwater. Model validation was performed by comparing simulated and observed soil matrix potentials and groundwater levels. Hedgerow transpiration was calculated from sap flow measurements of four trees. Water balance components calculated with a one‐dimensional water balance equation were compared with simulations. Simulation runs with and without tree root uptake underlined the effect of hedgerow transpiration, increasing capillary rise and decreasing drainage. Results demonstrated that the spatial and temporal variability of water balance components was related to the hedgerow presence as well as to the meteorological context. The relations between transpiration, groundwater proximity and soil‐water availability determined the way in which water balance components were affected. Increased capillary rise and decreased drainage near hedges were related to the high transpiration of trees identified in this study. Transpiration reached twice the potential evapotranspiration when groundwater level and precipitation amounts were high. Water balance analysis showed that transpiration was a substantial component, representing 40% of total water output. These results may offer support for improving hydrological models by including the effect of land use and land cover on hydrological processes. Copyright © 2012 John Wiley & Sons, Ltd.     

Nb–Ta fractionation induced by fluid‐rock interaction in subduction‐zones: constraints from UHP eclogite‐ and vein‐hosted rutile from the Dabie orogen,Central‐Eastern China

Niobium and Ta concentrations in ultrahigh‐pressure (UHP) eclogites and rutile from these eclogites and associated high pressure (HP) veins were used to study the behaviour of Nb–Ta during dehydration and fluid‐rock interaction. Samples were collected through a ～2 km profile at the Bixiling complex in the Dabie orogenic belt, Central‐Eastern China. All but one eclogite away from veins (EAVs) display nearly constant Nb/Ta ratios ranging from 16.1 to 19.2, with an average of 16.9 ± 0.8 (2 SE), similar to that of their gabbroic protolith from the Yangtze Block. Nb/Ta ratios of rutile from the EAVs range from 12.7 to 25.3 among different individual grains, with the average values close to those of the corresponding bulk rocks. These observations show that Nb and Ta were not significantly fractionated by prograde metamorphism up to eclogite facies when no significant fluid‐rock interaction occurs. In contrast, Nb/Ta ratios of rutile from eclogites close to veins (ECVs) are highly variable from 17.8 to 49.8, which are systematically higher (by up to 17) than those of rutile from the veins. These observations demonstrate that Nb and Ta were mobilized and fractionated during localized fluid flow and intensive fluid‐rock interaction. This is strongly supported by Nb/Ta zoning patterns in single rutile grains revealed by in situ LA‐ICP‐MS analysis. Ratios of Nb/Ta in the ECV‐hosted rutile decrease gradually from cores towards rims, whereas those in the EAV‐hosted rutile are nearly invariable. Furthermore, the vein rutile shows Nb/Ta zoning patterns that are complementary to those in rutile from their immediate hosts (ECVs), suggesting an internal origin for the vein‐forming fluids. The Nb/Ta ratios of such fluids evolved from low values at the early stage of subduction to higher values at later supercritical conditions with increased temperature and pressure. Quantitative modelling was conducted to constrain the compositional evolution of metamorphic fluids during dehydration and fluid‐rock interaction focusing on Nb–Ta distribution. The modelling results based on our proposed multistage fluid phase evolution path can essentially reproduce the natural observations reported in the present study.