SEARCH · Search NASA

Results for “Drainage”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Continental rift evolution and drainage reorganization along the Dead Sea rift since the Miocene

The Dead Sea fault is a section of the Arabian-African plate boundary. Widespread field relations indicate that three major drainage systems (stages) occupied the landscape west of the Dead Sea fault since its initiation at ca. 20 Ma. Specifically, (1) an early to middle Miocene drainage system, only minorly reconfigured by the fault (all sediments of this system belong to the Hazeva Formation); (2) a late Miocene to early Pleistocene fault-parallel drainage system named Paran-Neqarot (all sediments of this system belong to the Arava and Zehiha Formations; and (3) the early Pleistocene to present drainage configuration. The temporal and spatial frameworks of drainage stage 1 are generally constrained by radiometric dating of interfingering volcanic units, and the onset and temporal and spatial frameworks of drainage stage 3 are well constrained by cosmogenic 10 Be surface exposure ages. The timing and longevity of the rift-parallel drainage (stage 2) have until now been evasive to direct dating. The overall time gap between stages 1 and 3 is ~12–13 million years. Thus, an early age (within this time gap) of stage 2 would imply an immediate response of drainage reorganization to rift tectonics, while a later age of this drainage system would imply a delayed response. We present 11 10 Be- 26 Al cosmogenic burial ages of alluvial and colluvial units related to the fault-parallel drainage system (stage 2), which collectively constrain the time of deposition of the Arava Formation sediments in the central Negev to ca. 8 Ma. The general lack of stratigraphic order, together with the large dispersion of ages both across and within the sampling sites, attests to significant recycling of sediments from drainage stage 1 into the Arava Formation deposits. The termination of Arava and Zehiha Formation sediment deposition at ca. 1.8 Ma was determined previously using cosmogenic exposure ages of desert pavements that cover the formations. Combining the previously published data with our new data, we established the longevity and character of the Paran-Neqarot drainage system. In conlcusion, this framework highlights the temporal aspect of drainage system build-up and collapse as it responded to transform and extensional plate boundary tectonics during the Neogene.

58 GEOSCIENCES↗

DOE OSTI2021

New insights into the drainage of inundated ice-wedge polygons using fundamental hydrologic principles

Abstract. The pathways and timing of drainage from the inundated centers of ice-wedge polygons in a warming climate have important implications for carbon flushing, advective heat transport, and transitions from methane to carbon dioxide dominated emissions. Here, we expand on previous research using a recently developed analytical model of drainage from a low-centered polygon. Specifically, we perform (1) a calibration to field data identifying necessary model refinements and (2) a rigorous model sensitivity analysis that expands on previously published indications of polygon drainage characteristics. This research provides intuition on inundated polygon drainage by presenting the first in-depth analysis of drainage within a polygon based on hydrogeological first principles. We verify a recently developed analytical solution of polygon drainage through a calibration to a season of field measurements. Due to the parsimony of the model, providing the potential that it could fail, we identify the minimum necessary refinements that allow the model to match water levels measured in a low-centered polygon. We find that (1) the measured precipitation must be increased by a factor of around 2.2, and (2) the vertical soil hydraulic conductivity must decrease with increasing thaw depth. Model refinement (1) accounts for runoff from rims into the ice-wedge polygon pond during precipitation events and possible rain gauge undercatch, while refinement (2) accounts for the decreasing permeability of deeper soil layers. The calibration to field measurements supports the validity of the model, indicating that it is able to represent ice-wedge polygon drainage dynamics. We then use the analytical solution in non-dimensional form to provide a baseline for the effects of polygon aspect ratios (radius to thaw depth) and coefficient of hydraulic conductivity anisotropy (horizontal to vertical hydraulic conductivity) on drainage pathways and temporal depletion of ponded water from inundated ice-wedge polygon centers. By varying the polygon aspect ratio, we evaluate the relative effect of polygon size (width), inter-annual increases in active-layer thickness, and seasonal increases in thaw depth on drainage. The results of our sensitivity analysis rigorously confirm a previous analysis indicating that most drainage through the active layer occurs along an annular region of the polygon center near the rims. This has important implications for transport of nutrients (such as dissolved organic carbon) and advection of heat towards ice-wedge tops. We also provide a comprehensive investigation of the effect of polygon aspect ratio and anisotropy on drainage timing and patterns, expanding on previously published research. Our results indicate that polygons with large aspect ratios and high anisotropy will have the most distributed drainage, while polygons with large aspect ratios and low anisotropy will have their drainage most focused near their periphery and will drain most slowly. Polygons with small aspect ratios and high anisotropy will drain most quickly. These results, based on parametric investigation of idealized scenarios, provide a baseline for further research considering the geometric and hydraulic complexities of ice-wedge polygons.

Harp, Dylan R. (ORCID:0000000197778000)↗

DOE OSTI2024

Modeling the Effects of Artificial Drainage on Agriculture-Dominated Watersheds Using a Fully Distributed Integrated Hydrology Model

In agriculture-dominated watersheds where natural drainage is poor, agricultural ditches (narrow engineered channels) and tile drains (perforated pipes) are widely employed to enhance surface and subsurface drainage, respectively. Despite their relatively small scale, these features exert substantial control over the hydro-biogeochemical function of watersheds and their effects need to be represented in the models. We introduce a novel strategy to incorporate the effects of artificial agricultural drainage into a fully distributed basin-scale integrated surface-subsurface hydrology models. In our approach, narrow agriculture ditches for surface drainage are resolved efficiently using ditch-aligned computational meshes that are hydrologically conditioned to ensure connectivity in the stream/ditch network. For tile drainage in the subsurface, we use the physically based Hooghoudt's drainage equation as a subgrid model and route the water drained through tiles to the nearest ditch. Without site-specific calibration, this model reproduced observed streamflow in the Portage River Watershed (>1,000 km 2 ) as recorded by a USGS gauge with good accuracy (normalized KGE = 0.81) and outperformed a calibrated SWAT model (normalized KGE = 0.68). Numerical experiments confirm that artificial drainage reduces surface inundations and effectively controls the water table. At the watershed scale, artificial drainage increases baseflow but has little effect on watershed discharges above the 90th percentile. The strong physical underpinnings and reduced need for calibration allow us to study the impacts of artificial drainage on distributed hydrological response in terms of fluxes and states and provide a platform for investigating watershed-scale nutrient transport.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Data for Soil Oxygen Dynamics: A Key Mediator of Tile Drainage Impacts on Coupled Hydrological, Biogeochemical, and Crop Systems

Tile drainage removes excess water and is an essential, widely adopted management practice to enhance crop productivity in the US Midwest and throughout the world. Tile drainage has been shown to significantly change hydrological and biogeochemical cycles by lowering the water table and reducing the residence time of soil water, although examining the complex interactions and feedbacks in an integrated hydrology–biogeochemistry–crop system remains elusive. Oxygen dynamics are critical to unraveling these interactions and have been ignored or oversimplified in existing models. Understanding these impacts is essential, particularly so because tile drainage has been highlighted as an adaptation under projected wetter springs and drier summers in the changing climate in the US Midwest. We used the ecosys model that uniquely incorporates first-principle soil oxygen dynamics and crop oxygen uptake mechanisms to quantify the impacts of tile drainage on hydrological and biogeochemical cycles and crop growth in corn–soybean rotation fields. The model was validated with data from a multi-treatment, multi-year experiment in Washington, IA. The relative root mean square error (rRMSE) for the corn and soybean yield in validation is 5.66 % and 12.57 %, respectively. The Pearson coefficient (r) of the monthly tile flow during the growing season is 0.78. Plant oxygen stress turns out as an emergent property of the equilibrium between the soil oxygen supply and biological demand. The impact of tile drainage on the system is achieved through a series of coupled feedback mechanisms. The model results show that tile drainage reduces the soil water content and enhances soil oxygenation. It additionally increases the subsurface discharge and elevates inorganic nitrogen leaching, with seasonal variations influenced by climate and crop phenology. The improved aerobic condition alleviates crop oxygen stress during wet springs, thereby promoting crop root growth during the early growth stage. The development of greater root density, in turn, mitigates water stress during dry summers, leading to an overall increase in the crop yield by ∼6 %. These functions indicate the potential of tile drainage in bolstering crop resilience to climate change and the use of this modeling tool for large-scale assessments of tile drainage. The model reveals the underlying causal mechanisms that drive the agroecosystem response to drainage on the coupled hydrology, biogeochemistry, and crop system dynamics.

Modeling↗

DOE OSTI2025

Soil oxygen dynamics: a key mediator of tile drainage impacts on coupled hydrological, biogeochemical, and crop systems

60 APPLIED LIFE SCIENCES↗

DOE OSTI2025

Computationally efficient subglacial drainage modelling using Gaussian process emulators: GlaDS-GP v1.0

Subglacial drainage models represent water flow at the ice–bed interface through coupled distributed and channelized systems to determine water pressure, discharge, and drainage system geometry. While they are used to understand processes such as the relationship between surface melt and ice flow, the number of uncertain model parameters and the computational cost of running models makes it difficult to adequately explore the high-dimensional parameter space and evaluate uncertainty in model predictions. Here, we develop Gaussian process (GP) emulators that make fast predictions with associated uncertainty of subglacial drainage model outputs. Using a truncated principal component (PC) basis representation, we construct a GP emulator for diurnally averaged subglacial water pressure. We also explore emulation of scalar variables describing drainage efficiency and configuration. We train the emulators using ensembles of up to 512 simulations varying eight parameters of the Glacier Drainage System (GlaDS) model on a synthetic domain intended to represent an ice-sheet margin. The emulators make predictions ∼ 1000 times faster than GlaDS simulations, with errors <3 % for the water pressure field and ∼ 5 %–9 % for drainage efficiency and configuration. We apply the emulators to explore the eight-dimensional parameter space by computing variance-based parameter sensitivity indices, finding that three parameters (ice flow coefficient, bed bump aspect ratio, and the subglacial cavity system conductivity) explain 90 % of the variance in modelled water pressure in response to parameter changes. The GP emulator approach described here is well suited to integrating observational data with models to make calibrated, credible predictions of subglacial drainage.

58 GEOSCIENCES↗

DOE OSTI2021

Opportunities for Treatment and Reuse of Agricultural Drainage in the United States

Irrigation accounts for 42% of the total freshwater withdrawals in the United States. Climate change, the pressure of a growing population, degrading water quality, and increased competition from other sectors could constrain continuous supply to meet future agricultural water demand. This study presents an evaluation framework to assess the potential reuse of agricultural drainage water for crop irrigation. Using a regional approach, we review the current state of agricultural drainage treatment and reuse and the institutional, economic, and other barriers that can influence the reuse decision. In the 31 eastern states, agricultural drainage contains valuable nutrients that can be reused for irrigation with minimal treatment, while the 17 western states struggle with large volumes of saline drainage that can contain constituents of concern (e.g., selenium), preventing reuse without treatment. Using a new decision-support tool called WaterTAP3, a potential treatment train for saline agricultural drainage was analyzed to identify treatment challenges, research needs, and the potential implementation at a larger scale. As demonstrated by our case study, desalination of agricultural drainage is costly and energy intensive and will require sizable investments to fully develop and optimize technologies as well as manage the generated waste and brine.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2023

Coupled effects of future rainfall and land use on urban stormwater drainage system in Tampa, Florida (USA)

Future rainfall and land use are two important factors for flood management since they both directly and indirectly affect the functionality and performance of urban stormwater drainage infrastructures. Although some studies have already examined the coupled effects of future rainfall and land use conditions, they have not concluded how such effects on both surface runoff and flood areas would change with regional surface elevation variation and infrastructure conditions. This paper analyzed the coupled effects of future rainfall and land use on urban drainage systems in terms of surface runoff quantity and flood area changes using EPA SWMM. Future downscaled and bias-corrected precipitation projections for 2040–2060 and 2080–2099 from “Downscaled CMIP3 and CMIP5 Climate and Hydrology Projections” Archive and future land use for 2050 s and 2090 s from EPA ICLUS project were used for the City of Tampa in Florida (USA). It was shown that potential runoff volume and flood area changes ranged from –40% to 160%, and –40% to 400% due to rainfall change, while potential changes induced by land use change ranged from 0% to 3.5%, and 0% to 18% respectively. Additionally, this study examined the coupled effects of future rainfall and land use changes on flood area variations, considering the capacity of drainage infrastructure and elevation variation of study site. Results showed that flood area variations are not solely influenced by runoff quantity changes due to the two external drivers, but also by other factors. Specifically, in low-relief areas with inadequate drainage infrastructure and less elevation variation, these two drivers have no additive effects on the percentage of flood area changes. In contrast, in areas with adequate drainage infrastructure and greater surface elevation variation, there are additive effects from rainfall and land use changes on flood area expansions. These findings have important implications for stakeholders involved in city planning and flooding control, particularly for urban areas facing future precipitation and urbanization challenges.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Potential of water quality wetlands to mitigate habitat losses from agricultural drainage modernization

Given widespread biodiversity declines, a growing global human population, and demands to improve water quality, there is an immediate need to explore land management solutions that support multiple ecosystem services. Agricultural water quality wetlands designed to provide both water quality benefits and wetland and grassland habitat are an emerging restoration solution that may reverse habitat declines in intensive agricultural areas. Installation of water quality wetlands in the Upper Midwest, USA, when considered alongside the repair and modification of aging agricultural tile drainage infrastructure, is a likely scenario that may mitigate nutrient pollution exported from agricultural systems and improve crop yields. The capacity of water quality wetlands to provide habitat within the wetland pool and the surrounding grassland is not well-studied, particularly with respect to potential habitat changes resulting from drainage infrastructure upgrades. For the current study, we produced spatially explicit models of 37 catchments distributed throughout an important region for agriculture and biodiversity, the Des Moines Lobe of Iowa. Four scenarios were considered - with and without improved drainage and with and without water quality wetlands - to estimate the net potential habitat implications of these scenarios for amphibians, grassland birds, and wild bees. Model results indicate that drainage modification alone will likely result in moderate direct losses of suitable amphibian habitat and large declines in overall habitat quality. However, inclusion of water quality wetlands at the catchment scale may mitigate these amphibian habitat losses while also increasing grassland bird and pollinator habitat. In conclusion, the impacts of water quality wetlands and drainage modernization on waterfowl in the region require additional study.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Calibration of V-Notch and Compound Weirs for Subsurface Drainage Water Level Control Structures

Highlights Accurate discharge estimation is important when evaluating edge-of-field conservation practices. V-notch weir equations were developed for three sizes of subsurface drainage water level control structures. Compound weir equations were developed for subsurface drainage water level control structures. The compound weir equation accurately estimates discharge for flows within and overtopping the V-notch. Abstract.Numerous edge-of-field conservation practices use subsurface drainage water level control structures to monitor water levels and estimate discharge. In a control structure, procedures for calculating discharge when flow depth (head) exceeds the V-notch depth and overflows in the rectangular portion of the compound weir (CW) are ambiguous. In this study, we developed calibration equations for V-notch weirs in Agri Drain inline water level control structures of different sizes for flows within the V-notch and overtopping flow events. The discharge equation for overtopping events (Q CW , L s -1 ) was determined as: Q CW = a 1 (h b1 -h 1 b1 )+a 2 (W e -W v )h 1 b2 , where h and h 1 are heads above vertex/bottom and top of V-notch (cm), respectively, W is the effective crest width of rectangular weir (cm), W v is the top width of V-notch (cm), a 1 and b 1 are parameters for V-notch weir obtained by calibration, and a 2 and b 2 are calibration parameters for rectangular weir obtained from literature. Results were compared with a weir equation available in the literature (Q V+R ), which combines a V-notch equation with a head equal to V-depth and a rectangular weir equation for flow above V-depth. Discharge at overflow was estimated with high accuracy with Q CW, whereas Q V+R underestimated discharge (e.g., PBIAS of 0.67% vs. 17.82% for a 15.2 cm structure). An example using Q V+R resulted in a 14% lower annual estimation of nitrate-N load diverted to a saturated buffer than Q CW due to underestimation of drainage discharge during overflow events. Results suggest that the developed equation (Q CW ) accurately estimates discharge and will thus improve the estimated N load compared to Q V+R . Keywords: Compound weir, Flow monitoring, Subsurface drainage, V-notch weir, Water level control structure, Weir calibration.

Agriculture↗

DOE OSTI2025

Emulator-based Bayesian calibration of a subglacial drainage model

Subglacial drainage models, often motivated by the relationship between hydrology and ice flow, sensitively depend on numerous unconstrained parameters. We explore using borehole water-pressure time series to calibrate the uncertain parameters of a popular subglacial drainage model, taking a Bayesian perspective to quantify the uncertainty in parameter estimates and in the calibrated model predictions. To reduce the computation time associated with Markov Chain Monte Carlo sampling, we construct a fast Gaussian process emulator to stand in for the subglacial drainage model. We first carry out a calibration experiment using synthetic observations consisting of model simulations with hidden parameter values as a demonstration of the method. Using real borehole water pressures measured in western Greenland, we find meaningful constraints on four of the eight model parameters and a factor-of-three reduction in uncertainty of the calibrated model predictions. These experiments illustrate Gaussian process-based Bayesian inference as a useful tool for calibration and uncertainty quantification of complex glaciological models using field data. However, significant differences between the calibrated model and the borehole data suggest that structural limitations of the model, rather than poorly constrained parameters or computational cost, remain the most important constraint on subglacial drainage modelling.

58 GEOSCIENCES↗

DOE OSTI2024

Transforming Drainage Research Data (USDA-NIFA Award No. 2015-68007-23193)

This dataset contains research data compiled by the “Managing Water for Increased Resiliency of Drained Agricultural Landscapes” project a.k.a. Transforming Drainage. This project was funded from 2015-2021 by the United States Department of Agriculture, National Institute of Food and Agriculture (USDA-NIFA, Award No. 2015-68007-23193). Data are also available from a separate web-accessible application (drainagedata.org). At drainagedata.org, users can visualize the data with customized tools, query based on specific sites and measurements of interest, and access site photographs, maps, summaries, and publications. Additional data or edits made following the publication of this data here at USDA NAL Ag Data Commons will be posted under the Versions tab on drainagedata.org. These data began in 1996 and include plot- and field-level measurements for 39 experiments across the Midwest and North Carolina. Practices studied include controlled drainage, drainage water recycling, and saturated buffers. In total, 219 variables are reported and span 207 site-years for tile drainage, 154 for nitrate-N load, 181 for water quality, 92 for water table, and 201 for crop yield.

Modeling↗

DOE OSTI2024

Modeling the Effects of Artificial Drainage on Agriculture-dominated Watersheds using a Fully Distributed Integrated Hydrology Model: Datasets, scripts, model files

This model-data archive supports the research paper that demonstrates the integration of agricultural drainage features—specifically, narrow engineered ditches and tile drains—into a fully distributed, basin-scale integrated surface-subsurface hydrology model (ISSHM), Amanzi-ATS. The model employs innovative computational meshes aligned with agricultural ditches and incorporates the physically based Hooghoudt's drainage equation to simulate tile drainage, offering a novel strategy that enhances the accuracy of hydrological simulations.The archived dataset includes input parameters, model configurations, and select simulation outputs for the Amanzi-ATS model that successfully captured the streamflow patterns in the Portage River Watershed as validated by USGS gauge readings. Jupyter notebook for the preparation of model inputs and post-processing of outputs are also included. The model's predictive performance achieved a normalized Kling-Gupta Efficiency (KGE) of 0.81, surpassing SWAT without the necessity for site-specific calibration.The Amanzi-ATS model presented in this modeL-data archive allows for numerical experiments to explore the shifts in the flow structure under different drainage scenarios. As a tool for advancing the understanding of distributed hydrological responses and nutrient cycling, this archived model provides valuable insights for researchers, modelers, and decision-makers involved in watershed management and environmental modeling.The Watershed Workflow package is implemented in Python3. The Jupyter notebooks can be executed through multiple open-source tools, for example, Anaconda Jupyter Lab, VS Studio Code, etc. Other data files include CSV and HDF5 files, which can be read through Python scripts. The input files for the ATS model, open-source integrated hydrology, and transport model, are in XML format and can be edited in any commonly used text editors.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Recovery of Rare Earth Element from Acid Mine Drainage Using Organo-Phosphorus Extractants and Ionic Liquids

Acid mine drainage is a legacy environmental issue and one of the largest pollutants in many mining districts throughout the world. In prior work, the authors have developed a process for the recovery of critical materials, including the rare earth elements, from acid mine drainage using a preconcentration step followed by solvent extraction as a concentration and purification technology. As part of the downstream technology development efforts, we have synthesized a suite of ionic liquid extractants that facilitate greater separation factors leading to lower capital costs and reduced environmental impacts. This article provides a comparison of the conventional extractants D2EHPA, EHEHPA and C572 with their respective ionic liquids [c101][D2EHP,c101][EHEHP] and [c101][C572] for the recovery of rare earth elements from acid mine drainage. In the study, laboratory-scale, multi-contact solvent extraction tests were conducted at high and low extractant/dosages. The results show that the ionic liquids varied in performance, with [c101][D2EHP] and [c101][EHEHP] performing poorer than their conventional counterparts and [c101][c572] performing better. Recommendations for further study on [c101][c572] include stripping tests, continuous pilot testing, and techno-economic analysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Experimental study of three-dimensional CO 2 -water drainage and fracture-matrix interactions in fractured porous media

Estimation of CO 2 storage capacity in fractured porous reservoirs requires a better understanding of the CO 2 -water flow fundamentals in fracture-matrix systems. There are few core-flood experiments on three-dimensional CO 2 -water drainage in a fracture-matrix system, yet none have examined the impacts of the capillary continuity of fractures and fracture-matrix interactions. In this study, twelve drainage experiments were conducted in four fracture-matrix columns, each comprising a vertical stack of a cylindrical rock core, a filter paper serving as an analogue to a horizontal fracture, and a ceramic plate. Here, the core sample was surrounded by open space at the top and circumferential sides to model fractures, allowing for 3-D CO 2 -water drainage in the rock core and displaced water draining across the horizontal fracture. X-ray computed tomography was conducted to visualize the dynamic invasion/drainage processes in four rock core samples showing contrasts in anisotropy, permeability, and heterogeneity. Experimental results show (1) the equilibrium CO 2 saturations in the rock matrix vary from 0.10 to 0.60 at controlled capillary pressures up to 200 kPa, (2) the CO 2 saturations in the matrix increase with water saturation in the horizontal fracture resulting in better capillary continuity for water to drain across the fracture, and (3) the fracture water saturation can be enhanced by the non-uniform fracture capillary pressure and countercurrent flow of CO 2 and water across the fracture-matrix interface. In the case of high fracture water saturation, matrix CO 2 saturation largely depends on matrix anisotropy and heterogeneity. The core-scale experimental results contribute to understand the fracture-matrix interactions and CO 2 storage efficiency in fractured porous media.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Improved representation of laminar and turbulent sheet flow in subglacial drainage models

Subglacial hydrology models struggle to reproduce seasonal drainage patterns that are consistent with observed subglacial water pressures and surface velocities. We modify the standard sheet-flow parameterization within a coupled sheet–channel subglacial drainage model to smoothly transition between laminar and turbulent flow based on the locally computed Reynolds number in a physically consistent way (the ‘transition’ model). We compare the transition model to standard laminar and turbulent models to assess the role of the sheet-flow parameterization in reconciling observed and modelled water pressures under idealized and realistic forcing. Relative to the turbulent model, the laminar and transition models improve seasonal simulations by increasing winter water pressure and producing a more prominent late-summer water pressure minimum. In contrast to the laminar model, the transition model remains consistent with its own internal assumptions across all flow regimes. Based on the internal consistency of the transition model and its improved performance relative to the standard turbulent model, we recommend its use for transient simulations of subglacial drainage.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Pore-Scale Dynamics of Liquid CO2–Water Displacement in 2D Axisymmetric Porous Micromodels Under Strong Drainage and Weak Imbibition Conditions: High-Speed μPIV Measurements

Resolving pore-scale transient flow dynamics is crucial to understanding the physics underlying multiphase flow in porous media and informing large-scale predictive models. Surface properties of the porous matrix play an important role in controlling such physics, yet interfacial mechanisms remain poorly understood, in part due to a lack of direct observations. This study reports on an experimental investigation of the pore-scale flow dynamics of liquid CO 2 and water in two-dimensional (2D) circular porous micromodels with different surface characteristics employing high-speed microscopic particle image velocimetry (μPIV). The design of the micromodel minimized side boundary effects due to the limited size of the domain. The high-speed μPIV technique resolved the spatial and temporal dynamics of multiphase flow of CO 2 and water under reservoir-relevant conditions, for both drainage and imbibition scenarios. When CO 2 displaced water in a hydrophilic micromodel (i.e., drainage), unstable capillary fingering occurred and the pore flow was dominated by successive pore-scale burst events (i.e., Haines jumps). When the same experiment was repeated in a nearly neutral wetting micromodel (i.e., weak imbibition), flow instability and fluctuations were virtually eliminated, leading to a more compact displacement pattern. Energy balance analysis indicates that the conversion efficiency between surface energy and external work is less than 30%, and that kinetic energy is a disproportionately smaller contributor to the energy budget. This is true even during a Haines jump event, which induces velocities typically two orders of magnitude higher than the bulk velocity. These novel measurements further enabled direct observations of the meniscus displacement, revealing a significant alteration of the pore filling mechanisms during drainage and imbibition. While the former typically featured burst events, which often occur only at one of the several throats connecting a pore, the latter is typically dominated by a cooperative filling mechanism involving simultaneous invasion of a pore from multiple throats. This cooperative filling mechanism leads to merging of two interfaces and releases surface energy, causing instantaneous high-speed events that are similar, yet fundamentally different from, burst events. Finally, pore-scale velocity fields were statistically analyzed to provide a quantitative measure of the role of capillary effects in these pore flows.

Li, Yaofa↗

DOE OSTI2022

Metal (Fe, Cu, and As) transformation and association within secondary minerals in neutralized acid mine drainage characterized using X-ray absorption spectroscopy

The formation and transformation of secondary minerals in areas affected by neutralized acid mine drainage (AMD) determine the behavior of toxic elements. To better understand the binding form of toxic elements in the secondary minerals, we focused on the distribution of secondary precipitates and heavy metals in neutralized AMD near a typical Cu-polymetallic deposit located in central Tibet. Compared to background values, the average heavy metal (As, Cu, Co, Cu, Mo, Pb, and Zn) content is higher in sediments close to the neutralized AMD drain. X-ray absorption spectroscopy (XAS) characterization shows the Fe in the suspended particulate matter (SPM) first aggregate as small Fe(III) octahedra clusters and then precipitate as ferrihydrite along the river. The XAS results also show the by-products of neutralization-ferrihydrite and sulfate-are the primary forms of Fe and S in the sediments, respectively. Gypsum retains Cu in the sediments via structural fixation, surface adsorption, and surface coprecipitation. The difference of Cu-O coordination number between the SPM and the sediments reveals the Cu sulfate coprecipitation on the surface of gypsum accompanied by SPM transporting along the river. XAS results confirm that ferrihydrite act as an essential sink for heavy metals, especially As. As (III) mineral in the sediment near the neutralized AMD outfall is dissolved and oxidized to arsenate downstream, which then adsorb to ferrihydrite and/or form As(V)-Fe minerals along the river. Gypsum and ferrihydrite have the potential to release heavy metals with changes in water chemistry, which may pose a threat to the safety of drinking water for downstream residents. These results deepen the understanding of the influence of secondary minerals (especially gypsum and ferrihydrite) on the migration and transformation process of heavy metal pollutants in mine drainage system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