Search NASA⌕ Search

SEARCH · Search NASA

Results for “reactive flow simulations”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Plasma-assisted chemical-looping combustion: Mechanistic insights into low temperature methane oxidation with CuO

The low-temperature oxidation of CH 4 by CuO in a coaxial, fixed bed, double dielectric barrier discharge (DBD) reactor was investigated with time-dependent species measurements by an electron-ionization molecular beam mass spectrometer (EI-MBMS). In the experiment, 10% methane carried by noble gasses was flown at 50 sccm through 1 g CuO dispersed in quartz wool both under plasma and non-plasma conditions, while time-dependent gas-phase species profiles were collected. Plasma conditions were explored from 300 to 600 °C while the non-plasma conditions were set from 300 to 900 °C. Mechanistic insights into the oxidation of CH 4 by CuO with plasma discharge at lower temperatures (≤ 600 °C) were obtained by quantifying the fuel oxidation, intermediate species, and CO 2 production in comparison to the non-plasma conditions. Here, we observed significant enhancement of fuel oxidation from the plasma discharge between 400 and 500 °C. The CO 2 production at 500 °C with plasma was greater than that at 700 °C without plasma, reducing fuel oxidation temperature by 200+ °C. During tests, three distinct reaction stages were observed: a gas-phase transport limited stage, a surface reaction limited kinetic stage, and an oxygen ion diffusion limited stage. It was observed that plasma greatly improved the reactivity of the second stage at low temperature. In addition, no carbon deposits were observed on the resultant particles, even under the presence of plasma. M. species such as C 4 H 2 and C 6 H 6 not previously observed or predicted in CuO/CH 4 chemical looping were observed, with some species such as CH 3 OH only becoming detectable as total flowrate was increased from 50 to 1500sccm. A non-plasma reaction pathway for CH 4 based the observed species from the MBMS spectrum and previous predictions from reactive molecular dynamics simulations was created, providing a framework from which more complex plasma CuO mechanisms can be crafted from.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1991

Surface reaction and pore diffusion in flow-tube reactors

The interaction of gas diffusion with surface reaction in porous solids is discussed and applied specifically to heterogeneous rate measurements in flow-tube reactors. External diffusion to the outer surface of a reactive solid, internal diffusion within the pores, surface reaction, and laminar flow are considered. A procedure is developed to correct observed surface rate constants for the interaction of these processes. Measured surface areas and bulk densities are used to construct a semiempirical model for porous diffusion in vapor-formed HNO3-H2O ices which are used to simulate polar stratospheric cloud surfaces. The model is tested experimentally by varying the thickness of these ices from about 15 to 120 microns. The results are consistent with the model predictions and show that the HNO3-H2O ices used are highly porous, and the internal surface must be considered in calculating kinetic parameters from observed loss rates. The best fit of the data yields a tortuosity factor of 3.3 +/-1.1 for the ice substrates.

Keyser, Leon F.↗
DOE OSTI2024

Approaches for the Simulation of Coupled Processes in Evolving Fractured Porous Media Enabled by Exascale Computing

Models have historically represented fractured porous media with continuum descriptions that characterize the media using bulk parameters. The impact of small-scale features is not captured in these models, although they may be controlling the performance of subsurface applications. Pore-scale models can simulate processes in small-scale features by representing the pore space geometry explicitly but are computationally expensive for large domains. The alternative multiscale approach entails the combination of pore-scale and continuum-scale descriptions in a single framework. We use Chombo-Crunch, a computational capability that discretizes complex geometries with an adaptive, embedded boundary method to contrast these two approaches. Chombo-Crunch takes advantage of recent computational performance and memory bandwidth improvements resulting from the emergence of exascale computing resources. These combined improvements enable the efficient simulation of reactive transport in fractured media with a high degree of fidelity and the ability to capture the control small-scale processes exert on the overall medium evolution.

42 ENGINEERING↗
NASA NTRS1995

OMV/VTE variable thrust engine analysis

