Saturated Media Research Articles

Tracking microbial movement in saturated media with spectral induced polarization

Pathogenic microorganisms in the subsurface can contaminate soil and water supplies, potentially posing great danger to human health. Early contamination detection routines rely on sparse direct sampling which is spatiotemporally limited. Thus, the path of microorganisms in the subsurface remains ambiguous and this can cause delays in detection of biohazardous threats. The geophysical spectral induced polarization (SIP) technique, sensitive to microbes' presence and activity in porous media, is a promising method to monitor microbial transport pathways. Here we evaluated the efficiency of SIP in monitoring the chemotactic movement of Sporosarcina pasteurii in saturated porous media. A cylindrical sample holder was packed with Ottawa sand and saturated with sterile KCl solution. The sample holder was oriented vertically and S. pasteurii was introduced at the bottom, forcing the movement of the microbes against gravity, towards a carbon source available at the top of the column. Temporal SIP measurements were collected at 3 regions of the sample holder: bottom (microbial injection point), middle and top (carbon source). Both the real (σ′) and imaginary (σ″) conductivity parts of the SIP signal increased over time with the σ″ showing a peak signal magnitude following the upward movement of the microbes. We repeated the experiment excluding the carbon source in experiment 2 and omitting microbial injection in experiment 3. However, we did not observe any significant SIP signal changes in these two experiments. This is the first study to indicate the strong SIP signal correlation with microbial chemotaxis.

Error Estimates and Adaptivity of the Space-Time Discontinuous Galerkin Method for Solving the Richards Equation

We present a higher-order space-time adaptive method for the numerical solution of the Richards equation that describes a flow motion through variably saturated media. The discretization is based on the space-time discontinuous Galerkin method, which provides high stability and accuracy and can naturally handle varying meshes. We derive reliable and efficient a posteriori error estimates in the residual-based norm. The estimates use well-balanced spatial and temporal flux reconstructions which are constructed locally over space-time elements or space-time patches. The accuracy of the estimates is verified by numerical experiments. Moreover, we develop the hp-adaptive method and demonstrate its efficiency and usefulness on a practically relevant example.

Complex dynamics of a fractional-order epidemic model with saturated media effect

Complex dynamics of a fractional-order epidemic model with saturated media effect

Mode I sub- and super-shear rupture and forerunning in porous media

Numerical simulations demonstrate the occurrence of both sub- and supershear ruptures in dry and saturated porous media during mode I fracture. Initially, the Finite Element Method is applied to investigate a case of sub-shear velocity in a dry beam on an elastic foundation subjected to sinusoidal loading. Subsequently, the Extended Finite Element Method is employed for fluid injection in both an infinite saturated porous medium and a plate. The results reveal a transition from smooth behavior to a stepwise pattern and eventually to forerunning as the injection rate increases. Notably, this progression has been observed for the first time. Finally, a hybrid Finite Element/Peridynamic model is utilized to explore both dry and fully saturated media under mechanical loading. In the case of saturated media, fluid injection is also considered. Both methods document velocities surpassing the shear wave. The noteworthy aspect is that these results, obtained through entirely different numerical methods and constitutive relations, justify confidence in the assertion that, in mode I fracture, speeds exceeding the shear wave and even the dilatation wave speeds can be achieved through forerunning. The loading condition investigated may be relevant in geophysics.

Screening effect of pile rows for ground vibration induced by moving loads

This study aims to predict the screening efficiency of pile rows in transversely isotropic (TI) saturated soils for ground vibrations subjected to moving loads. First, based on the fundamental solutions of TI saturated media, the boundary element method (BEM) was adopted to derive the force-deformation equation of the soil-pile boundary. The finite element method (FEM) was then applied to establish the stiffness matrix equation of the Euler-based pile rows. Subsequently, by coupling the BEM and FEM equations, the screening efficiency of the pile rows was developed according to the barrier isolation theory. The proposed approach was validated using published numerical results. Finally, the effect of soil stratification on the screening effect of pile rows in TI saturated soils was investigated.

Performance analysis of multilayered transversely isotropic saturated media under temperature and horizontally circular loads

