Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Plumes”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Impact of wildfires on regional ozone and PM 2.5 : Considering the light absorption of Brown carbon

The influence of wildfire brown carbon (BrC) in moderating the formation of O 3 and fine particulate matter (PM 2.5 ) due to its absorption of ultraviolet (UV) radiation is, for the first time, investigated using a modified community multiscale air quality (CMAQ) model. The wavelength-dependent imaginary refractive index of the organic carbon from wildfires, which are needed by the CMAQ model for its inline photolysis rate calculation, are generated experimentally from wood burning aerosols from a combustion chamber. On high emission days of the Bastrop County Complex fire in Texas in early September 2011, BrC UV absorption reduces the daytime average NO 2 photolysis rate by up to 38% in the core region of the wildfire plume compared to the base scenario that does not consider BrC absorption. Consequently, O 3 production in the smoke plume is suppressed compared to the base scenario due to reduced HO x radical concentrations. Further, in the core region, the predicted O 3 increase due to wildfire reaches as high as 47–123 ppb without considering BrC absorption, the predicted increase of O 3 is 5–15% lower when BrC absorption is considered. Similarly, considering the BrC UV absorption leads to approximately 1% (or ~2–3 µgm -3 ) lower estimation of the wildfire emissions’ impact on total PM 2.5 . This change is small because secondary aerosols, which are the components affected by BrC absorption, only account for a small fraction of the total PM 2.5 in wildfire impacted regions in this study. In addition, our study shows that assumptions about aerosol mixing state (core-shell vs. homogeneous) in the inline photolysis rate calculation would not significantly affect out assessment of the impact of wildfire BrC light absorption on O 3 and PM 2.5 .

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Dust-planetary boundary layer interactions amplified by entrainment and advections

Mineral dust contributes to more than half of the global aerosol loading. However, the radiative impacts of dust aerosols on planetary boundary layer (PBL) structure have not been explored sufficiently. During a typical dust storm event over Tarim Basin, dust aerosols exhibit a well-mixed distribution during the daytime in spite of a shallow layer of dust particles accumulated at higher altitudes. By contrast, nocturnal dust plumes are located near the surface due to stable stratification. We demonstrate that these differentiated vertical distributions determine the spatial heterogeneity of dust loading, radiative fluxes and PBL height variations. Dust aerosols cause daytime PBL suppression and nighttime PBL promotion through modulating surface and atmospheric radiative budgets. Specifically, dust-induced cooling effect within PBL directly inhibits the daytime PBL development. PBL suppression effect is then amplified by entrainment processes resulting in excessively low PBL height, especially for dust particles below but near the PBL top. Dust plumes weaken updrafts from PBL and downdrafts of the free atmosphere, which further reduce the entrainment mixing through attenuating horizontal and vertical advection, and eventually amplify PBL suppression. At night, near-surface dust aerosols stimulate a warm and unstable lower atmosphere, generate warm advection heating and promote the PBL development. Our study highlights the importance of specifying entrainment parameters and nighttime advection activities in quantifying the dust-PBL interactions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Long-range transported North American wildfire aerosols observed in marine boundary layer of eastern North Atlantic

Wildfire is a major source of biomass burning aerosols, which greatly impact Earth climate. Tree species in North America (NA) boreal forests can support high-intensity crown fires, resulting in elevated injection height and longer lifetime (on the order of months) of the wildfire aerosols. Given the long lifetime, the properties of aged NA wildfire aerosols are required to understand and quantify their effects on radiation and climate. Here we present comprehensive characterization of climatically relevant properties, including optical properties and cloud condensation nuclei (CCN) activities of aged NA wildfire aerosols, emitted from the record-breaking Canadian wildfires in August 2017. Despite the extreme injection height of ~12 km, some of the wildfire plumes descended into the marine boundary layer in the eastern North Atlantic over a period of ~2 weeks, owing to the dry intrusions behind mid-latitude cyclones. The aged wildfire aerosols have high single scattering albedos at 529 nm (ω 529 ; 0.92–0.95) while low absorption Ångström exponents (Å abs ) at 464 nm/648 nm (0.7–0.9). In comparison, Å abs of fresh/slightly aged ones are typically 1.4–3.5. This low Å abs 529 and low Å abs ~110 nm diameter) with thick non-absorbing coatings. The accelerated descent of the wildfire plume also led to strong increase of CCN concentration at the supersaturation levels representative of marine low clouds. The hygroscopicity parameter, κ CCN , of the aged wildfire aerosols varies from 0.2 to 0.4, substantially lower than that of background marine boundary layer aerosols. However, the high fraction of particles with large diameter (i.e., within accumulation size ranges, ~100–250 nm) compensates for the low values of κ, and as a result, the aged NA wildfire aerosols contribute more efficiently to CCN population. These results provide direct evidence that the long-range transported NA wildfires can strongly influence CCN concentration in remote marine boundary layer, therefore the radiative properties of marine low clouds. Given the expected increases of NA wildfire intensity and frequency and regular occurrence of dry intrusion following mid-latitude cyclones, the influence of NA wildfire aerosols on CCN and clouds in remote marine environment need to be further examined.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Effect of soot concentration on radiative transport in a 2 m JP-8 pool fire

