A model for moisture and heat flow in fractured unsaturated porous media

The present study investigates propagation of a cohesive crack in non‐isothermal unsaturated porous medium under mode I conditions. Basic points of skeleton deformation, moisture, and heat transfer for unsaturated porous medium are presented. Boundary conditions on the crack surface that consist of mechanical interaction of the crack and the porous medium, water, and heat flows through the crack are taken into consideration. For spatial discretization, the extended finite element method is used. This method uses enriched shape functions in addition to ordinary shape functions for approximation of displacement, pressure, and temperature fields. The Heaviside step function and the distance function are exploited as enrichment functions for representing the crack surfaces displacement and the discontinuous vertical gradients of the pressure and temperature fields along the crack, respectively. For temporal discretization, backward finite difference scheme is applied. Problems solved from the literature show the validity of the model as well as the dependency of structural response on the material properties and loading. Copyright © 2016 John Wiley & Sons, Ltd.     

Synthesis of soil‐hydraulic properties and infiltration timescales in wildfire‐affected soils

We collected soil‐hydraulic property data from the literature for wildfire‐affected soils, ash, and unburned soils. These data were used to calculate metrics and timescales of hydrologic response related to infiltration and surface runoff generation. Sorptivity (S) and wetting front potential (Ψ_f) were significantly different (lower) in burned soils compared with unburned soils, whereas field‐saturated hydraulic conductivity (K_fs) was not significantly different. The magnitude and duration of the influence of capillarity during infiltration was greatly reduced in burned soils, causing faster ponding times in response to rainfall. Ash had large values of S and K_fs but moderate values of Ψ_f, compared with unburned and burned soils, indicating ash has long ponding times in response to rainfall. The ratio of S²/K_fs was nearly constant (~100 mm) for unburned soils but more variable in burned soils, suggesting that unburned soils have a balance between gravity and capillarity contributions to infiltration that may depend on soil organic matter, whereas in burned soils the gravity contribution to infiltration is greater. Changes in S and K_fs in burned soils act synergistically to reduce infiltration and accelerate and amplify surface runoff generation. Synthesis of these findings identifies three key areas for future research. First, short timescales of capillary influences on infiltration indicate the need for better measurements of infiltration at times less than 1 min to accurately characterize S in burned soils. Second, using parameter values, such as Ψ_f, from unburned areas could produce substantial errors in hydrologic modeling when used without adjustment for wildfire effects, causing parameter compensation and resulting underestimation of K_fs. Third, more thorough measurement campaigns that capture soil‐structural changes, organic matter impacts, quantitative water repellency trends, and soil‐water content along with soil‐hydraulic properties could drive the development of better techniques for numerically simulating infiltration in burned areas.     

Hydrologic response of engineered media in living roofs and bioretention to large rainfalls: experiments and modeling

The hydrologic response of engineered media plays an important role in determining a stormwater control measure's ability to reduce runoff volume, flow rate, timing, and pollutant loads. Five engineered media, typical of living roof and bioretention stormwater control measures, were investigated in laboratory column experiments for their hydrologic responses to steady, large inflow rates. The inflow, medium water content response, and outflow were all measured. The water flow mechanism (uniform flow vs. preferential flow) was investigated by analyzing medium water content response in terms of timing, magnitude, and sequence with depth. Modeling the hydrologic process was conducted in the HYDRUS‐1D software, applying the Richards equation for uniform flow modeling, and a mobile–immobile model for preferential flow modeling. Uniform flow existed in most cases, including all initially dry living roof media with bimodal pore size distributions and one bioretention medium with unimodal pore size distribution. The Richards equation can predict the outflow hydrograph reasonably well for uniform flow conditions when medium hydraulic properties are adequately represented by appropriate functions. Preferential flow was found in two media with bimodal pore size distributions. The occurrence of preferential flow is more likely due to the interaction between the bimodal pore structure and the initial water content rather than the large inflow rate.     

The effective stress concept and the evaluation of changes in pore air pressure under jacketed isotropic compression tests for dry sand based on the 2‐phase porous theory