The objective of the present work is to develop a predictive CFD based analytical tool for the Variable Thrust Engine (VTE) in the Orbital Maneuvering Vehicle (OMV). This objective is being accomplished within the framework of the Los Alamos KIVA computer code for chemically reactive flows with sprays. For the OMV application, the main structure of KIVA is to be retained while reformulating many of the phenomenological submodels, enhancing some of the numerics, and adding more features. The analytical model consists of the general conservation equations for two-phase reactive flows and of submodels for turbulence, chemical reactions, and bipropellant sprays. Tailoring this model to the OMV engine brings about the added complexities of combustion and flow processes that occur in a liquid hypergolic propellant rocket chamber. This report exposes the foundation upon which the analytical tool is being constructed and developed. Results from a cursory computational exercise involving the simulation of the flow and combustion processes in a hypothetical N2H4/N204 rocket engine thrust chamber is presented and discussed.

Larosillere, L.↗
NASA NTRS1989

OMV/VTE variable thrust engine analysis

The objective of the present work is to develop a predictive CFD based analytical tool for the variable thrust engine (VTE) in the orbital maneuvering vehicle (OMV). This objective is being accomplished within the framework of the Los Alamos KIVA computer code for chemically reactive flows with sprays. For the OMV application, the main structure of KIVA is to be retained while reformulating many of the phenomenological submodels, enhancing some of the numerics, and adding more features. The analytical model consists of the general conservation equations for two phase reactive flows and of submodels for turbulence, chemical reactions, and bipropellant sprays. Tailoring this model to the OMV engine brings about the added complexities of combustion and flow processes that occur in a liquid hypergolic propellant rocket chamber. This report exposes the foundation upon which the analytical tool is being constructed and developed. Results from a cursory computational exercise involving the simulation of the flow and combustion processes in a hypothetical N2H4/N2O4 rocket engine thrust chamber is presented and discussed.

Larosillere, Louis↗
NASA NTRS2003

Alabama NASA EPSCoR Preparation Grant Program: Grant No. NCC5-391

The funded research projects under the Experimental Program to Stimulate Cooperative Research (EPSCoR) grant program and the student fellowship awards are summarized in this report. The projects include: 1) Crystallization of Dehydratase/DcoH: A Target in Lung Disease; 2) Measuring Velocity Profiles in Liquid Metals using an Ultrasonic Doppler Velocimeter; 3) Synthesis, Structure, and Properties of New Thermoelectric Materials; 4) Computational Determination of Structures and Reactivity of Phenol-Formaldehyde Resins; 5) Synthesis of Microbial Polyesters in the NASA Bioreactor; 6) Visualization of Flow-Fields in Magnetocombustion; 7) Synthesis of Fluorescent Saccharide Derivatives. The student fellowship awards include: 1) Distributed Fusion of Satellite Images; 2) Study of the Relationship between Urban Development, Local Climate, and Water Quality for the Atlanta, Georgia Metrop; 3) Computer Simulation of the Effectiveness of a Spring-Loaded Exercise Device.

Gregory, John C.↗
DOE OSTI2025

Enabling BWR fuel rod analysis in the BISON fuel performance code

