Refined modeling and movement characteristics analyses of irregularly shaped particles

Irregularly shaped (IRS) particles widely exist in many engineering and industrial fields. The macro physical and mechanical properties of the particle system are governed by the interaction between the particles in the system. The interaction between IRS particles is more complicated because of their complex geometric shape with extremely irregular and co‐existed concave and convex surfaces. These particles may interlock each other, making the sliding and friction of IRS particles more complex than that of particles with regular shape. In order to study the interaction of IRS particles more efficiently, a refined method of constructing discrete element model based on computed tomography scanning of IRS particles is proposed. Three parameters were introduced to control the accuracy and the number of packing spheres. Subsequently, the inertia tensor of the IRS particle model was optimized. Finally, laboratory and numerical open bottom cylinder tests were carried out to verify the refined modeling method. The influence of particle shape, particle position, and mesoscopic friction coefficient on the interaction of particles was also simulated. It is noteworthy that with the increase of mesoscopic friction coefficient, the fluidity of IRS particle assembly decreases, and intermittent limit equilibrium state may appear. Copyright © 2014 John Wiley & Sons, Ltd.     

Modelling surface runoff to evaluate the effects of wildfires in multiple semi‐arid,shrubland‐dominated catchments

Wildfires change the infiltration properties of soil, reduce the amount of interception and result in increased runoff. A wildfire at Northeast Attica, Central Greece, in August 2009, destroyed approximately one third of a study area consisting of a mixture of shrublands, pastures and pines. The present study simultaneously models multiple semi‐arid, shrubland‐dominated Mediterranean catchments and assesses the hydrological response (mean annual and monthly runoff and runoff coefficients) during the first few years following wildfires. A physically based, hydrological model (MIKE SHE) was chosen. Calibration and validation results of mean monthly discharge presented very good agreement with the observed data for the pre‐wildfire and post‐wildfire period for two subcatchments (Nash–Sutcliffe Efficiency coefficient of 79.7%). The model was then used to assess the pre‐wildfire and post‐wildfire runoff responses for each of seven catchments in the study area. Mean annual surface runoff increased for the first year and after the second year following the wildfires increased by 112% and 166%, respectively. These values are within the range observed in similar cases of monitored sites. This modelling approach may provide a way of prioritizing catchment selection with respect to post‐fire remediation activities. Additionally, this modelling assessment methodology would be valuable to other semi‐arid areas because it provides an important means for comprehensively assessing post‐wildfire response over large regions and therefore attempts to address some of the scaled issues in the specific literature field of research. Copyright © 2015 John Wiley & Sons, Ltd.     

Discrete element method modeling of inherently anisotropic rocks under uniaxial compression loading

A new numerical approach is proposed in this study to model the mechanical behaviors of inherently anisotropic rocks in which the rock matrix is represented as bonded particle model, and the intrinsic anisotropy is imposed by replacing any parallel bonds dipping within a certain angle range with smooth‐joint contacts. A series of numerical models with β = 0°, 15°, 30°, 45°, 60°, 75°, and 90° are constructed and tested (β is defined as the angle between the normal of weak layers and the maximum principal stress direction). The effect of smooth‐joint parameters on the uniaxial compression strength and Young's modulus is investigated systematically. The simulation results reveal that the normal strength of smooth‐joint mainly affects the behaviors at high anisotropy angles (β > 45°), while the shear strength plays an important role at medium anisotropy angles (30°–75°). The normal stiffness controls the mechanical behaviors at low anisotropy angles. The angle range of parallel bonds being replaced plays an important role on defining the degree of anisotropy. Step‐by‐step procedures for the calibration of micro parameters are recommended. The numerical model is calibrated to reproduce the behaviors of different anisotropic rocks. Detailed analyses are conducted to investigate the brittle failure process by looking at stress‐strain behaviors, increment of micro cracks, initiation and propagation of fractures. Most of these responses agree well with previous experimental findings and can provide new insights into the micro mechanisms related to the anisotropic deformation and failure behaviors. The numerical approach is then applied to simulate the stress‐induced borehole breakouts in anisotropic rock formations at reduced scale. The effect of rock anisotropy and stress anisotropy can be captured. Copyright © 2015 John Wiley & Sons, Ltd.     

A multiphase approach for evaluating the horizontal and rocking impedances of pile group foundations

