Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “equilibrium modeling”

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 1,603 records · Page 89
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Active Thermochemical Tables: the thermophysical and thermochemical properties of methyl, CH 3 , and methylene, CH 2 , corrected for nonrigid rotor and anharmonic oscillator effects

The thermophysical properties (isobaric heat capacity, entropy, enthalpy increment) of two prominent radicals, methyl, CH 3 , and methylene, CH 2 , were computed using the Nonrigid Rotor Anharmonic Oscillator (NRRAO) approach and compared to their RRHO counterparts, demonstrating significant differences between the results from the two approaches. Methylene presents a typical case in which the NRRAO thermophysical properties have significantly higher values than their RRHO counterparts at higher temperatures. In the case of methyl, the positive anharmonicity of the umbrella motion causes an opposite effect, and the NRRAO corrected thermophysical properties have lower values than their RRHO counterparts. The NRRAO corrected thermophysical properties, in turn, affect the resulting thermochemical properties. Here, two reactions important in combustion modelling were tested: the recombination of methyl radical with hydrogen atoms to form methane, and the recombination of two methyl radicals to form ethane. The related NRRAO equilibrium constants differ significantly from their RRHO analogs, and the consequences for chemical modelling are discussed. Also reported are the most current ATcT enthalpies of formation for CH n (n = 4-0) species and for C 2 H 6 , together with the tightly related sequential bond dissociation enthalpies along with the CHn series.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1994

A rapid-pressure correlation representation consistent with the Taylor-Proudman theorem materially-frame-indifferent in the 2D limit

A nonlinear representation for the rapid-pressure correlation appearing in the Reynolds stress equations, consistent with the Taylor-Proudman theorem, is presented. The representation insures that the modeled second-order equations are frame-invariant with respect to rotation when the flow is two-dimensional in planes perpendicular to the axis of rotation. The representation satisfies realizability in a new way: a special ansatz is used to obtain analytically, the values of coefficients valid away from the realizability limit: the model coefficients are functions of the state of the turbulence that are valid for all states of the mechanical turbulence attaining their constant limiting values only when the limit state is achieved. Utilization of all the mathematical constraints are not enough to specify all the coefficients in the model. The unspecified coefficients appear as free parameters which are used to insure that the representation is asymptotically consistent with the known equilibrium states of a homogeneous sheared turbulence. This is done by insuring that the modeled evolution equations have the same fixed points as those obtained from computer and laboratory experiments for the homogeneous shear. Results of computations of the homogeneous shear, with and without rotation, and with stabilizing and destabilizing curvature, are shown. Results are consistently better, in a wide class of flows which the model not been calibrated, than those obtained with other nonlinear models.

Ristorcelli, J. R.↗
NASA NTRS2002

Laser Induced Aluminum Surface Breakdown Model

Laser powered propulsion systems involve complex fluid dynamics, thermodynamics and radiative transfer processes. Based on an unstructured grid, pressure-based computational aerothermodynamics; platform, several sub-models describing such underlying physics as laser ray tracing and focusing, thermal non-equilibrium, plasma radiation and air spark ignition have been developed. This proposed work shall extend the numerical platform and existing sub-models to include the aluminum wall surface Inverse Bremsstrahlung (IB) effect from which surface ablation and free-electron generation can be initiated without relying on the air spark ignition sub-model. The following tasks will be performed to accomplish the research objectives.

Chen, Yen-Sen↗
DOE OSTI2020

Experimental Investigation of Uranium Volatility during Vapor Condensation

