Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reactive Flow”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Machine learning to discover mineral trapping signatures due to CO 2 injection

Mineral trapping is pursued as a geological CO 2 sequestration (GCS) mechanism because it permanently stores CO 2 in solid phases or minerals. However, CO 2 mineral-trapping mechanisms are poorly understood due to (1) lack of sufficient field and laboratory data characterizing these complex processes, and (2) challenges to develop site-specific reactive-transport models coupling fluid flow and geochemical reactions occurring at various temporal (from milliseconds to years) and spatial (from pore (millimeters) to field (kilometers)) scales. Reactive transport with additional complexities such as heterogeneity can make the simulation outputs even more difficult to interpret because of complex nonlinearity and multi-scale interdependencies. Furthermore, the values of model outputs such as concentrations can vary by several orders of magnitude, making it harder to correlate and characterize the impact of the variables via traditional data interpretation techniques such as exploratory data analyses. Recently, machine learning (ML) has shown promise in feature discovery and in highlighting hidden mechanisms that cannot be obtained by existing data-analytics and statistical methods. In this study, we applied an unsupervised ML approach, non-negative matrix factorization with custom -means clustering (NMF) to the data generated by reactive-transport simulations of GCS. The reactive-transport data consisted of 19 attributes, including four physio-chemical variables (pH, porosity, aqueous CO 2 , and sequestered CO 2 ), six chemical species (K + , Na + , HCO, Ca 2+ , Mg 2+ , Fe 2+ ), and four carbonate minerals (calcite, dolomite, siderite, and ankerite), a feldspar mineral (albite), and four clay minerals (illite, clinochlore, kaolinite, and smectite) over a period of 200 years of simulation time. Furthermore, the simulation data used was for Morrow B sandstone at the Farnsworth hydrocarbon unit in Texas. Data are sampled at two locations within the model domain: (1) at the injection well and (2) 200 m west of the injection well. The injection was performed for a period of 10 years. Using NMF, we estimated the temporal interdependencies among the 19 attributes over a span of 200 years. We found that NMF was able to identify four reaction stages and their dominant attributes; these cannot be directly discerned through traditional visualization (e.g., line plots, Pareto analysis, Glyph-based visualization methods) or exploratory data analysis tools of the simulation data. The four stages were: reactions in the injection phase followed by short-, mid-, and long-term reactions. The NMF analysis also revealed that 10 among the 19 attributes are dominant. These dominant attributes for mineral trapping include calcite, dolomite at injection well, siderite at 200 m away from the injection well, clinochlore, kaolinite, Na + , K + , Ca 2+ , Mg 2+ , pH, and aqeuous CO 2 . Finally, at late times (65–200 years), our results showed that calcite plays a major role in mineral trapping with insignificant contribution from siderite, ankerite, and clay minerals. These findings make the proposed unsupervised ML-model attractive for reactive-transport sensing towards real-time GCS monitoring.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Rapid Mineral Precipitation During Shear Fracturing of Carbonate-Rich Shales

Target subsurface reservoirs for emerging low-carbon energy technologies and geologic carbon sequestration typically have low permeability and thus rely heavily on fluid transport through natural and induced fracture networks. Sustainable development of these systems requires deeper understanding of how geochemically mediated deformation impacts fracture microstructure and permeability evolution, particularly with respect to geochemical reactions between pore fluids and the host rock. Here, a series of triaxial direct shear experiments was designed to evaluate how fractures generated at subsurface conditions respond to penetration of reactive fluids with a focus on the role of mineral precipitation. Calcite-rich shale cores were directly sheared under 3.5 MPa confining pressure using BaCl 2 -rich solutions as a working fluid. Experiments were conducted within an X-ray computed tomography (xCT) scanner to capture 4-D evolution of fracture geometry and precipitate growth. Three shear tests evidenced nonuniform precipitation of barium carbonates (BaCO 3 ) along through-going fractures, where the extent of precipitation increased with increasing calcite content. Precipitates were strongly localized within fracture networks due to mineral, geochemical, and structural heterogeneities and generally concentrated in smaller apertures where rock:water ratios were highest. The combination of elevated fluid saturation and reactive surface area created in freshly activated fractures drove near-immediate mineral precipitation that led to an 80% permeability reduction and significant flow obstruction in the most reactive core. While most previous studies have focused on mixing-induced precipitation, this work demonstrates that fluid–rock interactions can trigger precipitation-induced permeability alterations that can initiate or mitigate risks associated with subsurface energy systems.

