SEARCH · Search NASA

Results for “nonequilibrium chemistry”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

NASA NTRS1995

Numerical Simulations of High Enthalpy Pulse Facilities

Axisymmetric flows within shock tubes and expansion tubes are simulated including the effects of finite rate chemistry and both laminar and turbulent boundary layers. The simulations demonstrate the usefulness of computational fluid dynamics for characterizing the flows in high enthalpy pulse facilities. The modeling and numerical requirements necessary to simulate these flows accurately are also discussed. Although there is a large body of analysis which explains and quantifies the boundary layer growth between the shock and the interface in a shock tube, there is a need for more detailed solutions. Phenomena such as thermochemical nonequilibrium. or turbulent transition behind the shock are excluded in the assumptions of Mirels' analysis. Additionally there is inadequate capability to predict the influence of the boundary layer on the expanded gas behind the interface. Quantifying the gas in this region is particularly important in expansion tubes because it is the location of the test gas. Unsteady simulations of the viscous flow in shock tubes are computationally expensive because they must follow features such as a shock wave over the length of the facility and simultaneously resolve the small length scales within the boundary layer. As a result, efficient numerical algorithms are required. The numerical approach of the present work is to solve the axisymmetric gas dynamic equations using an finite-volume formulation where the inviscid fluxes are computed with a upwind TVD scheme. Multiple species equations are included in the formulation so that finite-rate chemistry can be modeled. The simulations cluster grid points at the shock and interface and translate this clustered grid with these features to minimize numerical errors. The solutions are advanced at a CFL number of less than one based on the inviscid gas dynamics. To avoid limitations on the time step due to the viscous terms, these terms are treated implicitly. This requires a block tri-diagonal matrix inversion along each line of cells normal to the wall. The cost of this inversion is more than offset by the larger allowable time step. The source terms representing the finite-rate chemical kinetics are also treated implicitly. An algebraic turbulence model for compressible flow is used. The flow in a low pressure shock tube is computed and the results are compared with Mirels'analysis. The driven gas is nitrogen at 70 Pa, and the incident shock speed is approximately 2.9 km/sec so that there is little dissociation. The simulations include a laminar boundary layer and are run until the limiting flow regime is achieved. At this limit, the shock and interface travel at the same velocity because the amount of driven gas between these two features remains the same: the mass flow across the shock is equal to the mass of gas being entrained at the interface by the boundary layer. Simulations with several grids are presented to establish the grid independence of the solution, Good agreement is achieved between Mirels' correlations and the computations. This is expected since the flow conditions are chosen to be consistent with the assumptions used in Mirels' analysis. This comparison adds credibility to the numerical approach and highlights some of the differences between the theory and the detailed simulations. In addition, simulations of the HYPULSE expansion tube are presented for two operating conditions and the computations are compared to experimental data. The operating gas for both cases is nitrogen. One test condition is at a total enthalpy of 15.2 MJ/Kg and a relatively low pressure of 2 kPa. This case is characterized by a laminar boundary layer and significant chemical nonequilibrium. in the acceleration gas. The second test condition is at a total enthalpy of 10.2 MJ/Kg and a pressure of 38 kPa and is characterized by a turbulent boundary layer. The simulations compare well with experiment and reveal that the nonuniformity in pressure observed during the test time is related to variations in the boundary layer displacement thickness.

Wilson, Gregory J.↗

DOE OSTI2022

Fabricating Graphene Oxide/h-BN Metal Insulator Semiconductor Diodes by Nanosecond Laser Irradiation

To employ graphene’s rapid conduction in 2D devices, a heterostructure with a broad bandgap dielectric that is free of traps is required. Within this paradigm, h-BN is a good candidate because of its graphene-like structure and ultrawide bandgap. We show how to make such a heterostructure by irradiating alternating layers of a-C and a-BN film with a nanosecond excimer laser, melting and zone-refining constituent layers in the process. With Raman spectroscopy and ToF-SIMS analyses, we demonstrate this localized zone-refining into phase-pure h-BN and rGO films with distinct Raman vibrational modes and SIMS profile flattening after laser irradiation. Furthermore, in comparing laser-irradiated rGO-Si MS and rGO/h-BN/Si MIS diodes, the MIS diodes exhibit an increased turn-on voltage (4.4 V) and low leakage current. The MIS diode I-V characteristics reveal direct tunneling conduction under low bias and Fowler-Nordheim tunneling in the high-voltage regime, turning the MIS diode ON with improved rectification and current flow. This study sheds light on the nonequilibrium approaches to engineering h-BN and graphene heterostructures for ultrathin field effect transistor device development.