The effective stress concept for solid‐fluid 2‐phase media was revisited in this work. In particular, the effects of the compressibility of both the pore fluid and the soil particles were studied under 3 different conditions, i.e., undrained, drained, and unjacketed conditions based on a Biot‐type theory for 2‐phase porous media. It was confirmed that Terzaghi effective stress holds at the moment when soil grains are assumed to be incompressible and when the compressibility of the pore fluid is small enough compared to that of the soil skeleton. Then, isotropic compression tests for dry sand under undrained conditions were conducted within the triaxial apparatus in which the changes in the pore air pressure could be measured. The ratio of the increment in the cell pressure to the increment in the pore air pressure, m, corresponds to the inverse of the B value by Bishop and was obtained during the step loading of the cell pressure. In addition, the m values were evaluated by comparing them with theoretically obtained values based on the solid‐fluid 2‐phase mixture theory. The experimental m values were close to the theoretical values, as they were in the range of approximately 40 to 185, depending on the cell pressure. Finally, it was found that the soil material with a highly compressible pore fluid, such as air, must be analyzed with the multi‐phase porous mixture theory. However, Terzaghi effective stress is practically applicable when the compressibilities of both the soil particles and the pore fluid are small enough compared to that of the soil skeleton.     

Why asymptotic limit of the acid dissolution capacity can lead to a sharp dissolution front in chemical dissolution of porous rocks?

The use of the asymptotic limit can greatly simplify the theoretical analysis of chemical dissolution front instabilities in fluid‐saturated rocks and therefore make it possible to obtain mathematical solutions, which often play a crucial role in understanding the propagation behavior of chemical dissolution fronts in chemical dissolution systems. However, there has been a debate in recent years that the asymptotic limit of the acid dissolution capacity (i.e., the acid dissolution capacity number approaching zero) alone cannot lead to a sharp dissolution front of the Stefan type in the acidization dissolution system, in which the dissolvable minerals of carbonate rocks are chemically dissolved by the injected acid flow. The acid dissolution capacity number is commonly defined as the ratio of the volume of the carbonate rock dissolved by an acid to that of the acid. In this paper, we use four different proof methods, including (i) direct use of the fundamental concepts; (ii) use of the mathematical governing equations of an acidization dissolution system; (iii) use of the different time scaling approach; and (iv) use of a moving coordinate system approach, to demonstrate that the asymptotic limit of the acid dissolution capacity can indeed lead to sharp dissolution fronts of the Stefan type in acidization dissolution systems on a much larger time scale (than the dissolution time scale). Our new finding is that on the reaction time scale, the condition of the conventional time derivative of porosity approaching zero alone can ensure that the acidization dissolution front has a sharp shape of the Stefan type. Copyright © 2017 John Wiley & Sons, Ltd.     

Torsional dynamic response of a large‐diameter pipe pile in viscoelastic saturated soil

An analytical solution is developed in this paper to investigate the dynamic response of a large‐diameter end‐bearing pipe pile subjected to torsional loading in viscoelastic saturated soil. The wave propagation in saturated soil and pile are simulated by Biot's two‐phased linear theory and one‐dimensional elastic theory, respectively. The dynamic equilibrium equations of the outer soil, inner soil, and pile are established. The solutions for the outer and inner soils in frequency domain are obtained by Laplace transform technique and the separation of variables method. Then, the dynamic response of the pile is obtained on the basis of the perfect contacts between the pile and the outer soil as well as the inner soil. The results in this paper are compared with that of a solid pile in elastic saturated soil to verify the validity of the solution. Furthermore, the solution in this paper is compared with the classic plane strain solution to verify the solution further and check the accuracy of the plane strain solution. Numerical results are presented to analyze the vibration characteristics and illustrate the effect of the soil parameters and the geometry size of the pile on the complex impedance and velocity admittance of the pile head. Finally, the displacement of the soil at different depth and frequency is analyzed. Copyright © 2014 John Wiley & Sons, Ltd.     

HYDROGEOLOGY & ENGINEERING GEOLOGY

