Search NASA⌕ Search

SEARCH · Search NASA

Results for “Mixing thermodynamics”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Simulation of a rapid compression machine for evaluation of ignition chemistry and soot formation using gasoline/ethanol blends

Due to the projected decline of demand for gasoline in light duty engines and the advent of ethanol as a green fuel, the use of gasoline/ethanol blend fuels in heavy duty applications are being investigated as they are projected to have lower cost and lower lifecycle green house gas (GHG) emissions. In heavy duty engines, the primary mode of combustion is mixing controlled combustion where wide range of mixture conditions (equivalence ratio) exist. Soot emissions of these fuels in richer conditions are not well understood. The goal of this research is to evaluate some commercially available soot modeling codes for the particulate matter emissions from gasoline/ethanol fuel blends, especially at fuel rich conditions. A Rapid Compression Machine (RCM) is modeled in a three-dimensional numerical simulation using CONVERGE computational software using a reduced chemical kinetic mechanism with SAGE chemistry solver and a RANS k-ϵturbulence model with a sector model including the creviced piston. The creviced piston is used in the experimental setup to reduce boundary layer effects and to maintain a homogeneous core in the reaction cylinder. Computational fluid dynamics simulations are conducted for different gasoline-ethanol fuel blends from E10 (10% ethanol v/v) to E100. The fuel blend is modeled as a surrogate mixture of toluene, iso-octane, n-heptane for gasoline content, and ethanol. The computational results were validated against experimental results using pressure measurements and laser extinction diagnostics. Different soot models are investigated to evaluate their capability of predicting the sooting tendencies of fuel blends, especially in richer conditions experienced during mixing-controlled combustion. The experimental combustion characteristics such as the ignition delay of different blends of fuel are reasonably well predicted. The Particulate Size Mimic (PSM) model accurately predicts the soot generation characteristics of the different fuels, but the Hiroyasu-NSC model falls short in this regard. For accurate prediction of soot with the PSM model, the thermodynamic conditions during combustion must be accurately modeled. While the current computational modeling tools can produce accurate results for the prediction of particulate matter emissions, there is much work to be done in improving our understanding of the underlying fundamental processes.

Energy & Fuels↗
DOE OSTI2022

Octane Requirements of Lean Mixed-Mode Combustion in a Direct-Injection Spark-Ignition Engine

Here, this study investigates the octane requirements of a hybrid flame propagation and controlled autoignition mode referred to as mixed-mode combustion (MMC), which allows for strong control over combustion parameters via a spark-initiated deflagration phase. Due to the throughput limitations associated with both experiments and 3-D computational fluid dynamics calculations, a hybrid 0-D and 1-D modeling methodology was developed, supported by experimental validation data. This modeling approach relied on 1-D, two-zone engine simulations to predict bulk in-cylinder thermodynamic conditions over a range of engine speeds, compression ratios, intake pressures, trapped residual levels, fueling rates, and spark timings. Those predictions were then transferred to a 0-D chemical kinetic model, which was used to evaluate the autoignition behavior of fuels when subjected to temperature–pressure trajectories of interest. Finally, the predicted autoignition phasings were screened relative to the progress of the modeled deflagration-based combustion in order to determine if an operating condition was feasible or infeasible due to knock or stability limits. The combined modeling and experimental results reveal that MMC has an octane requirement similar to modern stoichiometric spark-ignition engines in that fuels with high research octane number (RON) and high octane sensitivity (S) enable higher loads. Experimental trends with varying RON and S were well predicted by the model for 1000 and 1400 rpm, confirming its utility in identifying the compatibility of a fuel’s autoignition behavior with an engine configuration and operating strategy. However, the model was not effective in predicting (nor designed to predict) operability limits due to cycle-to-cycle variations, which experimentally inhibited operation of some fuels at 2000 rpm. Putting the operable limits and efficiency from MMC in the context of a state-of-the-art engine, the MMC showed superior efficiencies over the range investigated, demonstrating the potential to further improve fuel economy.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Thermodynamic modeling of the Nb-Ni system with uncertainty quantification using PyCalphad and ESPEI