Nuclear fuel vendors around the world are pursuing approaches to sustain the existing nuclear reactor fleet consisting primarily of pressurized-water reactors (PWRs) and boiling-water reactors (BWRs). To support the industry's efforts, advanced modeling and simulation tools need to be capable of analyzing both legacy reactor concepts. BWR fuel rods are significantly different than those used in PWRs, which can affect fuel performance analysis. BWR fuel rods include: (1) an extensive use of gadolinia dopant as a burnable absorber, (2) an axial variation in fuel enrichment and gadolinia content, (3) the inclusion of a liner on the inner cladding surface to mitigate the impact of pellet-clad mechanical interaction (which impacts hydrogen and hydride distribution), (4) a lower initial fill gas pressure, (5) bottom-entry control rods, and (6) a lower coolant pressure that results in the two-phase flow boiling phenomenon. Although the primary focus of BISON has been in the area of PWR and advanced reactor fuel analyses, this paper presents the developments in BISON to support BWR fuel performance analysis. An overview of the models that account for the effects of gadolinia is highlighted. Internal mesh generation capabilities to include a liner is presented. Normal operating and transient (reactivity insertion accident) demonstration problems are presented to illustrate the impact of gadolinia, the hydrogen and hydride evolution due to the presence of the liner, and BISON's ability to simulate axially varying enrichments and dopant concentration. Bottom-entry control effects are captured by the axial power peaking factors present in the demonstration cases. Comparisons to integral experiments from the Halden IFA-681 experiments are discussed as initial validation. Reasonable comparisons are obtained for fuel centerline temperature and rod internal pressure as a function of time. In conclusion, simulations of additional experiments containing Gd 2 O 3 -bearing fuel are necessary to completely validate the code for BWR applications.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

CompLaB v1.0: a scalable pore-scale model for flow, biogeochemistry, microbial metabolism, and biofilm dynamics

Abstract. Microbial activity and chemical reactions in porous media depend on the local conditions at the pore scale and can involve complex feedback with fluid flow and mass transport. We present a modeling framework that quantitatively accounts for the interactions between the bio(geo)chemical and physical processes and that can integrate genome-scale microbial metabolic information into a dynamically changing, spatially explicit representation of environmental conditions. The model couples a lattice Boltzmann implementation of Navier–Stokes (flow) and advection–diffusion-reaction (mass conservation) equations. Reaction formulations can include both kinetic rate expressions and flux balance analysis, thereby integrating reactive transport modeling and systems biology. We also show that the use of surrogate models such as neural network representations of in silico cell models can speed up computations significantly, facilitating applications to complex environmental systems. Parallelization enables simulations that resolve heterogeneity at multiple scales, and a cellular automaton module provides additional capabilities to simulate biofilm dynamics. The code thus constitutes a platform suitable for a range of environmental, engineering and – potentially – medical applications, in particular ones that involve the simulation of microbial dynamics.

58 GEOSCIENCES↗
DOE OSTI2023

A New Proposal Generalized Predictive Control Algorithm With Polynomial Reference Tracking Applied for Sodium Fast Reactors

This paper proposes a generalized predictive control (GPC) with constraints and orthonormal Laguerre functions using the simplified model of the primary system (reactor core and intermediate heat exchanger (IHX)) of a prototypical sodium fast reactor (SFR). This paper develops a multiple-input multiple-output (MIMO) GPC with input constraints able to track polynomial references of any degree applied in coolant temperature difference across the core and fractional power. The manipulated variables of the GPC-SFR are the reactivity and the sodium flow rate of the primary and secondary pipes. Moreover, orthonormal Laguerre functions and step down condition number techniques were also applied to avoid the numerical ill-conditioning issue in quadratic programming of large systems. Thus, a GPC type-2 was designed to control fractional power, coolant temperature difference across the core and sodium tank temperature of the SFR primary system when temperature references change according to a linear ramp after reaching their steady-state operation, sustaining 100% power operation on the reactor. In order to analyze the load tracking capability of the GPC-SFR type-2, the load following from 100% fractional power (FP) to 60% FP at 0.8% FP/min rate is simulated. Constraints on the rate of coolant temperature difference across the core and reactivity were applied for the design safety. For comparison criteria, this paper compares the GPC-SFR type-2 with the GPC-SFR type-1, i.e., standard model predictive control (MPC), to verify the viability and superior performance of the proposal regarding: (a) ramp-tracking capability of temperature and load; (b) the rejections of a reactivity disturbance of -1 cent and a secondary sodium inlet temperature disturbance of +10°F; and (c) a simulation with uncertainty in reactor design. The simulations show that the GPC-SFR type-2 overcome the GPC-SFR type-1 robustness and performance.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Molecular fluctuations inhibit intermittency in compressible turbulence