58 GEOSCIENCES↗
DOE OSTI2016

Effects of Coolant Void Condition on Melt Fragmentation using Simulants in Initial Phase of Metal-fueled SFR Severe Accident

Metal fuel has inherent safety characteristics, especially in severe accident conditions. However, the metal-fueled SFR severe accident scenarios need to be categorized since overall severe accident mitigation strategies have not been prepared. In addition, these mitigation strategies should be classified in specific phases. In an initial phase of metal-fueled SFR severe accident scenarios, fuel melting, and cladding failures would occur. After the occurrence of the cladding breach, the molten metal fuel would be ejected into a coolant channel due to fission gas in fuel cavity or rapid evaporation of sodium bond. This is called Ex-pin phenomena. There are few studies considering the Ex-pin phenomena of metal fuel. Related representative studies are transient reactor test facility (TREAT) M series tests. The TREAT M series tests provided fundamental metal fuel performance in transient condition. They informed partial experimental data such as the amount of melt ejected into coolant channel after the cladding failure. However, these experimental data did not examine the melt behavior during the Ex-pin phenomena. In metal-fueled SFR, it is a known fact that the negative reactivity feedback occurs in the initial phase of severe accident. In an unprotected loss-of-flow (ULOF) accident, the void-reactivity feedback may be positive depending on where in the reactor core the sodium voiding occurs. However, molten metal fuel would be better dispersed upward compared to those of an oxide fuel and significant negative reactivity feedback could be induced through such behaviors. Thus, a key safety strategy to mitigate the accident progression at the initial phase of severe accident is to discharge melt upwards out of the core. In this paper, the upward discharge of melt is defined as melt levitation. Such melt levitation could ensure early termination of severe accident. There could be different melting behaviors inside the pin depending on reactor types or the type of fuel. In case of pressurized water reactor (PWR) oxide fuel-water system, the rapid power excursion would initiate melting of fuel at the surface of the pellets. In this paper, the ULOF accident for the metal-fueled SFR is considered as an accident scenario. To induce the melt levitation, the driving force should be clarified. Heo et al. observed a built-up vapor pressure was suggested as one of the driving forces for the upward dispersion of the molten materials. The following factors also have the potential to induce the melt levitation; the pressure induced by fission gas in the cavity or rapid boiling of sodium bond, natural circulation flow of the sodium, the sodium vapor flow generated by the difference between inlet and outlet plenum in the reactor core, the built-up vapor pressure due to the fuel-coolant interaction etc. These factors are balanced each other. They should be examined through the Ex-pin experiments. Before conducting a parametric study for each factor, the general behavior of the molten metal fuel should be investigated in the coolant channel. Since rapid sodium coolant boiling could occur in the initial phase of severe accident scenarios, the molten fuel would be ejected and fragmented in voided channel. Unlike previous studies concerning an un-voided channel, differences in melt behavior in the voided channel are expected. In present study, molten material behavior was evaluated to reveal their characteristics in voided channel. The Ex-pin experiments were performed in both un-voided and voided channel, thus such behavior differences were explicitly observed. In the lab-scale experiments, the un-voided and voided channel were reflected by coolant-filled and air-filled channel, respectively. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1992

Flux-split algorithms for hypersonic flows

This study reviews numerical techniques for the simulation of fluid flows spanning the range of reactive regimes, from local chemical equilibrium to full thermochemical non-equilibrium. In particular, characteristic-based algorithms are considered, of the flux-vector and the flux-difference type, which are becoming ever more popular due to their accurate rendition of physically complicated, shock-wave dominated flows. Consideration is given to the problems associated with modeling the thermo-chemical behavior of reactive mixtures of thermally perfect gases, and a few topics of current research are outlined.