AbstractThis study constructs a multilayered transversely isotropic saturated model under thermal and horizontally circular loads, and further investigates the model's thermo‐hydro‐mechanical coupling response. Firstly, the ordinary differential matrix equations of thermoelastic saturated media in the integral transformed domain are derived. Secondly, the solution for multilayered thermoelastic saturated media is developed using the extended precise integration method (EPIM), along with the boundary conditions at both ends of the foundation and the continuity conditions between adjacent layers. After that, the solution in the physical domain is further attained with the use of Laplace‐Hankel integral transform inversion. Finally, the accuracy of the proposed theory is confirmed through numerical examples, and the influences of anisotropic parameters, the soil's stratification and porosity on the thermo‐hydro‐mechanical coupling response of the media are studied.

Transversely isotropic effects on the coupled thermo‐hydro‐mechanical performance for layered saturated media

AbstractIn this manuscript, a novel transformed differential quadrature solution to the coupled thermo‐hydro‐mechanical (THM) problem of layered transversely isotropic (TI) saturated media is proposed, accompanied by a sensitivity analysis of pertinent parameters. Initially, the THM governing equations that encompass the transverse isotropy characteristics of thermal, permeable, and mechanical properties are introduced. Subsequently, utilizing the transformed differential quadrature method (TDQM), the discretized algebraic equations in the transformed domain are derived. Boundary conditions related to stresses, displacements, and temperatures are defined leveraging integral transform techniques and the stress‐strain relationship, while algebraic equations for continuity conditions across adjacent layers are incorporated. The global matrix equation is assembled and solved. After verifying the proposed TDQM solution, a sensitivity analysis on the complex transversely isotropic parameters is conducted to identify the principal factors influencing THM coupling behaviors, thereby simplifying the parameters for discussion in practical engineering applications.

Environmental colloid behaviors of humic acid - Cadmium nanoparticles in aquatic environments

Humic acid (HA), a principal constituent of natural organic matter (NOM), manifests ubiquitously across diverse ecosystems and can significantly influence the environmental behaviors of Cd(II) in aquatic systems. Previous studies on NOM-Cd(II) interactions have primarily focused on the immobilization of Cd(II) solids, but little is known about the colloidal stability of organically complexed Cd(II) particles in the environment. In this study, we investigated the formation of HA-Cd(II) colloids and quantified their aggregation, stability, and transport behaviors in a saturated porous media representative of typical subsurface conditions. Results from batch experiments indicated that the relative quantity of HA-Cd(II) colloids increased with increasing C/Cd molar ratio and that the carboxyl functional groups of HA dominated the stability of HA-Cd(II) colloids. The results of correlation analysis between particle size, critical aggregation concentration (CCC), and zeta potential indicated that both Derjaguin-Landau-Verwey-Overbeek (DLVO) and non-DLVO interactions contributed to the enhanced colloidal stability of HA-Cd(II) colloids. Column results further confirmed that the stable HA-Cd(II) colloid can transport fast in a saturated media composed of clean sand. Together, this study provides new knowledge of the colloidal behaviors of NOM-Cd(II) nanoparticles, which is important for better understanding the ultimate cycling of Cd(II) in aquatic systems.

An efficient fully Crouzeix-Raviart finite element model for coupled hydro-mechanical processes in variably saturated porous media

