Search NASA⌕ Search

SEARCH · Search NASA

Results for “Two-phase injection”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Laboratory evaluation of cyclic underground hydrogen storage in the Temblor sandstone of the San Joaquin Basin, California

Underground Hydrogen Storage (UHS) in depleted oil and gas reservoirs could provide a cost-effective solution to balance seasonal fluctuations in renewable energy generation. However, data and knowledge on UHS at subsurface conditions are limited so it is difficult to estimate how effective this type of storage could be. In this study, we perform high pressure experiment to measure the effectiveness of cyclic hydrogen (H 2 ) storage in a specimen of Temblor sandstone retrieved from the San Joaquin Basin of California. Our experiment mimics reservoir pressure conditions to measure H 2 -brine relative permeability and fluid-rock interactions over the course of ten charging and discharging cycles. Initial gas breakthrough occurred at 15 % to 25 % H2 saturation in the specimen with 3 % NaCl brine as the resident fluid. Continuing injecting to 4 pore volumes (PV) of H 2 yielded an asymptotic H 2 saturation of 38 % to 41 %, a level often referred to as the irreducible gas saturation based on two-phase flow. The boundary condition in this study mimics the near wellbore region, which experiences bi-directional H 2 flow. This bi-directional flow led to evaporative drying of the specimen resulting in 94 % H 2 saturation at the end of 10th cycle. This indicates that cyclic flow and evaporative drying can lead to more efficient reservoir storage where a larger fraction of the reservoir porosity is usable to store H 2 . The produced gas stream consisted of H 2 mixed with 8 % to 22 % H 2 O, indicating formation dry-out by evaporation. Meanwhile, produced water chemistry indicated calcite and silicate dissolution, with calcite sourced from fossil fragments. This led to a loss of cementation and weakened the rock sample. Combined, our results indicate dry-out, compaction, increased H 2 saturation, rock weakening, and permeability loss during cyclic UHS. Overall, we anticipate that the combined effects should lead to higher than anticipated UHS storage efficiency per volume of sandstone reservoir rock.

08 HYDROGEN↗
DOE OSTI2023

Motion of Air Bubbles in a Cement Slurry

The dynamics of air (gas) bubbles in a column of cement slurry is examined numerically. The air injected at the bottom of a laboratory-scale column through a porous distributor plate spatially distributes and migrates as a swarm of bubbles throughout the slurry toward the freeboard. The two-phase system of the cement slurry and the air bubbles is modeled using the conservation equations of mass and linear momentum in the framework of the volume-of-fluid (VOF) approach. The cement slurry is modeled using the Herschel–Bulkley and Bingham fluid models. Results show that the mean Sauter diameter and the mean rise velocity of the bubbles decrease with the gas flow rate. Meanwhile, it is found that the rising of the bubbles is controlled by breakup events, along with relatively weak path instabilities of the bubbles resulting in relatively straight trajectories, independent of the gas flow rate. The extent of the yielded region appears larger for the Herschel–Bulkley model compared to the Bingham fluid model (by approximately 10%).

36 MATERIALS SCIENCE↗
DOE OSTI2022

Two-Phase Flow Mechanisms Controlling CO2 Intrusion into Shaly Caprock

Abstract Geologic carbon storage in deep saline aquifers has emerged as a promising technique to mitigate climate change. CO 2 is buoyant at the storage conditions and tends to float over the resident brine jeopardizing long-term containment goals. Therefore, the caprock sealing capacity is of great importance and requires detailed assessment. We perform supercritical CO 2 injection experiments on shaly caprock samples (intact caprock and fault zone) under representative subsurface conditions. We numerically simulate the experiments, satisfactorily reproducing the observed evolution trends. Simulation results highlight the dynamics of CO 2 flow through the specimens with implications to CO 2 leakage risk assessment in field practices. The large injection-induced overpressure drives CO 2 in free phase into the caprock specimens. However, the relative permeability increase following the drainage path is insufficient to provoke an effective advancement of the free-phase CO 2 . As a result, the bulk CO 2 front becomes almost immobile. This implies that the caprock sealing capacity is unlikely to be compromised by a rapid capillary breakthrough and the injected CO 2 does not penetrate deep into the caprock. In the long term, the intrinsically slow molecular diffusion appears to dominate the migration of CO 2 dissolved into brine. Nonetheless, the inherently tortuous nature of shaly caprock further holds back the diffusive flow, favoring safe underground storage of CO 2 over geological time scales.

