A multi‐level parallelized substructuring‐frontal solution for coupled thermo/hydro/mechanical problems in unsaturated soil

This paper proposes a multi‐level parallelized substructuring–frontal combined algorithm for the analysis of the problem of thermo/hydraulic/mechanical behaviour of unsaturated soil. Temperature, displacement, pore water pressure and pore air pressure are treated as the primary variables in a non‐linear analysis. Details are given firstly of the substructuring–frontal combined approach. The incorporation of the algorithm in a multi‐level parallel strategy is then discussed. The parallel processing can thus be carried out at different substructural levels. The method thus developed impacts, in a positive way, on both computer storage requirement and execution time. Copyright © 2003 John Wiley & Sons, Ltd.     

Analytical solutions for advective–dispersive solute transport in double‐layered finite porous media

Analytical solutions for advection and dispersion of a conservative solute in a one‐dimensional double‐layered finite porous media are presented. The solutions are applicable to five scenarios that have various combinations of fixed concentration, fixed flux and zero concentration gradient conditions at the inlet and outlet boundaries that provide a wide number of options. Arbitrary initial solute concentration distributions throughout the media can be considered via explicit formulations or numerical integration. The analytical solutions presented have been verified against numerical solutions from a finite‐element‐based approach and an existing closed‐form solution for double‐layered media with an excellent correlation being found in both cases. A practical application pertaining to advective transport induced by consolidation of underlying sediment layers on contaminant movement within a capped contaminated sediment system is presented. Comparison of the calculated concentrations and fluxes with alternative approaches clearly illustrates the need to consider advection processes. Consideration of the different features of contaminant transport due to varying pore‐water velocity fields in primary consolidation and secondary consolidation stages is achieved via the use of non‐uniform initial concentration distributions within the proposed analytical solutions. Copyright © 2010 John Wiley & Sons, Ltd.     

Multi-dimensional chemo-osmotic consolidation of clays

A model of chemo-osmotic consolidation of clays in multi-dimensional domains is presented, with volumetric strains induced by both changes in the chemistry and osmotically driven pore water flow considered. Three fully coupled governing equations considering force equilibrium, pore water transport and chemical transport are presented and solved using the transient finite element method. The proposed approach is verified via consideration of chemo-osmotic consolidation of a compacted clay landfill liner and then applied to investigation of a hypothetical case with a local leachate leak in the compacted clay liner. An assessment of the impact of the two-dimensional nature of the system is made. The consolidation process is found to be dominated by osmotic consolidation in the early stages and subsequently by chemical consolidation. It is found that the surface settlement and the leachate concentration in the compacted clay liner may be highly overestimated by a one-dimensional analysis. Moreover, the peak negative excess pore pressure predicted by the two-dimensional solution remains in a shallow region under the leak while in the one-dimensional solution it moves progressively downwards.     

Consolidation of sensitive clays: a numerical investigation

Consolidation of sensitive clay layers, which have significant compressibility at different stress states, is investigated via a nonlinear one-dimensional consolidation approach with a piecewise linear e ~ log₁₀σ′ model. The behaviour of sensitive clays during consolidation and the limitations of conventional consolidation theory are addressed. It is shown that (1) the sensitive clay layer can be divided by the preconsolidation pressure into two zones, that is, high- and low-compressibility zones. The progressive destruction of particle cementation bonding through the soil layer is shown by the moving front of the interface between these two zones; (2) the excess pore pressure dissipation primarily takes place in the low-compressibility zone, which results in a small settlement during the early stages of consolidation; (3) conventional consolidation theory highly overestimates the settlement and gives a poor prediction of effective stress distribution.     

Hydro-geomechanical characterisation of a coastal urban aquifer using multiscalar time and frequency domain groundwater-level responses

Hydraulic properties of coastal, urban aquifers vary spatially and temporally with the complex dynamics of their hydrogeology and the heterogeneity of ocean-influenced hydraulic processes. Traditional aquifer characterisation methods are expensive, time-consuming and represent a snapshot in time. Tidal subsurface analysis (TSA) can passively characterise subsurface processes and establish hydro-geomechanical properties from groundwater head time-series but is typically applied to individual wells inland. Presented here, TSA is applied to a network of 116 groundwater boreholes to spatially characterise confinement and specific storage across a coastal aquifer at city-scale in Cardiff (UK) using a 23-year high-frequency time-series dataset. The dataset comprises Earth, atmospheric and oceanic signals, with the analysis conducted in the time domain, by calculating barometric response functions (BRFs), and in the frequency domain (TSA). By examining the damping and attenuation of groundwater response to ocean tides (OT) with distance from the coast/rivers, a multi-borehole comparison of TSA with BRF shows this combination of analyses facilitates disentangling the influence of tidal signals and estimation of spatially distributed aquifer properties for non-OT-influenced boreholes. The time-series analysed covers a period pre- and post-impoundment of Cardiff’s rivers by a barrage, revealing the consequent reduction in subsurface OT signal propagation post-construction. The results indicate that a much higher degree of confined conditions exist across the aquifer than previously thought (specific storage?=?2.3 × 10^?6 to 7.9 × 10^?5 m^?1), with implications for understanding aquifer recharge, and informing the best strategies for utilising groundwater and shallow geothermal resources.