Results for “Flow chemistry” — Page 29

DOE OSTI2025

Energy-Efficient Capacitive Deionization through Electrode Modification and Process Development

Electrochemical separation technologies, such as capacitive deionization (CDI), are promising for addressing global energy and water challenges. However, there is a need to improve the performance, better understand property-performance relationships, and evaluate the longevity of CDI electrodes. This study explores the chemical modification of electrodes and the adjustment of CDI operating parameters. Results indicate that nitric acid (HNO3) conditioning of activated carbon cloth (ACC) electrodes removes metal oxides, introduces oxygen and nitrogen functionalities, and increases the specific capacitance (16% at 1 mV/s). Moreover, these changes in electrode properties positively impact device-level CDI performance. Through HNO3-conditioning of the ACC and tuning of the operational parameters, this work demonstrates higher electrosorption capacity (4.0x), greater charge efficiency (90% vs 24%), and lower energy consumption (3.8x). Despite these enhancements, limitations of the HNO 3 -conditioned ACC include decreased desorption kinetics and a 32% loss in electrosorption capacity after 200 cycles. Overall, this work provides guidance on using oxidative pretreatment via HNO 3 to modify ACC electrodes for CDI and evaluates the trade-offs associated with varying operational parameters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Modeling Microwave-Enhanced Chemical Vapor Infiltration Process for Preventing Premature Pore Closure

The chemical vapor infiltration (CVI) process involves infiltrating a porous preform with reacting gases that undergo chemical transformation at high temperatures to deposit the ceramic phase within the pores, ultimately leading to a dense composite. The conventional CVI process in composite manufacturing needs to follow an isothermal approach to minimize temperature differences between the external and internal surfaces of the preform, ensuring that reactive gases infiltrate internal pores before external surfaces seal. Here, this study addresses the challenge of premature pore closure in CVI processes through microwave heating. A frequency-domain microwave solver is developed in OpenFOAM to investigate volumetric heating mechanisms within the preform. Through numerical studies, we demonstrate the capability of microwave heating of creating an inside-out temperature inversion. This inversion accelerates reactions proximal to the preform center, effectively mitigating the risk of premature external pore closure and ensuring uniform densification. The results reveal a significant enhancement in temperature inversion when high-permittivity reflectors are incorporated to generate resonant waves. This microwave heating strategy is then coupled with high-fidelity direct numerical simulation (DNS) of reacting flow, enabling the analysis of resulting densification processes. The DNS includes detailed chemistry and realistic diffusion coefficients. The numerical results can be used to estimate the impact of microwave-induced temperature inversion on densification in productions.

42 ENGINEERING↗

DOE OSTI2020

Revisiting the Terawatt Challenge

Richard E. Smalley, in 2003, defined the Terawatt (TW) Challenge as “Adapting our energy infrastructure to simultaneously address diminishing oil resources and rising levels of atmospheric CO 2 .” Smalley, best known for the discovery of C 60 , for which he received the 1996 Nobel Prize in Chemistry, continued to address the challenges of anthropomorphic and natural global energy flows until he passed away in 2005. Smalley challenged the world to transform the energy sector. He envisioned electricity transmitted by high-voltage direct current (DC) lines from massively deployed solar plants in sunny areas and remotely sited nuclear plants. He also envisioned using advanced batteries for local storage of energy. To meet the needs of ~10 people in a world with a dwindling oil supply, Smalley asserted that the world would need to transform its fossil-fuel-driven 14-TW (average power) energy used in 2003 to a largely renewable-energy-driven 30–60 TW (average power) in 2050. This would be possible only if solar-electricity costs could be drastically reduced. The challenges associated with this transition have been called the “Terawatt Challenge.” Fifteen years later, solar-module costs have been reduced by tenfold and annual deployment of solar photovoltaic (PV) modules has grown by a factor of 100,from ~1 gigawatt (GW) in 2004 to ~100 GW in 2018, with a total of 500 GW installed worldwide, producing 2% of the planet’s electricity. As global installed solar generating capacity approaches1 TW, we revisit Smalley’s TW challenge to identify what has changed and quantify the TW Challenge for a baseline scenario and for two scenarios designed as upper and lower bounds determined by the degree we implement electrification and storage. In this paper, we show that the energy choices we make today will dramatically affect the magnitude of future global energy requirements.

SOLAR ENERGY↗

DOE OSTI2020