58 GEOSCIENCES↗
DOE OSTI2024

, Experimental Investigation of Electrically Charged Water Droplets Two-Phase Cross-Flow Interactions with Humid

Electrically charged water droplets can capture water vapor molecules in surrounding moist air and promote vapor condensation. This concept unfolds an alternative to mechanical cooling-based air dehumidification systems that separate sensible and latent cooling loads in HVAC applications. The effect of dielectrophoresis, electro-hydrodynamic, and micro-diffusive flows due to highly electrically charged water droplets traveling through moist air were investigated as new methods to reduce the air’s water vapor content. This paper presents new experimental data on humidification and power consumption for a system consisting of highly energized electrosprays in combination with mist eliminators. Electrosprays produced highly charged water spray cones crossing the air perpendicularly and upwardly. The moist air stream was circulated in a horizontal duct at 22C, 80% relative humidity, and 0.24 ms-1 airspeed. Mist eliminators made of stainless wire clothes of 36 and 25 μm pore diameter were used to screen the droplets growing during their flight trajectory. The mist eliminators blocked some of the largest diameter droplets. Still, they could not prevent re-evaporation of the droplets into the bulk air from the droplets deposited on their wire clothes. The equilibrium between the dielectrophoresis condensation from the electrosprays and the droplets’ re-evaporation phenomena resulted in reduced humidification when compared to conventional-type spray evaporative cooling systems. However, to dehumidify the air, the droplets injected by the electrosprays must be separated entirely and promptly removed from the bulk airflow at the end of their flight trajectory.

Electrostatic droplets, Electrospray, humidificati↗
DOE OSTI2024

Breakup dynamics in a pressure-swirl injector for urea-water solution applications: A computational study

The co-optimization of in-cylinder combustion and after-treatment technology has become a major aspect in engine design and development, with the goal of meeting the increasingly restrictive emission regulations in the transportation industry. Selective Catalytic Reduction is a robust technology to control the emission of NO x , and the injection of urea in water solution is the exhaust tailpipe is a key aspect of its operation. The proposed work uses high-fidelity Computational Fluid Dynamics to characterize the atomization dynamics of the liquid jet in relevant cross-flow conditions. The study focuses on a commercial low-pressure (9 bar) pressure-swirl injector which is characterized in its internal geometry through high-resolution X-ray micro-computational tomography. The internal two-phase flow has been modeled according to the volume-of-fluid approach in a large eddy simulation framework and validated against near-nozzle X-ray radiography measurement. Moreover, characterizing the breakup dynamics for the swirling hollow cone formation, and assessing the influence of the cross-flow in the breakup dynamics was completed. The results have been reported proposing Re-Oh maps and probability density functions of the spray kinematics. Higher cross-flow momentum generates an increase in the jet intact length and a reduction of the liquid droplet diameters. The axial momentum of the jet is affected by the cross-flow already in the near-nozzle region, determining a relevant deviation of the spray velocities. In conclusion, this work aims to inform the initialization of Eulerian-Lagrangian spray models through the assignment of droplet kinematics and static one-way coupling between volume-of-fluid results and Lagrangian spray parcels, to be used for system-size domain simulations.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Probing multiscale dissolution dynamics in natural rocks through microfluidics and compositional analysis