The predictive models that describe the fate and transport of radioactive materials in the atmosphere following a nuclear incident (explosion or reactor accident) assume that uranium-bearing particulates would attain chemical equilibrium during vapor condensation. Here, we show that kinetically driven processes in a system of rapidly decreasing temperature can result in substantial deviations from chemical equilibrium. This can cause uranium to condense out in oxidation states (e.g., UO 3 vs UO 2 ) that have different vapor pressures, significantly affecting uranium transport. To demonstrate this, we synthesized uranium oxide nanoparticles using a flow reactor under controlled conditions of temperature, pressure, and oxygen concentration. The atomized chemical reactants passing through an inductively coupled plasma cool from ~5000 to 1000 K within milliseconds and form nanoparticles inside a flow reactor. The ex situ analysis of particulates by transmission electron microscopy revealed 2–10 nm crystallites of fcc-UO 2 or α-UO 3 depending on the amount of oxygen in the system. α-UO 3 is the least thermodynamically preferred polymorph of UO 3 . The absence of stable uranium oxides with intermediate stoichiometries (e.g., U 3 O 8 ) and sensitivity of the uranium oxidation states to local redox conditions highlight the importance of in situ measurements at high temperatures. Therefore, we developed a laser-based diagnostic to detect uranium oxide particles as they are formed inside the flow reactor. Our in situ measurements allowed us to quantify the changes in the number densities of the uranium oxide nanoparticles (e.g., UO 3 ) as a function of oxygen gas concentration. Our results indicate that uranium can prefer to be in metastable crystal forms (i.e., α-UO 3 ) that have higher vapor pressures than the refractory form (i.e., UO 2 ) depending on the oxygen abundance in the surrounding environment. This demonstrates that the equilibrium processes may not dominate during rapid condensation processes, and thus kinetic models are required to fully describe uranium transport subsequent to nuclear incidents.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1992

Stress observations and driving force models for the South American plate

A comparison of stress observations to predictions of driving force models for the South American plate is presented. Conditions for dynamical equilibrium are used to estimate the magnitudes of forces as a function of age, with ridge push and slab pull as active forces, and plate drag and slab resistance as passive forces. The results are consistent with earlier published models, and it is found that the primary balance for an oceanic plate is between slab pull and the resisting drag forces. The model is generalized for the possible forces driving the South American plate, and it is compared with observed stress measurements and kinematics. The plate balance is found to be largely between ridge push and plate drag. The magnitude of a possible trench suction force is estimated using a model with an eddy between the Nazca slab and the overlying South American plate.

Stefanick, Michael↗
NASA NTRS1991

Static shape and instability of an acoustically levitated liquid drop

Dynamical aspects of a drop drastically flattened by acoustic radiation stress are considered. Its static equilibrium has been studied, starting with a dislike shape and modeling the sound field and the associated radiation stress according to this geometry. It is suggested that, at low viscosity, the ripples are capillary waves generated by the parametric instability excited by the membrane vibration, which is driven by the sound pressure. Atomization occurs whenever the membrane becomes so thin that the vibration is sufficiently intense. Buckling occurs when an existent equilibrium is unstable to a radial oscillation of the membrane because of the Bernoulli effect. The radiation stress at the rim of the flattened drop is also destabilizing and leads to horizontal expansion and subsequent breakup.

Lee, C. P.↗
DOE OSTI2022

Assessment of SCALE and MELCOR for a generic pebble bed fluoride high-temperature reactor

This paper presents the development and application of SCALE and MELCOR models for a fluoride salt-cooled high-temperature reactor (FHR) based on publicly available specifications. SCALE version 6.3beta15 was used to generate power distributions and decay heat curves, and MELCOR version 2.2.18019 was used to calculate the thermal hydraulic response of an assumed FHR primary system. Here, an approach was developed for determining the equilibrium state of the core using depletion of a core slice model and blending of fuel compositions at different burnups to provide three-dimensional fuel composition in the core. Results compared between a core with entirely fresh fuel and one with an equilibrium fuel composition revealed that the equilibrium core led to lower steady-state temperatures but slower cooldown during a loss of flow accident (LOFA). A sensitivity study was then conducted to explore the transient response of our FHR system to variations in thermal hydraulic parameters of the system using the equilibrium core. The inlet temperature, graphite thermal conductivity, and SCRAM time all provided significant control over peak fuel and coolant temperatures. Uncertainties in radionuclide decay data in SCALE were used to perform a decay heat sensitivity study, and we found that uncertainties in decay heat led to negligible impact on peak temperatures during the course of the transient. In all cases evaluated, the observed peak fuel temperatures remained approximately 700 K below anticipated failure limits for the LOFA.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Modeling the electrochemical behavior and interfacial junction profiles of bipolar membranes at solar flux relevant operating current densities