Dynamic Simulation Technoeconomic Model for Power Generation

Sandia National Laboratories has built and successfully tested a dynamic simulation technoeconomic model of the Palo Verde Generating Station that is now being updated to help other US power plants improve operations. Palo Verde, located west of Phoenix, Arizona, is the largest electricity generator in the US at 4 GW. Palo Verde uses — 60 million gallons per day of treated wastewater from Phoenix to cool reactors, and disposes of blowdown in evaporation ponds. The model built for Palo Verde numerically evaluates the economic impact of changing, for example, alternative cooling technologies, water usage and treatment, and influent water chemistry, and is based on detailed accounting of mass, energy, and cash flows.

24 POWER TRANSMISSION AND DISTRIBUTION↗

NASA NTRS1975

An analysis of communications blackout for Pioneer Venus entry probes

Estimates of the maximum duration of communications blackout are presented for a proposed Pioneer Venus mission involving a multiprobe entry targeted to individual impact sites. Calculations of the inviscid flow about the probes have been made using both finite-rate chemistry and equilibrium chemistry computer codes. Equilibrium boundary-layer calculations including ablation are also presented both with and without the presence of alkali metal contaminants from the heat shield. Methods are described for estimating maximum electron density in the inviscid wake and the recirculation region at the base of the probes. The results presented provide a conservative basis for planning for storage and subsequent telemetry of data during the entry phase of the mission.

Grose, W. L.↗

NASA NTRS1976

Second-order closure analysis of turbulent reacting flows

A reacting shear layer program for a complete second-order closure of turbulent reacting flows is developed. The program uses a typical eddy box model and many of the fluid mechanical turbulence models developed in earlier studies. The computer program actually solves a total of 30 equations, where the numerical integration of the equations is performed by a forward-time-centered space quasi-implicit upwind finite-difference scheme. Results for nonreacting and simple isothermal reacting flows are presented along with preliminary results for propane-air flames. It is shown that the program predictions for simple nonreacting flows are in good agreement with available data. The calculations for isothermal reacting flows demonstrate the importance of including a mixedness correlation in finite rate chemistry calculations, where significant errors in the flowfield properties can result from neglect of the interaction between turbulence and the chemistry. Results for propane-air diffusion flames further support this conclusion.

Varma, A. K.↗

DOE OSTI2022

Lignin-First Biorefinery Development

Cost-effective biomass fractionation is an enabling, grand challenge for biorefining, especially when both carbohydrates and lignin are targeted for valorization. The advent of Reductive Catalytic Fractionation (RCF) - an active stabilization approach that solubilizes lignin from biomass and catalytically depolymerizes it into a narrow slate of monomers and oligomers - represents a potential step forward for this important goal. In the Lignin-First Biorefinery Development project, we are employing techno-economic analysis (TEA) and life cycle assessment (LCA) to guide bench-scale R&D efforts towards cost-effective RCF-based biorefining. This project is critical to lignin valorization efforts and could enable the use of woody feedstocks in a traditional biochemical conversion context. To date, we have conducted the first rigorous TEA and LCA study of the RCF process, developed a flow-through system to separate biomass and the chemical catalyst in RCF chemistry, developed models for solvolysis chemistry and transport phenomena for poplar, and co-led an international, authoritative perspective on guidelines for the research community on how to best practice lignin-first biorefining. The primary challenges for the lignin-first RCF process going forward are catalyst stability, the need to utilize, recover, and recycle high boiling point solvents to lower RCF reactor pressure, and the challenge of operating RCF continuously - all of which are being directly tackled by this project.

biomass↗

DOE OSTI2020

Engineering scale molten salt corrosion and chemistry code development

A new engineering scale transport code Mole is described. Mole was developed using the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework for solving species transport problems in flowing liquid. It is applied to leaching of Cr from alloys into molten salts part a of a multi-physics, multi-scale approach for coupling chemistry with mechanical behavior modeling. Demonstrations of the leaching phenomenon were benchmarked to experimental data provided by the Molten Salt Reactor campaign and the results are presented.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

NASA NTRS1972

The chemistry of neutral species in the D and E regions.