Mineral dissolution significantly impacts many geological systems. Carbon released by diagenesis, carbon sequestration, and acid injection are examples where geochemical reactions, fluid flow, and solute transport are strongly coupled. The complexity in these systems involves interplay between various mechanisms that operate at timescales ranging from microseconds to years. Current experimental techniques characterize dissolution processes using static images that are acquired with long measurement times and/or low spatial resolution. These limitations prevent direct observation of how dissolution reactions progress within an intact rock with spatially heterogeneous mineralogy and morphology. We utilize microfluidic cells embedded with thin rock samples to visualize dissolution with significant temporal resolution (100 ms) in a large observation window (3 × 3 mm). Here we injected acidic fluid into eight shale samples ranging from 8 to 86 wt % carbonate. The pre- and postreaction microstructures are characterized at the scale of pores (0.1 to 1 µm) and fractures (1 to 1,000 µm). We observe that nonreactive particle exposure, fracture morphology, and loss of rock strength are strongly dependent on both the relative volume of reactive grains and their distribution. Time-resolved images of the rock unveil the spatiotemporal dynamics of dissolution, including two-phase flow effects in real time and illustrate the changes in the fracture interface across the range of compositions. Moreover, the dynamical data provide an approach for characterizing reactivity parameters of natural heterogeneous samples when porous media effects are not negligible. The platform and workflow provide real-time characterization of geochemical reactions and inform various subsurface engineering processes.

58 GEOSCIENCES↗
DOE OSTI2021

Evaluating the Feasibility of Measuring and Analyzing Strain Caused by Changes in Barometric Pressure in the Vadose Zone

Changes in barometric pressure propagate into the subsurface where they can affect water level measurements in wells and cause vertical strain. Previous strain sensors were developed to measure vertical strain, but they were limited to measurements at a single depth, making it difficult to evaluate strains from a migrating pressure wave. Technology known as Coherence-length-gated Microwave Photonics Interferometry (CMPI) uses optical fiber sensors to measure strain at multiple locations. It has the potential to detect strain caused by variations in barometric pressure at multiple depths in the subsurface. The goal of this study is to evaluate the feasibility of using this technology to record the changes in vertical strain in the vadose zone from small fluctuations in air pressure. Propagation of these air pressure fluctuations is affected by hydrogeological properties such as water content and permeability. A second goal of this study is to evaluate how these properties affect the air pressure distribution, strain, and gas-phase diffusivity with a long-term goal of using strain to monitor the vadose zone. The focus of this study is on a suite of laboratory experiments that used a sand-filled column with an open head space. The CMPI fiber was packaged and embedded along the axis of the column to measure the strain at multiple locations while air pressure transducers were installed through the wall of the column. Small periodic fluctuations in air pressure, similar to barometric pressure fluctuations, were created in the head space of the column using an audio driven speaker driven at 4 Hz using a sinusoidal signal created by a function generator. Then, the fluctuations propagated along the column where they were measured by the transducers and recorded as functions of time. Initial tests were conducted using dry sand, but then additional tests were conducted after injecting water into the column, which changed the hydrologic and mechanical properties. Another experiment was conducted after creating a thin barrier to air flow, which was designed to simulate a thin layer of saturated ground during rainfall. A final experiment was conducted by changing the grain size of the material to create a heterogeneity in the upper half of the column. The suite of six experiments was conducted to highlight effects of variations in water content and heterogeneities on air pressure and vertical strain. Numerical simulations were also created to evaluate the laboratory data and to provide a baseline analysis for the results. The numerical simulations used the governing equations to linear poroelasticity and two-phase flow with boundary conditions representing the experiments. The conceptual model for strain caused by barometric pressure recognizes that a fluctuating barometric pressure causes fluctuations in the air pore pressure that decrease in amplitude and lag in time with increase depth. The pressure distribution results in two different loads that cause strain: 1.) a mechanical load as the barometric pressure acts on the ground surface; and 2.) pressure loading in the pore space. The results show that the total vertical strain is a contribution of both of these loads, and this is the strain that was observed in the laboratory. The lab experiments show that both air pressure and strain vary as sinusoids that propagate with depth. The amplitude of the air pressure and strain decrease as a function of depth in all the conditions that were evaluate. For example, in the dry sand, the amplitude of the air pressure decreases from 25 Pa at the head space to 5 Pa at a depth of 70 cm, whereas the amplitude of the strain decreases from 0.15 to 0.06 µε over the same depth interval. In many cases the amplitudes decrease is approximately negative exponential functions of depth. When water was injected into the sand and the saturation increased, the amplitude of the air pressure increased in the partially saturated sand but decreased sharply to zero where the soil was saturated. The strain decreased with depth, but the transition from partially to fully saturated conditions had little effect on the strain profile, even though it had a major effect on the air pressure. This apparently is because strain is caused by variations in both water and air pressure. The fluctuating air pressure caused the water pressure to fluctuate, which caused strain throughout the column. Phase delay of the air pressure and the strain both increase as linear to bi-linear functions of depth. The phase delays indicate that pressure propagates at a velocity of 7 to 8 m/s in dry or partially saturated sand, but it drops by more than an order of magnitude in the vicinity of the capillary fringe. The phase delays indicate the strain propagates faster than the pressure. The data indicate that the strain propagates from 1.2 to more than 4x faster than the pressure. Switching from sand to silt in the column caused the velocity of the pressure to drop by roughly half for both pressures, but curiously, it had little effect on the velocity of the strain.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Towards understanding the development and characteristics of under-expanded flash boiling jets