In the standard picture of fully developed turbulence, highly intermittent hydrodynamic fields are nonlinearly coupled across scales, where local energy cascades from large scales into dissipative vortices and large density gradients. Microscopically, however, constituent fluid molecules are in constant thermal (Brownian) motion, but the role of molecular fluctuations in large-scale turbulence is largely unknown, and with rare exceptions, it has historically been considered irrelevant at scales larger than the molecular mean free path. Recent theoretical and computational investigations have shown that molecular fluctuations can impact energy cascade at Kolmogorov length scales. Here, we show that molecular fluctuations not only modify energy spectrum at wavelengths larger than the Kolmogorov length in compressible turbulence, but also significantly inhibit spatio-temporal intermittency across the entire dissipation range. Using large-scale direct numerical simulations of computational fluctuating hydrodynamics, we demonstrate that the extreme intermittency characteristic of turbulence models is replaced by nearly Gaussian statistics in the dissipation range. These results demonstrate that the compressible Navier–Stokes equations should be augmented with molecular fluctuations to accurately predict turbulence statistics across the dissipation range. Our findings have significant consequences for turbulence modelling in applications such as astrophysics, reactive flows and hypersonic aerodynamics, where dissipation-range turbulence is approximated by closure models.

compressible turbulence↗
DOE OSTI2025

Direct Comparison of the Activity and Selectivity of Rh 1 Cu and Ni 1 Cu Single-Atom Alloy Sites for Ethanol Decomposition

Ethanol is an important source of clean hydrogen, acetaldehyde, acetic acid, acetate esters, and light hydrocarbons. Controlling the divergent reaction pathways to these products requires understanding how different active sites influence the elementary steps involved. Herein, we present a combined surface science, theory, and nanoparticle catalysis study demonstrating how two single-atom dopants (Rh and Ni) in a Cu host can distinctively alter the selectivity of alcohol conversion. Specifically, our model studies reveal that ethanol reacts on Ni 1 Cu single-atom alloys to selectively produce acetaldehyde, whereas methane and CO are also formed on Rh 1 Cu single-atom alloys. Interestingly, these different reactivities are in contrast to the behavior of the pure metals as Ni(111) and Rh(111) surfaces favor methane/CO and surface carbon/CO, respectively. DFT calculations of reaction pathways and simulated product desorption based on microkinetic analyses explain these reactivity differences, demonstrating that C–C cleavage leading to methane formation has a lower barrier on Rh single-atom sites. To test the catalytic relevance of these fundamental results we synthesized and characterized supported Ni 1 Cu and Rh 1 Cu single-atom alloy nanoparticles with dopant:Cu ratios of 1:200. Flow reactor results revealed that both Ni and Rh increased ethanol conversion over Cu and that Ni 1 Cu catalysts were >99.9% selective to acetaldehyde, while Rh 1 Cu also produced 0.6%–2.6% of equimolar methane and CO between 433 and 493 K, demonstrating that C–C bond cleavage is enabled by isolated Rh sites. Furthermore, these catalytic results bridge the pressure and materials gaps, and together, this study provides insights into how different isolated dopant sites promote different catalytic pathways.

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

Learning the factors controlling mineral dissolution in three-dimensional fracture networks: applications in geologic carbon sequestration