The photodissociation of molecular species such as oxygen, nitrogen, and water vapor in the upper atmosphere leads to the formation of a reservoir of chemical energy consisting of highly reactive atomic gases. The reactions of these with each other, and with the rest of the atmosphere, give a complex neutral chemistry involving many species and reactions. Local equilibrium can rarely be assumed for these reactions; transport by the horizontal and vertical mean flow, by eddy diffusion, and by molecular diffusion has quite different effects on the various constituents, depending on the time constants. The entire subject of the neutral chemistry in the D and E regions (between 50 and 160 km altitude) is reviewed, and related to existing experimental results. Needs for further research on the subject are indicated.

Bowhill, S. A.↗

NASA NTRS2026

Thermochemically-Closed Sonic-Flow Inversion for Enthalpy and Temperature in Multispecies Arc-Jet Flows

A thermochemically-closed sonic-flow inversion framework (TSIF) is developed to infer bulk enthalpy and total temperature upstream of a choked nozzle in arc-jet flows. The formulation recasts a pressure-rise total enthalpy quantification technique as an inverse problem in characteristic-velocity c * space using measured mass flow rate, upstream total pressure, gas composition, and nozzle throat geometry as inputs. Unlike calorimetric energy-balance approaches or optical diagnostics, the method relies primarily on routinely measured facility quantities combined with explicit thermochemical closure. Thermochemical states are obtained using NASA’s open-source Chemical Equilibrium with Applications (CEA) code, enabling construction of a chemistry-consistent relation between characteristic velocity, total enthalpy, and total temperature under equilibrium or frozen assumptions. A discharge coefficient is self-calibrated using cold-flow (arc-off) operation data and applied to hot-flow (arc-on) measurements, enabling upstream losses to be accounted for without empirical correlations. The framework is applied to air, N 2 , and CO 2 –N 2 arc-jet flows and demonstrates expected trends for the inferred thermochemical states as function of arc power, specific energy input, mass-flow, heater configuration, and test gas. In the air limit, under equilibrium assumptions, the method recovers the classical high-enthalpy asymptotic correlation of Winovich with a mean residual of 4.4%, demonstrating compatibility with established sonic-flow scaling, while extending applicability to arbitrary multi-species mixtures and non-equilibrium chemistry. The framework provides a mixture-flexible methodology for determining bulk thermochemical states in modern arc-jet environments using routine facility pressure, mass-flow, gas-composition, and nozzle-geometry information together with a cold-flow calibration.

stagnation heat flux↗

NASA NTRS2026

Thermochemically-Closed Sonic-Flow Inversion for Enthalpy and Temperature in Multispecies Arc-Jet Flows

A thermochemically-closed sonic-flow inversion framework (TSIF) is developed to infer bulk enthalpy and total temperature upstream of a choked nozzle in arc-jet flows. The formulation recasts a pressure-rise total enthalpy quantification technique as an inverse problem in characteristic-velocity c * space using measured mass flow rate, upstream total pressure, gas composition, and nozzle throat geometry as inputs. Unlike calorimetric energy-balance approaches or optical diagnostics, the method relies primarily on routinely measured facility quantities combined with explicit thermochemical closure. Thermochemical states are obtained using NASA’s open-source Chemical Equilibrium with Applications (CEA) code, enabling construction of a chemistry-consistent relation between characteristic velocity, total enthalpy, and total temperature under equilibrium or frozen assumptions. A discharge coefficient is self-calibrated using cold-flow (arc-off) operation data and applied to hot-flow (arc-on) measurements, enabling upstream losses to be accounted for without empirical correlations. The framework is applied to air, N 2 , and CO 2 –N 2 arc-jet flows and demonstrates expected trends for the inferred thermochemical states as function of arc power, specific energy input, mass-flow, heater configuration, and test gas. In the air limit, under equilibrium assumptions, the method recovers the classical high-enthalpy asymptotic correlation of Winovich with a mean residual of 4.4%, demonstrating compatibility with established sonic-flow scaling, while extending applicability to arbitrary multi-species mixtures and non-equilibrium chemistry. The framework provides a mixture-flexible methodology for determining bulk thermochemical states in modern arc-jet environments using routine facility pressure, mass-flow, gas-composition, and nozzle-geometry information together with a cold-flow calibration.

inviscid theory↗

DOE OSTI2024

Advancing a Watershed Hydro-Biogeochemical Theory: Linking Water Travel Time and Reaction Rates Under Changing Climate