2D materials↗

NASA NTRS1990

Comparison of three-dimensional nonequilibrium PNS codes

A comparison study has been conducted using four recently developed parabolized Navier-Stokes (PNS) codes which have the capability of predicting finite-rate, chemically reacting flows over three-dimensional bodies. These are the (1) UPS code, (2) the STUFF code, (3) the TONIC code, and (4) the VRA-PNS code. All of the codes use the same seven-species, single-temperature air chemistry model, but otherwise they are unique, with different capabilities and characteristics. The differences include upwinding vs central differencing, strongly-coupled vs weakly-coupled chemistry, shock capturing vs shock fitting, finite volume vs finite difference, and full PNS vs thin-layer PNS equations. Three test cases were utilized to compare the codes. The comparisons presented indicate a good agreement among the codes tested.

Buelow, Philip E.↗

NASA NTRS1989

Computing viscous reactive flow over blunt winged vehicles

A hypersonic flow simulation methodology is developed for analyzing thermochemical nonequilibrium and viscous phenomena associated with space-flight vehicles at realistic conditions. The computer code consists of grid generation, Navier-Stokes and species models, and equation solvers for a complete configuration at arbitrary angles of attack. Provisions are implemented in the code to consider chemically frozen or equilibrium chemistry, inviscid model, conic configurations and simplifications of physical parameters for verifying solution consistency and accuracy. Implications of physical modeling and input parameters on a blunt-cone flow are investigated in details. Applications to USA Orbiter-like vehicles have emphasized the leeside canopy flow as well as the bow-wing shock interaction problem at Mach 22 and 40 deg angle of attack. Comparison of windward centerline pressures and heating rates with available wind-tunnel and flight data is in progress.

Li, C. P.↗

NASA NTRS1993

A coupled implicit method for chemical non-equilibrium flows at all speeds

The present time-accurate coupled-solution procedure addresses the chemical nonequilibrium Navier-Stokes equations over a wide Mach-number range uses, in conjunction with the strong conservation form of the governing equations, five unknown primitive variables. The numerical tests undertaken address steady convergent-divergent nozzle flows with air dissociation/recombination, dump combustor flows with n-pentane/air chemistry, and unsteady nonreacting cavity flows.

Shuen, Jian-Shun↗

NASA NTRS1992

Nonequilibrium synthesis of NbAl3 and Nb-Al-V alloys by laser cladding. II - Oxidation behavior

Isothermal oxidation behaviors of NbAl3 alloy synthesized by laser cladding were investigated at temperatures between 800 and 1400 C, and the effect of vanadium microalloying on the oxidation of the laser-clad alloy was examined. The oxidation kinetics of the two alloys were monitored using thermal gravimetric weight gain data, and the bulk and surface chemistries were analyzed using XRD and XPS, respectively. It was found that NbAl3 did not form an exclusive layer of protective Al2O3. The oxidation products at 800 C were found to be a mixture of Nb2O5 and Al2O3. At 1200 C, a mixture of NbAlO4, Nb2O5, and Al2O3 formed; and at 1400 C, a mixture of NbAlO4, Al2O3, NbO2, NbO(2.432), and Nb2O5 formed. The addition of V led to a dramatic increase of the oxidation rate, which may be related to the formation of (Nb, V)2O5 and VO2, which grows in preference to protective Al2O3.

Haasch, R. T.↗

DOE OSTI2020

Sulfonyl PIM–1: A diverse separation membrane with dilation resistance

A sulfone modified variant (soPIM–1) of the first polymer of intrinsic microporosity has been studied through molecular simulations to analyze its applicability for adsorption–based separations of diverse nonpolar, quadrupolar, and dipolar adsorbate species. Single component gas phase adsorption isotherms of each adsorbate are provided. Here, the adsorptive properties of soPIM–1 are reported from the application of two methods: (a) traditional grand canonical Monte Carlo (GCMC) simulations and (b) a combination of Monte Carlo and molecular dynamics (GCMC/MD) techniques, which accounts for sorption–induced polymer dynamics. The calculated isosteric heats of adsorption are compared to those from the parent PIM–1 structure and reveals increased CO2 affinity with relatively no change in hydrocarbon affinity. Moreover, to quantify soPIM–1's resistance to dilation the evolution of the microporous structure across the adsorption isotherm loading pressures has been evaluated. Relatedly, nonequilibrium MD simulations have been utilized to determine that soPIM–1 is approximately twice as stiff as PIM–1.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Trapping and imaging dynamic battery nanointerfaces via electrified cryo-EM