Hubbard et al. (2024) investigated conjugate heat transfer in a 7.9 m JP-8 pool fire, utilizing both heavily sooting and non-sooting fuel models. Their findings indicated that all models produced pool radiative heat fluxes within the measurement uncertainty for large-scale outdoor experiments. This study aims to further investigate the influence of soot concentration on radiative transport. To this end, Large Eddy Simulations (LES) of a 2 m diameter JP-8 pool fire were conducted, employing four turbulent combustion chemistry models: the Eddy Dissipation Concept (EDC) and three variants of flamelet model. Additionally, three variations of the one-equation subgrid scale turbulence kinetic energy model, and multiple mesh resolutions were utilized. The primary validation metrics were pool radiative heat flux and plume radiative fraction, with simulation data compared to experimental results published by Blanchat et al. (2010). The EDC model produced high soot concentrations near the pool surface, resulting in a significant shielding effect. In contrast, the flamelet combustion model has different kinetic limiters for soot production, leading to soot formation higher in the flame, where it can be more readily oxidized. Flamelet models incorporating both acetylene and benzene contributions to soot growth demonstrated better agreement with experimental data. To assess the sensitivity of pool radiative heat flux to soot mass concentration, we used a one-dimensional approximation over a wide range of soot concentrations and compared it to three-dimensional solutions obtained using Discrete Ordinates. Furthermore, this analysis illustrates the impacts of soot shielding and Turbulence Radiation Interactions (TRI), with TRI emission enhancement factors predicted to be in the range of 3–5 near the edges of the plume.

Computational fluid dynamics↗
DOE OSTI2021

Dimensionality effects on multicomponent ionic transport and surface complexation in porous media

Coulombic interactions between charged species in pore water and at surface/solution interfaces are of pivotal importance for multicomponent reactive transport in porous media. In this study, we investigate the impact of domain dimensionality on electrostatically coupled dispersion and surface-solution reactions during transport of acidic plumes and major ions in porous media. Column and quasi two-dimensional flow-through experiments were performed, with identical silica porous media and under the same advection-dominated conditions. Equal mass fluxes of different electrolyte solutions (i.e., HCl - pH ~ 2.8, NaBr - 100 mM, HCl - pH ~ 2.8 plus NaBr - 100 mM) were continuously injected in the 1-D and 2-D experiments and breakthrough curves of pH and major ions were measured at the outlet of the domains. The presence of pronounced ionic strength gradients in the transverse direction in the 2-D setup caused distinct retardation and transport behaviors of protons and major ions which were not observed in the one-dimensional column experiments. Furthermore, in the cases of salt electrolytes injection, considerably enhanced release of H + (>61%) from the quartz surface was observed in the multidimensional system compared to the one-dimensional setup. Reactive transport modeling was performed to reproduce the experimental outcomes and to analyse the coupling between transport processes, based on the Nernst-Planck formulation of diffusive/dispersive fluxes and on surface complexation reactions at the solid-solution interface. Electrostatic interactions between Na + , B – , and H + , and deprotonation of the quartz surface upon the formation of sodium outer-sphere complexes, are the primary controllers of the spatial and temporal features displayed by the pH and major ions measurements. The reactive transport simulations allowed us to interpret the experimental observations, to visualize the distribution and spatio-temporal evolution of dissolved and solid species, to identify a spatially heterogeneous zonation of Coulombic interactions with distinct behavior at the fringe and core of the injected plumes in the multidimensional setup, and to quantify the different components of the Nerst-Planck fluxes of the charged solutes. This study demonstrates that the domain dimensionality directly affects electrostatic interactions between charged aqueous species in the pore water and surface complexation reactions at the solid-solution interface. The non-trivial effects of dimensionality on multicomponent ionic transport result in a significantly different behavior in 1-D and 2-D systems.