Here, the Nb–Ni system is remodeled with uncertainty quantification (UQ) using software tools of PyCalphad and ESPEI (the Extensible, Self-optimizing Phase Equilibria Infrastructure) with the presently implemented capability of modeling site fraction based on Wyckoff positions. The five- and three-sublattice models are used to model the topologically close pack (TCP) μ-Nb 7 Ni 6 and δ-NbNi 3 phases according to their Wyckoff positions. The inputs for CALPHAD-based thermodynamic modeling include the thermochemical data as a function of temperature predicted by first-principles and phonon calculations based on density functional theory (DFT), ab initio molecular dynamics (AIMD) simulations, together with phase equilibrium and site fraction data in the literature. In addition to phase diagram and thermodynamic properties, the CALPHAD-based predictions of site fractions of Nb in μ-Nb 7 Ni 6 agree well with experimental data. Furthermore, the UQ estimation using the Markov Chain Monte Carlo (MCMC) method as implemented in ESPEI is applied to study the uncertainty of site fraction in μ-Nb 7 Ni 6 and enthalpy of mixing (ΔH mix ) in liquid.

36 MATERIALS SCIENCE↗
DOE OSTI2026

The Factors Governing Metal Dependence of an Emergent Superfamily of Bimetallic Oxygenases

Metalloenzyme superfamilies are typically defined by their protein scaffolds and active sites. Owing to the high tunability of protein structures, members of a single superfamily can catalyze diverse reactions with the same metallocofactor. Some superfamilies, such as amidohydrolase-related dinuclear oxygenases (AROs), display further versatility by utilizing multiple metallocofactors. We have shown that certain AROs catalyze monooxygenation reactions with diiron, dimanganese, and/or mixed manganese−iron cofactors, but the molecular factors governing the selection of a particular cofactor remain unknown, and the extent of this superfamily in biology is unclear. Here, we report bioinformatic analyses that expand the ARO superfamily to approximately 17,000 unique UniProt sequences, far exceeding the number of previously characterized enzymes. Through the integration of structural, spectroscopic, and thermodynamic analyses of representative proteins with a bioinformatic pipeline that identifies key secondary- and tertiary-sphere residues, we can predict in silico the metal preference for the majority of reported ARO sequences. These annotations were validated via the characterization of multiple new AROs, including ones implicated in key oxidative steps of natural product biosyntheses. This study establishes the key structure−function relationships governing metal preferences in AROs and highlights their vastly underappreciated role in myriad biological processes.