The electrified interface between a liquid and a solid underpins diverse phenomena, from ion-transfer during battery operation to action potentials enabling biological communication. However, conventional tools are blind to the nanoscale dynamics of this metastable interface. Here, we leverage electrified cryo–electron microscopy (eCryo-EM), a technique that rapidly freezes and kinetically traps these dynamic, nonequilibrium states during battery operation for nanoscale characterization. Collective snapshots of the electrified interface at controlled time intervals quantifies early-stage growth kinetics of the solid electrolyte interphase (SEI), a passivation film that governs electron and ion transport. Unexpectedly, the diffusivity of charged species of the two SEI films with differing chemistry and performance are estimated to be within 10% of the other, indicated by the slope of their diffusion-limited SEI growth regimes. Instead, the slope of the reaction-limited SEI growth regimes differs by a factor of 3, suggesting that lowered reactivity of the high-performance electrolyte is largely responsible for its high coulombic efficiency.

Science & Technology - Other Topics↗

DOE OSTI2020

Surfactant-Induced Interfacial Aggregation of Porphyrins for Structuring Color-Tunable Liquids

Locking nonequilibrium shapes of liquids into targeted architectures by interfacial jamming of nanoparticles is an emerging area in material science. 5,10,15,20-tetrakis(4-sulfonatophenyl) porphyrin (H 6 TPPS) shows three different aggregation states that present an absorption imaging platform to monitor the assembly and jamming of supramolecular polymer surfactants (SPSs) at the liquid/liquid interface. The interfacial interconversion of H 6 TPPS, specifically H 4 TPPS 2- dissolved in water, from J- to an H-aggregation was induced by strong electrostatic interactions with amine-terminated polystyrene dissolved in toluene at the water/toluene interface. We find that this resulted in color-tunable liquids due to interfacial jamming of the SPSs formed between H 4 TPPS 2- and amine-terminated polystyrene. However, the formed SPSs cannot lock in nonequilibrium shapes of liquids. In addition, a self-wrinkling behavior was observed when amphiphilic triblock copolymers of PS-block-poly(2-vinylpyridine)-block-poly(ethylene oxide) were used to interact with H 4 TPPS 2- . Subsequently, the SPSs formed can lock in nonequilibrium shapes of liquids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Polynitrogen High Energy Density Materials Synthesized by Nonequilibrium Plasma

Metal pentazolate compounds that contain 5-nitrogen singly charged anion rings are an interesting class of materials from a thermodynamics perspective. These compounds are believed to be equilibrium phases only at extremely high pressures, and thus, they constitute exotic states of matter at ambient pressure and temperature. Moreover, the energy released by these pentazolate compounds upon relaxation to the equilibrium state at ambient conditions is in the range 1 to 10 MJ per kg, rivaling that of combustion reactions. The synthesis of pentazolates is not an easy task, and there are only two known methods, each of which has significant challenges. Here, the development of new scalable synthesis routes could enable the production of these exotic materials in sufficient quantity to explore their properties more widely beyond the context of explosives and propellants. In this Perspective, it is argued that nonequilibrium plasma is a promising reaction medium for the synthesis of metal pentazolate compounds, for example from the corresponding metal azide and N 2 gas activated by the discharge at low background temperature.

Chemical reactions↗

DOE OSTI2023

Learning heterogeneous reaction kinetics from X-ray videos pixel by pixel

Reaction rates at spatially heterogeneous, unstable interfaces are notoriously difficult to quantify, yet are essential in engineering many chemical systems, such as batteries and electrocatalysts. Experimental characterizations of such materials by operando microscopy produce rich image datasets, but data-driven methods to learn physics from these images are still lacking because of the complex coupling of reaction kinetics, surface chemistry and phase separation. Here we show that heterogeneous reaction kinetics can be learned from in situ scanning transmission X-ray microscopy (STXM) images of carbon-coated lithium iron phosphate (LFP) nanoparticles. Combining a large dataset of STXM images with a thermodynamically consistent electrochemical phase-field model, partial differential equation (PDE)-constrained optimization and uncertainty quantification, we extract the free-energy landscape and reaction kinetics and verify their consistency with theoretical models. We also simultaneously learn the spatial heterogeneity of the reaction rate, which closely matches the carbon-coating thickness profiles obtained through Auger electron microscopy (AEM). Across 180,000 image pixels, the mean discrepancy with the learned model is remarkably small (<7%) and comparable with experimental noise. Our results open the possibility of learning nonequilibrium material properties beyond the reach of traditional experimental methods and offer a new non-destructive technique for characterizing and optimizing heterogeneous reactive surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Computing Nonequilibrium Responses with Score-Shifted Stochastic Differential Equations