58 GEOSCIENCES↗
DOE OSTI2023

Molybdenum in basalt-hosted seafloor hydrothermal systems: Experimental, theoretical, and field sampling approaches

Here, seafloor hydrothermal vents represent potential sources of Mo and other biologically relevant transition metals to the global ocean, complementing continental runoff. Here, we use a combination of experimental, theoretical, and field-sampling approaches to investigate the behavior of Mo in basalt-hosted seafloor hydrothermal systems to provide insight into the processes controlling Mo concentrations in hydrothermal fluids and to derive estimates of vent fluid Mo concentrations and fluxes. Results of this study demonstrate that reaction fluids generated from 350°C, 500 bar hydrothermal basalt alteration experiments contain 775–801 nmol/kg Mo and are thus comparable to a recently collected time series of natural seafloor vent fluids that contained 200–220 nmol/kg Mo at 302°C and 29–30 nmol/kg at 281–282°C (Evans et al., 2023). Synchrotron-based analyses of experimentally altered basalt produced in this study and additional natural samples of altered oceanic crust originally collected from Pito Deep Rift reveal the presence of Mo-rich particles consistent with trace molybdenite. Comparisons of Mo:Cu ratios in natural vent fluids and near-vent sediment trap samples from Main Endeavour Field indicate that vent fluid Mo is readily incorporated into buoyant plume particles and advected out of the near-vent field, analogous to previous mass balance studies of Cu in this region. Thermodynamic calculations of molybdenite solubility in the context of mineral-buffered hydrothermal fluids and comparisons with natural and experimental hydrothermal fluids suggest that high-temperature vent fluids contain 30–1500 nmol/kg Mo. While a minor component of the modern Mo budget, hydrothermal Mo fluxes are estimated to have constituted 0.3–200× the contemporaneous continental weathering fluxes prior to the ~2.4 Ga ago “Great Oxidation Event” and widespread oxidative continental weathering. Overall, identification of hydrothermal vents as a source of Mo-rich plume particles with potential for dispersal into the wider marine environment has significant implications for hypotheses regarding the co-evolution of Life and Earth’s environments, specifically the form and availability of Mo in anoxic Archean-Eon oceans, where Mo-dependent enzymatic pathways are thought to have emerged and subsequently evolved.

58 GEOSCIENCES↗
DOE OSTI2023

Activity model for 36 elements in Fe-Ni-Si-S-C liquids with application to terrestrial planet accretion and mantle geochemistry: New data for Ru, Re, Pt, Os, Ti, Nb, and Ta

Understanding siderophile element partitioning between metal and silicate melts under diverse conditions can be used to place important constraints on the materials and conditions of planetary accretion and core formation, as well as post core formation processes. However, the effects of Si on the partitioning and activity coefficients for these elements are not well known, despite Si likely being one of the dominant light elements in Earth’s core. To address this gap in understanding, we have undertaken a systematic study of the highly siderophile elements Re, Pt, Os, and Ru, and the refractory lithophile elements Nb, Ta and Ti at 1600 °C and 1 GPa, to derive epsilon interaction parameters for these elements in FeSi metallic liquids. Positive epsilon interaction parameters were measured for Nb, Ta, Ti, Ru, Re, Pt, and Os, indicating that dissolved Si in Fe liquids causes a decrease in their metal/silicate partition coefficients (or ‘silicophobic’ behavior). Furthermore, ε$_{Re, Os, or Ru}^{Si}$ > ε$_{Re, Os, or Ru}^{S}$ which means Si causes a larger decrease in D(metal/silicate) than S, and the chalcophile behavior expected from some elements will be completely masked by the presence of Si in a metallic liquid. The new parameters are used to update an activity model that now includes 36 siderophile elements in Fe-Ni-Si-S-C liquids (27 trace elements considered here). Systematic assessment of these 27 elements shows which have the strongest affinity for Si, C, and S, and also how activity coefficients for these elements would vary during accretion and core formation in Earth, Mars, and Mercury of widely differing fO 2 and core compositional conditions. The activity model is combined with new partitioning expressions for Mo, W, Cr, Re Ru, Pt, and Os and applied to aspects of post core formation mantle geochemistry of Earth, Mars, and Mercury. Our updated expressions show that the BSE Mo/W ratio can easily be achieved with metal/silicate partitioning during growth of the Earth, whereas Re, Os and Ru become lower than and highly fractionated compared with BSE values during core formation and accretion, and thus nearly 99% of their BSE abundances are likely contributed by late accretion. Ru isotopes should be a very good indicator of the source material for the late accretion. The high Pt/Os and Re/Os developed in a deepening magma ocean during the growth of the Earth, indicates 186 Os and 187 Os isotopes could be coupled if this ancient material remained isolated and subsequently became entrained in mantle plumes and measured in surficial lavas. The extent to which this occurred will be limited by the low Os content of this ancient material, thus requiring mixing as a major component in plume sources. Martian mantle Hf/W ratio stays low during accretion and core formation modelling, suggesting that W isotope anomalies are more likely due to solid/liquid silicate fractionation than to core formation. Finally, Ti contents measured by MESSENGER at Mercury’s surface can be explained by segregation of either a metallic core (IW-6 to -8) or metallic core + sulfide (IW-4 to -7.5) followed by mantle melting.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2020

