Low-order stabilized finite element for the full Biot formulation in soil mechanics at finite strain

This article presents a novel finite element formulation for the Biot equation using low-order elements. Additionally, an extra degree of freedom is introduced to treat the volumetric locking steaming from the effective response of the medium; its balance equation is also stabilized. The accuracy of the proposed formulation is demonstrated by means of numerical analyses.     

Implicit integration under mixed controls of a breakage model for unsaturated crushable soils

This paper discusses a series of stress point algorithms for a breakage model for unsaturated granular soils. Such model is characterized by highly nonlinear coupling terms introduced by breakage‐dependent hydro‐mechanical energy potentials. To integrate accurately and efficiently its constitutive equations, specific algorithms have been formulated using a backward Euler scheme. In particular, because implementation and verification of unsaturated soil models often require the use of mixed controls, the incorporation of various hydro‐mechanical conditions has been tackled. First, it is shown that the degree of saturation can be replaced with suction in the constitutive equations through a partial Legendre transformation of the energy potentials, thus changing the thermomechanical state variables and enabling a straightforward implementation of a different control mode. Then, to accommodate more complex control scenarios without redefining the energy potentials, a hybrid strategy has been used, combining the return mapping scheme with linearized constraints. It is shown that this linearization strategy guarantees similar levels of accuracy compared with a conventional strain–suction‐controlled implicit integration. In addition, it is shown that the use of linearized constraints offers the possibility to use the same framework to integrate a variety of control conditions (e.g., net stress and/or water‐content control). The convergence profiles indicate that both schemes preserve the advantages of implicit integration, that is, asymptotic quadratic convergence and unconditional stability. Finally, the performance of the two implicit schemes has been compared with that of an explicit algorithm with automatic sub‐stepping and error control, showing that for the selected breakage model, implicit integration leads to a significant reduction of the computational cost. Such features support the use of the proposed hybrid scheme also in other modeling contexts, especially when strongly nonlinear models have to be implemented and/or validated by using non‐standard hydro‐mechanical control conditions. Copyright © 2015 John Wiley & Sons, Ltd.     

Chemo‐mechanics of cemented granular solids subjected to precipitation and dissolution of mineral species

This paper studies the chemo‐mechanics of cemented granular solids in the context of continuum thermodynamics for fluid‐saturated porous media. For this purpose, an existing constitutive model formulated in the frame of the Breakage Mechanics theory is augmented to cope with reactive processes. Chemical state variables accounting for the reactions between the solid constituents and the solutes in the pore fluid are introduced to enrich the interactions among the microstructural units simulated by the model (i.e., grains and cement bonds). Two different reactive processes are studied (i.e., grain dissolution and cement precipitation), using the chemical variables to describe the progression of the reactions and track changes in the size of grains and bonds. Finally, a homogenization strategy is used to derive the energy potentials of the solid mixture, adopting probability density functions that depend on both mechanical and chemical indices. It is shown that the connection between the statistics of the micro‐scale attributes and the continuum properties of the solid enables the mathematical capture of numerous mechanical effects of lithification and chemical deterioration, such as changes in stiffness, expansion/contraction of the elastic domain, and development of inelastic strains during reaction. In particular, the model offers an interpretation of the plastic strains generated by aggressive environments, which are here interpreted as an outcome of chemically driven debonding and comminution. As a result, the model explains widely observed macroscopic signatures of geomaterial degradation by reconciling the energetics of the deformation/reaction processes with the evolving geometry of the microstructural attributes. Copyright © 2015 John Wiley & Sons, Ltd.     

Response of suspended sediment particle size distributions to changes in water chemistry at an Andean mountain stream confluence receiving arsenic rich acid drainage

Acid drainage is an important water quality issue in Andean watersheds, affecting the sustainability of urban, agricultural, and industrial activities. Mixing zones receiving acid drainage are critical sites where changes in pH and chemical environment promote the formation and dissolution of iron and aluminum oxy/hydroxides. These particles can significantly change the speciation of toxic metals and metalloids throughout drainage networks via sorption, desorption, and settling processes. However, little is known about the behavior of particle size distributions (PSDs) in streams affected by acid drainage and their relationship to metal speciation. This work studied: (a) the PSDs for a wide range of mixing ratios found at a fluvial confluence affected by acid drainage, and (b) the response of PSDs and arsenic speciation to environmental changes found when the particles approach complete mixing conditions. The confluence between the Azufre River (pH ~ 2, high concentration of dissolved metals) and Caracarani River (pH = 8.6, low concentration of dissolved metals) was used as a representative model for study. Field measurements show a bimodal PSD with modal diameters of ~50 and 300 μm. At shorter distances from the junction, the smaller modes with smaller particle volumes were dominant across the stream cross‐sections. A systematic shift towards larger particle sizes and larger particle volumes occurred downstream. The analysis of laboratory PSDs for Azufre/Caracarani mixing ratios between 0.01 and 0.5 (pH from 6.2 to 2.3) showed a bimodal trend with ~15 and 50 μm characteristic diameters; larger particles formed at pH>4. When particle suspensions were transferred in laboratory experiments from very low pH to full mixing conditions (pH ~ 2.8 and mixing ratio ~ 0.25) particle sizes varied, and the dissolved arsenic concentration decreased. The observed reaction kinetics were slow compared to the time scale of advective transport, creating opportunities for engineered controls for arsenic. This work contributes to a better understanding of the chemical‐hydrodynamic interactions in watersheds affected by mining, and identifying opportunities to improve water quality at points of use.     

A GLUE‐based uncertainty assessment framework for tritium‐inferred transit time estimations under baseflow conditions

The last decade has seen major technical and scientific improvements in the study of water transfer time through catchments. Nevertheless, it has been argued that most of these developments used conservative tracers that may disregard the oldest component of water transfer, which often has transit times greater than 5 years. Indeed, although the analytical reproducibility of tracers limits the detection of the older flow components associated with the most dampened seasonal fluctuations, this is very rarely taken into account in modelling applications. Tritium is the only environmental tracer at hand to investigate transfer times in the 5‐ to 50‐year range in surface waters, as dissolved gases are not suitable due to the degassing process. Water dating with tritium has often been difficult because of the complex history of its atmospheric concentration, but its current stabilization together with recent analytical improvements open promising perspectives. In this context, the innovative contribution of this study lies in the development of a generalized likelihood uncertainty estimation‐based approach for analysing the uncertainties associated with the modelling of transit time due to both parameter identification and tracer analytical precision issues. A coupled resampling procedure allows assessment of the statistical significance of the transfer time differences found in diverse waters. This approach was developed for tritium and the exponential‐piston model but can be implemented for virtually any tracer and model. Stream baseflow, spring and shallow aquifer waters from the Vallcebre research catchments, analysed for tritium in different years with different analytical precisions, were investigated by using this approach and taking into account other sources of uncertainty. The results showed three groups of waters of different mean transit times, with all the stream baseflow and spring waters older than the 5‐year threshold needing tritium. Low sensitivity of the results to the model structure was also demonstrated. Dual solutions were found for the waters sampled in 2013, but these results may be disambiguated when additional analyses will be made in a few years. Copyright © 2016 John Wiley & Sons, Ltd.