A Pore-Scale Model for Dispersion and Mass Transfer during Acoustically Assisted Miscible Displacements in Porous Media

Break-up and mobilization of DNAPL by acoustic excitation: Experimental evidence and pore network modeling

Dense non-aqueous phase liquids (DNAPLs) are long-term groundwater contaminants due to their high toxicity and slight solubility in water. The use of acoustic waves to remobilize trapped ganglia in subsurface porous systems have some advantages over pre-existing solutions including eliminating the bypassing effect and new environmental hazards. Designing an effective acoustically assisted remediation method for such purposes relies on understanding the underlying mechanisms and developing validated models. In this work, pore-scale microfluidic experiments were run to investigate the interplay between break-up and remobilization under sonication at different levels of flow rate and wettability conditions. Based on the experimental observation and pore-scale physical characteristics, a pore network model was developed and verified against the experimental results. Such a model was developed based on a two-dimensional network and scaled up to three-dimensional networks. In the experiments, processing of two-dimensional images showed that acoustic waves can remobilize trapped ganglia. The other observed effect of vibration is to break up blobs and reduce the mean ganglia size. Recovery enhancements were greater in hydrophilic micromodels as compared to hydrophobic system. A strong correlation was found between the remobilization and breakup indicating that the trapped ganglia are breaking up due to acoustic stimulation firstly and then a background viscous force may get them flowing under the new generated fluid distribution. In modeling, the simulation results of residual saturation reasonably matched with experimental observations. The differences between the prediction by the model and the experimental data at verification points is less than 2% for data before and after the acoustic excitation. The transitions from three-dimensional simulations were used to propose a modified capillary number. This study gives a better understanding of the mechanisms behind the effect of acoustic waves in porous media and provides a predictive tool for evaluating enhancement in fluid displacement.

Effects of freeze-thaw cycles on methanogenic hydrocarbon degradation: Experiment and modeling

Cold regions are warming much faster than the global average, resulting in more frequent and intense freeze-thaw cycles (FTCs) in soils. In hydrocarbon-contaminated soils, FTCs modify the biogeochemical and physical processes controlling petroleum hydrocarbon (PHC) biodegradation and the associated generation of methane (CH₄) and carbon dioxide (CO₂). Thus, understanding the effects of FTCs on the biodegradation of PHCs is critical for environmental risk assessment and the design of remediation strategies for contaminated soils in cold regions. In this study, we developed a diffusion-reaction model that accounts for the effects of FTCs on toluene biodegradation, including methanogenic biodegradation. The model is verified against data generated in a 215 day-long batch experiment with soil collected from a PHC contaminated site in Ontario, Canada. The fully saturated soil incubations with six different treatments were exposed to successive 4-week FTCs, with temperatures oscillating between -10 °C and +15 °C, under anoxic conditions to stimulate methanogenic biodegradation. We measured the headspace concentrations and ¹³C isotope compositions of CH₄ and CO₂ and analyzed the porewater for pH, acetate, dissolved organic and inorganic carbon, and toluene. The numerical model represents solute diffusion, volatilization, sorption, as well as a reaction network of 13 biogeochemical processes. The model successfully simulates the soil porewater and headspace concentration time series data by representing the temperature dependencies of microbial reaction and gas diffusion rates during FTCs. According to the model results, the observed increases in the headspace concentrations of CH₄ and CO₂ by 87% and 136%, respectively, following toluene addition are explained by toluene fermentation and subsequent methanogenesis reactions. The experiment and the numerical simulation show that methanogenic degradation is the primary toluene attenuation mechanism under the electron acceptor-limited conditions experienced by the soil samples, representing 74% of the attenuation, with sorption contributing to 11%, and evaporation contributing to 15%. Also, the model-predicted contribution of acetate-based methanogenesis to total produced CH₄ agrees with that derived from the ¹³C isotope data. The freezing-induced soil matrix organic carbon release is considered as an important process causing DOC increase following each freezing period according to the calculations of carbon balance and SUVA index. The simulation results of a no FTC scenario indicate that, in the absence of FTCs, CO₂ and CH₄ generation would decrease by 29% and 26%, respectively, and that toluene would be biodegraded 23% faster than in the FTC scenario. Because our modeling approach represents the dominant processes controlling PHC biodegradation and the associated CH₄ and CO₂ fluxes, it can be used to analyze the sensitivity of these processes to FTC frequency and duration driven by temperature fluctuations.

Dissolution and remobilization of NAPL in surfactant-enhanced aquifer remediation from microscopic scale simulations

In this paper, the dissolution and mobilization of non-aqueous phase liquid (NAPL) blobs in the Surfactant-Enhanced Aquifer Remediation (SEAR) process were upscaled using dynamic pore network modeling (PNM) of three-dimensional and unstructured networks. We considered corner flow and micro-flow mechanisms including snap-off and piston-like movement for two-phase flow. Moreover, NAPL entrapment and remobilization were evaluated using force analysis to develop the capillary desaturation curve (CDC) and predict the onset of remobilization. The corner diffusion mechanism was also applied in the modeling of interphase mass transfer to represent NAPL dissolution as the dominant mass transfer process. In addition, the effect of pore-scale heterogeneity on mass transfer rate coefficient and recovered residual NAPL was considered in the simulations. Sodium dodecyl sulfate (SDS) and Triton X-100 were used as the surfactant for the SEAR process. The results indicate that although surfactants enhance NAPL recovery during two-phase flow, surfactant-enhanced remediation of residual NAPL through dissolution is highly dependent on surfactant type. When SDS ─as a surfactant with high critical micelle concentration (CMC) and low micelle partition coefficient (K_m)─ was injected into a NAPL contaminated site, the mass transfer rate coefficient decreased (due to considerable changes in interface chemical potentials) which leads to a significant reduction in NAPL recovery after the end of two-phase flow. In contrast, Triton X-100 (with low CMC and high K_m) improved NAPL recovery, by enhancing solubility at surfactant concentrations greater than CMC which overcompensates the interphase mass transfer reduction.