Here in this research, the n-hexane jets discharged from a straight channel were studied numerically to reveal the development and characteristics of under-expanded flashing jets. Firstly, the development of flash boiling jets was investigated and compared to that of two-phase jets without phase transition. It was shown that the two-phase jets without phase transition can be under-expanded only if the flow is choked. In contrast, flash boiling jets can be under-expanded whether or not the choked conditions are verified. Secondly, the factors affecting the flashing jet behaviors were discussed. The relaxation time in phase transition, the expansion ratio, and the injection pressure have been identified as three key factors determining the jet characteristics. With high enough expansion ratio and low enough relaxation time, the expansion was so strong that the flow became transonic and shock structures similar to those of under-expanded gaseous jets appeared near the orifice exit. Finally, two types of choking have been identified for flash boiling jets, i.e., vaporization-induced choking and expansion-induced choking. The former one is mainly ascribed to the vapor formation in the orifice which reduces the local speed of sound, while the latter one is due to the expansion of under-expanded two-phase flow and the resultant flow acceleration.

42 ENGINEERING↗
DOE OSTI2022

Evaluation of subsurface transport processes of delayed gas signatures applicable to underground nuclear explosions

Radioactive gas signatures from underground nuclear explosions (UNEs) result from gas-migration processes occurring in the subsurface. The processes considered in this study either drive or retard upward migration of gases from the detonation cavity. The relative importance of these processes is evaluated by simulating subsurface transport in a dual-permeability medium for the multi-tracer Noble Gas Migration Experiment (NGME) originally intended to study some aspects of transport from a UNE. For this experiment, relevant driving processes include weak two-phase convection driven by the geothermal gradient, over pressuring of the detonation cavity, and barometric pumping while gas sorption, dissolution, radioactive decay, and usually diffusion represent retarding processes. From deterministic simulations we found that over-pressuring of the post-detonation chimney coupled with barometric pumping produced a synergistic effect amplifying the tracer-gas reaching the surface. Bounding simulations indicated that the sorption and dissolution of gases, tending to retard transport, were much smaller than anticipated by earlier laboratory studies. The NGME observations themselves show that differences in gas diffusivity have a larger effect on influencing upward transport than do the combined effects of tracer-gas sorption and dissolution, which is consistent with a Sobol’ sensitivity analysis. Both deterministic simulations and those considering parametric uncertainties of transport-related properties predict that the excess in concentration of SF 6 compared to 127 Xe as might be captured in small volumetric samples should be much smaller than the order-of-magnitude contrast found in the large-volume gas samples taken at the site. While extraction of large-volume subsurface gas samples is shown to be capable of distorting in situ gas compositions, the highly variable injection rate of SF 6 into the detonation cavity relative to that of 127 Xe at the start of the field experiment is the most likely explanation for the large difference in observed concentrations.

