Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distributed model predictive control”

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.

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

At least 451 records · Page 25

Temporal evolution of soil moisture statistical fractal and controls by soil texture and regional groundwater flow

Soil moisture statistical fractal is an important tool for downscaling remotely-sensed observations and has the potential to play a key role in multi-scale hydrologic modeling. The fractal was first introduced two decades ago, but relatively little is known regarding how its scaling exponents evolve in time in response to climatic forcings. Previous studies have neglected the process of moisture re-distribution due to regional groundwater flow. Here, we used a physically-based surface-subsurface processes model and numerical experiments to elucidate the patterns and controls of fractal temporal evolution in two U.S. Midwest basins. Groundwater flow was found to introduce large-scale spatial structure, thereby reducing the scaling exponents (τ), which has implications for the transferability of calibrated parameters to predict τ. However, the groundwater effects depend on complex interactions with other physical controls such as soil texture and land use. The fractal scaling exponents, while in general showing a seasonal mode that correlates with mean moisture content, display hysteresis after storm events that can be divided into three phases, consistent with literature findings: (a) wetting, (b) re-organizing, and (c) dry-down. Modeling experiments clearly show that the hysteresis is attributed to soil texture, whose "patchiness" is the primary contributing factor. We generalized phenomenological rules for the impacts of rainfall, soil texture, groundwater flow, and land use on τevolution. Grid resolution has a mild influence on the results and there is a strong correlation between predictions of τ from different resolutions. Overall, our results suggest that groundwater flow should be given more consideration in studies of the soil moisture statistical fractal, especially in regions with a shallow water table.

54 ENVIRONMENTAL SCIENCES↗

Generation of Subsurface Voids, Incubation Effect, and Formation of Nanoparticles in Short Pulse Laser Interactions with Bulk Metal Targets in Liquid: Molecular Dynamics Study

The ability of short pulse laser ablation in liquids to produce clean colloidal nanoparticles and unusual surface morphology has been employed in a broad range of practical applications. In this paper, we report the results of large-scale molecular dynamics simulations aimed at revealing the key processes that control the surface morphology and nanoparticle size distributions by pulsed laser ablation in liquids. The simulations of bulk Ag targets irradiated in water are performed with an advanced computational model combining a coarse-grained representation of liquid environment and an atomistic description of laser interaction with metal targets. For the irradiation conditions that correspond to the spallation regime in vacuum, the simulations predict that the water environment can prevent the complete separation of the spalled layer from the target, leading to the formation of large subsurface voids stabilized by rapid cooling and solidification. The subsequent irradiation of the laser-modified surface is found to result in a more efficient ablation and nanoparticle generation, thus suggesting the possibility of the incubation effect in multipulse laser ablation in liquids. The simulations performed at higher laser fluences that correspond to the phase explosion regime in vacuum reveal the accumulation of the ablation plume at the interface with the water environment and the formation of a hot metal layer. The water in contact with the metal layer is brought to the supercritical state and provides an environment suitable for nucleation and growth of small metal nanoparticles from metal atoms emitted from the hot metal layer. The metal layer itself has limited stability and can readily disintegrate into large (tens of nanometers) nanoparticles. The layer disintegration is facilitated by the Rayleigh–Taylor instability of the interface between the higher density metal layer decelerated by the pressure from the lighter supercritical water. The nanoparticles emerging from the layer disintegration are rapidly cooled and solidified due to the interaction with water environment, with a cooling rate of ~2 × 10 12 K/s observed in the simulations. The computational prediction of two distinct mechanisms of nanoparticle formation yielding nanoparticles with different characteristic sizes provides a plausible explanation for the experimental observations of bimodal nanoparticle size distributions in laser ablation in liquids. The ultrahigh cooling and solidification rates suggest the possibility for generation of nanoparticles featuring metastable phases and highly nonequilibrium structures.

Chemistry↗

Correlated Anion Disorder in Heteroanionic Cubic TiOF 2

