Three dimensional numerical simulation of residential building on shrink–swell soils in response to climatic conditions

Shrink–swell soils can cause distresses in buildings, and every year, the economic loss associated with this problem is huge. This paper presents a comprehensive system for simulating the soil–foundation–building system and its response to daily weather conditions. Weather data include rainfall, solar radiation, air temperature, relative humidity, and wind speed, all of which are readily available from a local weather station or the Internet. These data are used to determine simulation flux boundary conditions. Different methods are proposed to simulate different boundary conditions: bare soil, trees, and vegetation. A coupled hydro‐mechanical stress analysis is used to simulate the volume change of shrink–swell soils due to both mechanical stress and water content variations. Coupled hydro‐mechanical stress‐jointed elements are used to simulate the interaction between the soil and the slab, and general shell elements are used to simulate structural behavior. All the models are combined into one finite element program to predict the entire system's behavior. This paper first described the theory for the simulations. A site in Arlington, Texas, is then selected to demonstrate the application of the proposed system. Simulation results are shown, and a comparison between measured and predicted movements for four footings in Arlington, Texas, over a 2‐year period is presented. Finally, a three‐dimensional simulation is made for a virtual residential building on shrink–swell soils to identify the influence of various factors. Copyright © 2015 John Wiley & Sons, Ltd.     

Computational model coupling mode II discrete fracture propagation with continuum damage zone evolution

We propose a numerical method that couples a cohesive zone model (CZM) and a finite element‐based continuum damage mechanics (CDM) model. The CZM represents a mode II macro‐fracture, and CDM finite elements (FE) represent the damage zone of the CZM. The coupled CZM/CDM model can capture the flow of energy that takes place between the bulk material that forms the matrix and the macroscopic fracture surfaces. The CDM model, which does not account for micro‐crack interaction, is calibrated against triaxial compression tests performed on Bakken shale, so as to reproduce the stress/strain curve before the failure peak. Based on a comparison with Kachanov's micro‐mechanical model, we confirm that the critical micro‐crack density value equal to 0.3 reflects the point at which crack interaction cannot be neglected. The CZM is assigned a pure mode II cohesive law that accounts for the dependence of the shear strength and energy release rate on confining pressure. The cohesive shear strength of the CZM is calibrated by calculating the shear stress necessary to reach a CDM damage of 0.3 during a direct shear test. We find that the shear cohesive strength of the CZM depends linearly on the confining pressure. Triaxial compression tests are simulated, in which the shale sample is modeled as an FE CDM continuum that contains a predefined thin cohesive zone representing the idealized shear fracture plane. The shear energy release rate of the CZM is fitted in order to match to the post‐peak stress/strain curves obtained during experimental tests performed on Bakken shale. We find that the energy release rate depends linearly on the shear cohesive strength. We then use the calibrated shale rheology to simulate the propagation of a meter‐scale mode II fracture. Under low confining pressure, the macroscopic crack (CZM) and its damaged zone (CDM) propagate simultaneously (i.e., during the same loading increments). Under high confining pressure, the fracture propagates in slip‐friction, that is, the debonding of the cohesive zone alternates with the propagation of continuum damage. The computational method is applicable to a range of geological injection problems including hydraulic fracturing and fluid storage and should be further enhanced by the addition of mode I and mixed mode (I+II+III) propagation. Copyright © 2016 John Wiley & Sons, Ltd.     

A nonstationary extreme value distribution for analysing the cessation of karst spring discharge

The effects of climate change and population growth in recent decades are leading us to consider their combined and potentially extreme consequences, particularly regarding hydrological processes, which can be modeled using a generalized extreme value (GEV) distribution. Most of the GEV models were based on a stationary assumption for hydrological processes, in contrast to the nonstationary reality due to climate change and human activities. In this paper, we present the nonstationary generalized extreme value (NSGEV) distribution and use it to investigate the risk of Niangziguan Springs discharge decreasing to zero. Rather than assuming the location, scale, and shape parameters to be constant as one might do for a stationary GEV distribution analysis, the NSGEV approach can reflect the dynamic processes by defining the GEV parameters as functions of time. Because most of the GEV model is designed to evaluate maxima (e.g. flooding, represented by positive numbers), and spring discharge cessation is a ?minima’, we deduced an NSGEV model for minima by applying opposite numbers, i.e. negative instead of positive numbers. The results of the model application to Niangziguan Springs showed that the probability of zero discharge at Niangziguan Springs will be 1/80 in 2025, and 1/10 in 2030. After 2025, the rate of decrease in spring discharge will accelerate, and the probability that Niangziguan Springs will cease flowing will dramatically increase. The NSGEV model is a robust method for analysing karst spring discharge. Copyright © 2013 John Wiley & Sons, Ltd.     

A new standardized Palmer drought index for hydro‐meteorological use

Accepting the concept of standardization introduced by the standardized precipitation index, similar methodologies have been developed to construct some other standardized drought indices such as the standardized precipitation evapotranspiration index (SPEI). In this study, the authors provided deep insight into the SPEI and recognized potential deficiencies/limitations in relating to the climatic water balance it used. By coupling another well‐known Palmer drought severity index (PDSI), we proposed a new standardized Palmer drought index (SPDI) through a moisture departure probabilistic approach, which allows multi‐scalar calculation for accurate temporal and spatial comparison of the hydro‐meteorological conditions of different locations. Using datasets of monthly precipitation, temperature and soil available water capacity, the moisture deficit/surplus was calculated at multiple temporal scales, and a couple of techniques were adopted to adjust corresponding time series to a generalized extreme value distribution out of several candidates. Results of the historical records (1900–2012) for diverse climates by multiple indices showed that the SPDI was highly consistent and correlated with the SPEI and self‐calibrated PDSI at most analysed time scales. Furthermore, a simple experiment of hypothetical temperature and/or precipitation change scenarios also verified the effectiveness of this newly derived SPDI in response to climate change impacts. Being more robust and preferable in spatial consistency and comparability as well as combining the simplicity of calculation with sufficient accounting of the physical nature of water supply and demand relating to droughts, the SPDI is promising to serve as a competent reference and an alternative for drought assessment and monitoring. Copyright © 2013 John Wiley & Sons, Ltd.     