Coupled hydro-mechanical processes in variably saturated porous media can be encountered in several applications in geosciences and soil mechanics. They are ruled by highly nonlinear and strongly coupled equations of fluid flow in unsaturated porous media and quasi-static mechanical deformation. In this work, we develop an efficient and robust numerical scheme for unsaturated poroelasticity that allows for balancing computational requirements and accuracy. The spatial discretization is based on the Crouzeix-Raviart (CR) finite element method for both the displacement field and the hydraulic head. The CR method is locally mass conservative and uses low-order approximation, which alleviates the computational burden. The time discretization of the coupled nonlinear hydraulic mechanical equations is performed using high order time integration methods and efficient time stepping schemes via the method of lines (MOL).Four test cases are investigated to highlight the efficiency and accuracy of the developed numerical scheme. The first test case deals with the settlement induced by gravity in the case of saturated and unsaturated soils. A new semi-analytical solution is developed for the settlement in the unsaturated case and used for the validation of the numerical scheme. The second test case is the Cantilever bracket problem, extended here to variably saturated media and studied to investigate the nonphysical oscillations of the pressure variable with different numerical schemes. The results are compared to finite element solutions obtained with COMSOL with either linear or quadratic approximation of the displacement field. The third test case deals with the deformation of an unsaturated soil under partial infiltration and surface load. This test case is simulated to investigate the accuracy and efficiency of the CR scheme for the challenging problem of infiltration into dry soils. The last test case is simulated to show the applicability of the developed model to a realistic problem dealing with the stability of a hillslope under intense rainfall periods.Results of the numerical experiments show that the CR scheme coupled with high order time integration methods allows for significant gain in both computational requirements and accuracy.

Tomato Seedling Performance in Commercial Organic Growing Media

Healthy transplants are critical to productivity in the field. For certified organic production in the United States, seedlings must be grown in media that meet the standards of the US Department of Agriculture’s National Organic Program. Many commercial organic media options are available, they vary substantially in composition, and it is unknown to what extent this influences seedling performance. This project compared tomato (Solanum lycopersicum L.) seedling emergence and growth in seven commercially available media for organic production and evaluated posttransplant performance. Tomato seedlings were grown in greenhouses at Wanatah, West Lafayette, and Vincennes, IN, USA. Chemical characteristics of the media measured in saturated media extract ranged as follows: pH 5.2–7.5; electrical conductivity (EC) 0.79–4.68 dS·m−1; 1–332 ppm nitrate-nitrogen, 5–69 ppm phosphorus, 41–451 ppm potassium, 78–714 ppm calcium, and 25–121 ppm magnesium. Higher media EC was associated with slower and less uniform seedling emergence and reduced total emergence. Seedling aboveground dry weights were significantly greater in media that contained compost. Relative performance in media containing compost varied across trials. The aboveground dry weight of tomato seedlings 4 weeks after transplanting did not differ for seedlings started in the five compost-based media, and those plants were significantly larger than plants started in the two media without compost. Larger plants tended to flower and set fruit earlier. Media testing protocols that predict nutrient supply over the production cycle could likely improve management in organic transplant production.

Threshold dynamics of a stochastic infectious disease model with vaccination age under saturated media coverage

Threshold dynamics of a stochastic infectious disease model with vaccination age under saturated media coverage

A two-stage based model for partially-buried pile group response in cross-anisotropic saturated media adjacent to moving harmonic loads

The pile foundations are frequently affected by adjacent traffic loads except for general active loads. However, the dynamic responses of the pile groups under moving loads have been rarely reported before. In addition, the influence of material anisotropy is often neglected. In this paper, a three-dimensional (3D) analytical model for the dynamic analysis of partially buried pile groups in stratified saturated cross-anisotropic media under adjacent moving harmonic loadings is developed. Specifically, a 3D Bernoulli–Euler beam theory is adopted to simulate the monopile and then superimposed into the finite element (FE) formulation for the partially embedded pile groups. Subsequently, the basic solutions of layered saturated cross-anisotropic soils are derived by the analytical layer-element approach (ALEA). By combination with ALEA, the two-stage based boundary element (BE) equations for the pile-soil interface are obtained. Finally, by coupling the FE and BE formulations, the dynamical response equation of pile-soil is established. Based on the verification of the accuracy of the proposed methodology by comparisons with existing solutions and FE results from ABAQUS, the impacts of material anisotropy, pile buried length, loading velocity, and pile stiffness on the time-domain dynamic behavior of pile groups are analyzed. The results show that increasing the pile-buried ratio or pile stiffness reduces the magnitude of the dynamic response. Moreover, with the increase of material anisotropic parameters, the peak pile response decreases. And the peak dynamic response first increases and then reduces as the loading velocity increases.

Numerical modeling of drying from saturated to unsaturated conditions with two-phase partitioning boundaries at the medium surface and interface-dependent phase change