Resolving anion configurations in heteroanionic materials is crucial for understanding and controlling their properties. For anion-disordered oxyfluorides, conventional Bragg diffraction cannot fully resolve the anionic structure, necessitating alternative structure determination methods. We have investigated the anionic structure of anion-disordered cubic (ReO 3 -type) TiOF 2 using X-ray pair distribution function (PDF), 19 F MAS NMR analysis, density functional theory (DFT), cluster expansion modeling, and genetic-algorithm structure prediction. Our computational data predict short-range anion ordering in TiOF 2 , characterized by predominant cis-[O 2 F 4 ] titanium coordination, resulting in correlated anion disorder at longer ranges. To validate our predictions, we generated partially disordered supercells using genetic-algorithm structure prediction and computed simulated X-ray PDF data and 19 F MAS NMR spectra, which we compared directly to experimental data. To construct our simulated 19 F NMR spectra, we derived new transformation functions for mapping calculated magnetic shieldings to predicted magnetic chemical shifts in titanium (oxy)fluorides, obtained by fitting DFT-calculated magnetic shieldings to previously published experimental chemical shift data for TiF 4 . We find good agreement between our simulated and experimental data, which supports our computationally predicted structural model and demonstrates the effectiveness of complementary experimental and computational techniques in resolving anionic structure in anion-disordered oxyfluorides. From additional DFT calculations, we predict that increasing anion disorder makes lithium intercalation more favorable by, on average, up to 2 eV, highlighting the significant effect of variations in short-range order on the intercalation properties of anion-disordered materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Accounting for hydrological and microbial processes on greenhouse gas budgets from river systems (Final Report)

River systems have a disproportionately large impact on carbon and nitrogen cycling on this planet. For instance, despite their much smaller size in area, it is estimated that global carbon dioxide emission from river systems exceed the emission from lakes and reservoirs, and are on par with the estimated oceanic and land carbon sinks. Closing carbon and nitrogen budgets in river systems requires knowledge of end member greenhouse gas (GHG) products (e.g. carbon dioxide, nitrous oxide, methane), which typically have not been measured in river systems. Moreover, while the production of these gasses is thought to be microbially derived, little is known about microbial activities in sediments of large rivers, despite their potentially high influence on biogeochemical budgets. Typically, when we think of rivers we visualize the channels that carve through the earth, but this only represents a small part of the hydrologic exchange that a river participates in. Rivers are more than conduits; they are dynamically exchanging water and materials with surrounding terrestrial and subsurface environments. In fact, 96% of the carbon dioxide produced within river systems occurs in hyporheic zones, or sediment regions beneath and alongside a river where groundwater and river water mix. In fact, it is widely recognized that microbial processes in hyporheic zones, also known as the river’s liver, strongly influence the fate of carbon, nutrients, and contaminants transported through watersheds. Despite this general knowledge, observational studies linking nitrogen and carbon gas production to microbial activity, and studies of overall gas emissions, are limited across river systems, but especially absent from larger rivers. To identify microbial processes responsible for GHG emissions from large rivers, and to quantify their environmental drivers and dependencies, we selected the Columbia River as a representative field site. This is the largest river west of the North American Continental Divide, and produces a drainage and discharge roughly double of the Nile in Egypt. This river is also an important source for agricultural irrigation and generates more hydroelectric power than any other North American river. Dams manage the release of water in the river, with these changes in river discharge patterns strongly impacting microbial activities in hyporheic zones. Specifically, we target a region down-river from the Priest Rapids Dam (WA) associated with the PNNL SFA, where water discharge from the dam is highly controlled, resulting in 1.5 m fluctuations in river stage daily. This dynamic river stage creates an ideal research opportunity for investigating hydrological perturbations in GHG production, uptake, and emission from large river systems. In this exploratory 1-year project, we conducted data-driven field research to provide observations needed to support predictions about GHG flux from river ecosystems. Using porewater dialysis and surface flux measurements, we observed where and when GHGs were produced and is the emission impacted by river stage. Using meta-omics tools, we determined what microbial processes control production and consumption of GHG in these systems. Finally, we provided key observations to parameterize models that could predict the effect of environmental conditions on GHG emission in the Columbia River near shore environment, at the scale that the modeled processes operate. Closing knowledge gaps in the linkages between river stage, microbial processes, and GHG production and emission, as well as the spatial distribution and ecosystem dependencies of these processes, is necessary to better parameterize land surface models. We provide some of the first, detailed GHG measurements paired microbial activity from large river systems. Our findings have direct relevance to many other large rivers with flows altered by hydroelectric energy, reservoirs, snowmelt, tides, or agriculture.

54 ENVIRONMENTAL SCIENCES↗

Prediction of wing aeroelastic effects on aircraft lift and pitching moment characteristics

