Search NASA⌕ Search

SEARCH · Search NASA

Results for “number of injections”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Measurement of accreted helium ion velocity using Doppler spectroscopy in the CTIX compact-toroid plasma accelerator

A compact-toroid (CT) accelerator is a device for the production of plasma bodies which contain a relatively strong internal magnetic field, which permits the plasma body to withstand large forces of acceleration while retaining its magnetic field structure. Typically, the CT occupies an annular region between inner and outer radii, defined by coaxial inner and outer electrodes, across which a large electric potential is applied. The internal magnetic fields of the CT contain both a toroidal and poloidal component, as shown in Fig. 1(a). These fields are typically of comparable magnitude. Because of the effective rigidity provided by internal magnetic fields, the CT can be accelerated to high velocities in short distances. Since internal fields increase in proportion to accelerating magnetic field, there is no strength at which the accelerating magnetic fields overcome the internal field (blow-by). Because of this, the kinetic energy density of CTs can be made to exceed the energy density of magnetic fields typical of plasma containment devices such as tokamaks. CTs have the potential to be launched into the interior of such devices for either delivering fuel (typically deuterium or hydrogen), or to deliver particles of high atomic number (high-Z). An important application of high-Z CT injection into tokamaks is to aid in plasma quenching by radiation, to mitigate the effects of plasma disruptions and prevent damage to the plasma confinement device. In addition, high-Z CTs may be used as sources of kinetic energy density for plasma-compression experiments. In this case, the use of plasma species of high atomic number is also desirable, since kinetic energy density is proportional to atomic number for given CT velocity and electron density gas puffing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Time-Lapse Vertical Seismic Profiling (VSP) for CO 2 Storage in a Depleted Oil Field in Northern Michigan

The Midwest Regional Carbon Sequestration Partnership (MRCSP) was founded in 2003 as part of the U.S. Department of Energy’s (DOE’s) Regional Carbon Sequestration Partnership initiative. Since its founding, MRCSP has made significant strides toward making CCUS a viable option for states in the region. The public/private consortium, funded through the DOE Regional Carbon Sequestration Initiative, brings together nearly 40 industry partners and 10 states. Battelle, as the project lead, oversees research, development and operations and coordinates activities among the partners. The incremental, phased approach has built a valuable knowledge base for the industry and paved the way for commercial-scale adoption of CCUS technologies. From 2008 to 2020, MRCSP Phase III focused on the development of large-scale injection projects. This report is part of a series of reports prepared under the Midwestern Regional Carbon Sequestration Partnership (MRCSP) Phase III (Development Phase). These reports summarize and detail the findings of the work conducted under the Phase III project. This report describes the Dover 33 VSP study to test the effectiveness of time-lapse Vertical Seismic Profile (VSP) for detecting and delineating a plume of more than 271,000 tonnes of CO2 injected into the Brown Niagaran and A-1 Carbonate formations within the Dover 33 reef between March 2013 and September 2016. Five 2D walkaway VSP (WVSP) source lines were acquired by SIGMA3 in September 2016 to investigate the possible time-lapse response in both P-wave and PS-wave seismic data. The data was compared to the same survey geometry acquired in March 2013 by SR2020. The seismic energy was recorded into an 80-level, three-component geophone array deployed into the Dover 1-33 well and placed just above the reef. P-wave and PS-wave reflection images were produced for each of the three source lines from both the 2013 baseline and 2016 monitor surveys. The images were then compared to look for changes in the reflectivity at and around the injection location that might indicate how the CO 2 has moved over this time period. This technique revealed several localized areas with sizable impedance differences inside the reef where CO 2 would be expected; however, a large number of similar impedance “hotspots” were also detected outside the reef in areas where injected CO 2 would not be expected. Therefore, a second analysis that involved calculating P-wave and S-wave travel time differences between the 2013 and 2016 VSPs was conducted to look for a change that could be caused by the CO 2 plume.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2022

Metal fuel relocation experiments with pressure injection