Liu, Chang [University of California, Berkeley, CA↗
DOE OSTI2025

Electrified vapour deposition at ultrahigh temperature and atmospheric pressure for nanomaterials synthesis

Vapour-phase synthesis methods have shown promise for the scalable synthesis of nanomaterials and coatings. However, the vaporization of different precursors for the synthesis of a broad nanomaterial space, particularly at atmospheric pressure, while maintaining compositional and structural control of the final product is challenging. Here we report the generation of an ultrahigh-temperature atomic vapour at atmospheric pressure based on electrified heating, for the growth of multi-elemental nanomaterials and thin films. This process relies on a reactor design whereby solid-state precursors are vaporized within a semi-confined space beneath an electrified heater that can reach ~3,000 K. The proximity of the heater rapidly breaks down the bonds of metal salt precursors and decomposes them into an atomic vapour that expands into a high-temperature (>2,000 K), highly reactive and high-flux vapour (10 21 –10 22 atoms per cm 2 per second) that travels upwards in a directional flow. When mixed with entrained ambient gases, the highly reactive atomic species rapidly nucleate and grow into the desired final products, including alloys, oxides, sulfides and thin films, which can be deposited on a low-temperature substrate. This EVD approach can synthesize a broad range of functional nanomaterials at atmospheric pressure, including single-phase multi-elemental nanomaterials formed under thermodynamically non-equilibrium conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

REDOTHERM (Redox Countercurrent Thermodynamic Limits Model) [SWR-24-88]

REDOTHERM is an open-source, MATLAB-based thermodynamic modeling framework developed to evaluate the performance of redox-active materials for water (H2O) and carbon dioxide (CO2) splitting. It includes models of all major unit operations and supports comparative analysis of different redox-active material candidates. The model is tailored for systems of moving oxide under a parallel/cocurrent flow (PF) and countercurrent flow (CF) configurations. Unvalidated mixed flow reactor (MFR, also known as CSTR) model is also included as an optional addition.

Lidor, Alon [National Renewable Energy Laboratory ↗
DOE OSTI2024

Thermodynamic modeling of CsF with LiF-NaF-KF for molten fluoride-fueled reactors

Gibbs energy models were developed to describe the thermochemical behavior of CsF in molten FLiNaK (46.5LiF-11.5NaF-42KF mol%), a proposed molten salt reactor (MSR) fuel solvent and coolant, as cesium is of concern due to its high radiotoxicity and volatility. Initially, it was necessary to obtain a more accurate Gibbs energy function for CsF which required fitting parameters to reported vapor pressures over condensed phase CsF. The pseudo-binary systems CsF-LiF, CsF-NaF and CsF-KF were then evaluated utilizing phase equilibria and enthalpy of mixing (Δ mix H) values, together with original differential scanning calorimetry (DSC) measurements performed for the CsF-LiF and CsF-KF systems. The CsF-LiF-NaF, CsF-LiF-KF and CsF-NaF-KF pseudo-ternary system representations were obtained by interpolation of the constituent pseudo-binary systems, with DSC measurements performed for the CsF-LiF-NaF system to corroborate the calculated liquidus temperature. Ultimately, the pseudo-ternary systems were interpolated to obtain Gibbs energy models for the pseudo-quaternary CsF-LiF-NaF-KF system, supported by DSC measurements at low CsF compositions (1–10 mol%), yielding computed equilibria and cesium-containing vapor pressures that compare favorably with reported values. In conclusion, the Molten Salt Thermal Properties Database – Thermochemical (MSTDB-TC) was subsequently expanded to include these Gibbs energy models allowing description of the thermochemical behavior of the CsF-LiF-NaF-KF system.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Stabilization of Miscible Aqueous Phases via Diffusion‐Controlled Multifunctional Nanoparticle‐Ligand Complexation

Liquid‐in‐liquid structuring by harnessing miscible aqueous domains can be achieved by inducing thermodynamically defined phase separation using aqueous solutions of incompatible solutes, where immiscibility is dictated by the concentrations of the solutes. For instance, in aqueous two‐phase systems (ATPSs), the interfacial tension between the two different aqueous solutions is so small that it stabilizes the phase‐separated domains, but not their shape, thus falling short of designing tunable, robust structures. Here, we introduce a diffusion‐controlled strategy that enables liquid‐in‐liquid compartmentalization absent a thermodynamically defined interface or bulk liquid‐liquid phase separation using a barrier formed in situ at the initial contact boundary of an aqueous dispersion of multifunctional nanoparticles and a solution of multifunctional ligands in pure water or a water/alcohol mixture. The nanoparticles and ligands are entropically driven to disperse uniformly but interact at the initial contact boundary more rapidly than bulk mixing, forming a layer that establishes a barrier akin to a semi‐permeable membrane and allows the formation of tubular compartments. A wealth of opportunities is enabled by the proposed technique, where examples in separations, soft electronic devices, and ultralight porous electromagnetic shielding materials are shown.

cryogels↗
DOE OSTI2022

Comparison of Na 2 SO 4 , K 2 SO 4 and Na 2 SO 4 -K 2 SO 4 deposit induced hot corrosion of a β-NiAl coating

Here, hot corrosion behaviour of a β-NiAl coating induced by Na 2 SO 4 , K 2 SO 4 , and Na 2 SO 4 -K 2 SO 4 deposits was investigated at 700 °C in air + 300 ppm SO 2 . To provide an insight into the corrosion mechanisms, phase equilibria of the Na 2 SO 4 -K 2 SO 4 -NiO-SO 3 system were calculated using an in-house developed thermodynamic database. The attack rates and morphologies could be correlated with the deposit compositions and their reactions with oxides. K 2 SO 4 and Na 2 SO 4 -K 2 SO 4 deposits are more corrosive than Na 2 SO 4 . This is attributed to a reduction of the melting temperature of the deposit and the formation of Ni-containing mixed sulphate phases in the presence of K 2 SO 4 .

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

Methods and apparatus to facilitate alkali metal transport during battery cycling, and batteries incorporating same

An anode includes a mixed ionic-electronic conductor (MIEC) with an open pore structure. The open pore structure includes open pores to facilitate motion of an alkali metal into and/or out of the MIEC. The open pore structure thus provides open space to relieve the stresses generated by the alkali metal when charging/discharging a battery. The MIEC is formed from a material that is thermodynamically and electrochemically stable against the alkali metal to prevent the formation of solid-electrolyte interphase (SEI) debris and the formation of dead alkali metal. The MIEC may also be passive (the MIEC does not store or release alkali metal). In one example, the open pore structure may be an array of substantially aligned tubules with a width less than about 300 nm, a wall thickness between about 1 nm to about 30 nm, and a height of at least 10 um arranged as a honeycomb.

Li, Ju↗
DOE OSTI2021

Adsorbate Free Energies from DFT-Derived Translational Energy Landscapes

Adsorption free energies are fundamental to surface chemistry and catalysis. Standard models combine some assumed analytical form of the translational potential energy surface, often parameterized against density functional theory (DFT) calculations, with an analytical expression for the resultant translational densities of states (DOS), free energy, and entropy. Here we compare the performance of such models against numerical evaluations of the DOS and thermodynamic functions derived from solutions to the translational Schrodinger equation. We compare results for a translational potential energy surface (PES) derived from nudged eleastic band calculations with those obtained from adsorbate rastering across a series of monatomic (O, S, C, N, H) and polyatomic (NH x ) adsorbates on (100) Pt and Au facets. We find that analytical models as commonly parameterized have mixed performance for describing the translational PES and that the consequences for computed free energies are modest but potentially significant in microkinetic models. Here, numerical solutions are possible for modest to no additional computational cost over analytical models and thus should be considered when reliable free energy estimates are needed or translational potential energy surfaces are available.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Thermodynamic controls on the synthesis, structure and reactivity of materials for energy (Final Report)

The major objectives have been: a) to advance and use unique calorimetric capabilities to determine the energetics of oxide and other materials having possible energy applications, with an emphasis on new materials for batteries, fuel cells and solar energy conversion, b) to understand, in terms of both macroscopic energetics and microscopic structure and bonding, the interplay of defect chemistry, oxidation-reduction, and size effects at the bulk and nanoscale in determining materials properties. This project emphasizes mixed ionic-electronic conduction materials containing transition metals of variable valence, especially perovskites and lithium containing phases, with M = Mn, Fe, Co, Ni. The thermochemistry of new families of complex halide, chalcogenide, nitride and hybrid (organic-inorganic) perovskite materials were explored. The major and unique technique employed is high temperature oxide melt solution calorimetry, with careful attention given to materials synthesis and characterization. The thermochemical studies will be coupled with structural studies, especially to characterize order-disorder phenomena by in situ high temperature diffraction techniques.