58 GEOSCIENCES↗
DOE OSTI2022

Using Direct Numerical Simulation of Pore-Level Events to Improve Pore-Network Models for Prediction of Residual Trapping of CO2

Direct numerical simulation and pore-network modeling are common approaches to study the physics of two-phase flow through natural rocks. For assessment of the long-term performance of geological sequestration of CO 2 , it is important to model the full drainage-imbibition cycle to provide an accurate estimate of the trapped CO 2 . While direct numerical simulation using pore geometry from micro-CT rock images accurately models two-phase flow physics, it is computationally prohibitive for large rock volumes. On the other hand, pore-network modeling on networks extracted from micro-CT rock images is computationally efficient but utilizes simplified physics in idealized geometric pore elements. This study uses the lattice-Boltzmann method for direct numerical simulation of CO 2 -brine flow in idealized pore elements to develop a new set of pore-level flow models for the pore-body filling and snap-off events in pore-network modeling of imbibition. Lattice-Boltzmann simulations are conducted on typical idealized pore-network configurations, and the interface evolution and local capillary pressure are evaluated to develop modified equations of local threshold capillary pressure of pore elements as a function of shape factor and other geometrical parameters. The modified equations are then incorporated into a quasi-static pore-network flow solver. The modified model is applied on extracted pore-network of sandstone samples, and saturation of residual trapped CO 2 is computed for a drainage-imbibition cycle. The modified model yields different statistics of pore-level events compared with the original model; in particular, the occurrence of snap-off in pore-throats is reduced resulting in a more frontal displacement pattern along the main injection direction. Compared to the original model, the modified model is in closer agreement with the residual trapped CO 2 obtained from core flow experiments and direct numerical simulation.

Kohanpur, Amir H.↗
DOE OSTI2020

Combined visualisation of cavitation and vortical structures in a real-size optical diesel injector

A high-speed flow visualisation set-up comprising of combined diffuse backlight illumination (DBI) and schlieren imaging has been developed to illustrate the highly transient, two-phase flow arising in a real-size optical fuel injector. The different illumination nature of the two techniques, diffuse and parallel light respectively, allows for the capturing of refractive-index gradients due to the presence of both interfaces and density gradients within the orifice. Hence, the onset of cavitation and secondary-flow motion within the sac and injector hole can be concurrently visualised. Experiments were conducted utilising a diesel injector fitted with a single-hole transparent tip (ECN spray D) at injection pressures of 700–900 bar and ambient pressures in the range of 1–20 bar. High-speed DBI images obtained at 100,000 fps revealed that the orifice, due to its tapered layout, is mildly cavitating with relatively constant cavity sheets arising mainly in regions of manufacturing imperfections. Nevertheless, schlieren images obtained at the same frame rate demonstrated that a multitude of vortices with short lifetimes arise at different scales in the sac and nozzle regions during the entire duration of the injection cycle but the vortices do not necessarily result in phase change. The magnitude and exact location of coherent vortical structures have a measurable influence on the dynamics of the spray emerging downstream the injector outlet, leading to distinct differences in the variation of its cone angle depending on the injection and ambient pressures examined.

42 ENGINEERING↗
DOE OSTI2021

Implementation of two-phase gas transport into VERA for molten salt reactor analysis