A pool-type sodium-cooled fast reactor (SFR) using metal fuels has a number of inherent safety features that can support benign consequences for design basis accidents (DBAs). Even in postulated severe accident conditions, the core is designed to remain under sub-critical condition in a passive coolable geometry. In case of the postulated severe accidents in SFRs, fuel relocation in the core region along the coolant channel is an important negative reactivity feedback factor that lowers the reactor power level and consequently eliminates the possibility of recriticality. Therefore, understanding the relocation behavior of fuels and the coolability of the relocated fuels in the postulated severe accident is one of the most important factors in the safety assessment of SFRs. In the present study, the relocation behavior of the metal fuel in a pin bundle geometry was investigated by injecting the metal fuels into the coolant channels with pressure. The first metal fuel relocation with pressure injection (RPI-1) experiment showed that many of the metal fuels levitated to upper plenum. In case of RPI-2 experiment, the Real Time X-ray Video System (RTXVS) was constructed and used to obtain a real time X-ray video of the experiment. This real time X-ray video clearly showed the relocation behavior of the pressure injected metallic uranium in the sodium coolant channel. Through this video, it was confirmed that part of the fuel was dispersed downward and part of the fuel was dispersed upward, and the relocation behavior of the injected fuel proceeded simultaneously in the downward and upward directions. So, it can be concluded that in case of pressure injection of the metallic fuel, there is a possibility that the metallic fuel can be quickly removed from the core region, resulting in a negative reactivity feedback effect.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Understanding syringeability and injectability of high molecular weight PEO solution through time-dependent force-distance profiles

Opioid medications play a vital role in treating moderate to severe pain. Unfortunately, many drug misusers and abusers attempt to alter the formulations or properties of these drugs by manipulation, (e.g., crushing, chewing, smoking, snorting, injecting). The intravenous (IV) route is most dangerous to abusers, as the drugs directly enter the circulatory system and produce intense euphoria. To obtain a full understanding of the impact of syringe factors (e.g., needle gauge size, needle length, syringe barrel size), on the ease of injection, we undertook a comprehensive assessment of syringeability and injectability of manipulated abuse-deterrent formulations (ADFs). A texture analyzer-based testing method was developed for the measurement of the resistance force of pulling, holding, and pushing phases of injections. Results showed that the finer needle gauge sizes required higher injection force to withdraw drug solutions. In addition, the syringed liquid volume was highly dependent on needle gauge size, holding time, and sample viscosity. In most cases, a lower needle gauge number and a longer holding time increased the syringed volume. Needle length was highly correlated to injection force (R 2 = 0.99). Using longer needles to inject drug solution requires greater force. Furthermore, large barrel size was correlated to pushing force (R 2 = 0.99); thus, increasing the difficulty of pushing the plunger of a large syringe with one hand. Lastly, relationships between injection force, sample viscosity, and testing conditions were elucidated using a mathematical model, which could be used in the future to assess and predict injection force of solution samples.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Fireball antinucleosynthesis

The tentative identification of approximately ten relativistic antihelium ( He ¯ ) cosmic-ray events at AMS-02 would, if confirmed, challenge our understanding of the astrophysical synthesis of heavy antinuclei. We propose a novel scenario for the enhanced production of such antinuclei that is triggered by isolated, catastrophic injections of large quantities of energetic Standard Model (SM) antiquarks in our galaxy by physics beyond the Standard Model (BSM). We demonstrate that SM antinucleosynthetic processes that occur in the resulting rapidly expanding, thermalized fireballs of SM plasma can, for a reasonable range of parameters, produce the reported tentative ∼ 2 : 1 ratio of He ¯ 3 to He ¯ 4 events at AMS-02, as well as their relativistic boosts. Moreover, we show that this can be achieved without violating antideuterium or antiproton flux constraints for the appropriate antihelium fluxes. A plausible BSM paradigm for the catastrophic injections is the collision of macroscopic composite dark-matter objects carrying large net antibaryon number. Such a scenario would require these objects to be cosmologically stable, but to destabilize upon collision, promptly releasing a fraction of their mass energy into SM antiparticles within a tiny volume. We show that, in principle, the injection rate needed to attain the necessary antihelium fluxes and the energetic conditions required to seed the fireballs appear possible to obtain in such a paradigm. We leave open the question of constructing a BSM particle physics model to realize this, but we suggest two concrete scenarios as promising targets for further investigation. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

