Numerical evaluation of soft inter‐slab joint in concrete‐faced rockfill dam with dual mortar finite element method

Recently constructed concrete‐faced rockfill dams (CFRDs) often use soft inter‐slab joints to prevent axial compression‐induced extrusion damage in the concrete face. Due to the complexity of the multibody contact and the lack of information on the actual behavior of soft joints, it is highly challenging to numerically assess the effect of soft joints in CFRDs. In this paper, we present a numerical approach for the three‐dimensional modeling of CFRDs with hard and soft joints. A dual mortar finite element method with Lagrange multiplier is developed to treat the multibody contact in hard joints with impenetrability condition. The soft joint slab‐filler‐slab contact system is modeled using an equivalent contact interface approach, where the soft contact constraints are imposed using a perturbed Lagrange formulation. Through a series of laboratory tests, the mechanical behavior of soft joint is investigated. An extrusion model for the soft joint is presented and implemented in the dual mortar finite element method. The proposed numerical method is applied to the three‐dimensional analysis of Tianshengqiao‐1 CFRD. Despite the complex multibody contact and strong material and geometry nonlinearities in the CFRD, the proposed method is stable and capable of capturing salient characteristics of the CFRD. Numerical results show that in Tianshengqiao‐1, the employment of soft joints can effectively reduce the axial compression stress, thus greatly alleviating the risk of extrusion damage in the concrete face.     

The role of fog,orography, and seasonality on precipitation in a semiarid,tropical island

Isotopes of water (²H/¹H and ¹⁸O/¹⁶O) are commonly used to trace hydrological processes such as moisture recycling, evaporation loss, and moisture source region and often vary temporally in a given region. This study provides a first‐ever characterization of temporally variable precipitation mechanisms of San Cristóbal Island, Galápagos. We collected fog, rain, and throughfall samples over three field seasons to understand the mechanisms driving seasonal‐ and event‐based variability in the isotopic composition of precipitation in Galápagos. We establish that fog is a common phenomenon in San Cristóbal, especially during the dry season, and we found that fog, compared with cocollected rainfall, is consistently enriched. We further suggest that the relative contribution of fog formed via different mechanisms (orographic, advective, radiation) varied seasonally. We found that the source region is the most dominant control of the isotopic composition of rainfall in the Galápagos at both the seasonal and event scales, but subcloud evaporative processes (the nontraditional manifestation of the amount effect) became a dominant control on the isotopic composition of rainfall during the dry season. Overall, our findings suggest that understanding seasonally variable water‐generating mechanisms is required for effective water resource management in San Cristóbal Island and other semiarid island ecosystems under current and future regimes of climate change.     

Changes in run‐off and sediment load in the three parts of the Yellow River basin,in response to climate change and human activities

Hydrological regimes in the Yellow River have changed significantly because of climate change and intensive human interventions. These changes present severe challenges to water resource utilization and ecological development. Variation of run‐off, suspended sediment load (SSL), and eight precipitation indices (P1: 0–12 mm·day^?1, P12: 12–25 mm·day^?1, P25: 25–50 mm·day^?1, P50: P ≥ 50 mm·day^?1 and corresponding rainfall day: Pd1, Pd12, Pd25, Pd50 day year^?1) in three critical parts of the Yellow River basin (source region: SRYRB, upper reaches: URYRB, middle reaches: MRYRB) were investigated for the period from 1960 to 2015. The results show that run‐off and SSL significantly decreased (P < 0.01) in the URYRB and the MRYRB, whereas their decline in the SRYRB was insignificant (P > 0.05). Moreover, run‐off in the URYRB had one change point in 1987, and SSL in the URYRB as well as run‐off and SSL in the MRYRB had two change points (in the 1970s and the 1990s). Over the same period, only Pd1 and Pd12 in the SRYRB showed significant increasing trends, and an abrupt change appeared in 1981. The optimal precipitation indices for assessing the effects of precipitation on run‐off and SSL in the URYRB and MRYRB were Pd50 and P12, respectively. A double‐mass curve analysis showed that precipitation and human activities contributed to approximately 20% and 80% of the reduction in run‐off, respectively, for both the SRYRB and the MRYRB. However, the contribution rate of precipitation and human activities on SSL reduction was approximately 40% and 60% in the URYRB and 5% and 95% in the MRYRB, respectively. Human activities, primarily soil and water conservation measures and water extraction (diversion), were the main factors (>50%) that reduced the run‐off. However, the dominant driving factors for SSL reduction were soil and water conservation measures and reservoir interception, for which the contribution rate was higher than 70% in the MRYRB. This work strengthens the understanding of hydrological responses to precipitation change and provides a useful reference for regional water resource utilization.     