The distribution of flight loads on an aircraft structure determine the lift and pitching moment characteristics of the aircraft. When the load distribution changes due to the aeroelastic response of the structure, the lift and pitching moment characteristics also change. An estimate of the effect of aeroelasticity on stability and control characteristics is often required for the development of aircraft simulation models of evaluation of flight characteristics. This presentation outlines a procedure for incorporating calculated linear aeroelastic effects into measured nonlinear lift and pitching moment data from wind tunnel tests. Results are presented which were obtained from applying this procedure to data for an aircraft with a very flexible transport type research wing. The procedure described is generally applicable to all types of aircraft.

Eckstrom, Clinton V.↗

Prediction of wing aeroelastic effects on aircraft life and pitching moment characteristics

The distribution of flight loads on an aircraft structure determine the lift and pitching moment characteristics of the aircraft. When the load distribution changes due to the aeroelastic response of the structure, the lift and pitching moment characteristics also change. An estimate of the effect of aeroelasticity on stability and control characteristics is often required for the development of aircraft simulation models of evaluation of flight characteristics. This presentation outlines a procedure for incorporating calculated linear aeroelastic effects into measured nonlinear lift and pitching moment data from wind tunnel tests. Results are presented which were obtained from applying this procedure to data for an aircraft with a very flexible transport type research wing. The procedure described is generally applicable to all types of aircraft.

Eckstrom, Clinton V.↗

Joint Geophysical and Numerical Insights of the Coupled Thermal‐Hydro‐Mechanical Processes During Heating in Salt

Abstract Salt offers an optimal medium for the permanent isolation of heat‐producing radioactive waste due to its impermeability, high thermal conductivity, and ability to close fractures through creep. A thorough understanding of the thermal‐hydrological‐mechanical (THM) processes, encompassing brine migration, is fundamental for secure radioactive waste disposal within salt formations. At the Waste Isolation Pilot Plant (WIPP), we conducted joint in situ geophysical monitoring experiments during active heating to investigate brine migration near excavations. This experiment incorporated electrical resistivity tomography (ERT) alongside high‐resolution fiber‐optic‐based distributed temperature sensing within a controlled heating experiment. Additionally, discrete element model (DEM) based numerical simulations were conducted to simulate THM processes during heating, providing a more mechanistic understanding of the coupled processes leading to the observed changes in geophysical measurements. During heating, resistivity shifts near the heater were reasonably explained by temperature effects. However, in more distant, cooler regions, the resistivity decrease exceeded predictions based solely on temperature. DEM simulations highlighted brine migration, propelled by pore pressure gradients, as the likely primary factor contributing to the additional resistivity decline beyond temperature effects. The comparison between the predicted ERT responses and observations was much improved when considering the effects of brine migration based on the DEM simulations. These geophysical and simulation findings shed light on brine migration in response to salt heating, enhancing our understanding of the coupled THM processes in salt for safe radioactive waste disposal.

58 GEOSCIENCES↗

Simulated Phase-dependent Spectra of Terrestrial Aquaplanets in M Dwarf Systems

Orbital phase-dependent variations in thermal emission and reflected stellar energy spectra can provide meaningful constraints on the climate states of terrestrial extrasolar planets orbiting M dwarf stars. Spatial distributions of water vapor, clouds, and surface ice are controlled by climate. In turn, water, in each of its thermodynamic phases, imposes significant modulations to thermal and reflected planetary spectra. Here we explore these characteristic spectral signals, based on 3D climate simulations of Earth-sized aquaplanets orbiting M dwarf stars near the habitable zone. By using 3D models, we can self-consistently predict surface temperatures and the location ofwater vapor, clouds, and surface ice in the climate system. Habitable zone planets in M dwarf systems are expected to be in synchronous rotation with their host star and thus present distinct differences in emitted and reflected energy fluxes depending on the observed hemisphere. Here we illustrate that icy, temperate, and incipient runaway greenhouse climate states exhibit phase-dependent spectral signals that enable their characterization.

Wolf, E. T.↗

Coupling Big Data Analytics and Reactive Transport Modeling for Cost-effective Groundwater Monitoring - 17163