A 1-dimensional, multi-physics model that accounts for the migration and diffusion of solution species, electrostatics, and chemical reactions, in particular water dissociation (WD), at bipolar membrane (BPM) interfaces was developed to study the electrochemical behavior of bipolar membranes (BPMs) at solar flux relevant operating current densities (tens of mA cm -2 ). Significant partial current densities for WD were observed at BPM voltages much less than the equilibrium voltage, e.g., 59 mV × ΔpH from both experiments and modeling. The co-ion leakage across the BPM at pH differentials accounted for the early presence of the partial current density for WD. Two distinctive electric field dependent WD pathways, the un-catalyzed pathway and the catalyzed pathway, were quantitatively and parametrically studied to improve the turn-on potential of the BPM. The catalyzed pathway accounted for the majority of the partial current density for WD at low voltages, while the un-catalyzed pathway dominated the WD at relatively high voltages. Significant WD was observed only within the interfacial CL (<5 nm), in which a large electric field was present. To improve the electrochemical behavior and the turn-on potential of BPMs, the impacts of the pKa of the immobilized WD catalysts, the electric-field dependent rate constant, the thickness of the catalyst layer and fixed charge density in BPMs on the partial current densities for WD were studied systematically. In addition, the electrochemical behavior and concentration profiles of BPMs in a buffered electrolyte were studied and contrasted with those in an un-buffered electrolyte from both modeling and experiments.

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

Machine learning nonequilibrium electron forces for spin dynamics of itinerant magnets

Abstract We present a generalized potential theory for conservative as well as nonconservative forces for the Landau-Lifshitz magnetization dynamics. Importantly, this formulation makes possible an elegant generalization of the Behler-Parrinello machine learning (ML) approach, which is a cornerstone of ML-based quantum molecular dynamics methods, to the modeling of force fields in adiabatic spin dynamics of out-of-equilibrium itinerant magnetic systems. We demonstrate our approach by developing a deep-learning neural network that successfully learns the electron-mediated exchange fields in a driven s-d model computed from the nonequilibrium Green’s function method. We show that dynamical simulations with forces predicted from the neural network accurately reproduce the voltage-driven domain-wall propagation. Our work also lays the foundation for ML modeling of spin transfer torques and opens a avenue for ML-based multi-scale modeling of nonequilibrium dynamical phenomena in itinerant magnets and spintronics.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Modified Energy Span Analysis of Catalytic Parallel Pathways and Selectivity

Mechanistic modeling provides vital insights into catalytic reactions. To analyze complex reaction networks with parallel pathways, we leverage the graph theory approach of the Energy Span Model (ESM) to develop a modified energy span analysis (MESA). A new method of cycle plots is proposed to perform reaction pathways analysis visually. We demonstrate this method on two published models: one describing carbon monoxide oxidation and the other simulating ethylene conversion to propanal via hydroformylation or ethane via hydrogenation. Fundamental insights explain kinetic observables, such as a reactant’s negative reaction order. General principles are revealed, such as rate-determining surface species being outside the primary reaction flux cycle and pathway selectivity being a purely kinetic property when reaction conditions are not near equilibrium. Lastly, we demonstrate MESA’s consistency with published microkinetic modeling results, highlighting this technique’s extension of the ESM to heterogeneous catalysts using collision theory to describe adsorption steps and concentration effects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1981

The thermal emission spectrum of Io and a determination of the heat flux from its hot spots

Observations of thermal emission from Io in the near infrared made during an eclipse were combined with unpublished 8- to 13-micron intermediate band photometry and a 16- to 22-micron spectrum to specify Io's emission spectrum from 2.2 to 22 microns. Models were calculated having 'hot spots' at several different temperatures superposed on a surface, the major part of which is assumed to be at the solar equilibrium temperature. It was possible to fit the entire composite spectrum with this model. It is argued that the total emission from the hot spots can be equated to the nonsolar energy input into Io. The disk-averaged heat radiated by the hot spots is found to be 180 + or 60 microwatts/sq cm = 43 + or - 14 microcalories/sq cm-sec. A possible bimodal temperature distribution of the hot spots is discussed.

Sinton, W. M.↗
NASA NTRS1989

Properties of grains derived from IRAS observations of dust