Dissolved organic matter in surface runoff in the Loess Plateau of China: The role of rainfall events and land-use

Exploring the chemical characterization of dissolved organic matter (DOM) is important for understanding the fate of laterally transported organic matter in watersheds. We hypothesized that differences in water-extractable organic matter (WEOM) in soils of varying land uses and rainfall events may significantly affect the quality and the quantity of stream DOM. To test our hypotheses, characteristics of rainfall-runoff DOM and WEOM of source materials (topsoil from different land uses and gullies, as well as typical vegetation) were investigated at two adjacent catchments in the Loess Plateau of China, using ultraviolet–visible absorbance and excitation emission matrix fluorescence with parallel factor analysis (PARAFAC). Results indicated that land-use types may significantly affect the chemical composition of soil WEOM, including its aromaticity, molecular weight, and degree of humification. The PARAFAC analysis demonstrated that the soils and stream water were dominated by terrestrial/allochthonous humic-like substances and microbial transformable humic-like fluorophores. Shifts in the fluorescence properties of stream DOM suggested a pronounced change in the relative proportion of allochthonous versus autochthonous material under different rainfall patterns and land uses. For example, high proportions of forestland could provide more allochthonous DOM input. This study highlights the relevance of soils and hydrological dynamics on the composition and fluxes of DOM issuing from watersheds. The composition of DOM in soils was influenced by land-use type. Precipitation patterns influenced the proportion of terrestrial versus microbial origins of DOM in surface runoff. Contributions of allochthonous, terrestrially derived DOM inputs were highest from forested landscapes.     

Landform‐oriented flow‐routing algorithm for the dual‐structure loess terrain based on digital elevation models

The loess landform in the Loess Plateau of China is with typical dual structure, namely, the upper smooth positive terrain and the lower cliffy negative terrain (P–N terrain for short). Obvious differences in their morphological feature, geomorphological mechanism, and hydrological process could be found in the both areas. Based on the differences, a flow‐routing algorithm that separately addresses the dual‐structure terrain would be necessary to encompass this spatial variation in their hydrological behaviour. This paper proposes a mixed flow‐routing algorithm to address aforementioned problems. First, the loess landform surface is divided into P–N terrains based on digital elevation models. Then, specific catchment area is calculated with the new algorithm to simulate the water flows in both positive and negative terrain areas. The mixed algorithm consists of the multiple flow‐routing algorithm (multiple‐flow direction) for positive areas and the D8 algorithm for negative areas, respectively. The approach is validated in two typical geomorphologic areas with low hills and dense gullies in the northern Shaanxi Loess Plateau. Four indices are used to examine the results, which show that the new algorithm is more suitable for loess terrain in simulating the spatial distribution of water accumulation, as well as in modeling the flow characteristics of the true surface by considering the morphological structures of the terrain. Copyright © 2013 John Wiley & Sons, Ltd.     

Parallelized combined finite-discrete element (FDEM) procedure using multi-GPU with CUDA

This paper focuses on the efficiency of finite discrete element method (FDEM) algorithmic procedures in massive computers and analyzes the time-consuming part of contact detection and interaction computations in the numerical solution. A detailed operable GPU parallel procedure was designed for the element node force calculation, contact detection, and contact interaction with thread allocation and data access based on the CUDA computing. The emphasis is on the parallel optimization of time-consuming contact detection based on load balance and GPU architecture. A CUDA FDEM parallel program was developed with the overall speedup ratio over 53 times after the fracture from the efficiency and fidelity performance test of models of in situ stress, UCS, and BD simulations in Intel i7-7700K CPU and the NVIDIA TITAN Z GPU. The CUDA FDEM parallel computing improves the computational efficiency significantly compared with the CPU-based ones with the same reliability, providing conditions for achieving larger-scale simulations of fracture.     

Frequency analysis of nonstationary annual maximum flood series using the time‐varying two‐component mixture distributions

The most popular practice for analysing nonstationarity of flood series is to use a fixed single‐type probability distribution incorporated with the time‐varying moments. However, the type of probability distribution could be both complex because of distinct flood populations and time‐varying under changing environments. To allow the investigation of this complex nature, the time‐varying two‐component mixture distributions (TTMD) method is proposed in this study by considering the time variations of not only the moments of its component distributions but also the weighting coefficients. Having identified the existence of mixed flood populations based on circular statistics, the proposed TTMD was applied to model the annual maximum flood series of two stations in the Weihe River basin, with the model parameters calibrated by the meta‐heuristic maximum likelihood method. The performance of TTMD was evaluated by different diagnostic plots and indexes and compared with stationary single‐type distributions, stationary mixture distributions and time‐varying single‐type distributions. The results highlighted the advantages of TTMD with physically‐based covariates for both stations. Besides, the optimal TTMD models were considered to be capable of settling the issue of nonstationarity and capturing the mixed flood populations satisfactorily. Copyright © 2016 John Wiley & Sons, Ltd.