This study presents an innovative approach for sustainable and cost-effective groundwater monitoring. This approach takes advantage of recent advances in various technologies: (1) in situ autonomous sensors, (2) big data analytics, and (3) parallel high-performance computing for flow and reactive transport modeling. In situ sensors are used to periodically measure the key variables (such as pH, redox potential, electrical conductivity, and groundwater level), which control contaminant mobility and the plume spatial and temporal distribution. Based on a limited number of groundwater sampling, the data analytics methods-data mining and machine learning-allow us to identify and quantify the correlations between the in situ-measured variables and contaminant concentrations, and also to detect significant changes associated with the plume mobility. In addition, a state-of-the-art parallel numerical flow and reactive transport simulator Amanzi, and uncertainty quantification software Agni are used to provide an improved physical and mechanistic understanding of the contaminated groundwater system behavior, and to predict the long-term plume distribution for optimizing and adapting the monitoring strategy. Amanzi and Agni were developed as part of the Advanced Simulation Capability for Environmental Management (ASCEM) program of the DOE Office of Environmental Management. Such modeling is critical, particularly, for assessing the impact of climate change and associated hydrological shifts. The developed approach is expected to significantly reduce the groundwater sampling frequency and associated cost. In addition, the real-time information on plume mobility serves as an early warning system, improving the resiliency of contaminated or potentially contaminated sites. We demonstrate this approach using as an example the Savannah River Site (SRS) F-Area, where groundwater is contaminated by various radionuclides, including uranium, tritium and technetium. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

An Investigation into the Effects of Swirl on the Performance and Emissions of an Opposed-Piston Two-Stroke Engine using Large Eddy Simulations

Opposed-piston two-stroke (OP-2S) engines have the potential to achieve higher thermal efficiency than a conventional four-stroke diesel engine. However, the uniflow scavenging process is difficult to control over a wider range of speed and loads due to its sensitivity to pressure dynamics, port timings, and port design. Specifically, the angle of the intake ports can be used to generate swirl which has implications for open and closed cycle effects. This study proposes an analysis of the effects of port angle on the in-cylinder flow distribution and combustion performance of an OP-2S using computational fluid dynamics engine. Large Eddy Simulation (LES) was used to model turbulence given its ability to predict in-cylinder mixing and cyclic variability. A three-cylinder model was validated to experimental data collected by Achates Power and the grid was verified using an LES quality approach from the literature. The model was used to simulate port angles from 12 to 29 degrees at constant pressure and temperature boundary conditions. Results indicated that the higher bulk swirl ratio generated by larger port angles tends to trap more internal residuals. This effect on the scavenging performance, combined with the larger trapped swirl ratio, also has a significant impact on the combustion performance in a two-stroke engine. It was concluded that there exists a tradeoff of efficiency and emissions that must be considered when increasing the port angle in a uniflow two-stroke engine. 1. Introduction Two-stroke engines were developed in the 1800s and are commonly used in practice today for lightweight power applications such as motorcycles or handheld power tools, as well as large-bore low speed engines for stationary power and marine applications. Two-stroke engines provide the distinct advantage of higher torque/power output but present challenges due to its coupled intake and exhaust process. Four-stroke engines have the benefit of using two full piston strokes to induct intake charge and to expel exhaust gases. Controlling the trapped charge composition in a twostroke engine via a scavenging process is crucial for reliable and efficient engine operation.

33 ADVANCED PROPULSION SYSTEMS↗

Numerical analysis of experimental studies of methane hydrate dissociation induced by depressurization in a sandy porous medium

Methane Hydrates (MHs) are a promising energy source abundantly available in nature. Understanding the complex processes of MH formation and dissociation is critical for the development of safe and efficient technologies for energy recovery. Many laboratory and numerical studies have investigated these processes using synthesized MH-bearing sediments. A near-universal issue encountered in these studies is the spatial heterogeneous hydrate distribution in the testing apparatus. In the absence of direct observations (e.g. using X-ray computed tomography) coupled with real time production data, the common assumption made in almost all numerical studies is a homogeneous distribution of the various phases. In an earlier study (Yin et al., 2018) that involved the numerical description of a set of experiments on MH-formation in sandy medium using the excess water method, we showed that spatially heterogeneous phase distribution is inevitable and significant. In the present study, we use as a starting point the results and observations at the end of the MH formation and seek to numerically reproduce the laboratory experiments of depressurization-induced dissociation of the spatially-heterogeneous MH distribution. This numerical study faithfully reproduces the geometry of the laboratory apparatus, the initial and boundary conditions of the system, and the parameters of the dissociation stimulus, capturing accurately all stages of the experimental process. Using inverse modelling (history-matching) that minimized deviations between the experimental observations and numerical predictions, we determined the values of all the important flow, thermal, and kinetic parameters that control the system behaviour, which yielded simulation results that were in excellent agreement with the measurements of key monitored variables, i.e. pressure, temperature, cumulative production of gas and water over time. We determined that at the onset of depressurization (when the pressure drop – the driving force of dissociation – is at its maximum), the rate of MH dissociation approaches that of an equilibrium reaction and is limited by the heat transfer from the system surroundings. Furthermore, as the effect of depressurization declines over time, the dissociation reaction becomes kinetically limited despite significant heat inflows from the boundaries, which lead to localized temperature increases in the reactor.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermal stress in degraded forests in the Brazilian Amazon Arc of Deforestation