The authors used the results of Infrared Astronomy Satellite (IRAS) observations of diffuse medium dust to develop a theoretical model of the infrared properties of grains. Recent models based entirely on traditional observations of extinction and polarization include only particles whose equilibrium temperatures do not exceed 20 K in the diffuse interstellar medium. These classical grains, for which the authors have adopted the multipopulation model developed by Hong and Greenberg (1980), can explain only the emission in the IRAS 100 micron band. The measurements at shorter wavelengths (12, 25 and 60 microns) require two new particle populations. Vibrational fluorescence from aromatic molecules provides the most likely explanation for the emission observed at 12 microns, with polycyclic aeromatic hydrocarbons (PAHs) containing about 10 percent of cosmic carbon. A simplified model of the emission process shows that PAH molecules can also explain most of the emission measured by IRAS at 25 microns. The authors identified the warm particles responsible for the excess 60 microns emission with small (a approx. equals 0.01 microns) iron grains. A compilation of the available data on the optical properties of iron indicates that the diffuse medium temperature of small iron particles should be close to 50 K and implies that a large, possibly dominant, fraction of cosmic iron must be locked up in metallic particles in order to match the observed 60 microns intensities. The model matches the infrared fluxes typically observed by IRAS in the diffuse medium and can also reproduce the infrared surface brightness distribution in individual clouds. In particular, the combination of iron and classical cool grains can explain the surprising observations of the 60/100 microns flux ratio in clouds, which is either constant or increases slightly towards higher opacities. The presence of metallic grains has significant implications for the physics of the interstellar medium, including catalytic H2 formation, for which iron grains could be the main site; differences in depletion patterns between iron and other refractory elements (Mg, Si); and superparamagnetic behavior of large grains with embedded iron clusters giving rise to the observed high degree of alignment by the galactic magnetic field.

Wesselius, P. R.↗
DOE OSTI2023

Low-Temperature Plasma Assisted Kinetics Study of Ethanol

Non-equilibrium plasma assisted kinetics of ethanol at 0.5 atm over a temperature range of 523-1203 K were studied at stoichiometric conditions by using a plasma flow reactor (PFR). The effects of plasma chemistry on fuel and its influence on neutral chemistry is decoupled with both experimental and modelling efforts. Enhanced reactivity was observed experimentally in both plasma-assisted oxidation and pyrolysis cases with fuel consumption beginning at temperatures as low as 523 K. Absolute fuel consumption was observed at temperatures nearly 200 K earlier than pure thermal reactions highlighting the effectiveness of non-equilibrium plasma to promote reactivity. Furthermore, modelling results illustrated the effect of electron impact reactions towards accelerating fuel consumption and promoting low temperature chemistry by generating excited state N2 and reactive O and OH radicals.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

A dissipation-induced superradiant transition in a strontium cavity-QED system

Driven-dissipative many-body systems are ubiquitous in nature and a fundamental resource for quantum technologies. However, they are also complex and hard to model because they cannot be described by the standard tools in equilibrium statistical mechanics. Probing nonequilibrium critical phenomena in pristine setups can illuminate fresh perspectives on these systems. Here, we use an ensemble of cold 88 Sr atoms coupled to a driven high-finesse cavity to study the cooperative resonance fluorescence (CRF) model, a classic driven-dissipative model describing coherently driven dipoles superradiantly emitting light. We observe its nonequilibrium phase diagram characterized by a second-order phase transition. Below a critical drive strength, the atoms quickly reach the so-called superradiant steady state featuring a macroscopic dipole moment; above the critical point, the atoms undergo persistent Rabi-like oscillations. At longer times, spontaneous emission transforms the second-order transition into a discontinuous first-order transition. Our observations pave the way for harnessing robust entangled states and exploring boundary time crystals in driven-dissipative systems.

Science & Technology - Other Topics↗
DOE OSTI2022

On the effect of mixing-driven vaporization in a homogeneous relaxation modeling framework

The homogeneous relaxation model (HRM) is one of the most widely used models to describe the liquid–gas phase transition in multiphase flows due to the occurrence of cavitation. However, in its original formulation, the HRM does not account for the presence of ambient gas species, which generally limits its applicability to the injector's internal flow where ambient gases are negligible. In this work, a mixing-driven vaporization (MDV) model was developed to extend the capability of the HRM in handling the mixing effect in the regions external to the nozzle, where vapor–liquid equilibrium for multi-species mixtures of fuel and ambient gas is considered. Herein, to assess the model performance, simulations of the Engine Combustion Network's Spray G injector were performed with the HRM and the MDV model under both flash-boiling and evaporating conditions. It was found that the MDV model led to a better match against x-ray measurements of fuel density in the near-nozzle region. In contrast to the HRM, the MDV model was able to reproduce the vaporization process in the mixing zone at the edge of the fuel jet, which aligns with the expected physics. This resulted in substantial differences in the prediction of other flow characteristics such as mixture temperature and pressure. Furthermore, this work demonstrates that evaporation timescales have a considerable effect on the MDV model's predictions, as shown by a parametric study in which a time factor was introduced to mimic the effect of different timescales due to different phase change mechanisms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1994