We perform a set of high-fidelity simulations of geochemical reactions within three-dimensional discrete fracture networks (DFN) and use various machine learning techniques to determine the primary factors controlling mineral dissolution. The DFN are partially filled with quartz that gradually dissolves until quasi-steady state conditions are reached. At this point, we measure the quartz remaining in each fracture within the domain as our primary quantity of interest. We observe that a primary sub-network of fractures exists, where the quartz has been fully dissolved out. This reduction in resistance to flow leads to increased flow channelization and reduced solute travel times. However, depending on the DFN topology and the rate of dissolution, we observe substantial variability in the volume of quartz remaining within fractures outside of the primary subnetwork. This variability indicates an interplay between the fracture network structure and geochemical reactions. We characterize the features controlling these processes by developing a machine learning framework to extract their relevant impact. Specifically, we use a combination of high-fidelity simulations with a graph-based approach to study geochemical reactive transport in a complex fracture network to determine the key features that control dissolution. We consider topological, geometric and hydrological features of the fracture network to predict the remaining quartz in quasi-steady state. We found that the dissolution reaction rate constant of quartz and the distance to the primary sub-network in the fracture network are the two most important features controlling the amount of quartz remaining. This study is a first step towards characterizing the parameters that control carbon mineralization using an approach with integrates computational physics and machine learning.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Solid-State Transformer and Hybrid Transformer With Integrated Energy Storage in Active Distribution Grids: Technical and Economic Comparison, Dispatch, and Control

Solid-state transformer (SST) and hybrid transformer (HT) are promising alternatives to the line-frequency transformer (LFT) in smart grids. The SST features medium-frequency isolation, full controllability for voltage regulation, reactive power compensation, and the capability of battery energy storage system (BESS) integration with multiport configuration. The HT has a partially-rated converter for fractional controllability and can integrate a small BESS. Fast grid-edge voltage fluctuations from increased solar photovoltaic (PV) and electric vehicle (EV) penetration are difficult to manage for mechanical load tap changers. Hence, along with the trend towards more BESS in the grid, the controllability and the storage integration capability of the SST and HT are of strong interest. However, a review of literature shows existing SST and HT research is mostly at converter level, while system-level assessments are scarce. Assessing technical and economic impacts is critical to understanding the benefits and role of the SST and HT to guide future research, which is presented for the first time in this article. Experimental results from medium-voltage (MV) SST and MV HT prototypes are shown to confirm equipment-level feasibility, where the voltage controllability waveforms of a MV HT prototype are reported for the first time. Comparative simulations are performed on a modified IEEE 34-bus system. Here, a grid-model-less decentralized grid-edge voltage control method and a day-ahead BESS dispatch method are proposed for the SST and HT. The simulations show that the SST and HT with integrated storage can host more PV, achieve peak shaving, mitigate voltage fluctuation and reverse power flow, and support energy arbitrage for operational cost reduction, as compared to the LFT. Moreover, comprehensive analyses of net present value (NPV) and internal rate of return (IRR) are performed under different installed PV capacities, HT’s partial converter ratings, and BESS capacities. Sensitivities to future cost reductions of the PV and BESS are studied. Although the NPV and IRR are currently negative, 60% capital cost reduction or 150% revenue increase will make the SST and HT economically viable in the use case studied.

14 SOLAR ENERGY↗
DOE OSTI2022

The influence of particle size on the fluid dynamics of a laser-induced plasma

The interaction of a laser-induced shock wave with nanoparticles and microparticles of aluminum oxide is investigated through experiments and modeling. The chemistry and physics of the interaction between the particles and plasma generated from laser ablation shows similarities and discrete differences for the two particle sizes. For both particle sizes, early stage (< 10 μs) ionization was dominant and evidenced by higher concentrations of Al II. While both sizes exhibit ionization over the same duration, the intensity of emission was greater for nanoparticles indicating greater concentrations of ionized species. Moreover, the dispersion of species was notably more elongated for microparticles while radial dispersion was more pronounced for nanoparticles with elevated drag forces. At later stages (i.e., > 10 μs), oxidation reactions were dominant for both particle sizes, but the same distinctions in flow field were observed and attributed to particle drag. In all stages of interaction, microparticles expand axially with less drag that suppresses their radial expansion. As a result, the dispersion of reactive species was mapped over an up to 80% larger area for nanoparticles relative to microparticles. Results shown here can be applied toward advancing experimental diagnostics and particle-shock wave modeling and simulation efforts for energetic materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A Field Enhancement Integration Design Featuring Misalignment Tolerance for Wireless EV Charging Using LCL Topology