Understanding thermal stress in tropical forests has taken on new urgency in light of accelerating climate change and expansion of deforestation and forest degradation. Degraded tropical forests in particular may be approaching critical temperature thresholds even more rapidly than intact forests, with implications for tree survival and ecosystem recovery. We investigate thermal stress in degraded tropical forests within the Brazilian Amazon Arc of Deforestation. Using land surface temperature data from the ECOsystem Spaceborne Thermal Radiometer Experiment on the international Space Station (ECOSTRESS), we compared canopy temperatures of intact, selectively logged, and burned forests in Feliz Natal, Mato Grosso, Brazil. Upper canopy temperatures in previously burned forests were 4.1% higher (mean = 36.5 °C) and 50.9% more variable compared to intact and logged forests, which showed remarkably similar temperature distributions (means of 34.9 °C and 35.1 °C, respectively). Modeled leaf temperature distributions based on canopy temperature measurements from one of the warmest days in a 2 year dry season record indicated that 87% of leaves in the warmest burned forest patches exceeded the temperature threshold where respiration surpasses photosynthesis, compared to approximately 72%–74% in intact and logged forests, respectively, in the same time period. After controlling for environmental factors, burned forests were predicted to be 2.6 °C warmer on average [1.39 °C–3.96 °C, 95% credible interval] than intact forests across a 5–40 m canopy height range. Burned forests showed modest thermal recovery over time, with temperatures decreasing by approximately 1.2 °C over a 30 year recovery period. In contrast, logged forests showed minimal thermal differences from intact forests (−0.353 on average, [−1.125–0.274, 95% credible interval]) and negligible change (0.15 °C) over the same timeframe. While the absolute probabilities of exceeding damaging thermal thresholds remain low across all forest types under current climate conditions, the probability of leaves exceeding temperatures that cause permanent leaf damage was ten times higher in burned forests, with implications for the future of burned forest regeneration in water-limited regions of the Amazon basin. In particular, the combination of higher mean temperatures, greater temperature variability, and more frequent exposure to damaging thermal thresholds implies that burned tropical forests will experience substantially higher mortality rates and slower biomass recovery compared to intact and selectively logged forests, especially in water-limited regions where trees cannot rely on evaporative cooling to moderate canopy temperatures.

Brazilian Amazon Arc of Deforestation↗

World Ocean Circulation Experiment

The oceans are an equal partner with the atmosphere in the global climate system. The World Ocean Circulation Experiment is presently being implemented to improve ocean models that are useful for climate prediction both by encouraging more model development but more importantly by providing quality data sets that can be used to force or to validate such models. WOCE is the first oceanographic experiment that plans to generate and to use multiparameter global ocean data sets. In order for WOCE to succeed, oceanographers must establish and learn to use more effective methods of assembling, quality controlling, manipulating and distributing oceanographic data.

Clarke, R. Allyn↗

Combining Generative Modeling and Advanced Control for Building Scenario Generation

Buildings make up a large portion of energy consumption in the U.S. today. Understanding their energy consumption patterns can improve their efficiency, but requires detailed models that rely on incomplete or unknown information. Previous work has shown that artificial intelligence (AI) can be used to predict missing information and even suggest upgrades to improve building efficiency. However, building upgrades may require undesirable upfront costs. Oppositely, advanced control could improve building efficiency with negligible upfront cost. To explore the tradeoffs between these two approaches, in this work we propose a workflow to compute optimal temperature setpoint schedules to minimize energy consumption and operational cost. Results show that modifying the temperature setpoints in a building using model predictive control (MPC) can effectively reduce its energy consumption and operational cost. This optimal operation cannot fully meet a desired goal. However, we show that by considering MPC in addition to component upgrades, a desired goal can be met with significantly less upfront costs.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Measurement of free chlorine levels in water using potentiometric responses of biofilms and applications for monitoring and managing the quality of potable water