Rill flow resistance law under equilibrium bed‐load transport conditions

In this paper, a recently deduced flow resistance equation for open channel flow was tested under equilibrium bed‐load transport conditions in a rill. First, the flow resistance equation was deduced applying dimensional analysis and the incomplete self‐similarity condition for the flow velocity distribution. Then, the following steps were carried out for developing the analysis: (a) a relationship (Equation  13 ) between the Γ function of the velocity profile, the rill slope, and the Froude number was calibrated by the available measurements by Jiang et al.; (b) a relationship (Equation  17 ) between the Γ function, the rill slope, the Shields number, and the Froude number was calibrated by the same measurements; and (c) the Darcy–Weisbach friction factor values measured by Jiang et al. were compared with those calculated by the rill flow resistance equation with Γ estimated by Equations  13 and 17 . This last comparison demonstrated that the rill flow resistance equation, in which slope and Shields number, representative of sediment transport effects, are introduced, is characterized by the lowest values of the estimate errors.     

Responses of chemically active and naturally fractured shale under time‐dependent mechanical loading and ionic solution exposure

In this paper, the analytical dual‐porosity dual‐permeability poromechanics solution for saturated cylinders is extended to account for electrokinetic effects and material transverse isotropy, which simulate the responses of chemically active naturally fractured shale under time‐dependent mechanical loading and ionic solution exposure. The solution addresses the stresses, fracture pore pressure, matrix pore pressure, fluid fluxes, ion concentration evolution, and displacements due to the applied stress, pore pressure, and solute concentration difference between the sample and the circulation fluid. The presented solution will not only help validate numerical simulations but also assist in calibrating and interpreting laboratory results on dual‐porosity dual‐permeability shale. It is recommended that the analytical solutions of radial and axial displacements be used to match the corresponding laboratory‐recorded data to determine shale dual permeability and chemo‐electrical parameters including membrane coefficient, ions diffusion coefficients, and electro‐osmotic permeability.     

Stability analysis of shotcrete supported crown of NATM tunnels with discontinuity layout optimization

For tunnel constructed by New Austrian Tunnelling Method, the crown is the upper part of tunnel section, constructed during excavation process and supported by shotcrete. The stability of the crown has great influence on the safety of tunnel itself and the buildings above, which correlates, among others, with geometrical setup of tunnel and material properties of shotcrete and soil/rock. In this paper, aiming at analyzing the stability of shotcrete supported crown, a recently presented numerical method discontinuity layout optimization is adopted, which introduces a great amount of potential discontinuities cross over one another and provides a wide search space for efficient upper limit analysis. In the analysis, a well‐established hydration model of cementitious material is implemented for accounting the hydration of shotcrete. Then assumptions based on convergence‐confinement method are used for accounting the 3‐dimensional effect in 2‐dimensional analysis, finally providing time‐space–dependent assessments of stability of shotcrete supported crown.     

Comprehensive study of the drying behavior of Boom clay: Experimental investigation and numerical modeling

This paper presents a thermo‐hydro‐mechanical framework to model the drying behavior of Boom clay. First, the experimental campaign conducted Noémie Prime is briefly presented because it is used to validate the model. The data acquisition and processing is emphasized because of the use of X‐ray microtomography to be able to more accurately compare experimental and numerical strain fields. The different submodels are introduced. Numerical simulations are performed to illustrate the capability of the proposed model to reproduce the observed behavior. Finally, a comprehensive sensitivity study on several key model parameters associated with the water retention curve, and the permeability of the medium, is performed to get a better understanding of the physics behind the coupled model.     