This paper advocates the use of a multiphase model, already developed for static or quasi‐static geotechnical engineering problems, for simulating the behaviour of piled raft foundations subject to horizontal as well as rocking dynamic solicitations. It is shown that such a model, implemented in a FEM code, yields appropriate predictions for the foundation impedance characteristics, provided that shear and bending effects in the piles are taken into account, thus corroborating the findings of the asymptotic homogenization theory. Besides, it is notably pointed out that such a multiphase‐based computational tool makes it possible to assess the dynamic behaviour of pile groups in a much quicker way than when using direct numerical simulations, which may face oversized problems when large pile groups are concerned. Copyright © 2016 John Wiley & Sons, Ltd.     

Semi‐analytical solution to one‐dimensional consolidation for unsaturated soils with semi‐permeable drainage boundary under time‐dependent loading

This paper presents semi‐analytical solutions to Fredlund and Hasan's one‐dimensional consolidation of unsaturated soils with semi‐permeable drainage boundary under time‐dependent loadings. Two variables are introduced to transform two coupled governing equations of pore‐water and pore‐air pressures into an equivalent set of partial differential equations, which are easily solved by the Laplace transform. The pore‐water pressure, pore‐air pressure and settlement are obtained in the Laplace domain. Crump's method is adopted to perform the inverse Laplace transform in order to obtain semi‐analytical solutions in time domain. It is shown that the present solutions are more general and have a good agreement with the existing solutions from literatures. Furthermore, the current solutions can also be degenerated into conventional solutions to one‐dimensional consolidation of unsaturated soils with homogeneous boundaries. Finally, several numerical examples are provided to illustrate consolidation behavior of unsaturated soils under four types of time‐dependent loadings, including instantaneous loading, ramp loading, exponential loading and sinusoidal loading. Parametric studies are illustrated by variations of pore‐air pressure, pore‐water pressure and settlement at different values of the ratio of air–water permeability coefficient, depth and loading parameters. Copyright © 2017 John Wiley & Sons, Ltd.     

Comparison of runoff characteristics of two adjacent basins in a tropical rainforest using a modified hydrologic cycle model with outflow

We propose a new runoff model including an outflow process that was applied to two adjacent basins (CL, TL) located in Lambir Hills National Park in north‐central Sarawak, Malaysia. Rainfall, runoff, topography, and soil layer thickness were observed. About 19% of annual runoff was observed in the CL basin (21.97 ha), whereas about 46% was observed in the TL basin (23.25 ha). It was inferred that the CL basin has an outflow because of low base flow, small runoff peak, and excessive water loss. By incorporating the outflow process into the HYdrological CYcle MODEL, good agreement between the data generated by the model and that observed was shown, with the exception of the data from the rainless period. Then, the fitting parameters for each basin were exchanged, except for the outflow parameter, and the characteristics of each basin were compared by calculating virtual runoff. As a result, the low base flow of the CL basin was estimated by the movement of the rainwater that escaped from the basin as deep percolation or lateral flow (11% of rainfall). The potential of the CL basin for mitigating flood and drought appeared to be higher than that of the TL basin. This is consistent with the topographic characteristics of the CL basin, which has a gentler slope than the TL basin. Copyright © 2012 John Wiley & Sons, Ltd.     

Development of a soil moisture‐based distributed hydrologic model for determining hydrologically based critical source areas

A simple grid cell‐based distributed hydrologic model was developed to provide spatial information on hydrologic components for determining hydrologically based critical source areas. The model represents the critical process (soil moisture variation) to run‐off generation accounting for both local and global water balance. In this way, it simulates both infiltration excess run‐off and saturation excess run‐off. The model was tested by multisite and multivariable evaluation on the 50‐km² Little River Experimental Watershed I in Georgia, U.S. and 2 smaller nested subwatersheds. Water balance, hydrograph, and soil moisture were simulated and compared to observed data. For streamflow calibration, the daily Nash‐Sutcliffe coefficient was 0.78 at the watershed outlet and 0.56 and 0.75 at the 2 nested subwatersheds. For the validation period, the Nash‐Sutcliffe coefficients were 0.79 at the watershed outlet and 0.85 and 0.83 at the 2 subwatersheds. The per cent bias was less than 15% for all sites. For soil moisture, the model also predicted the rising and declining trends at 4 of the 5 measurement sites. The spatial distribution of surface run‐off simulated by the model was mainly controlled by local characteristics (precipitation, soil properties, and land cover) on dry days and by global watershed characteristics (relative position within the watershed and hydrologic connectivity) on wet days when saturation excess run‐off was simulated. The spatial details of run‐off generation and travel time along flow paths provided by the model are helpful for watershed managers to further identify critical source areas of non‐point source pollution and develop best management practices.     