Cinnella, P.↗
DOE OSTI2017

Dynamic reduced-order models of integrated physics-specific systems for carbon sequestration

In this study, CO 2 leakage from a deep storage reservoir into a shallow groundwater aquifer through a wellbore is simulated using physics-specific models conceptualized as nonisothermal two-phase flow in a deep storage reservoir and through the wellbore, and reactive transport in the aquifer. Supercritical CO 2 is injected into the reservoir subsystem with uncertain permeabilities of reservoir and injection strategies (rate and duration as decision variables). The simulated pressure and CO 2 saturation at an uncertain location are connected to the wellbore model as a boundary condition to output CO 2 flux to the reactive transport model in the shallow aquifer. Uncertainties are propagated from the deep reservoir model, to the wellbore leakage model, and eventually to the reactive transport model, thus contributing to system performance. Finally, to minimize the risk associated with the CO 2 leakage, we develop dynamic reduced order models for quantifying and managing the uncertainty propagation.

58 GEOSCIENCES↗
DOE OSTI2022

Dataset documenting reaction-induced changes to five fractured foamed wellbore cement cores during CO2 fluid flow

This dataset consists of reconstructed CT scans for four fractured foamed cement cores, generated in the lab using Portland Class H cement, and associated processed data. The cements were exposed to reactive brines, and CT scans were take at different time points through the experiment. The CT imaging allows for the tracking of dissolution-induced changes to fracture morphology, as well as alteration and dissolution in the cement matrix.

Computed Tomography↗
DOE OSTI2024

Limit cycle oscillations in the zonal-flow-catalyzed interactions of ion-temperature-gradient turbulence

Limit-cycle oscillations are studied for ion temperature gradient turbulence, which, in the absence of large diamagnetic (mean) shear flows, saturates through energy transfer from unstable modes to large-scale stable modes via zonal-flow intermediary modes. Oscillations of zonal flow and turbulence levels are strongly constrained by the reactive, largely non-dissipative character of the zonal flows. Since existing predator–prey models for observed oscillations in experiments do not include energy transfer through zonal flows to stable modes, low-order fluid models with this physics are constructed and investigated. A simple three-wave truncation produces low-amplitude zonal flows that slowly oscillate around a zero mean, with turbulence oscillations between coupled wavenumbers that exceed linear frequencies by orders of magnitude. This inconsistency with experimental observations is caused by the weak non-linear drive of zonal flows in three-wave systems and the lack of multiple-wavenumber turbulent interactions. A more comprehensive model that preserves multiple wavenumber interactions within the context of conservative zonal-flow-mediated energy transfer to stable modes accurately reflects observed dynamics when the phase between stable and unstable modes is occasionally randomized.

Li, P. -Y. (ORCID:0000000295254171)↗
DOE OSTI2016

Modeling the influence of preferential flow on the spatial variability and time-dependence of mineral weathering rates

Inferences of weathering rates from laboratory and field observations suggest significant scale and time-dependence. Preferential flow induced by heterogeneity (manifest as permeability variations or discrete fractures) has been suggested as one potential mechanism causing scale/time-dependence. In this paper, we present a quantitative evaluation of the influence of preferential flow on weathering rates using reactive transport modeling. Simulations were performed in discrete fracture networks (DFNs) and correlated random permeability fields (CRPFs), and compared to simulations in homogeneous permeability fields. The simulations reveal spatial variability in the weathering rate, multidimensional distribution of reactions zones, and the formation of rough weathering interfaces and corestones due to preferential flow. In the homogeneous fields and CRPFs, the domain-averaged weathering rate is initially constant as long as the weathering front is contained within the domain, reflecting equilibrium-controlled behavior. The behavior in the CRPFs was influenced by macrodispersion, with more spread-out weathering profiles, an earlier departure from the initial constant rate and longer persistence of weathering. DFN simulations exhibited a sustained time-dependence resulting from the formation of diffusion-controlled weathering fronts in matrix blocks, which is consistent with the shrinking core mechanism. A significant decrease in the domain-averaged weathering rate is evident despite high remaining mineral volume fractions, but the decline does not follow a math formula dependence, characteristic of diffusion, due to network scale effects and advection-controlled behavior near the inflow boundary. Finally, the DFN simulations also reveal relatively constant horizontally averaged weathering rates over a significant depth range, challenging the very notion of a weathering front.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Evolution of planar fractures in limestone: The role of flow rate, mineral heterogeneity and local transport processes