14 SOLAR ENERGY↗
DOE OSTI2025

Insights into the Oxidation Mechanism of Vivianite to Metavivianite from First-Principles Calculations

Vivianite, a hydrous ferrous iron-phosphate mineral (Fe 3 (PO 4 ) 2 ·8H 2 O), readily oxidizes in contact with air yielding the less hydrous mixed-valent iron-phosphate mineral metavivianite. This topotactic transformation nominally occurs by oxidative dehydrogenation in which outgoing electrons from the iron sublattice are charge compensated by hydrolysis of structural water. However, the details of this internal charge balancing mechanism that allows the structure to remain electrostatically stable remain unknown. Here, in this study, we use density functional theory (DFT) calculations and ab initio thermodynamics (AIT) to evaluate the energetics of this process in terms of hydrogen release as a function of environmental variables such as partial pressure and temperature. The results show a thermodynamic driving force for vivianite phase transformation as oxidation progresses that is triggered by its rigid structure that has a limited accommodation for hydrogen vacancies. In contrast, metavivianite has a more flexible lattice and hydrogen bond network that stabilizes these hydrogen defects. Metavivianite is shown to be a stable intermediate for 66% residual Fe 2+ down to 33% Fe 2+ , below which other phases/structures such as santabarbarite should be more thermodynamically favorable. Our study provides a basis for experimental tests of our mechanistic findings and helps fill a basic knowledge gap about the solid-state process that defines how vivianite interacts with its surrounding environment.