This article proposes a magnetic integration design for EV wireless power transfer (WPT) systems, where the compensation and transmitting coils overlap one on top of the other to share the ferrite layer without any decoupling consideration. The magnetic field generated by both the compensation and transmitting coils are exploited to transfer power. To this end, a compensation method is proposed to enable magnetic field enhancement without any reactive power flow between the compensation and the transmitting coils, and to achieve input zero phase angle and constant current output. In addition, an efficient finite element analysis-based coil optimization algorithm is proposed to improve the coil misalignment tolerance on the horizontal plane, in which a reversely connected inner coil is used to stabilize the system output under misalignment conditions. Analytical and simulation results confirm transmission flux density enhancement and reduced leakage field characteristics of the proposed coil design. Finally, a scaled-down WPT prototype is built and tested to verify the performance and effectiveness of the design. The proposed design achieves 91.17% efficiency and misalignment tolerance up to 200 mm in any XY-direction while maintaining the power transfer and its efficiency. The results of this work provide insights into magnetic integration design of high-order compensation topologies featuring higher compactness, less ferrite usage, magnetic field enhancement, and misalignment tolerance for WPT systems.

42 ENGINEERING↗
DOE OSTI2023

Temperature-based reactive flow model for TATB plastic bonded explosives

A new reactive flow model is presented for TATB-based plastic bonded explosives, applicable to shock initiation and detonation problems of differing initial temperature. Temperature disequilibrium is assumed between unreacted explosive, material in the vicinity of compressed defects (called hot spots), and detonation products. The model incorporates temperature-dependent decomposition reaction rates. Arrhenius model parameters were derived from subscale quantum-based molecular dynamics simulations of shock-induced TATB decomposition. Further, a model of detonation carbon aggregation is incorporated, describing the late-time slow release of energy inherent to detonation in TATB-based materials. The predictive ability of the model in the shock initiation regime is tested against recent thin pulse experiments on PBX 9502. Here, the model is found to perform equally well in predicting the size effect curve of ambient, cold, and hot PBX 9502 rate sticks.

Chemistry - Chemical explosives↗
DOE OSTI2021

TOUGHREACT V4.12

TOUGHREACT V4.12 is a 3-D reactive transport code that handles to non-isothermal, multiphase/multicomponent fluid flow, heat transport, aqueous and gaseous species advection-diffusion, and equilibrium/kinetic water-gas-rock-biological reactions. TOUGHREACT V4.12 is based on TOUGHREACT V3.3, with several new features including: an option to compute activity coefficients using the Pitzer ion interaction model; implementation of ECO2N V2.0 for temperatures up to 300 C; options to simulate forward and reverse osmosis, including 2D Poiseuille flow; addition of EOS8 and revisions for EOS5 and EWASG with H2(g); addition of gas species decay, various new output formats and time step control options; many bug fixes.

Spycher, Nicolas↗
DOE OSTI2023

Electrochemical Residence Time Distribution as a Diagnostic Tool for Redox Flow Batteries

The fluid dynamic and electrochemical performance of redox flow batteries (RFBs) stems from the relationship between the flow field and the porous electrode, whose interplay determines how active species move and react during device operation. While characterization techniques, such as residence time distribution, offer insights into species mobility within a reactive volume for a traditional chemical reactor, electrochemical reactors also enable simultaneous measurement of the redox reactions, unlocking another dimension of analysis. Herein, we demonstrate how potentiodynamic measurements, using injections of electrolyte examined through moment analysis, can provide electrode-specific performance scaling relationships across a matrix of carbon paper and cloth electrodes with flow through and interdigitated flow fields. We further combine experimental campaigns with multiphysics simulations to demonstrate how electrode surface area can be estimated with this technique, which we then validate with activated and unactivated commercial carbon cloth electrodes. These studies reveal the multiscale observations that potentiodynamic measurements afford, augmenting existing electrochemical techniques for holistic electrochemical reactor diagnostics.

Electrochemistry↗