Molten salt reactors (MSRs) are a class of next-generation nuclear reactors that have received recent industrial and research interest. A generalized species transport solver was implemented in the Virtual Environment for Reactor Applications (VERA) computing suite to extend this tool to analyze liquid-fueled MSRs. This core simulator has been extended to model the transport of fission product gases into a collection of circulating gas bubbles with the purpose of removing the gases. This paper presents the governing species transport equation, along with various nuclear source terms. Development of the source term for phase migration is discussed, along with a simplified interfacial area tracking method. Finally, a case study on a simplified MSR loop is presented in which modeling parameters were varied to assess their impact on gas removal. The steady state results show that parameters such as bubble diameter, gas injection rate and mass transfer coefficient have a low to moderate effect on the fraction of xenon in the core region. Removal efficiency has the greatest effect on the fraction in the core region. After the pump bowl, bubble diameter has a minor effect on the fraction of xenon in the gas void. Finally, these results point out that increasing parameters such as mass transfer coefficient, gas injection rate, and removal efficiency drives the xenon into the circulating gas void, while decreasing bubble diameter also drives xenon into the gas void by increasing interfacial area.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Flow visualisation in real-size optical injectors of conventional, additised, and renewable gasoline blends

Research on renewable and alternative fuels is crucial for improving the energy and environmental efficiency of modern gasoline internal combustion engines. To highlight the influence of fuel rheological and thermodynamic properties on phase change and atomisation processes, three types of gasoline blends were tested. More specifically, the campaign comprised a reference gasoline, an ethanol/gasoline blend (10% v/v) representative of renewable fuels, and an additised gasoline sample treated with viscoelasticity-inducing agents. High-speed imaging of the transient two-phase flow field arising in the internal geometry and the near-nozzle spray region of gasoline injectors was performed employing Diffuse Backlight Illumination. The metallic body of a commercial injector was modified to fit transparent tips realising two nozzle layouts, namely a two-hole real size model resembling the Engine Combustion Network spray G injector and an enraged replica with an offset hole. Experiments were conducted at realistic operating conditions comprising an injection pressure of 100 bar and ambient pressures in the range of 0.1–6.0 bar to cover the entire range of chamber pressures prevailing in Gasoline Direct Injection engines. The action of viscoelastic additives was verified to have a suppressive effect on in-nozzle cavitation (6% reduction in cavitation extent) , while also enhancing spray atomisation at flash-boing conditions, in a manner resembling the more volatile gasoline/ethanol blends. Finally, persisting liquid ligaments were found to form after the end of injection for the additised sample, owing to the surfactant nature of the additives.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2026

Determining the extent of potential fugitive fluid migration from geologic carbon storage in hydrocarbon-bearing reservoirs: Insights from one-dimensional numerical modeling

Numerical modeling of Geologic Carbon Sequestration in permeable reservoirs initially containing hydrocarbons is conducted using the multi-phase, multi-component thermohydrologic simulator TOGA (TOUGH Oil, Gas, Aqueous; TOUGH stands for Transport Of Unsaturated Groundwater and Heat), to determine how phase and composition of the original fluids influence the extent of the zone where upward fugitive fluid migration could potentially occur, denoted R f . The area within R f comprises regions of substantially elevated pressure and free-phase CO 2 saturation, where a breach in reservoir sealing capacity would lead to upward fugitive fluid migration. The model examines the conditions within the storage reservoir that could lead to fugitive flow, but does not model the fugitive flow itself. A one-dimensional radial model of the storage reservoir is used, and three initial phase conditions are considered: single-phase aqueous, two-phase gas-aqueous, and three-phase oil-gas-aqueous. Components that may be present are H 2 O, CO 2 , CH 4 , C 4 H 10 , and C 10 H 22 . The most important factors controlling Rf are (1) the initial gas-phase saturation within the reservoir, and (2) the lateral extent of multi-phase initial conditions, particularly CO 2 . The composition of liquid and gas phases has a secondary effect. The impact of reservoir depth, thickness, injection rate, and hydrologic properties are also briefly examined, with thickness (or equivalently injection rate) having the biggest effect. These results can help to understand important trends in potential response of CO 2 -EOR fields being considered for dedicated CO 2 storage.

CO₂ plume migration↗
DOE OSTI2021