In this study, the processes that affect the evolution of fracture geometry and reactive transport in limestone fractures were evaluated using an integrated experimental and modeling approach. Flow-through experiments were performed to investigate the alteration of three artificially made planar limestone fractures, which were exposed to HCl-acidified water flow (pH 3.8 ± 0.1) at three different rates. Reactive infiltration instability developed in the fractures because of small fracture aperture variations and naturally occurring heterogeneity in the rock, i.e., the presence of clay spots of higher porosity. The dissolution regime changed from conical to dominant and ramified wormhole with increasing flow rate, i.e., Péclet number. Maps of fracture aperture obtained from X-ray micro-tomography imaging, together with measurements of permeability and analyses of fluid chemistry, provided experimental constraints for modeling. Numerical simulations were performed using a 2.5D continuum reactive transport model to examine the development of flow patterns and the reactivity and permeability evolution over time, and to interrogate the impacts of fluid reactivity, mineral heterogeneity and diffusive boundary layer. The modeling results showed that fracture alteration is very sensitive to initial fluid reactivity and mineralogical heterogeneity, in addition to flow rate. Finally, the changes of local transport processes during fracture opening also affected the reactivity of the fractured samples.

58 GEOSCIENCES↗
DOE OSTI2023

Sensitivity Analysis of Transient Critical Heat Flux by RIA Under High-Pressure Flow Boiling Conditions in TRTL

A reactivity-initiated accident (RIA) is a design-basis accident under which critical heat flux (CHF) is likely to be exceeded. The operational margin for RIAs is currently determined using steady-state CHF lookup tables, which provide conservative estimates relative to transient CHF phenomena. The Transient Reactor Test Loop (TRTL) facility at Oregon State University is capable of performing out-of-pile rapid heating experiments representative of a RIA at conditions representative of a pressurized water reactor (PWR). Here, to further our understanding of and ability to predict transient CHF under PWR conditions, we performed a sensitivity analysis on a RELAP5-3D model of the TRTL facility coupled to the RAVEN code framework to define a proposed experimental test matrix to be performed at the TRTL facility. We then implemented a flow boiling CHF correlation into RELAP5-3D and performed a secondary sensitivity analysis inspecting the impact of the built-in RELAP5-3D CHF and heat transfer multipliers on both the prediction of CHF and key safety parameters, such as peak cladding temperature and heat flux. The results show that the multiplier with the highest influence toward the prediction of CHF occurrence and the safety parameters is the transient CHF multiplier. Operational performance envelopes have been developed for each of the test matrix cases and will be used for validation once the experiments are performed. The TRTL facility is currently performing shakedown testing to verify system performance prior to proceeding with the experimental campaign. Restart testing results include pump curve restart testing, pressure tests, and heater rod thermocouple transients.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Impact of Boron-Containing Lubricant Additive on Three-Way Catalyst Reactivity and Physicochemical Properties

Boron-containing compounds are one of the lubricant additive options due to their suitable properties for additives and have been used as commercial lubricant additives. In the present study, the impact of a boron-containing lubricant oil additive, AR9100 (BR), on Pd/Rh-based three-way catalyst (TWC) performance is investigated, and the results are compared with the baseline no-additive (NA) case and the industry standard zinc dialkyl-dithiophosphate (ZDDP) results. Accelerated engine aging is performed using a genset to expose the catalysts to lubricant additives at high temperatures. All aged TWC samples are investigated for reactivity in a bench-flow reactor and characterized using a variety of analytical techniques. Compared with the no-additive case, the temperatures of 90% conversion (T 90 ) of NO, CO, C 3 H 6 , and C 3 H 8 for the ZDDP-aged TWC sample increased by 34, 30, 37, and 48 °C. However, the T 90 of all gas species for the BR-aged TWC sample are similar to the NA-aged TWC sample. Additionally, a significant decrease in water–gas shift reactivity and oxygen storage capacity is observed in the ZDDP-aged sample, but not in the BR-aged sample. Inductively coupled plasma-optical emission spectrometry (ICP-OES) analysis and electron probe microanalysis (EPMA) maps of accelerated engine aging samples show the presence of phosphorus and boron in ZDDP- and BR-aged TWC samples, respectively. However, no boron-related peaks are observed in the X-ray diffraction (XRD) pattern of the BR-aged TWC sample, which may exist in the form of an amorphous phase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Mechanistic understanding of carbon mineralization in fracture systems using microfluidics