Sassi, Michel [Pacific Northwest National Laborato↗
DOE OSTI2022

Local cation order and ferrimagnetism in compositionally complex spinel ferrites

We present an exploration of a family of compositionally complex cubic spinel ferrites featuring combinations of Mg, Fe, Co, Ni, Cu, Mn, and Zn cations, systematically investigating the average and local atomic structures, chemical short-range order, magnetic spin configurations, and magnetic properties. All compositions result in ferrimagnetic average structures with extremely similar local bonding environments; however, the samples display varying degrees of cation inversion and, therefore, differing apparent bulk magnetization. Additionally, first-order reversal curve analysis of the magnetic reversal behavior indicates varying degrees of magnetic ordering and interactions, including potentially local frustration. Finally, reverse Monte Carlo modeling of the spin orientation demonstrates a relationship between the degree of cation inversion and the spin collinearity. Collectively, these observations correlate with differences in synthesis procedures. This work provides a framework for understanding magnetic behavior reported for “high-entropy spinels,” revealing many are likely compositionally complex oxides with differing degrees of chemical short-range order—not meeting the community established criteria for high or medium entropy compounds. Moreover, this work highlights the importance of reporting complete sample processing histories and investigating local to long-range atomic arrangements when evaluating potential entropic mixing effects and assumed property correlations in high entropy materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

High-pressure phase diagram of beryllium from ab initio free-energy calculations

In this report we use first-principles molecular dynamics simulations coupled with the thermodynamic integration method to study the hexagonal close-packed (hcp) to body-centered cubic (bcc) transition and melting of beryllium up to a pressure of 1600 GPa. We derive the melting line by equating solid and liquid Gibbs free energies and represent it by a Simon-Glatzel fit T m = 1564 K [1 + P/(15.6032 GPa)] 0.383 , which is in good agreement with previous two-phase simulations <6000 K. We also derive the hcp-bcc solid-solid phase boundary and show that the quasiharmonic approximation underestimates the stability of the hcp structure, predicting lower transition pressures between hcp and bcc phases. Our results are consistent with the stability regime predicted by the phonon quasiparticle method. We also predict that the hcp-bcc-liquid triple point is located at 164.7 GPa and 4314 K. In addition, we compute the shock Hugoniot curve and show that it is in good agreement with experiments, intersecting our derived melting curve at ~235 GPa and 4900 K. Finally, we make predictions for future ramp compression experiments. Starting with an isentropic compression of the liquid, we predict the path to intersect the melting line at low pressure and temperature, then to continue along the melting line over a large temperature interval of 7000 K as the sample remains in the mixed solid-liquid state before it enters the solid phase.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Vortex Structure of Head Bubble in Convective Cloud Starting Plume

Abstract A developing cumulus cloud (Cu) was modeled, and dynamic, thermodynamic, and microphysical properties of an ascending head bubble reproducing the upper part of a developing Cu were investigated. The data for analysis are taken from 10-m-resolution LES of trade wind Cu under BOMEX conditions. The detection of a rising bubble is carried out using wavelet filtering of the velocity fields and microphysical fields, while a low-frequency signal of the filtering is associated with the convective-scale structure of cloud. We substantiate and discuss the representation of the bubble as a vortex ring, and estimate the parameters of this vortex ring. The simplest Hill’s vortex was chosen as a model of a vortex ring inside cloud. Analytical approximations of the radial profiles of the vertical velocity and of conservative quantities (such as total water mixing ratio and liquid water potential temperature inside and outside the bubble) are obtained. The spatial structure of these quantities is investigated using analytical expressions. Analytical models for spatial distributions of liquid water content (LWC) and adiabatic fraction (AF) are also designed and analyzed. The results demonstrate the existence of a cloud core with high values of LWC and AF up to the height of 1800 m. The horizontally averaged value of the adiabatic fraction, calculated analytically using the Hill’s vortex concept, is evaluated as 0.39, which is the typical AF value in the upper parts of such Cu. The vertical profiles of different important quantities characterizing cloud structure are presented. The analysis performed in this study allows us to conclude that a rising vortex ring plays the dominating role in formation of the thermodynamic and microphysical structure of developing Cu. Significance Statement 1) Dynamic and thermodynamic fields of a developing cumulus cloud simulated by high-resolution LES with spectral bin microphysics are separated into convective and turbulent components by means of the wavelet technique. 2) The analysis of convective component of the cloud revealed the existence of a vortex ring at the developing stage of the cloud and evaluate its parameters. 3) The analysis performed in this study allows us to conclude that a rising vortex ring plays the important role in formation of the thermodynamic and microphysical structure of developing Cu. 4) The study provides a novel insight into the cloud–environment interaction. 5) The approximating equations describing the vortex ring can be usefully applied for developing new schemes of convective parameterization.