NRAP-Open-IAM Analytical Reservoir Model: Development and Testing

Geological carbon sequestration (GCS) is a key technology for reducing global carbon dioxide (CO 2 ) emissions. Over the last decade, the U.S. Department of Energy has invested in understanding the science base, developing practical implementation methods, and demonstrating secure GCS technologies to mitigate the environmental impacts associated with the atmospheric release of CO 2 . As part of the National Risk Assessment Partnership, a systems-level risk assessment tool, called the NRAP-Open-IAM, has been developed to conduct risk assessment and enable safe operations at a GCS site. The current NRAP-Open-IAM contains a simple reservoir model component that calculates the evolution of CO 2 saturation and fluid pressure in a storage reservoir during CO 2 injection operations. This report presents the development and testing of a new analytical reservoir reduced-order model (ROM), which is extended from an existing semi-analytical model for estimation of CO 2 and brine leakage along legacy wells, and enhances the capability of the NRAP-Open-IAM to simulate more types of reservoir conditions. The developed model is validated against three reference studies, and the results indicate that the new ROM predicts the behavior of the two-phase fluids (brine and injected CO 2 ) well and is applicable to different reservoir simulation boundary conditions (i.e., constant pressure boundary and infinite-acting boundary) without a priori user specification of the boundary type. Sensitivity analysis for a set of model parameters is performed using 4,000 synthetic cases prepared via a fully automated process and using machine-learning-based feature selection. The stochastic analysis identifies gravitational number (i.e., ratio of gravitational forces to viscous force) and distance between the injection well and observation location as the most impactful parameters for matching the pressure and CO 2 saturation, respectively, between the numerical simulations and the ROM. This report details the possible ROM uncertainties and serves as a guide for users to understand the use and limitations of this ROM. The code implementation of the model will be released as a module within the NRAP-Open-IAM.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Model-Free Voltage Control of Active Distribution System with PVs Using Surrogate Model-Based Deep Reinforcement Learning

Accurate knowledge of the distribution system topology and parameters is required to achieve good voltage control performance, but this is difficult to obtain in practice. This paper proposes a physical-model-free voltage control method based on a surrogate-model-enabled deep reinforcement learning approach. Specifically, a surrogate model is trained in a supervised manner using the recorded limited number of historical data to learn the relationship between the power injections and voltage fluctuations of each node. Then, the deep reinforcement learning algorithm is applied to learn an optimal control strategy from the experiences obtained by continuous interactions with the surrogate model. The proposed method can achieve physical-model-free control of unbalanced distribution network and inform real-time decisions to deal with fast voltage fluctuations caused by the rapid variation of PV generation. Simulation results on an unbalance IEEE 123-bus system show that the proposed method can achieve similar performance as that of perfect physical-model-based approaches while being advantageous over other traditional methods.

active distribution network↗
DOE OSTI2022

Enhancing $\phi$-sensitivity of ignition delay times through dilution of fuel-air mixture