The role of climate change and human interventions in affecting watershed runoff responses

Identifying the role of the two main driving factors—climate change and human interventions—in influencing runoff processes is essential for sustainable water resources management. For this purpose, runoff regime change detection methods were used to divide the available hydroclimatic variables into a baseline and a disturbed period. We applied hydrological modelling and the climate elasticity of runoff method to determine the contribution of climate change and human interventions to changes in runoff. The hydrological model, SWAT, was calibrated during the baseline period and used to simulate the naturalized runoff pattern for the disturbed period. Significant changes in runoff in the study watershed were detected from 1982, suggesting that human interventions play a dominant role in influencing runoff. The combined effects of climate change and human interventions resulted in a 41.3 mm (23.9%) decrease in runoff during the disturbed period, contributing about 40% and 60% to the total runoff change, respectively. Furthermore, analysis of changes in land cover dynamics in the watershed over the past four decades supported these changes in runoff. Contrary to other decades, the discrepancy between naturalized and observed runoff was small in the 2010s, likely due to increased baseflow as a result of storage and/or release of excess water during the dry season. This study contributes to our understanding of how climate change and human interventions affect hydrological responses of watersheds, which is important for future sustainable water management and drought adaptation.     

Using raw regional climate model outputs for quantifying climate change impacts on hydrology

General circulation model outputs are rarely used directly for quantifying climate change impacts on hydrology, due to their coarse resolution and inherent bias. Bias correction methods are usually applied to correct the statistical deviations of climate model outputs from the observed data. However, the use of bias correction methods for impact studies is often disputable, due to the lack of physical basis and the bias nonstationarity of climate model outputs. With the improvement in model resolution and reliability, it is now possible to investigate the direct use of regional climate model (RCM) outputs for impact studies. This study proposes an approach to use RCM simulations directly for quantifying the hydrological impacts of climate change over North America. With this method, a hydrological model (HSAMI) is specifically calibrated using the RCM simulations at the recent past period. The change in hydrological regimes for a future period (2041–2065) over the reference (1971–1995), simulated using bias‐corrected and nonbias‐corrected simulations, is compared using mean flow, spring high flow, and summer–autumn low flow as indicators. Three RCMs driven by three different general circulation models are used to investigate the uncertainty of hydrological simulations associated with the choice of a bias‐corrected or nonbias‐corrected RCM simulation. The results indicate that the uncertainty envelope is generally watershed and indicator dependent. It is difficult to draw a firm conclusion about whether one method is better than the other. In other words, the bias correction method could bring further uncertainty to future hydrological simulations, in addition to uncertainty related to the choice of a bias correction method. This implies that the nonbias‐corrected results should be provided to end users along with the bias‐corrected ones, along with a detailed explanation of the bias correction procedure. This information would be especially helpful to assist end users in making the most informed decisions.     

Transient seepage analysis in zoned anisotropic soils based on the scaled boundary finite‐element method

The scaled boundary finite‐element method, a semi‐analytical computational scheme primarily developed for dynamic stiffness of unbounded domains, is applied to the analysis of unsteady seepage flow problems. This method is based on the finite‐element technology and gains the advantages of the boundary element method as well. Only boundary of the domain is discretized, no fundamental solution is required and singularity problems can be modeled rigorously. Anisotropic and non‐homogeneous materials satisfying similarity are modeled with no additional efforts. In this study, firstly, formulation of the method for the transient seepage flow problems is derived followed by its solution procedures. The accuracy, simplicity and applicability of the method are demonstrated via four numerical examples of transient seepage flow – three of them are available in the literature. Homogenous, non‐homogenous, isotropic and anisotropic material properties are considered to show the versatility of the technique. Excellent agreement with the finite‐element method is observed. The method out‐performs the finite‐element method in modeling singularity points. Copyright © 2014 John Wiley & Sons, Ltd.