正1.HYDROGEOLOGY20142452Chang Chengcao(Institute of Coalfield Geological Survey and Design of Jilin Province,Changchun 130062,China)Hydrogeological Type Division of Coal Mine No.2Mine in Yitong Manchu Autonomous County,Jilin Province(Jilin Geology,ISSN1001-2427,CN22-1099/P,32(4),2013,p.129-131,3refs.)     

Seasonal variations in infiltrability of moss‐dominated biocrusts on aeolian sand and loess soil in the Chinese Loess Plateau

Biocrust effects on soil infiltration have attracted increasing attention in dryland ecosystems, but their seasonal variations in infiltrability have not yet been well understood. On the Chinese Loess Plateau, soil infiltrability indicated by saturated hydraulic conductivity (K_s) of biocrusts and bare soil, both on aeolian sand and loess soil, was determined by disc infiltrometer in late spring (SPR), midsummer (SUM), and early fall (FAL). Then their correlations with soil biological and physiochemical properties and water repellency index (RI) were analysed. The results showed that the biocrusts significantly decreased K_s both on sand during SPR, SUM, and FAL (by 43%, 66%, and 35%, respectively; P < .05) and on loess (by 42%, 92%, and 10%, respectively; P <.05). As compared with the bare soil, the decreased K_s in the biocrusted surfaces was mostly attributed to the microorganism biomass and also to the increasing content of fine particles and organic matter. Most importantly, both the biocrusts and bare soil exhibited significant (F ≥ 11.89, P ≤ .003) seasonal variations in K_s, but their patterns were quite different. Specifically, the K_s of bare soil gradually decreased from SPR to SUM (32% and 42% for sand and loess, respectively) and FAL (29% and 39%); the K_s of biocrusts also decreased from SPR to SUM (59% and 92%) but then increased in FAL (36% and 588%). Whereas the seasonal variations in K_s of the biocrusts were closely correlated with the seasonal variations in RI, the RI values were not high enough to point at hydrophobicity. Instead of that, the seasonal variations of K_s were principally explained by the changes in the crust biomass and possibly by the microbial exopolysaccharides. We conclude that the biocrusts significantly decreased soil infiltrability and exhibited a different seasonal variation pattern, which should be carefully considered in future analyses of hydropedological processes.     

The exponential decline in saturated hydraulic conductivity with depth: a novel method for exploring its effect on water flow paths and transit time distribution

The strong vertical gradient in soil and subsoil saturated hydraulic conductivity is characteristic feature of the hydrology of catchments. Despite the potential importance of these strong gradients, they have proven difficult to model using robust physically based schemes. This has hampered the testing of hypotheses about the implications of such vertical gradients for subsurface flow paths, residence times and transit time distribution. Here we present a general semi‐analytical solution for the simulation of 2D steady‐state saturated‐unsaturated flow in hillslopes with saturated hydraulic conductivity that declines exponentially with depth. The grid‐free solution satisfies mass balance exactly over the entire saturated and unsaturated zones. The new method provides continuous solutions for head, flow and velocity in both saturated and unsaturated zones without any interpolation process as is common in discrete numerical schemes. This solution efficiently generates flow pathlines and transit time distributions in hillslopes with the assumption of depth‐varying saturated hydraulic conductivity. The model outputs reveal the pronounced effect that changing the strength of the exponential decline in saturated hydraulic conductivity has on the flow pathlines, residence time and transit time distribution. This new steady‐state model may be useful to others for posing hypotheses about how different depth functions for hydraulic conductivity influence catchment hydrological response. Copyright © 2016 John Wiley & Sons, Ltd.     

ISA model: A constitutive model for soils with yield surface in the intergranular strain space

In this article, a new constitutive model for soils is proposed. It is formulated by means of plasticity, but in contrast to the precedent works, it presents a yield function describing a surface within the intergranular strain space. This latter is a state variable providing information of the recent strain history. An expression for the plastic strain rate has been proposed to guarantee the stress rate continuity. Under the application of medium or large strain amplitudes, the constitutive equation becomes independent of the intergranular strain and delivers a mathematical structure similar to some Karlsruhe hypoplastic models. Some simulations of monotonic and cyclic triaxial test are provided to evaluate and analyze the model performance. Copyright © 2015 John Wiley & Sons, Ltd.