The high $\phi$-sensitivity (η) of ignition delay time (τ IDT ) is one of the desirable fuel properties for the high-load extension of advanced compression ignition engine. Recent studies revealed the effectiveness of a high dilution rate (x D ) for enhancing η at the engine-relevant conditions. This study aims to quantify the effect of dilution on the η of isooctane using a combined experiment-simulation approach. The ignition delay of the isooctane/air mixture was measured with an Advanced Fuel Ignition Delay Analyzer (AFIDA) over the temperature range of 623 - 923 K at 10 bar pressure with global $\phi$ = 0.3 - 0.6, with and without 28.6% of additional N 2 dilution. For precise evaluation of experimental η, the facility effect of the AFIDA experiment was characterized with three-dimensional computational fluid dynamics (3-D CFD) simulation. The temperature in the combustion chamber from 3-D CFD indicated a substantial temporal dependency, varying up to ~52 K by charge-cooling of fuel injection and heat transfer from the wall. We introduced the dimensionless number θ(t) for characterizing the temporal profile of chamber temperature. Consideration of the facility effect using θ(t) resulted in better agreement between the experimental ..eta.. and zero-dimensional (0-D) kinetics simulation. The refined η were then further utilized to quantify the effectiveness of dilution to η The extent of η enhancement with dilution strategy was maximized at the low-temperature chemistry regime, increasing η by 77% with a 28.6% dilution rate. Further analysis on the dilution effect was carried out using 0-D kinetics simulation, revealing the critical dimensionless numbers relevant to the effectiveness of dilution to η enhancement. Here, this study is the first experiment-simulation combined research to quantify the effect of dilution on η, facilitating the kinetics model refinement for better reproduction of $\phi$-sensitivity.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Simultaneous velocity and density measurements of fully developed Rayleigh-Taylor mixing

Here, the dynamics of molecular mixing and the energy transfer process in the Rayleigh-Taylor instability (RTI) are studied through the collection of simultaneous density-velocity measurements. These experiments provide simultaneous density-velocity field measurements, in contrast to previous point measurements. Statistically stationary experiments are performed in a “convective-type” gas tunnel facility, with density contrast achieved through the injection of helium into the bottom stream. Three experiments at Atwood number ≈ 0.1 are captured at three outer scale Reynolds numbers Re = 520, 2260, and 4050. Particle image velocimetry and laser induced fluorescence are employed simultaneously. Statistics of the density and velocity show self-similar collapse of RTI profiles at large Reynolds number Re >2000. Flat velocity profiles indicate homogeneous turbulence characteristics in the core of the mixing region. Significant anisotropy develops in the flow, with horizontal velocity fluctuations being only 60% of the vertical velocity fluctuations. The turbulent mass flux is found to be asymmetric about the centerline, with increased peak towards the spike. Measurements of the molecular mixing show that mixing is maximized at the core of the flow and increases with increased Reynolds number. The transport equation of density-specific-volume correlation b shows that it is mostly produced in the core of the mixing region, and that the spatial evolution of its profile is the result of transport by bulk motion of the bubble and spike. Energy transfer from gravitational potential energy to turbulent kinetic energy and viscous dissipation is observed to occur in the experiment with a ratio of dissipated energy to potential energy released of 38%. The analysis of the turbulent kinetic energy transport equation budget reveals that production is the primary mechanism towards the growth of turbulent kinetic energy in the core of the flow, and is asymmetrically slightly skewed towards the spike. However, it is through the transport that the strong advection at the edges of the mixing region is maintained.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

The Potential Impact of Offshore Wind Energy on a Future Power System in the U.S. Northeast

This study aims to understand and quantify the potential impact of offshore wind on a future electricity system in the U.S. Northeast. In this analysis, a detailed representation of the Northeast power system is adopted, using a generation portfolio for 2024 paired with offshore wind nameplate capacities of 0 gigawatts (GW), 2 GW, and 7 GW. The analysis identifies points of offshore wind interconnection in the ISO-New England and New York Independent System Operator control areas and uses hourly wind profiles from the Wind Integration National Dataset Toolkit. Hourly simulations of the 2024 power system operations show an ability to accommodate the prescribed offshore wind capacities by adapting the system's generation dispatch. Curtailment levels of offshore wind range between 4% and 5%. Offshore wind generation displaces primarily natural-gas combined-cycle generation; however, requires increased flexibility from combined cycles through more frequent start-ups. The number of hours with transmission congestion increases because of offshore wind injection, with varying impact on a subregional level. Offshore wind's capacity credit was found to be 14.5%-28.3% and is lower than estimated in other large-scale power system studies, in part because this study considered a different weather year and a power system representation with higher shares of solar PV and onshore wind. The reliability contribution of resources is a topic area identified for more research and collaboration. The 7-GW scenario shows a reduction in locational marginal price of 11%, with production cost savings of up to 18% compared to the 0-GW scenario.