Carbon mineralization in mafic and ultramafic rocks presents an opportunity for permanent carbon storage in the Earth's subsurface. However, due to their lower permeability, pre-existing fracture networks are key for mineralization to occur. Therefore, to fully develop this technology, a mechanistic understanding of the mineralization behavior in fractures with the consideration of hydrodynamic components is required. We use high-pressure microfluidics to investigate key mechanisms influencing dissolution–precipitation in a fracture network. The experiments were conducted in micromodels made of natural rocks with a comb-shaped flow channel to mimic a fracture network. This enabled studying the effect of injection rate on coupled dissolution–precipitation in advection and diffusion-dominated flow paths. We used gypsum carbonation as an analog reaction to allow for realistic experimental time frames due to its rapid reaction kinetics. The experimental work is coupled with high-fidelity numerical simulations to enhance our understanding of the parameters affecting the mineralization reaction. Our results demonstrate the importance of flow rate on the rate and nature of the gypsum carbonation reaction revealing that higher flow rates enable deeper penetration of the mineral precipitation front into the dead-end channels. This is an important finding since for sustained mineralization in a fracture network, precipitation in dead-ends while still allowing for flowing fractures is critical. Detailed characterization of the precipitates showed that lower flow rates led to porous and loose precipitates in the form of aragonite while higher flow rates mimicked supersaturation behavior leading to the formation of calcite. The reactive transport simulations further demonstrated the significance of flow velocity in advection-dominated channels to influence the efficiency of carbon mineralization in diffusion-dominated channels, potentially clogging of dead-end channels. These findings highlight the need for coupling chemical, mechanical, and hydrodynamic processes to evaluate the nature and extent of carbon mineralization in fractured media critical for permanent storage in mafic and ultramafic formations. This research further highlights the need for more investigation in potential subsurface fracture generation techniques to aid carbon mineralization.

25 ENERGY STORAGE↗
DOE OSTI2021

Visualization and analysis of coupling between plasmas self-organization and plasma-induced fluid circulation in 1 atm DC glows with liquid anode

The physical processes prevailing when plasma contacts liquid water are poorly understood; however, it is this very interaction that is the basis of new, exciting technologies for water purification, the treatment of cancer and disease and the production of new, high value chemical products. A chief impediment to widespread application of plasmas for the aforementioned applications is scale up and cost. Here we investigate plasma induced flows which facilitate transport of reactivity from the plasma to liquid water at rates much faster than diffusion. We find that this plasma induced motion varies depending on whether or not the plasma is self organized. The insight from this work paves the way to a better understanding of transport and ultimately informing plasma treatment technologies on how to best optimize the rate at which processes such as plasma based water purification occurs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

GPU-based simulation of the two-dimensional unstable structure of gaseous oblique detonations

In this paper, the two-dimensional structure of unstable oblique detonations induced by the wedge from a supersonic combustible gas flow is simulated using the reactive Euler equations with a one-step Arrhenius chemistry model. A wide range of activation energy of the combustible mixture is considered. Computations are performed on the Graphical Processing Unit (GPU) to reduce the simulation runtimes. A large computational domain covered by a uniform mesh with high grid resolution is used to properly capture the development of instabilities and the formation of different transverse wave structures. After the initiation point, where the oblique shock transits into a detonation, an instability begins to manifest and in all cases, the left-running transverse waves first appear, followed by the subsequent emergence of right-running transverse waves forming the dual-head triple point structure. This study shows that for low activation energies, a long computational length must be carefully considered to reveal the unstable surface due to the slow growth rate of the instability. For high activation energies, the flow behind the unstable oblique detonation features the formation of unburnt gas pockets and strong vortex-pressure wave interaction resulting in a chaotic-like vortical structure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