The major goal of the project is to answer the question how and to what degree does a warming climate alter hydrological functioning and (biogeochemical and chemical weathering) reaction rates at the watershed scale? In particular, how does hydrological functioning (water storage, flow paths, and travel time) change responding to a warming climate? What are the corresponding alterations in biogeochemical and chemical weathering rates at the watershed scale? What is the general watershed-scale reaction rate law that causally links reaction rates to metrics of hydrological functioning that are determined by external hydroclimatic conditions and internal watershed structure characteristics? To answer these questions, we have used field measurements and reactive transport modeling. We have found that: 1) stream chemistry has altered tremendously in warmer years; 2) In this water-limited environment, hydrology is the dominant driver for soil respiration and carbon lateral export at the watershed scale; 3) deeper carbon respiration plays an essential role in regulating the patterns of stream carbon in low and high flow conditions, and deeper groundwater will become increasingly important as warming intensifies. This indicates in mountain streams, a warming climate will modify not only the timing and magnitude of streamflow as expected, but also the chemistry and quality of stream water.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Monitoring the long-term performance of organic redox flow battery by a distribution of relaxation time analysis

Organic redox flow batteries hold great promise as an energy storage technology, but their intricate chemistry makes them vulnerable to various degradation mechanisms. Monitoring this degradation is essential for identifying the limiting processes within the cells. Electrochemical impedance spectroscopy (EIS) offers a straightforward, in-situ method for measuring the total resistance of an operating cell. However, to pinpoint the limiting processes during long-term cycling, EIS data must be complemented by other techniques. Distribution of relaxation time (DRT) analysis is particularly effective for differentiating resistance components. Here, in this study, we perform a comprehensive analysis of resistance evolution and the separation of anode and cathode contributions during long-term cycling of a full cell employing 7,8-dihydroxyphenazine-2-sulfonic acid (DHPS) as the anolyte. Separate analyses of the DHPS anolyte and ferri-/ferrocyanide catholyte were conducted using a symmetric cell setup. The relaxation times derived from symmetric cells facilitate the identification of peaks in the DRT profiles from the full cell. Importantly, the DRT profiles indicate a correlation between the evolution of charge transfer resistance and the chemical degradation of DHPS. The methodologies and results outlined in this study offer significant insights for developing diagnostic tools applicable to other types of redox flow batteries.

Distribution of relaxation time↗

NASA NTRS1995

Computations of Axisymmetric Flows in Hypersonic Shock Tubes

A time-accurate two-dimensional fluid code is used to compute test times in shock tubes operated at supersonic speeds. Unlike previous studies, this investigation resolves the finer temporal details of the shock-tube flow by making use of modern supercomputers and state-of-the-art computational fluid dynamic solution techniques. The code, besides solving the time-dependent fluid equations, also accounts for the finite rate chemistry in the hypersonic environment. The flowfield solutions are used to estimate relevant shock-tube parameters for laminar flow, such as test times, and to predict density and velocity profiles. Boundary-layer parameters such as bar-delta(sub u), bar-delta(sup *), and bar-tau(sub w), and test time parameters such as bar-tau and particle time of flight t(sub f), are computed and compared with those evaluated by using Mirels' correlations. This article then discusses in detail the effects of flow nonuniformities on particle time-of-flight behind the normal shock and, consequently, on the interpretation of shock-tube data. This article concludes that for accurate interpretation of shock-tube data, a detailed analysis of flowfield parameters, using a computer code such as used in this study, must be performed.

Sharma, Surendra P.↗

DOE OSTI2025

A physics-constrained neural ordinary differential equations approach for robust learning of stiff chemical kinetics

The high computational cost associated with solving for detailed chemistry poses a significant challenge for predictive computational fluid dynamics (CFD) simulations of turbulent reacting flows. While deep learning techniques have been explored to develop faster surrogate models, they often fail to integrate reliably with CFD solvers. This instability arises because traditional deep learning approaches optimize for training error without ensuring compatibility with ordinary differential equation (ODE) solvers, resulting in accumulation of errors over time. Recently, neuralODE (NODE) based approaches have been shown to be a promising technique to emulate and accelerate detailed chemistry computations. Here, in the present work, we extend this NODE framework for stiff chemical kinetics by incorporating mass conservation constraints directly into the loss function during training. This ensures that the total mass as well as the individual elemental species masses are conserved in an a-posteriori manner. Proof-of-concept studies are performed with the novel physics-constrained NODE (PC-NODE) approach for homogeneous autoignition of hydrogen-air mixture over a range of composition and thermodynamic conditions. It is demonstrated that the PC-NODE framework not only improves the physical consistency of the resulting data-driven model with respect to mass conservation criteria, but also improves training efficiency. PC-NODE is shown to achieve 2–100× speedup relative to the hydrogen-air detailed chemical mechanism depending on the type of the ODE solver (implicit or explicit) used during autoregressive inference tests. Lastly, a-posteriori studies are performed wherein the trained PC-NODE model is coupled with a CFD solver. It is shown that higher accuracy is achieved with PC-NODE relative to the purely data-driven NODE approach. Moreover, PC-NODE also exhibits robustness and generalizability to unseen initial conditions from within (interpolative capability) as well as outside (extrapolative capability) the training regime.