17 WIND ENERGY↗
DOE OSTI2024

Surface Brightness Fluctuations in Two SPT Clusters: A Pilot Study

Studies of surface brightness fluctuations in the intracluster medium present an indirect probe of turbulent properties such as the turbulent velocities, injection scales, and the slope of the power spectrum of fluctuations toward smaller scales. With the advancement of Sunyaev–Zel’dovich (SZ) studies and surveys relative to X-ray observations, we seek to investigate surface brightness fluctuations in a sample of South Pole Telescope (SPT)-SZ clusters which also have archival XMM-Newton data. Here we present a pilot study of two typical clusters in that sample: SPT-CLJ0232-4421 and SPT-CLJ0638-5358. We infer injection scales larger than 500 kpc in both clusters and Mach numbers ≈ 0.5 in SPT-CLJ0232-4421 and Mach numbers ≈ 0.6–1.6 in SPT-CLJ0638-5358, which has a known shock. We find hydrostatic bias values for M$_{500}$ less than 0.2 for SPT-CLJ0232-4421 and less than 0.1 for SPT-CLJ0638-5358. These results show the importance to assess quantitative values via a detailed multiwavelength approach and suggest that the drivers of turbulence may occur at quite large scales.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Global Sensitivity Analysis of Large Distribution System with PVs using Deep Gaussian Process

Global sensitivity analysis (GSA) of the voltage to uncertain power injection variations plays an important role for appropriate Volt-VAR optimization. This paper proposes a data-driven GSA method for large-scale distribution systems with a large number of uncertain sources. Specifically, the deep Gaussian process is used to identify the mapping relationship between uncertain power injections and voltages. This allows us to resort to the analysis of variance framework to calculate the Sobol indices for GSA. Unlike the existing polynomial chaos expansion and Gaussian process-based approaches, our proposed method has much better scalability. Test results on the EPRI 1747-node K1 circuit with a different number and different probability distributions of uncertain sources demonstrate that the proposed method can achieve accurate GSA under various conditions.

14 SOLAR ENERGY↗
DOE OSTI2024

Initial design concepts for solid boron injection in ITER

As part of ITER’s consideration to change its first wall material from beryllium to tungsten, the ITER Organization has proposed studying the feasibility of real-time solid boron injection (SBI) into the plasma to coat the walls and divertor to supplement glow discharge boronization (GDB). Boron deposits getter oxygen and reduce sputtering of tungsten from plasma facing components (PFCs). Particularly in areas with significant plasma wall interactions, boron coatings are expected to be short-lived under high performance plasma conditions. The proposed SBI system aims to maintain boron layers in these areas to avoid excessive radiation from tungsten in the plasma as a risk mitigation to ensure ITER will be able to reach and sustain Q = 10 conditions. The system will be used sparingly, as redeposition of boron can lead to significant tritium retention, which must be minimized in ITER to comply with nuclear safety concerns. SBI is proposed to limit and precisely control the amount of boron injected in real-time during plasma operation. Here, some of the design requirements and initial concepts for an SBI system in ITER are presented based on previous results carried out with SBI systems on a number of tokamaks and stellarators around the world. Previous results using SBI systems installed by PPPL have injected boron particles from 5 µm–2 mm diameter at calibrated rates of 2–200 mg/s in real-time during plasma operation on AUG, DIII-D, EAST, KSTAR, LHD, TFTR, WEST, and W7-X, leading to improved wall conditions with reduced plasma impurity concentrations and radiated power and improved plasma performance. The boron is ionized in the plasma edge and then deposited on plasma-wetted surfaces. On AUG, EAST and WEST reduced tungsten sputtering sources were observed following several discharges with SBI. Extrapolation of these SBI results are presented to estimate the amount of boron needed for wall conditioning in ITER. Real-time SBI control requirements and plasma operation scenarios for ITER are also described.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Investigating the potential of a higher reactivity fuel to achieve faster heat-up of aftertreatment systems