We report residual free chlorine is not monitored continuously at scale in drinking water distribution systems because existing real-time sensor technologies require frequent maintenance, cleaning, and calibration, which makes these products too costly to be used throughout a distribution system. As a result, current measurement approaches require manual sampling, which is not feasible for the consistent monitoring of free chlorine because chlorine concentrations vary significantly throughout pipeline distribution and over time and space. This research presents an alternative and cost-effective method of predicting free chlorine levels in drinking water using graphite electrodes coated with naturally grown microbial biofilms. This Microbial Potentiometric Sensor (MPS) array was installed in a Continuously Mixed Batch Reactor (CMBR), and drinking water containing variable free chlorine concentrations. The chlorine concentrations were introduced in a controlled manner, and the MPS signals were monitored over time. MPS signals were measured from the change in Open Circuit Potential (OCP) across the MPS array in real-time. An empirically derived relationship between the normalized change in OCP and free chlorine was established by fitting individual and average MPS data to a decaying exponential growth function in order to predict free chlorine levels. The results show that free chlorine can be predicted with reasonable accuracy, with model validation showing an average absolute error of ±0.09 ppm below 1.1 ppm and ±0.30 ppm between 1.1 and 2.7 ppm. However, the accuracy of predictions was reduced at higher free chlorine levels. The researchers conclude that MPS systems may benefit drinking water distribution systems by measuring free chlorine. These advantages of the MPS are especially pronounced in the developing world because this system is inexpensive and does not require routine maintenance or cleaning. The system relies on a naturally forming and regenerating biofilm and an inexpensive potentiometric meter to produce stable measurements.

54 ENVIRONMENTAL SCIENCES↗

A numerical study of the meso-structure variability in the compaction process of prepreg platelet molded composites

Compaction deformations in a thermoplastic composite system molded from chopped prepreg platelets were analyzed to obtain statistical distributions of local distortions that define the final consolidated meso-structure. This study provides an insight into and allows quantification of the intrinsic geometric variability of the composite meso-structure that controls the variability of the effective composite mechanical properties. A two-phase viscoelastic finite element model was utilized to treat the behavior of molten prepreg platelets, which were assumed to be incompressible and inextensible in the fiber direction during processing. A two-step simulation procedure with discretely modeled platelets and platelet-to-platelet contacts is proposed. The analysis begins with a drop simulation of the uncontrolled material deposition followed by a 3D finite element analysis of the preform compaction and bulk material flow. The proposed model provides a method for generating a composite meso-structure ab initio starting from a statistical sampling of the initial conditions for platelet deposition. The predicted platelet thickness distribution, intra-platelet fiber waviness and out-of-plane platelet undulations were found to be in good agreement with experimental measurements. An estimate of the volume of resin pockets was inferred from the compaction simulation.

36 MATERIALS SCIENCE↗

Visualization and prediction of supercritical CO 2 distribution in sandstones during drainage: An in situ synchrotron X-ray micro-computed tomography study

Pore-scale distribution of supercritical CO 2 (scCO 2 ) exerts significant control on a variety of key hydrologic as well as geochemical processes, including residual trapping and dissolution. Despite such importance, only a small number of experiments have directly characterized the three-dimensional distribution of scCO 2 in geologic materials during the invasion (drainage) process. Here, we present a study which couples dynamic high-resolution synchrotron X-ray micro-computed tomography imaging of a scCO 2 /brine system at in situ pressure/temperature conditions with quantitative pore-scale modeling to allow direct validation of a pore-scale description of scCO2 distribution. The experiment combines high-speed synchrotron radiography with tomography to characterize the brine saturated sample, the scCO 2 breakthrough process, and the partially saturated state of a sandstone sample from the Domengine Formation, a regionally extensive unit within the Sacramento Basin (California, USA). The availability of a 3D dataset allowed us to examine correlations between grains and pores morphometric parameters and the actual distribution of scCO 2 in the sample, including the examination of the role of small-scale sedimentary structure on CO2 distribution. The segmented scCO 2 /brine volume was also used to validate a simple computational model based on the local thickness concept, able to accurately simulate the distribution of scCO 2 after drainage. The same method was also used to simulate Hg capillary pressure curves with satisfactory results when compared to the measured ones. Finally, this predictive approach, requiring only a tomographic scan of the dry sample, proved to be an effective route for studying processes related to CO 2 invasion structure in geological samples at the pore scale.

3D quantitative morphometric analysis↗