computational combustion↗

NASA NTRS1987

An analysis of artificial viscosity effects on reacting flows using a spectral multi-domain technique

Standard techniques used to model chemically-reacting flows require an artificial viscosity for stability in the presence of strong shocks. The resulting shock is smeared over at least three computational cells, so that the thickness of the shock is dictated by the structure of the overall mesh and not the shock physics. A gas passing through a strong shock is thrown into a nonequilibrium state and subsequently relaxes down over some finite distance to an equilibrium end state. The artificial smearing of the shock envelops this relaxation zone which causes the chemical kinetics of the flow to be altered. A method is presented which can investigate these issues by following the chemical kinetics and flow kinetics of a gas passing through a fully resolved shock wave at hypersonic Mach numbers. A nonequilibrium chemistry model for air is incorporated into a spectral multidomain Navier-Stokes solution method. Since no artificial viscosity is needed for stability of the multidomain technique, the precise effect of this artifice on the chemical kinetics and relevant flow features can be determined.

Macaraeg, M. G.↗

NASA NTRS1989

Analysis of artificial viscosity effects on reacting flows using a spectral multidomain technique

Standard techniques used to model chemically-reacting flows require an artificial viscosity for stability in the presence of strong shocks. The resulting shock is smeared over at least three computational cells, so that the thickness of the shock is dictated by the structure of the overall mesh and not the shock physics. A gas passing through a strong shock is thrown into a nonequilibrium state and subsequently relaxes down over some finite distance to an equilibrium end state. The artificial smearing of the shock envelops this relaxation zone which causes the chemical kinetics of the flow to be altered. A method is presented which can investigate these issues by following the chemical kinetics and flow kinetics of a gas passing through a fully resolved shock wave at hypersonic Mach numbers. A nonequilibrium chemistry model for air is incorporated into a spectral multidomain Navier-Stokes solution method. Since no artificial viscosity is needed for stability of the multidomain technique, the precise effect of this artifice on the chemical kinetics and relevant flow features can be determined.

Macaraeg, Michele G.↗

NASA NTRS1993

Study and modeling of finite rate chemistry effects in turbulent non-premixed flames

The development of numerical models that reflect some of the most important features of turbulent reacting flows requires information about the behavior of key quantities in well defined combustion regimes. In turbulent flames, the coupling between turbulent and chemical processes is so strong that it is extremely difficult to isolate the role played by one individual physical phenomenon. Direct numerical simulation (hereafter DNS) allows us to study in detail the turbulence-chemistry interaction in some restricted but completely defined situations. Globally, non-premixed flames are controlled by two limiting regimes: the fast chemistry case, where the turbulent flame can be pictured as a random distribution of local chemical equilibrium problems; and the slow chemistry case, where the chemistry integrates in time the turbulent fluctuations. The Damkoehler number, ratio of a mechanical time scale to chemical time scale, is used to distinguish between these regimes. Today most of the industrial computer codes are able to perform predictions in the hypothesis of local equilibrium chemistry using a presumed shape for the probability density function (pdt) of the conserved scalar. However, the finite rate chemistry situation is of great interest because industrial burners usually generate regimes in which, at some points, the flame is undergoing local extinction or at least non-equilibrium situations. Moreover, this variety of situations strongly influences the production of pollutants. To quantify finite rate chemistry effect, the interaction between a non-premixed flame and a free decaying turbulence is studied using DNS. The attention is focused on the dynamic of extinction, and an attempt is made to quantify the effect of the reaction on the small scale mixing process. The unequal diffusivity effect is also addressed. Finally, a simple turbulent combustion model based on the DNS observations and tractable in real flow configurations is proposed.

Vervisch, Luc↗