Meteorology & Atmospheric Sciences↗
DOE OSTI2025

Atomic-scale understanding of oxide growth and dissolution kinetics of Ni-Cr alloys

Aqueous corrosion of metals is governed by formation and dissolution of a passivating, multi-component surface oxide. Unfortunately, a detailed atomistic description is challenging due to the compositional complexity and the need to consider multiple kinetic factors simultaneously. To this end, we combine experiments with a first-principles-derived, multiscale computational framework that transcends thermodynamic descriptions to explicitly simulate the kinetic evolution of surface oxides of Ni-Cr alloys as a function of composition, temperature, pH, and applied voltage. In the absence of pitting, we identify three distinct voltage regimes, which are kinetically dominated by oxide growth, dissolution, and competitive dissolution and reprecipitation. Evolving compositional gradients and oxide thickness are revealed, including a transition between a metastable Ni-Cr mixed oxide and a thick, porous Ni-dominated oxide. Beyond elucidating the underlying physics, we highlight the need for competing kinetics in models to properly predict the transition from passivation to corrosion. Our results provide a key step towards co-design of alloy composition alongside environmental conditions for sustainable use across a variety of critical energy and infrastructure applications.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Solventless, Ambient-Pressure Production of Bio-Based Lactones Over Earth-Abundant, Mixed Metal Oxide Catalysts for Circular Polyesters

Transitioning to a circular plastics economy will require use of renewable feedstocks, energy-efficient processes, and closed-loop recyclable polymers, such as polyesters. A key challenge lies in sustainably sourcing monomers used to make recyclable polyesters. This work presents a catalytic platform utilizing earth-abundant Cu(x)Ca(1-x)O mixed metal oxides for the oxidative dehydrocyclization of bio-based diols to lactones, which are advantaged for energy-efficient ring-opening polymerization. Operating below 200 degrees C, at ambient pressure, and without solvent, the process uses air as the sole oxidant, achieving high yields of lactones across a broad substrate scope of C4-8 diols in the liquid phase. The oxidative dehydrocyclization reaction is thermodynamically downhill due to water formation and energy-efficient compared to incumbent, non-redox pathways utilized in fossil carbon-based industrial processes for lactone production. Mechanistic studies reveal facile redox cycling of Cu2+-O(Ca2+)-Cu2+ interfacial sites unique to the developed catalyst. Techno-economic analysis and life cycle assessment estimate 40% lower energy demand and 15% lower GHG intensity per mass of butyrolactone produced compared to the fossil carbon-based route. Liquid-phase oxidative dehydrocyclization offers a promising approach for scalable lactone production from renewable, bio-based diols to enable circular polyesters.

36 MATERIALS SCIENCE↗