Using equilibrium fluctuations to understand the response of a physical system to an externally imposed perturbation is the basis for linear response theory, which is widely used to interpret experiments and shed light on microscopic dynamics. For nonequilibrium systems, perturbations cannot be interpreted simply by monitoring fluctuations in a conjugate observable and general response results rely on path ensemble averaging. Furthermore, these techniques do not apply to perturbations that affect the diffusion tensor in a stochastic system. Here, we introduce an “effective” physical process that represents the diffusion perturbed dynamics and enables accurate calculations of responses to a change in the diffusion. Interestingly, the effective dynamics contain an additional drift involving the instantaneous “score” of the system, and we leverage score matching algorithms to carry out nonequilibrium response calculations on systems for which the exact stationary distribution is unknown.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2012

Modeling of Laser Vaporization and Plume Chemistry in a Boron Nitride Nanotube Production Rig

Flow in a pressurized, vapor condensation (PVC) boron nitride nanotube (BNNT) production rig is modeled. A laser provides a thermal energy source to the tip of a boron ber bundle in a high pressure nitrogen chamber causing a plume of boron-rich gas to rise. The buoyancy driven flow is modeled as a mixture of thermally perfect gases (B, B2, N, N2, BN) in either thermochemical equilibrium or chemical nonequilibrium assuming steady-state melt and vaporization from a 1 mm radius spot at the axis of an axisymmetric chamber. The simulation is intended to define the macroscopic thermochemical environment from which boron-rich species, including nanotubes, condense out of the plume. Simulations indicate a high temperature environment (T > 4400K) for elevated pressures within 1 mm of the surface sufficient to dissociate molecular nitrogen and form BN at the base of the plume. Modifications to Program LAURA, a finite-volume based solver for hypersonic flows including coupled radiation and ablation, are described to enable this simulation. Simulations indicate that high pressure synthesis conditions enable formation of BN vapor in the plume that may serve to enhance formation of exceptionally long nanotubes in the PVC process.

Gnoffo, Peter A.↗

NASA NTRS1993

A viscous shock-layer analysis of 2-D and axisymmetric flows

Results are obtained for cylindrical leading edges of proposed transatmospheric vehicles by employing a 2D viscous shock-layer code for nonequilibrium and perfect gas flows. The accuracy and efficiency of the planar code is verified through detailed comparisons with other predictions. It is found to be as accurate and robust as its axisymmetric counterpart. This study includes results for nose radii ranging from 0.01 to 2.0 ft and half-angles of 5 and 6 deg for both cylindrically-blunted wedges and spherically blunted cones (included for comparison). Some results are presented as a ratio of the noncatalytic to the corresponding fully catalytic heating value to illustrate the maximum potential for a heating reduction in dissociated nonequilibrium flows. This ratio and the individual heating rates are smaller for cylindrically-blunted wedges with small nose radii, relative to the spherically-blunted cones (at the same radius). Therefore, a larger potential exists for heating reduction in cylindrically-blunted as compared with the spherically-blunted surfaces for finite-rate chemistry.

Gupta, Roop N.↗

NASA NTRS2014

Study of Aerothermodynamic Modeling Issues Relevant to High-Speed Sample Return Vehicles

This paper examines the application of state-of-the-art coupled ablation and radiation simulations to highspeed sample return vehicles, such as those returning from Mars or an asteroid. A defining characteristic of these entries is that the surface recession rates and temperatures are driven by nonequilibrium convective and radiative heating through a boundary layer with significant surface blowing and ablation products. Measurements relevant to validating the simulation of these phenomena are reviewed and the Stardust entry is identified as providing the best relevant measurements. A coupled ablation and radiation flowfield analysis is presented that implements a finite-rate surface chemistry model. Comparisons between this finite-rate model and a equilibrium ablation model show that, while good agreement is seen for diffusion-limited oxidation cases, the finite-rate model predicts up to 50% lower char rates than the equilibrium model at sublimation conditions. Both the equilibrium and finite rate models predict significant negative mass flux at the surface due to sublimation of atomic carbon. A sensitivity analysis to flowfield and surface chemistry rates show that, for a sample return capsule at 10, 12, and 14 km/s, the sublimation rates for C and C3 provide the largest changes to the convective flux, radiative flux, and char rate. A parametric uncertainty analysis of the radiative heating due to radiation modeling parameters indicates uncertainties ranging from 27% at 10 km/s to 36% at 14 km/s. Applying the developed coupled analysis to the Stardust entry results in temperatures within 10% of those inferred from observations, and final recession values within 20% of measurements, which improves upon the 60% over-prediction at the stagnation point obtained through an uncoupled analysis. Emission from CN Violet is shown to be over-predicted by nearly and order-of-magnitude, which is consistent with the results of previous independent analyses. Finally, the coupled analysis is applied to a 14 km/s Earth entry representative of a Mars sample return. Although the radiative heating provides a larger fraction of the total heating, the influence of ablation and radiation on the flowfield are shown to be similar to Stardust.

