Pore-Scale Mechanisms of Solid Phase Emergence During DNAPL Remediation by Chemical Oxidation.

In situ chemical oxidation (ISCO) has proven successful in the remediation of aquifers contaminated with dense nonaqueous phase liquids (DNAPLs). However, the treatment efficiency can often be hampered by the formation of solids or gas, reducing the contact between remediation agents and residual DNAPLs. To further improve the efficiency of ISCO, fundamental knowledge is needed about the complex multiphase flow and reactive transport processes as new solid and fluid phases emerge at the microscale. Here, via microfluidic experiments, we study the pore-scale dynamics of trichloroethylene degradation by permanganate. We visualize how the remediation evolves under the influence of solid phase emergence and explore the roles of injection rate, oxidant concentration, and stabilization supplement. Combining image processing, pressure analysis, and stoichiometry calculations, we provide comprehensive descriptions of the oxidant concentration-dependent growth patterns of the solid phase and their impact on the remediation efficiency. We further corroborate the stabilization mechanism provided by phosphate supplement, which is effective in inhibiting solid phase generation and thus highly beneficial for the oxidation remediation. This work elucidates the pore-scale mechanisms during remediation of chlorinated solvents with a particular context in the solid phase production and the associated effects, which is of general significance to understanding various processes in natural and engineered systems involving solid phase emergence or aggregation phenomena, such as groundwater and soil remediation.

Dissolution of dense non-aqueous phase liquids in vertical fractures: effect of finger residuals and dead-end pools.

Understanding the dissolution behavior of dense non-aqueous phase liquids (DNAPLs) in rock fractures under different entrapment conditions is important for remediation activities and any related predictive modeling. This study investigates DNAPL dissolution in variable aperture fractures under two important entrapment configurations, namely, entrapped residual blobs from gravity fingering and pooling in a dead-end fracture. We performed a physical dissolution experiment of residual DNAPL blobs in a vertical analog fracture using light transmission techniques. A high-resolution mechanistic (physically-based) numerical model has been developed which is shown to excellently reproduce the experimentally observed DNAPL dissolution. We subsequently applied the model to simulate dissolution of the residual blobs under different water flushing velocities. The simulated relationship between the Sherwood number Sh and Peclet number Pe could be well fitted with a simple power-law function (Sh=1.43Pe°·4³). To investigate mass transfer from dead-end pools, another type of trapping in rock fractures, entrapment and dissolution of DNAPL in a vertical dead-end fracture was simulated. As the entrapped pool dissolves, the depth of the interface between the DNAPL and the flowing water increases linearly with decreasing DNAPL saturation. The interfacial area remains more or less constant as DNAPL saturation decreases, unlike in the case of residual DNAPL blobs. The decreasing depth of the contact interface changes the flow field and causes decreasing water flow velocity above the top of the DNAPL pool, suggesting the dependence of the mass transfer rate on the depth of the interface, or alternatively, the remaining mass percentage in the fracture. Simulation results show that the resultant Sherwood number Sh is significantly smaller than in the case of residual blobs for any given Peclet number, indicating slower mass transfer. The results also show that the Sh can be well fitted with a power-law function of Pe and remaining mass percentage. The obtained relationships of dimensionless groups concerning the mass transfer characteristics at the level of individual fractures can be further used in predictive modeling of dissolution at a larger (fracture network) scale.

Effects of single-fracture aperture statistics on entrapment, dissolution and source depletion behavior of dense non-aqueous phase liquids.

Understanding of the entrapment and dissolution behavior of dense non-aqueous phase liquids (DNAPLs) in single fractures is important for modeling contaminant flux generation from fractured sites. Here a systematic numerical study is presented to investigate the effect of fracture aperture statistics on DNAPL migration, entrapment and dissolution within individual, variable-aperture fractures. Both fractures with open and closed bottom boundaries were considered. For the simulation a continuum-based two-phase model was used with a capillary pressure function which calculates the entry pressure based on the local aperture. Prior to application the model was compared against the invasion percolation approach and found more suitable for the present study, in particular as it allows a more versatile presentation of boundary conditions. The results showed that increasing aperture standard deviation and/or decreasing correlation length lead to larger amounts of entrapped DNAPL (due to the fact that larger standard deviation produces more distinct contrast between small and large aperture regions and the fact that longer correlation length provides more possible channels through the fracture) as well as larger maximum and average sizes of DNAPL blobs, and subsequently lead to longer times for complete dissolution. To understand the relationship between the solute flux and the remaining mass, a simplified source depletion function which links the outflow concentration to the DNAPL saturation was found adequate to describe the dissolution process for the case where the bottom boundary is open for DNAPL migration and thus the DNAPL does not accumulate to form a pool. The parameters in this function were not very sensitive to variations in correlation length but were sensitive to aperture standard deviation. The same average entrapped DNAPL saturation produced considerably smaller solute concentrations in cases with larger aperture variability due to the larger average size of DNAPL blobs (i.e., smaller contact area for DNAPL dissolution). Boundary conditions had a significant impact on DNAPL entrapment and dissolution. A closed boundary at the bottom led to DNAPL pooling (i.e., large continuous blobs) which causes significant tailing in the dissolution breakthrough curve due to water bypassing.