Impacts of Pore Network-Scale Wettability Heterogeneity on Immiscible Fluid Displacement: A Micromodel Study

The mixed-wet nature of reservoir formations imposes a wide range of rock wettability from strong resident-fluid wetting to strong invading-fluid wetting. The characteristics of two-phase flow in porous media composed of mixed-wetting surfaces remain poorly understood. In this study, we investigated the displacement of resident ethylene glycol (EG) by hexane in two mixed-wet micromodels of identical 2.5-D geometry heterogeneity, with uniformly or heterogeneously distributed patches strongly wetting to hexane. These patches are mixed among pores with unaltered EG-wetting surfaces. Along with control tests in the originally EG-wet micromodel, here we show the classic fingering and transitions in flow regimes at $logCa$ (capillary number) from -7.2 to -3.9. Moreover, pore-scale distributions of wettability and their spatial correlation influence displacement efficiency. In the two mixed-wet micromodels, we found (a) an increase of steady-state hexane saturation at the end of experiments by up to 0.12 in the capillary fingering regime and a decrease of at most by 0.06 in the viscous fingering regime, compared to the EG-wet micromodel, and (b) dispersed and fragmented hexane distribution after displacement. Brine drainage during supercritical CO 2 (scCO 2 ) injections in these micromodels occurs with lower wettability contrasts, and under similar viscosity ratios and interfacial tensions resulted in higher displacement efficiency relative to displacement of EG by hexane. While mixed-wettability can enhance displacement efficiency compared to uniform wettability, the dynamics of immiscible fluids in strong mixed-wet reservoirs are expected to be less pronounced in contributing to the efficiency of geological CO 2 sequestration, oil recovery, and remediation of hydrocarbon-contaminated aquifers.

2.5-D micromodel↗
DOE OSTI2023

SILCC – VII. Gas kinematics and multiphase outflows of the simulated ISM at high gas surface densities

ABSTRACT We present magnetohydrodynamic (MHD) simulations of the star-forming multiphase interstellar medium (ISM) in stratified galactic patches with gas surface densities Σgas = 10, 30, 50, and 100 $\mathrm{M_\odot \, pc^{-2}}$. The silcc project simulation framework accounts for non-equilibrium thermal and chemical processes in the warm and cold ISM. The sink-based star formation and feedback model includes stellar winds, hydrogen-ionizing UV radiation, core-collapse supernovae, and cosmic ray (CR) injection and diffusion. The simulations follow the observed relation between Σgas and the star formation rate surface density ΣSFR. CRs qualitatively change the outflow phase structure. Without CRs, the outflows transition from a two-phase (warm and hot at 1 kpc) to a single-phase (hot at 2 kpc) structure. With CRs, the outflow always has three phases (cold, warm, and hot), dominated in mass by the warm phase. The impact of CRs on mass loading decreases for higher Σgas and the mass loading factors of the CR-supported outflows are of order unity independent of ΣSFR. Similar to observations, vertical velocity dispersions of the warm ionized medium (WIM) and the cold neutral medium (CNM) correlate with the star formation rate as $\sigma _\mathrm{z} \propto \Sigma _\mathrm{SFR}^a$, with a ∼ 0.20. In the absence of stellar feedback, we find no correlation. The velocity dispersion of the WIM is a factor ∼2.2 higher than that of the CNM, in agreement with local observations. For $\Sigma _\mathrm{SFR} \gtrsim 1.5 \times 10^{-2}\, \mathrm{M}_\odot \, \mathrm{yr}^{-1}\, \mathrm{kpc}^{-2}$ the WIM motions become supersonic.

Astronomy & Astrophysics↗
DOE OSTI2025

Dark Matter Search Results from 4.2 Tonne−Years of Exposure of the LUX-ZEPLIN (LZ) Experiment