ERSAP: Toward Better NP Data-Stream Analytics With Flow-Based Programming

This paper presents an reactive, actor-model and FBP paradigm based framework that we develop to design data-stream processing applications for HEP and NP. This framework encourages a functional decomposition of the overall data processing application into small mono-functional artifacts. Artifacts that are easy to understand, develop, deploy and debug. The fact that these artifacts (actors) are programmatically independent they can be scaled and optimized independently, which is impossible to do for components of the monolithic application. One of the important advantages of this approach is fault tolerance where independent actors can come and go on the data-stream without forcing the entire application to crash. Furthermore, it also makes it is easy to locate the faulty actor in the data pipeline. Due the fact that the actors are loosely coupled, and that the data (inevitably) carries the context, they can run on heterogeneous environments, utilizing different accelerators. This paper describes the main design concepts of the framework and presents a ?proof of concept? application design and deployment results obtain processing on-beam calorimeter streaming data.

Gyurjyan, V.↗
DOE OSTI2024

Impacts of irregularly-distributed acidified brine flow on geo-chemo-mechanical alteration in an artificial shale fracture under differential stress

The efficacy of geological carbon sequestration is reliant on the integrity of the caprock and its resistance to physical and chemical alteration. Caprocks with high abundance of reactive carbonates like calcite are susceptible to acid-promoted dissolution and can result in structural weakening. This work investigates the effect of acidified brine flow through an artificially fractured, high-carbonate (30 % by XRD) shale under differential compressive stress. Cylindrical samples were cut in half vertically and milled to create an artificial fracture with interlocking asperities and open channels. Samples were sheared with a single applied stress in a custom flow cell housed within an industrial CT scanner. Further, either acidic (pH 4) or reservoir-simulated (pH 9.5) brine was flowed through the artificial fracture for 7–8 days under reservoir pressure and room temperature. Model simulations indicate flow mainly occurred in open channels, with limited flow between overlapping asperities. Analysis of fracture surfaces by optical and scanning electron microscopy show increased surface alteration and roughness after exposure to pH 4 versus pH 9.5 brine indicating mineral dissolution/loss, and this effect is greater in areas that receive the highest brine flows. Similarly, surface analysis by scratch testing shows fracture toughness decreases more after exposure to acidic versus reservoir-simulated brine, with the greatest alteration in areas of highest acidic brine flows. Despite weakening, no shear slip occurred. Overall, the results indicate that acidified brine can result in significant physical and geomechanical alteration of irregular fracture surfaces in shale caprock, with greatest effects in preferential flow regions.

58 GEOSCIENCES↗
DOE OSTI2013

Upscaling of Solute Transport in Heterogeneous Media with Non-uniform Flow and Dispersion Fields

An analytical and computational model for non-reactive solute transport in periodic heterogeneous media with arbitrary non-uniform flow and dispersion fields within the unit cell of length ε is described. The model lumps the effect of non-uniform flow and dispersion into an effective advection velocity Ve and an effective dispersion coefficient De. It is shown that both Ve and De are scale-dependent (dependent on the length scale of the microscopic heterogeneity, ε), dependent on the Péclet number Pe, and on a dimensionless parameter α that represents the effects of microscopic heterogeneity. The parameter α, confined to the range of [-0.5, 0.5] for the numerical example presented, depends on the flow direction and non-uniform flow and dispersion fields. Effective advection velocity Ve and dispersion coefficient De can be derived for any given flow and dispersion fields, and . Homogenized solutions describing the macroscopic variations can be obtained from the effective model. Solutions with sub-unit-cell accuracy can be constructed by homogenized solutions and its spatial derivatives. A numerical implementation of the model compared with direct numerical solutions using a fine grid, demonstrated that the new method was in good agreement with direct solutions, but with significant computational savings.

advection↗