Assessment of 2DH and pseudo‐3D modelling platforms in a large saline aquatic system: Lake Urmia,Iran

The main objective of this paper is to provide comparative quantitative examinations on the capabilities of two‐dimensional horizontal and pseudo‐three‐dimensional (3D) modelling approaches for simulating spatial and temporal variability of the flow and salinity in Lake Urmia, Iran. The water quality in the lake has been an environmentally important subject partly because this shallow hypersaline aquatic ecosystem is considered to be one of the largest natural habitats of a unique multicellular organism, Artemia urmiana. This brine shrimp is the major food source for many of the protected and rare shorebirds that visit the lake. A. urmiana can grow and survive in certain ranges of salinity, and their disappearance could lead to an alteration of existing equilibria. The lake has also experienced considerable man‐made changes during the past three decades. A newly built crossing embankment almost divided the lake into two northern and southern halves. A relatively small opening of 1.25 km in the new embankment provides water connections between the two halves. As a result, the flow and salinity regimes have been significantly changed. This might have had adverse serious impacts on the lake ecosystem. In the current study, the two‐dimensional horizontal hydrodynamic model has been found to provide reasonable predictions for the flow regime in the lake, whereas its salinity predictions have not been consistent with the field observations. The pseudo‐3D model has produced results fairly close to the salinity measurements and its temporal and spatial variations. The pseudo‐3D model has been used for evaluating the embankment effects on the lake hydrodynamics and on the salinity conditions. The effectiveness of introducing a different number or length of openings in the embankment for restoring the pre‐embankment conditions has also been examined. These remedy options have been found not to offer substantial improvements to the lake's existing ecosystem. Copyright © 2013 John Wiley & Sons, Ltd.     

Modelling the internal stress field in argillaceous rocks under humidification/desiccation

This paper deals with analytical and numerical modelling of the internal stress generated in argillaceous rocks during humidification/desiccation processes, which is an essential issue for damage study. This local stress field arises from two mechanisms: (i) complex interactions between free swelling/shrinking clay matrix and non‐strained inclusions of carbonate and quartz and (ii) a self‐restraint effect induced by the moisture gradient during the transient moisture exchange process. The inclusion–matrix interaction is investigated in different cases. Firstly, the analytical solution of the stress around a cylindrical inclusion embedded in an infinite swelling matrix is derived: The inclusion would suffer tension (compression) under humidification (desiccation), and the resulting cracking patterns are discussed. Then, the problem of two inclusions with different distances in an infinite swelling matrix is considered, and it is shown that the local stress around an inclusion will be perturbed and amplified by neighbouring inclusions. Finally, an inclusion outcropping at the free surface of a swelling matrix is modelled as to investigate the effect of free surface: The inclusion–matrix interface undergoes shear stresses of which the maximum is found at the free surface. In addition to the inclusion–matrix interaction, the self‐restraint effect is investigated: The induced stress is maximal at the beginning of humidification/desiccation processes and vanishes gradually with time. The quantity of the self‐restraint stress is strongly controlled by the hydric loading rate. Copyright © 2014 John Wiley & Sons, Ltd.     

Study of the interaction between dry granular flows and rigid barriers with an SPH model

This study uses an incompressible smoothed‐particle hydrodynamics model to investigate the interaction between dry granular material flows and rigid barriers. The primary aim is to summarise some practical guidelines for the design of debris‐resisting barriers. The granular materials are modelled as a rigid‐perfectly plastic material where the plastic flow corresponds to the critical state. The coupled continuity equation and momentum equation are solved by a semi‐implicit algorithm. Compared with flows in controlled flume experiments, the model adequately reproduces both the kinetic of the flows and the impact force under various conditions. Then the numerical simulations are used to study the detailed interaction process. It is illustrated quantitatively that the interaction force consists of two parts, ie, the earth pressure force caused by the weight of the soil and a dynamic force caused by the internal deformation (flowing mass on top of a dead zone). For the estimation of impact load, this study suggests that an increased earth pressure coefficient depending on the Froude number should be incorporated into the hydrostatic model.