To accurately predict the drying of a porous medium from saturated to unsaturated conditions over time, this paper proposes two-phase partitioning boundaries for liquid water evaporation and water vapor transport at the medium surface. A non-isothermal two-phase (liquid water and gas mixtures), two-component (dry air and water vapor) transport model with nonequilibrium and interface-dependent phase change is then developed. The governing equations are then numerical solved using COMSOL Multiphysics software and applied to describe a series of laboratory soil drying tests as examples. The performance of the proposed boundaries at the surface is illustrated by comparison with other three common used boundaries. The results show that a specific flux of vapor at the surface with a zero flux of liquid water cannot be adopted to describe the drying process of initially saturated media due to the lack of interconnected gas-filled pores. The merged liquid-vapor boundary applied in the Darcy–Richards model significantly prolongs the constant rate stage when the liquid water evaporation is dominant. Then, the effects of air/gas–water interfacial area and water retention parameters on the drying process are investigated. A faster reduction in the air–water relative interfacial area at the medium surface will shorten the drying process with earlier dominance of the vapor transport. An increased entry capillary pressure and a decreased pore size distribution index for the water retention curve will result in slower drying, more loss by vapor transport and a more homogeneous distribution of pore water respect to depth. The current study can promote the understanding of the evolution of drying mechanisms of a porous medium from saturated to unsaturated conditions.

Approximation of van Genuchten Parameter Ranges from Hydraulic Conductivity Data.

The use of retention function and relative conductivity function is essential for the calculation of flow in a variably saturated media using the Richards equation. A widely used mathematical model for this is the Mualem-van Genuchten model which requires the shape parameters and . These, however, are difficult to obtain. When data is scarce, and are often taken from literature and may deviate largely from actual values. The current study presents a novel mathematical model for the approximation of and and for the further estimation of realistic value ranges, which may be used as parameter space, for example, for the calibration of a numerical model. The model was developed for cases where data is scarce and only values of saturated hydraulic conductivity are available. It is based on a large data set from literature and it is demonstrated that the model estimates mean values from that data set with a good accuracy. In order to show the applicability of the model, a second data set has been compiled anew (provided as Supporting Information). The model is incorporated into the current version of the freeware computer program HYPAGS, which enables easy usage.

On the effects of Br−, ClO4−, NO3−anions on the oxygen reduction reaction electrocatalysis at Fe phthalocyanine modified electrodes

Fe phthalocyanine with carbon nanotubes (FePc-CNT) has been used as catalysts for ORR using the following supporting electrolyte: 0.1 M NaOH, 1 M NaOH, 1 M HNO3, 1 M HClO4, and 1 M HBr. The thermodynamic and kinetic information was acquired from cyclic voltammetry in N2 saturated media, and then polarization curves saturating with O2. The rate-determining step is modulated by the formal potential (E0’Fe(III)/(II)) in basic media, coherent with the inner sphere electron transfer mechanism. However, in strong acid electrolytes, the Tafel slope value is ∼ -0.090 V dec−1 at low overpotentials and ∼ -0.120 V dec−1 at high overpotentials indicating that a poisoning process is hampering the ORR, and causing an outer sphere electron transfer mechanism. Impedance spectroscopy and poisoning test (recovery of the activity of the catalyst in NaOH after initial use in acid) support the thesis of ion adsorption/poisoning for the low electrocatalytic activity of FePc based catalysts in acid media.

Effectiveness of self-sealing after gas transport in Boom Clay