Time-Lapse Gravity Monitoring of CO2 Migration Based on Numerical Modeling of a Faulted Storage Complex

In this study, the performance of both surface and borehole time-lapse gravity monitoring to detect CO2 leakage from a carbon storage site is evaluated. Several hypothetical scenarios of CO2 migration in a leaky fault, and thief zones at different depths at the Kimberlina site (California, USA) constitute the basis of the approach. The CO2 displacement is simulated using the TOUGH2 simulator applied to a detailed geological model of the site. The gravity responses to these CO2 plumes are simulated using forward modeling with sensors at ground surface and in vertical boreholes. Results of inversion on one scenario are also presented. The surface-based gravity responses obtained for the different leakage scenarios demonstrate that leakage can be detected at the surface in all the scenarios but the time to detection is highly variable (10 to 40 years) and dependent on the detection threshold considered. Borehole measurements of the vertical component of gravity provide excellent constrains in depth when they are located in proximity of the density anomaly associated with the presence of CO2, thus discriminating multiple leaks in different thief zones. Joint inversion of surface and borehole data can bring valuable information of the occurrence of leakages and their importance by providing a reasonable estimate of mass of displaced fluids. This study demonstrates the importance of combining multiphase flow simulations with gravity modeling in order to define if and when gravity monitoring would be applicable at a given storage site.

Time-lapse gravity monitoring, leak detection, CCS↗
DOE OSTI2022

Development of lean, efficient, and fast physics-framed deep-learning-based proxy models for subsurface carbon storage

In this work, we present deep-learning-based surrogate models for CCUS developed with four different algorithms and a physics-framed two-phase flow problem involving displacement of water by CO 2 . The deep-learning models were trained using 3D datasets describing the pressure plume, CO 2 saturation plume, and water extraction rate generated by numerical simulation. The hyperparameters defining the architecture of the neural networks were optimized to determine the slimmest network size and training parameters that give the most efficient performance at the least training cost. To develop a robust model that closely mimics the governing physical laws, the discretized form of the two-phase fluid transport equation was used to formulate the supervised deep-learning task. The algorithms investigated in this study predicted the data to above 95% accuracy, with the multi-layer perceptron model demonstrating the best performance by balancing training speed, prediction time, and prediction accuracy with lean network capacity. Furthermore, the surrogate models simultaneously predict reservoir pressure and CO 2 saturation in every grid block, including the surface well extraction rate and bottomhole pressure, at all simulation times for a given static model realization in just a few seconds on a standard desktop computer. A key outcome of this study is that limits can be placed on network design parameters to avoid over designing neural networks, with associated efficiencies in training and prediction times. This is very useful because large volumes of data may be generated in CCUS projects and over-design of neural network architectures imposes penalties that are antithetical to the goal of near-real time forecasting.

58 GEOSCIENCES↗
DOE OSTI2022

Simulated CO 2 storage efficiency factors for saline formations of various lithologies and depositional environments using new experimental relative permeability data