Shortening the aftertreatment system heat-up period will be crucial to meet upcoming emissions regulations, which mandate significant reductions in tail-pipe emissions for on-road heavy-duty vehicles. This work studies the impact of increasing fuel reactivity on the combustion of post injections during the expansion stroke. The work demonstrates an operating strategy using a blend of di-butyl ether (a potential low-carbon biofuel component) and #2 diesel that achieves higher exhaust enthalpy, while maintaining similar engine-out emissions (NOx and combustion efficiency) in comparison to the baseline #2 diesel fuel (cetane number 41). On average, a 12% increase in exhaust enthalpy was observed for the blend (cetane number 60) relative to the baseline fuel. The ability to operate the higher reactivity fuel at a more retarded post-injection timing, where combustion with #2 diesel was found to be more unstable, and less complete was key to achieving higher exhaust enthalpy. However, a small fuel penalty (3%–5%) was also associated with the higher reactivity fuel with this operating strategy. The formation of fuel lean regions, due to the long ignition delays at late injection conditions, was predicted to be the leading contributor to the increased emissions of products of partial combustion. Here, chemical kinetic simulations using surrogates for the fuels used in the experiments, demonstrated the ability of a higher cetane fuel to outperform the baseline fuel under fuel lean conditions with equivalence ratios less than 0.7.

09 BIOMASS FUELS↗
DOE OSTI2020

A novel active optimization approach for rapid and efficient design space exploration using ensemble machine learning

In this work, a novel design optimization technique based on active learning, which involves dynamic exploration and exploitation of the design space of interest using an ensemble of machine learning algorithms, is presented. In this approach, a hybrid methodology incorporating an explorative weak learner (regularized basis function model) which fits high-level information about the response surface, and an exploitative strong learner (based on committee machine) that fits finer details around promising regions identified by the weak learner, is employed. For each design iteration, an aristocratic approach is used to select a set of nominees, where points that meet a threshold merit value as predicted by the weak learner are selected to be evaluated using expensive function evaluation. In addition to these points, the global optimum as predicted by the strong learner is also evaluated to enable rapid convergence to the actual global optimum once the most promising region has been identified by the optimizer. This methodology is first tested by applying it to the optimization of a two-dimensional multi-modal surface. The performance of the new active learning approach is compared with traditional global optimization methods, namely micro-genetic algorithm (pGA) and particle swarm optimization (PSO). It is demonstrated that the new optimizer is able to reach the global optimum much faster with a significantly fewer number of function evaluations. Subsequently, the new optimizer is also applied to a complex internal combustion (IC) engine combustion optimization case with nine control parameters related to fuel injection, initial thermodynamic conditions, and in-cylinder flow. It is again found that the new approach significantly lowers the number of function evaluations that are needed to reach the optimum design configuration (by up to 80%) when compared to particle swarm and genetic algorithm-based optimization techniques.

Machine learning, active learning, design optimiza↗
DOE OSTI2026

Haze–cloud correlations mediated by supersaturation fluctuations