Active chlorine and nitric oxide formation from chemical rocket plume afterburning

Chlorine and oxides of nitrogen (NO(x)) released into the atmosphere contribute to acid rain (ground level or low-altitude sources) and ozone depletion from the stratosphere (high-altitude sources). Rocket engines have the potential for forming or activating these pollutants in the rocket plume. For instance, H2/O2 rockets can produce thermal NO(x) in their plumes. Emphasis, in the past, has been placed on determining the impact of chlorine release on the stratosphere. To date, very little, if any, information is available to understand what contribution NO(x) emissions from ground-based engine testing and actual rocket launches have on the atmosphere. The goal of this work is to estimate the afterburning emissions from chemical rocket plumes and determine their local stratospheric impact. Our study focuses on the space shuttle rocket motors, which include both the solid rocket boosters (SRB's) and the liquid propellant main engines (SSME's). Rocket plume afterburning is modeled employing a one-dimensional model incorporating two chemical kinetic systems: chemical and thermal equilibria with overlayed nitric oxide chemical kinetics (semi equilibrium) and full finite-rate chemical kinetics. Additionally, the local atmospheric impact immediately following a launch is modeled as the emissions diffuse and chemically react in the stratosphere.

Leone, D. M.↗
NASA NTRS1987

Investigation of turbulent flow in highly curved ducts with application to turbomachinery components

Numerical predictions were performed using a semi-elliptic calculation procedure for the case of turbulent flow in passage through a 90 deg bend of square cross section. Two versions of the isotropic turbulent viscosity two equation k-epsilon model were used. The first, the wall function model (WFM), employs the logarithmic law-of-the-wall relation and the notion of equilibrium flow to set all the necessary boundary conditions at the first grid point adjacent to a solid wall. The second model, the Van Driest model (VDM), employs Prandtl's original mixing length formulation, in conjunction with Van Driest's semi-empirical relation for the mixing length, to calculate the turbulent viscosity in the near wall regions of the flow. In this case, boundary conditions for k and epsilon, required to calculate these quantities in the core of the flow, are obtained by matching the mixing length and Reynolds number model formulations in an overlapping region of the flow near the walls. In both cases the results obtained show an improvement over earlier calculations using an elliptic numerical procedure. This is attributed to the finer grids possible in the present work. Of the two models, the VDM formulation shows better overall conformity with the mean flow measurements. Neither model reproduces well the details of the stress distribution as a result of the implied isotropic turbulent viscosity.

Humphrey, J. A. C.↗
DOE OSTI2026

Non-Equilibrium Actinide Radiation Chemistry and the Nuclear Fuel Cycle

Actinides are inherently unstable elements that frequently coexist with other radioisotopes, generating intense ionizing radiation fields that drive the formation of non equilibrium oxidation states. These transient species exert a profound mechanistic influence on the radiation response of actinide containing systems due to their unique redox chemistry. Despite their importance, they remain poorly understood, yet such insight is essential for advancing actinide science and accurately predicting radiation driven behavior. Actinide separations—critical for nuclear energy technologies, strategic deterrence, space exploration, and nuclear medicine—depend on precise control of actinide oxidation states to recover targeted elements from complex matrices such as used nuclear fuel. However, during these processes, actinides, their coordination complexes, and the separation media are all exposed to intense, multicomponent (alpha, beta, gamma, etc.) radiation fields that can alter process efficiency, selectivity, and chemical stability. Understanding, controlling, and mitigating radiation induced reactions is therefore key to innovating and optimizing next generation separation technologies. This seminar will provide an overview of the nuclear fuel cycle and non equilibrium actinide radiation chemistry in the context of recovering actinides from used nuclear fuel, with a particular emphasis on direct dissolution–based reprocessing strategies. We will explore time resolved electron pulse radiolysis and alpha and gamma dose accumulation studies, integrated with multiscale computational modeling, to elucidate the molecular level roles of radiation driven, non equilibrium actinide species in process performance and in the radiolytic stability of organic ligands used for actinide recovery. These insights offer new pathways for designing advanced separation methods and next generation solvent systems, with broad implications for the future of the nuclear fuel cycle.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