Saline formations are attractive geologic reservoirs for permanent carbon dioxide (CO 2 ) storage. Here, the U.S. Department of Energy's National Energy Technology Laboratory (DOE-NETL) has worked to develop and refine methods and tools for the calculation of CO 2 storage potential in subsurface reservoirs. DOE-NETL's CO 2 -SCREEN provides an online tool for executing these storage methods. CO 2 storage efficiency terms are input parameters in DOE-NETL's methods and equations embedded in the CO 2 -SCREEN, which assesses pore space available for CO 2 storage. In this work, a modeling workflow was initiated to refine two CO 2 storage efficiency terms - volumetric displacement (E V ) and microscopic displacement (E d ). The models are based on new experimental relative permeability data that are specific to homogenous lithology and depositional environments of key subsurface saline formations targeted for CO 2 storage. In future work, heterogenous features will be added to this initial modeling effort to update efficiency factors as described in DOE-NETL's methods and CO 2 -SCREEN tool. E V accounts for the volume utilized in the reservoir under the areal plume, while E d accounts for saturation values in the plume to assess efficiency of CO 2 storage at the pore scale. The results of this work are significant in that prior values were based on a limited geologically non-specific relative permeability data set that were collected prior to 2009. Specifically, we applied numerical simulations using TOUGH3 models to update CO 2 storage efficiency values for supercritical CO 2 injection into brine-saturated reservoirs for three lithologies (clastics, limestone, dolomite) and six depositional environments (Marginal Marine, Strand Plain, Deltaic Complex Fluvial, Aeolian, Shallow Marine, and Reef) that have a high potential for geologic CO 2 storage. Experimental relative permeability data in cores from these environments were utilized in the models with corresponding rock type/sedimentary environment. Results of this study showed that dolomite followed by limestone generated higher ranges of storage efficiency compared to clastics. The updated values provided a tighter efficiency range for clastics, lower P 10 but higher P 90 range for limestone, and higher P 10 and P 90 for dolomite. In general, tighter reservoirs with relatively low permeability and porosity were associated with higher E V and E d , showing efficient reservoir and pore utilization in these scenarios. High reservoir pressure and temperature associated with increasing depth increased the E V , and high CO 2 injection rates resulted in increases in E V and E d , while the impact of permeability anisotropy was minimal after the 30-year injection period.

03 NATURAL GAS↗
DOE OSTI2022

Seismic monitoring of CO 2 geosequestration using multi-well 4D DAS VSP: Stage 3 of the CO2CRC Otway project

An important part of any CO 2 geosequestration project is to ensure CO 2 containment and conformance in the subsurface. This is generally done by implementing a comprehensive, risk-based Measurement, Monitoring and Verification plan, a key element of which is active time-lapse seismic monitoring. However, high cost and environmental impact of the standard surface seismic monitoring dictate the need for a cost-effective and environmentally friendly alternative. An opportunity to develop such method emerges with advances in distributed acoustic sensing (DAS) technology, which turns an optical fibre into a seismic sensor with dense spatial sampling. DAS can be permanently deployed in multiple wells across the geosequestration site providing a robust and non-intrusive network of seismic receivers. This approach was developed and tested in the CO2CRC Otway project, where injection of 15 kt of CO 2 at 1.5 km depth was monitored with a 4D vertical seismic profiling (VSP) using five borehole DAS arrays and mobile vibroseis sources. The 4D DAS VSP in each of the five wells provides broadly consistent images of the CO 2 plume with some differences due to different illumination of the target horizon, lateral variation of velocities, and seismic anisotropy. When the newly injected CO 2 reaches a CO 2 plume created as a result of an earlier injection into the same formation ~600 m updip, 4D DAS VSP shows a change in reflectivity in that area and beyond. Furthermore, this shows a potential of 4D DAS VSP for monitoring gas injection into gas-saturated reservoirs.

58 GEOSCIENCES↗
DOE OSTI2022

Fetch-trap Pairs: Exploring definition of carbon storage prospects to increase capacity and flexibility in areas with competing uses

As CCS grows beyond isolated first-mover projects, the question of how to identify potential storage sites (prospects) is key to siting projects amidst the constraints of neighboring storage projects, hydrocarbon fields and other considerations. To date, prospects have been defined either by the extent of a buoyant trap or by the predicted maximum extent of the CO 2 or pressure plume. Both approaches have worked, but they may struggle in places like the Gulf of Mexico where complex geologic structure can invalidate simple assumptions of radial plume spread and buoyant traps are commonly unattractive due to numerous legacy well penetrations. Here, we propose identifying prospective storage sites by fetch area (drainage cell) rather than by buoyant closure alone. Doing so offers 1) greater freedom in siting injectors to avoid surface and subsurface constraints, including legacy wells that tend to cluster on structural high; and 2) a coherent flow regime in which buoyancy drives all injected CO 2 toward a common high. This strategy gives space to dissipate injection pressure and minimize the number of legacy wells needing review. It may also offer 1) better injectivity by bringing synclines and associated channel axes into play; and 2) potentially improved capacity by tapping a larger rock volume and taking advantage of migration losses.