Atmospheric aerosol particles that contain water-soluble components can absorb water vapor in humid environments and form either haze particles or cloud droplets, depending on supersaturation conditions. Laboratory and in situ measurements have shown that haze particles and cloud droplets often coexist and compete for available water vapor in shallow clouds and fogs, especially under polluted conditions. It is expected that more aerosol particles can form more haze particles and cloud droplets, so that the haze and cloud number concentrations are positively correlated. However, recent large-eddy simulations show that haze and cloud number concentration can be negatively correlated under extremely polluted conditions. Haze–cloud interactions across different environmental settings remain poorly understood. In this study, experiments in a convection cloud chamber with the same aerosol injection rate show that, as supersaturation forcing increases (i.e., changing from polluted to clean conditions), the covariance between haze and cloud droplet number concentration changes from negative to positive and finally to zero. Large-eddy simulations (LES) of the cloud chamber with a fixed supersaturation forcing but varying aerosol injection rates show a similar result for the haze–cloud correlation: near zero in clean, mean-supersaturation-dominated activation conditions; positive in moderate, supersaturation-fluctuation-influenced activation conditions; and negative in polluted, supersaturation-fluctuation-dominated activation conditions. A theoretical covariance framework was developed to interpret this behavior based on the relative magnitude and signs of the correlations between supersaturation and the populations of haze and cloud droplets. Significantly, experiments, LES, and theory all yield the same three-regime behavior for the sign of the haze–cloud covariance. Furthermore, our results show that the haze–cloud covariance remains robust and easily measurable, thereby providing a useful metric for regime identification in the atmosphere, improving regime-aware parameterizations, and informing aerosol interventions such as fog dispersion, rainfall enhancement, and albedo modification.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Simulation of gas mixture dynamics in a pipeline network using explicit staggered-grid discretization

Here we develop an explicit staggered finite difference discretization scheme for simulating the transport of highly heterogeneous gas mixtures through pipeline networks. This study is motivated by the proposed blending of hydrogen into natural gas pipelines to reduce end use carbon emissions while using existing pipeline systems throughout their planned lifetimes. Our computational method accommodates an arbitrary number of constituent gases with very different physical properties that may be injected into a network with significant spatiotemporal variation. In this setting, the gas flow physics are highly location- and time- dependent, so that local composition and nodal mixing must be accounted for. The resulting conservation laws are formulated in terms of pressure, partial densities and flows, and volumetric and mass fractions of the constituents. We include non-ideal equations of state that employ linear approximations of gas compressibility factors, so that the pressure dynamics propagate locally according to a variable wave speed that depends on mixture composition and density. We derive compatibility relationships for network edge boundary values that are more complex than for a homogeneous gas. The simulation method is evaluated on initial boundary value problems for a single pipe and a small network, is cross-validated with a lumped element simulation, and used to demonstrate a local monitoring and control policy for maintaining allowable concentration levels.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Assessing reservoir performance for geologic carbon sequestration in offshore saline reservoirs

Carbon capture and storage (CCS) is a technological strategy to reduce CO 2 emissions from hydrocarbon-based energy production and non-energy industrial emitters, like fertilizer and cement production. To expand the geographic opportunities for CCS, offshore geologic settings are increasingly considered for saline CO 2 storage, particularly where legacy oil and gas infrastructure may be repurposed for CO 2 transport and injection. In this context, the Gulf of Mexico (GoM), United States (U.S.), may offer tremendous opportunity for CCS because the surrounding states are among the largest CO 2 producers in the U.S. and there has been extensive oil and gas development over the last half century. Nevertheless, there remains significant uncertainty in how GoM reservoir/seal systems will respond to industrial-scale CO 2 injections, particularly in the context of meso-scale (1–10s of meters) permeability variations that are difficult to identify prior to injection. Such permeability variations can have significant impacts on CO 2 migration and fluid pressure propagation, both of which govern the overall effectiveness of the CCS project. This research uses ensemble simulation methods to quantify the influence of spatially variable and a priori unknown permeability fields on CO 2 plume development and fluid pressure propagation when CO 2 is injected at 1 MMt per year for ten years. Results show a number of characteristic patterns that reveal the substantial influence of both near- and far-field permeability on reservoir performance. By focusing on basin characteristics that are common among offshore basins worldwide, results from this study suggest that offshore CCS may be a feasible carbon management strategy in regions where offshore basins are proximal to industrial CO 2 sources.

42 ENGINEERING↗