We report results of a search for nuclear recoils induced by weakly interacting massive particle (WIMP) dark matter using the LUX-ZEPLIN (LZ) two-phase xenon time projection chamber. This analysis uses a total exposure of 4.2 ±0.1 tonne-years from 280 live days of LZ operation, of which 3.3 ± 0.1 tonne-years and 220 live days are new. A technique to actively tag background electronic recoils from 214 Pb 𝛽 decays is featured for the first time. Enhanced electron-ion recombination is observed in two-neutrino double electron capture decays of 124 Xe, representing a noteworthy new background. After removal of artificial signal-like events injected into the dataset to mitigate analyzer bias, we find no evidence for an excess over expected backgrounds. World-leading constraints are placed on spin-independent (SI) and spin-dependent WIMP-nucleon cross sections for masses ≥9 GeV/𝑐 2 . The strongest SI exclusion set is 2.2×10 −48 cm 2 at the 90% confidence level and the best SI median sensitivity achieved is 5.1 ×10 −48 cm 2 , both for a mass of 40 GeV/𝑐 2 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Simulation of Particulate Transport for Delivery of Solid Amendments into the Subsurface: FY24 Status Report

For particulate-based amendments to be viable for field-scale remediation at the Hanford Site (e.g., 200 DV-1 Operable Unit), particles need to be delivered a sufficient radial distance from an injection well and retained at concentrations high enough for effective treatment. An accurate description of the particle radius of influence (ROI) is critical for developing an overall remediation strategy. However, field-scale particle simulations are currently limited due to insufficient simulation capabilities and a lack of experimental data to validate and parameterize particle transport models. To help build toward field-scale deployment, this fiscal year (FY) we have (1) developed a pre screening tool to estimate particle transport, (2) implemented particle transport models within PFLOTRAN, and (3) conducted preliminary estimations of particle ROI. While field-scale numerical simulations will ultimately be necessary before remedy design and field implementation, we have developed a pre-screening tool that offers valuable estimations of expected particle injectability and ROI in a 1-D system. The advantage of the tool is that it does not require extensive laboratory experiments and instead makes predictions based solely on routine laboratory measurements. This tool can assist in down-selection and decision-making by identifying which particle amendment systems are worth pursuing in future laboratory experiments, such as 1-D column tests and beyond. With any system, scaling up from the lab to the field presents challenges. Currently, there is no field data available for model calibration or validation. However, the theoretical particle models being developed herein are the best tools available to guide progress toward field deployment. To help bridge this gap and verify model predictions, larger-scale lab experiments are being proposed. To advance simulation capabilities, six particle transport models are being integrated into the reactive transport simulator PFLOTRAN. These include colloid filtration theory (CFT) and five additional particle transport models (M1-M5). Each model, from M1 to M5, progressively incorporates additional particle transport and retention processes. Ultimately, the simplest model capable of accurately describing 1-D column data will be selected and parameterized. During FY24, the CFT and M1 model have been fully implemented within PFLOTRAN. Using an existing 1 D column experiment, the two currently implemented particle transport models (CFT and M1), and associated parameters, were fit to this experiment. While simpler model formulations are helpful for estimations, these formulations could not fully describe particle transport and retention behavior in the previous 1-D column experiment. Thus, additional complexities will need to be considered, which will be accounted for in the M2-M5 model formulations. Additionally, because a viscous, shear thinning fluid was required to keep particles in suspension, considerations for flow will also need to also be accounted for. Therefore, a new immiscible two-phase flow mode is currently being implemented in PFLOTRAN. With some modifications, this new flow module could also support simulation of non-Newtonian liquid amendments, foams, and emulsions. We also estimated the expected ROI of solid amendments using 1-D simulations. The average predicted ROI was approximately 15 ft for micron-sized zero valent iron (mZVI) suspended in xanthan gum (XG). Using the pre screening tool and ROI estimates, additional amendment-delivery laboratory characterization and experiments are proposed. The results from additional experiments can be used to validate and parametrize particulate transport model formulations, which will ultimately provide predictive capabilities for field amendment-delivery systems.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