58 GEOSCIENCES↗
DOE OSTI2026

Elastic-wave sensitivity-guided adaptive seismic survey design for cost-effective monitoring of geological carbon storage

Effective seismic monitoring is essential for verifying CO₂ containment, detecting potential leakage, and optimizing operational decisions in geologic carbon storage. Here, this study presents a time-adaptive, elastic-wave sensitivity-guided framework for designing cost-effective seismic monitoring layouts for tracking CO₂ plume migration. The method is based on elastic-wave sensitivity analysis, which quantifies how variations in subsurface properties impact seismic wavefields. Two complementary design strategies are developed: one based on selecting a fixed number of seismic sources (Method A), and the other based on selecting source–receiver pairs contributing to a fixed fraction of cumulative elastic-wave sensitivity energy (Method B). The optimization workflow to identify source–receiver configurations with the highest detection potential is demonstrated using a hypothetical GCS scenario at the Kimberlina site in California using simulations of elastic-wave sensitivity data at multiple post-injection timesteps. Results show that both strategies adapt to evolving plume geometries and wavefield sensitivities, with Method B offering broader spatial coverage and Method A ensuring simpler deployment. This framework enables site-specific, cost-effective, and risk-informed seismic survey designs, enhancing the ability to monitor CO₂ migration over time in evolving geological environments

58 GEOSCIENCES↗
DOE OSTI2025

Modeling supercritical CO2 injection induced rupture of a minor fault embedded in a poroelastic layered reservoir-caprock system

CO2 injection for geologic carbon sequestration involves hydromechanical processes that lead to changes in fluid pressure and stresses that can activate existing faults. This paper presents a new method and workflow of modeling fault activation considering more complex three-dimensional geometry of natural faults using the TOUGH-FLAC multiphase fluid flow and geomechanical simulator. In this method and workflow, FLAC3D mechanical interfaces and TOUGH3 finite volume elements are discretized using computer aided design and gridding software along with a tailored mesh translation routine. The method and workflow are demonstrated with a model of a curved minor fault embedded in a poro-elastic layered reservoir-caprock system. The model is used for a comprehensive sensitivity analysis of fault responses to fault length, injection mass rate, injection schedule, well-fault distance, and well locations versus fault location. Four metrics (CO2 plume, shear state of fault, pressure and stress path at fault monitoring points) are selected to assess CO2 migration, pressure change, and the reactivation of faults. The results reveal that CO2 can bypass around the tip of the minor impermeable fault, building up pressure and poro-elastic stress on both sides that tends to impede fault rupture. Our study shows the benefit of carefully designing the injection to achieve the targeted final storage volume, starting at a relatively low rate for considerable time, and then ramping up the injection rate to the full rate of injection. The initial low injection has two distinct benefits: (1) it allows for the formation of an extensive CO2 plume with a much higher mobility through a low viscosity that will result in a lower pressure for a given injection rate, and (2) it allows for gradual build-up of horizontal poro-elastic stress within the reservoir that will tend to impede activation of steeply dipping faults. The injection scenario starting at a low injection rate, denoted here as conservative injection, can significantly reduce the risk of fault activation as high fluid mobility and reservoir strengthening poro-elastic stress has been established long before reaching the peak injection rates. Moreover, simultaneous injection in two injection wells on both sides of fault can provide further reservoir strengthening through poro-elastic stress buildup acting on a fault under normal faulting stress regime. The findings presented in the paper can provide practical and effective guidance on long-term, safe, and reliable geological CO2 storage.

Cao, Meng↗
DOE OSTI2023

Computational study of the ECN spray C via one-way coupling of internal nozzle flow and ensuing spray