The study of gas transport in low permeable materials is becoming a significant focus in energy-related geotechnics, particularly for managing deep geological disposal of long-lived and heat-emitting radioactive waste. Argillaceous rocks, studied to host the disposal, may present induced fractures caused by the excavation activities that increase the permeability to liquid and gas. The generation of gases can also lead to an excessive pressure build-up in these saturated media resulting in the development or reactivation of fractures/fissures creating preferential pathways for the gas flow [1]. Nevertheless, these rocks present the advantage of self-seal (via swelling of clay minerals due to re-saturation, consolidation or creep), reducing fissure permeability and potentially restoring the barrier function [2].
 This study focuses on Boom Clay, a Cenozoic clay candidate for hosting the repository in the Belgian programme, aiming at characterising the effectiveness of its self-sealing process after gas injection/dissipation tests due to the swelling of clay minerals during re-saturation. This poorly indurated rock presents sedimentary bedding planes that can act as preferential pathways during the gas invasion. In a previous experimental campaign, samples at two bedding orientations (parallel and orthogonal to the flow) were tested under oedometer conditions [3, 4]. The results indicate that gas transport induced slight expansion of the samples in both orientations, increasing their intrinsic permeability, which suggested the opening of preferential paths. Furthermore, the analyses of the pore network using mercury intrusion porosimetry (MIP) and micro-focus X-ray computed tomography (µ-CT) confirmed the development of gas pathways following the bedding direction or interconnecting bedding planes [4, 5].
 To assess the self-sealing capacity, oedometer tests were carried out at two bedding directions using an oedometer cell with lateral stress measurement to comprehensively understand the stress state and ensure that gas flow occurred through the sample rather than between the sample-ring interface. Prior to the gas injection stage, water permeability was measured. Subsequently, a gas injection/dissipation stage was performed at constant vertical stress. Immediately after, the sample was placed in contact with synthetic water to allow re-saturation. Finally, the water permeability was measured again. The self-sealing capacity of the clay was assessed by comparing the water permeability at both stages. If the obtained values are similar, the barrier function has been restored. Some tests also included a second gas injection to evaluate the potential pathway reopening. MIP and µ-CT data allowed for evaluating microstructural changes due to these processes.
 The experimental results demonstrate that gas can flow at pressures lower than the minor lateral stress, and the computed effective gas permeability is always higher than the initial intrinsic water permeability, pointing to the opening of preferential pathways (Figure 1a). Water permeability after re-saturation showed values comparable to the initial intrinsic permeability, recovering the hydraulic barrier function thanks to the clay minerals’ swelling (Figure 1a). This indicates the good self-sealing capacity of the Boom Clay. The gas permeability calculated during a second injection/dissipation stage is slightly higher than the initial one, suggesting some memory of the previously opened path (Figure 1a). The self-sealing effect was also observed in the CT images performed under unstressed conditions after water-undrained unloading of the tested samples. The fissures detected after the gas injection (Figure 1b left) are no longer visible after re-saturation within the technique resolution (fissures > 40µm) (Figure 1b middle). However, a small proportion of large and disconnected pores were identified in CT images, likely due to some gas exsolution. Regarding the subsequent gas injection, it again led to the development of fissures following the bedding direction detected with microstructural techniques (Figure 1b right).

Insight into Impact of Phosphate on the Cotransport and Corelease of Eu(III) with Bentonite Colloids in Saturated Quartz Columns

Understanding the fate and transport of radionuclides in porous media reduces the risk of contaminating soils and groundwater systems. While the cotransport of bentonite colloids (BC) with radionuclides in saturated media is well documented, the role of phosphate (P) in the colloid-driven transport of radionuclides in saturated porous media is still unaddressed; in particular, phosphate increases the mobilities of radionuclides in porous media, which should be subjected to an environmental risk assessment and model construction. In this work, the effects of phosphate on the transport and release of Eu(III) in different colloid systems (P-Eu(III), P-BC, P-BC-Eu(III)) was investigated with a fundamental colloid chemistry approach and a range of characterization techniques. The results showed that intrinsic europium colloids with size of 685 nm were formed by precipitation with phosphate, which affected the mobility of Eu(III) due to colloid stability and physical straining. Phosphate enhanced BC and BC-Eu(III) transport, and a high phosphate concentration promoted BC transport by eliminating physical straining and enhancing the electrostatic repulsions. The crystal structure of EuPO4 was not destroyed by the subsequent introduction of BC, which carried EuPO4 for further migration. However, when phosphate, bentonite and Eu(III) coexisted in a colloid suspension, the phosphate promoted Eu(III) transport by preferentially interacting with the BC to form ternary BC-P-Eu(III) pseudo-colloids rather than forming the intrinsic EuPO4 colloids. The synergetic role of P and BC on Eu(III) transport involved a relatively complex process and was not a simply additive effect. The findings in this work highlight the significance of phosphate in controlling the fate and transport of Ln(III)/Am(III) radionuclides in the presence of intrinsic colloids and pseudo-colloids in P-rich colloid-bearing environments.