Johnston, Christopher O.↗

NASA NTRS1992

Nonequilibrium synthesis of NbAl3 and Nb-Al-V alloys by laser cladding. I - Microstructure evolution

The evolution of the microstructure in NbAl3 synthesized by a laser cladding technique (a rapid solidification process, with cooling rates up to 10 exp 6 C/sec) is investigated, and the phases are identified using convergent beam electron diffraction. Two new metastable phases were identified and characterized in detail. The effect of adding V on the final microstructure was also investigated, and the various phase chemistries and the partitioning of different elements into different phases were studied.

Sircar, S.↗

DOE OSTI2022

Time-resolved measurements of HO 2 radical in a heated plasma flow reactor

Time-resolved, absolute HO 2 number density in diluted H 2 –O 2 –Ar, CH 4 –O 2 -Ar, and C 2 H 4 –O 2 –Ar mixtures excited by a repetitive ns pulse discharge in a heated plasma flow reactor is measured by Cavity Ringdown Spectroscopy (CRDS). The experimental results are obtained at $\textit{T}$ = 300-600 K and $\textit{P}$ = 130 Torr, both during the discharge pulse burst and in the afterglow. In this work, the HO 2 number density is inferred from the CRDS data using a spectral model exhibiting good agreement with previous measurements of absolute HO 2 absorption cross sections. In the room-temperature H 2 –O 2 mixture, as well as in CH 4 –O 2 and C 2 H 4 –O 2 mixtures over the entire temperature range studied, HO2 is generated only during the discharge burst and decays in the afterglow. However, in the H 2 –O 2 mixture at elevated temperatures, $\textit{T}$ = 400-600 K, HO 2 persists in the afterglow up to 10 ms after the discharge burst, comparable with the flow residence time in the reactor. Comparison with kinetic modeling shows that the sustained reactivity after the source of radicals is turned off is due to a chain propagation / hydrogen oxidation process, which dominates the radical recombination reactions. The kinetic modeling predictions are in good agreement with the relative HO 2 number density measured in all three mixtures, although the model underpredicts the absolute number densities in H 2 –O 2 at $\textit{T}$ = 400-600 K by up to a factor of two. Detection of the sustained low-temperature reactivity in H 2 –O 2 , initiated by the radical generation in the plasma, suggests that the plasma excitation may also affect kinetics of oxidation and reforming of fuels exhibiting low-temperature chemistry below hot ignition point.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Nonequilibrium Solvent Effects during Photodissociation in Liquids: Dynamical Energy Surfaces, Caging and Chemical Identity

In the gas phase, potential energy surfaces can be used to provide insight into the details of photochemical reaction dynamics. In solution, however, it is unclear what potential energy surfaces, if any, can be used to describe even simple chemical reactions such as the photodissociation of a diatomic solute. In this paper, we use mixed quantum/classical (MQC) molecular dynamics (MD) to study the photodissociation of $Na$ $^+_2$ in both liquid Ar and liquid tetrahydrofuran (THF). We examine both the gas-phase potential surfaces and potentials of mean force (PMF), which assume that the solvent remains at equilibrium with the solute throughout the photodissociation process and show that neither resemble a nonequilibrium dynamical energy surface that is generated by taking the time integral of work. For the photodissociation of $Na$ $^+_2$ in liquid Ar, the dynamical energy surface shows clear signatures of solvent caging, and the degree of caging is directly related to the mass of the solvent atoms. For $Na$ $^+_2$ in liquid THF, local specific interactions between the solute and solvent lead to changes in chemical identity that create a kinetic trap that effectively prevents the molecule from dissociating. Here, the results show that nonequilibrium effects play an important role even in simple solution-phase reactions, requiring the use of dynamical energy surface to understand such chemical events.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