Here, in this study, computational fluid dynamics simulations of the Engine Combustion Network Spray C injector were performed to investigate the effect of the internal nozzle flow on the development of the ensuing spray in automotive fuel injectors. A set of simulation best practices was first developed to achieve an accurate representation of the internal nozzle flow. Detailed, spatiotemporally defined information was extracted at the injector’s orifice exit and stored in maps that were subsequently employed to initialize a Lagrangian spray by means of a static, one-way coupling (OWC) approach. To highlight the impact of the in-nozzle flow on the spray, the OWC approach was compared with the conventional rate-of-injection (ROI) approach, for which the spray parcels are commonly initialized using a blob injection model. The results indicated that, using the ROI approach, considerable evaporation was observed within the plume in the near-nozzle region. On the other hand, the OWC approach showed evidence of evaporation farther downstream and at the plume periphery. The different behavior observed with the two approaches was mostly due to substantial differences in fuel mass distribution. The liquid fuel concentration predicted by the OWC approach matched well with X-ray tomography measurements of projected fuel mass per unit volume, in terms of both magnitude and distribution. In contrast, an unphysical distribution of fuel mass characterized by an extremely low liquid concentration was predicted by the ROI approach. In the far-field, the effect of the in-nozzle flow was found to be weaker as the spray penetrated and interacted with the ambient gas. In this region, the OWC and ROI approaches provided similar predictions of vapor penetration, spray morphology, and mixing-related characteristics.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2026

Finite-element boundary-integral simulation of thin wires and inhomogeneous penetrable bodies in subsurface multilayered anisotropic media

With the prevailing presence of drilling wells near the subsurface in mature oil and gas fields, the application of electromagnetic methods can be particularly challenging where the electromagnetic field is affected by the steel casing. In the past decades, borehole-to-surface and crosswell electromagnetic methods have been utilized for monitoring of reservoir and underground CO 2 storage. This paper presents a unified finite-element boundary-integral (FEBI) method capable of simultaneously modeling the complex electromagnetic interactions between thin metallic wires (representing steel casings) with 3D trajectory and arbitrary 3D inhomogeneous penetrable bodies (such as CO 2 plumes or hydrocarbon reservoirs) within anisotropic multilayered subsurface environments. Unlike existing approaches that treat these components separately or require dense discretization, or are limited to vertical wells, our unified formulation preserves flexible electromagnetic coupling while delivering improved computational efficiency. Assuming the background formation is multilayered anisotropic media, the surface integral equation method is applied to model the thin wires and boundaries of the inhomogeneous bodies. Meanwhile, the finite element method is applied to model the volume of inhomogeneous bodies. Here, the performance of the proposed FEBI method is assessed through comparison with reference numerical results and its practical significance is demonstrated through CO 2 plume monitoring scenarios.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Temporal flow variations interact with spatial physical heterogeneity to impact solute transport in managed river corridors

The interactions between surface water and groundwater in river corridors lead to temporal fluctuations in subsurface water fluxes which have a critical role on solute transport dynamics. Here, we develop a framework to analyze the relative impacts of different temporal frequencies of the flow field in a spatially heterogeneous aquifer on solute transport. Our analysis indicates that the advection-dispersion equation behaves as a low-pass filter by wiping out the effect of high-frequency velocity fluctuations on the first two spatial moments of the solute plume, namely its center of mass and spreading. The concepts discussed in the theoretical analysis are then applied to understand solute transport dynamics at the 300 Area of the Hanford site (USA) adjacent to the Columbia River. We examine the temporal behavior of the solute plume's spatial moments for different temporal frequencies utilizing geostatistical parameters estimated in the 300 Area. Due to the proximity to the Columbia river, groundwater fluxes at the Hanford site are highly dynamic resulting in a large range of characteristic temporal frequencies. Nonetheless, similar to the theoretical analysis, our results show that the effect of high-frequency fluctuations is filtered, with most of the solute transport dynamics being controlled by fluctuations characterized by a large characteristic period.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Conversion of hydrogen gas to water in soil: Implications to fusion contributions

With the continuing development of fusion energy, it is reasonable to be concerned about contamination from the potential release of tritiated water (HTO) plumes, as well as the oxidation of atmospheric tritiated hydrogen gas (HT) releases in nearby soils. We build on past studies by examining oxidation rates of HT to HTO at two sites representing dry and moist soil. We also examine diurnal changes to better constrain rates of microbial HT oxidation and subsequent deposition into soils. We estimate the potential effects of an HTO plume resulting from an oxidation HT. In conclusion, this methodology framework is expected to be useful in the determination of potential health and environmental impacts from HTO as a result of HT releases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