On the Possibility of Detecting Pore Pressure Changes in Marine Sediments Using Bottom Seismometer Data

This paper is devoted to the detection and analysis of overpressure zones in unconsolidated seafloor sediments using an ocean-bottom seismometer. The methodological aspects of creating a system of anomalous pore pressure zone detection in marine sediments are studied. The aim of this study is to establish the requirements for a pore pressure monitoring system necessary to successfully detect overpressure zones based on seismic response, and to analyze temporal changes in pore pressure distribution. Data from a certain offshore field are used as a basis from which to construct synthetic models of overpressure distribution in marine sediments. Synthetic models are constructed using specially developed rock physics models for unconsolidated saturated media. Seismic responses are calculated for these synthetic models to represent data that otherwise would be obtained from bottom seismometers placed on the seafloor. Resultant seismic responses are studied with respect to the detection of overpressure zones. Possibilities and limitations of bottom seismometer data are discussed. Requirements for the frequency bands of bottom seismometers are formulated based on the results that are obtained.

Experimental investigation on the efficiency of the phase change materials for enhancing the thermal performance of energy piles in sandy soils

The application of phase change materials (PCMs) is an efficient technique for enhancing the thermal properties of energy piles. In most previously applied studies, paraffin-based materials with low melting temperatures were implemented as PCM to improve energy piles in cold climate conditions. This study explored the efficiency of using lauric acid as a PCM in energy piles by performing small-scale physical model tests under warm climate conditions. Two concrete energy piles with the same diameter of 58 mm and length of 830 mm with PCM and without it were placed in the center of a tank filled with silica sand. The energy piles went through a series of heating tests with the loading temperature (corresponding to the initial temperature) of +10 to +40 °C. To analyze the effect of implementing lauric acid on temperature, and pore water pressure (PWP) distributions in sandy soil, digital thermometers and pore pressure transducers were used in different parts of the inlet side of the model. Results were compared between the two piles with and without PCM in dry and saturated media. The achieved results demonstrate that by adding 1.5% of lauric acid to the energy pile, temperature changes along the pile can be reduced up to 13% in the dry and 10% in the saturated condition. Also, the temperature changes in the soil media could be decreased by up to 5% in dry media. Additionally, the induced excess pore water pressure around the modified energy pile with lauric acid was about 2.14 times lower than the conventional pile. Besides the reduction in the temperature and PWP, the extracted power had an upward trend by increasing the thermal loading. The results revealed that due to the phase transition process, applying a small amount of lauric acid in the energy piles reduces the temperature fluctuations-related responses and could improve the thermal performance of the heat exchange system.

A 3D coupled ALEM-BEM method for assessing pile mitigation performance to moving-harmonic-load induced vibrations in cross-anisotropic porous saturated media

The present investigation aims at proposing a three-dimensional ALEM-BEM coupling analytical approach to evaluate the isolation effectiveness of pile barriers in stratified cross-anisotropic porous saturated media under moving harmonic loads. The governing equations of Biot’s porous elastic theory are decoupled through an analytical layer-element method (ALEM). The pile-soil boundary equation is derived using BEM combined with a two-stage theory. Two-node beam cells are adopted to discretize the piles based on Timoshenko beam theory. The pile barrier mitigation effectiveness is obtained depending upon a vibration screening theory. After carefully validating the presented solutions by various reduced instances, the impacts of two-phase compressibility, geomaterial cross anisotropy, and load frequency on the mitigation effectiveness are analyzed. Results illustrate that the isolation effectiveness is more sensitive to the velocity in the medium with high two-phase compressibility or low cross-anisotropic coefficient. The pile screening effectiveness influenced by a high-speed load is more sensitive to the frequency.