Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic Modelling”

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,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Theoretical and experimental study on the O( 3 P) + 2,5-dimethylfuran reaction in the gas phase

In this work we report a joint experimental and computational study on the 2,5-dimethylfuran oxidation reaction in the gas phase initiated by atomic oxygen O( 3 P). The experiments have been performed by using vacuum-ultraviolet synchrotron radiation at the Advanced Light Source (ALS) of the Lawrence Berkeley National Laboratory (LBNL), at a temperature of 550 K and a pressure of 8 Torr. The experimental data were supported by quantum-chemical calculations along with a kinetic model, also taking into account the possible involvement of different magnetic states, performed in the framework of the RRKM theory. Propyne, acetaldehyde, methylglyoxal, dimethylglyoxal, 3-penten-2-one, 2,5-dimethylfuran-3(2 H )-one, and 1,2-diacetyl ethylene have been identified as the main primary products arising under the conditions of the experiment. Our computational model suggests that these species can be formed at the concentration and branching ratio experimentally observed only in the presence of a non-negligible fraction of non-thermalized intermediates.

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

Identification of the acetaldehyde oxide Criegee intermediate reaction network in the ozone-assisted low-temperature oxidation of trans -2-butene

Uni- and bi-molecular reactions involving Criegee intermediates (CIs) have been the focus of many studies due to the role these molecules play in atmospheric chemistry. The reactivity of CIs is known to strongly depend on their structure. The reaction network of the second simplest CI, acetaldehyde oxide (CH 3 CHOO), is investigated in this work in an atmospheric pressure jet-stirred reactor (JSR) during the ozonolysis of trans-2-butene to explore the kinetic pathways relevant to atmospheric chemistry and low-temperature combustion. The mole fraction profiles of reactants, intermediates, and final products are determined by means of molecular-beam mass spectrometry in conjunction with single-photon ionization employing tunable synchrotron-generated vacuum ultraviolet radiation. A network of CI reactions is identified in the temperature region below 600 K, characterized by CI addition to trans-2-butene, water, formaldehyde, formic acid, and methanol. No sequential additions of the CH 3 CHOO CI are observed, in contrast with the reactivity of the simplest CI (H 2 COO) and the earlier observation of an extensive reaction network with up to four H 2 COO sequential additions (Phys. Chem. Chem. Phys., 2019, 21, 7341–7357). Experimental photoionization efficiency scans recorded at 300 K and 425 K and ab initio threshold energy calculations lead to the identification and quantification of previously elusive intermediates, such as ketohydroperoxide and hydroperoxide species. Specifically, the C 4 H 8 + O 3 adduct is identified as a ketohydroperoxide (KHP, 3-hydroperoxybutan-2-one, CH 3 C([double bond, length as m-dash]O)CH(CH 3 )OOH), while hydroxyacetaldehyde (glycolaldehyde, HCOCH 2 OH) formation is attributed to unimolecular isomerization of the CIs. Other hydroperoxide species such as methyl hydroperoxide (CH 3 OOH), ethyl hydroperoxide (C 2 H 5 OOH), butyl hydroperoxide (OOH), hydroperoxyl acetaldehyde (HOOCH 2 CHO), hydroxyethyl hydroperoxide (CH 3 CH(OH)OOH), but-1-enyl-3-hydroperoxide, and 4-hydroxy-3-methylpentan-2-one (HOCH(CH 3 )CH(CH 3 )C([double bond, length as m-dash]O)CH 3 ) are also identified. Detection of additional oxygenated species such as methanol, ethanol, ketene, and aldehydes suggests multiple active oxidation routes. Finally, these results provide additional evidence that CIs are key intermediates of the ozone-unsaturated hydrocarbon reactions providing critical inputs for improved kinetics models.

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

Unconventional gas-phase preparation of the prototype polycyclic aromatic hydrocarbon naphthalene (C 10 H 8 ) via the reaction of benzyl (C 7 H 7 ) and propargyl (C 3 H 3 ) radicals coupled with hydrogen-atom assisted isomerization

Polycyclic aromatic hydrocarbons (PAHs) are ubiquitous in the interstellar medium and in meteorites such as Murchison and Allende and signify the missing link between resonantly stabilized free radicals and carbonaceous nanoparticles (soot particles, interstellar grains). However, the predicted lifetime of interstellar PAHs of some 108 years imply that PAHs should not exist in extraterrestrial environments suggesting that key mechanisms of their formation are elusive. Exploiting a microchemical reactor and coupling these data with computational fluid dynamics (CFD) simulations and kinetic modeling, we reveal through an isomer selective product detection that the reaction of the resonantly stabilized benzyl (C 7 H 7 ) and the propargyl (C 3 H 3 ) synthesizes the simplest representative of PAHs – the 10π Hückel aromatic naphthalene (C 10 H 8 ) molecule – via the novel Propargyl Addition–BenzAnnulation (PABA) mechanism. The gas-phase preparation of naphthalene affords a versatile concept of the reaction of combustion and astronomically abundant propargyl radicals with aromatic radicals carrying the radical center at the methylene moiety (aromatic-CH 2 ) as a previously passed over source of aromatics in high temperature environments thus bringing us closer to an understanding of the aromatic universe we live in.

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

Iodide oxidation by ozone at the surface of aqueous microdroplets

The kinetics of iodide oxidation by ozone are studied in single microdroplets to explore chemical reactivity at the air–water interface. Molecular simulations and kinetic modeling are employed to understand the underlying multiphase mechanism.

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

Tracking the reaction networks of acetaldehyde oxide and glyoxal oxide Criegee intermediates in the ozone-assisted oxidation reaction of crotonaldehyde

The reaction of unsaturated compounds with ozone (O 3 ) is recognized to lead to the formation of Criegee intermediates (CIs), which play a key role in controlling the atmospheric budget of hydroxyl radicals and secondary organic aerosols. The reaction network of two CIs with different functionality, i.e. acetaldehyde oxide (CH 3 CHOO) and glyoxal oxide (CHOCHOO) formed in the ozone-assisted oxidation reaction of crotanaldehyde (CA), is investigated over a temperature range between 390 K and 840 K in an atmospheric pressure jet-stirred reactor (JSR) at a residence time of 1.3 s, stoichiometry of 0.5 with a mixture of 1% crotonaldehyde, 10% O 2 , at an fixed ozone concentration of 1000 ppm and 89% Ar dilution. Molecular-beam mass spectrometry in conjunction with single photon tunable synchrotron vacuum-ultraviolet (VUV) radiation is used to identify elusive intermediates by means of experimental photoionization energy scans and ab initio threshold energy calculations for isomer identification. Addition of ozone (1000 ppm) is observed to trigger the oxidation of CA already at 390 K, which is below the temperature where the oxidation reaction of CA was observed in the absence of ozone. The observed CA + O 3 product, C 4 H 6 O 4 , is found to be linked to a ketohydroperoxide (2-hydroperoxy-3-oxobutanal) resulting from the isomerization of the primary ozonide. Products corresponding to the CIs uni- and bi-molecular reactions were observed and identified. A network of CI reactions is identified in the temperature region below 600 K, characterized by CIs bimolecular reactions with species like aldehydes, i.e., formaldehyde, acetaldehyde, and crotonaldehyde and alkenes, i.e., ethene and propene. The region below 600 K is also characterized by the formation of important amounts of typical low-temperature oxidation products, such as hydrogen peroxide (H 2 O 2 ), methyl hydroperoxide (CH 3 OOH), and ethyl hydroperoxide (C 2 H 5 OOH). Detection of additional oxygenated species such as alcohols, ketene, and aldehydes are indicative of multiple active oxidation routes. This study provides important information about the initial step involved in the CIs assisted oligomerization reactions in complex reactive environments where CIs with different functionalities are reacting simultaneously. It provides new mechanistic insights into ozone-assisted oxidation reactions of unsaturated aldehydes, which is critical for the development of improved atmospheric and combustion kinetics models.

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

Radiating electron source generation in ultraintense laser-foil interactions

A radiating electron source is shown to be created by a laser pulse (with intensity of 10{sup 23 }W/cm{sup 2} and duration equal to 30 fs) interacting with a near-critical density plasma. It is shown that the back radiation reaction resulting from high energy synchrotron radiation tends to counteract the action of the ponderomotive force. This enhances the collective dynamics of the radiating electrons in the highest field areas, resulting in the production of a compact radiation source (containing 80% of the synchrotron radiation emission), with an energy on the order of tens of MeV over the laser pulse duration. These phenomena are investigated using a QED-particle-in-cell code, and compared with a kinetic model accounting for the radiation reaction force in the electron distribution function. The results shed new light on electron-photon sources at ultra-high laser intensities and could be tested on future laser facilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Plasma filamentation and shock wave enhancement in microwave rockets by combining low-frequency microwaves with external magnetic field

A filamentary plasma is reproduced based on a fully kinetic model of electron and ion transports coupled with electromagnetic wave propagation. The discharge plasma transits from discrete to diffusive patterns at a 110-GHz breakdown, with decrease in the ambient pressure, because of the rapid electron diffusion that occurs during an increase in the propagation speed of the ionization front. A discrete plasma is obtained at low pressures when a low-frequency microwave is irradiated because the ionization process becomes more dominant than the electron diffusion, when the electrons are effectively heated by the low-frequency microwave. The propagation speed of the plasma increases with decrease in the incident microwave frequency because of the higher ionization frequency and faster plasma diffusion resulting from the increase in the energy-absorption rate. An external magnetic field is applied to the breakdown volume, which induces plasma filamentation at lower pressures because the electron diffusion is suppressed by the magnetic field. The thrust performance of a microwave rocket is improved by the magnetic fields corresponding to the electron cyclotron resonance (ECR) and its higher-harmonic heating, because slower propagation of the ionization front and larger energy-absorption rates are obtained at lower pressures. It would be advantageous if the fundamental mode of ECR heating is coupled with a lower frequency microwave instead of combining the higher-harmonic ECR heating with the higher frequency microwave. This can improve the thrust performance with smaller magnetic fields even if the propagation speed increases because of the decrease in the incident microwave frequency.

42 ENGINEERING↗
DOE OSTI2017

Self-consistent perturbed equilibrium with neoclassical toroidal torque in tokamaks

Toroidal torque is one of the most important consequences of non-axisymmetric fields in tokamaks. The well-known neoclassical toroidal viscosity (NTV) is due to the second-order toroidal force from anisotropic pressure tensor in the presence of these asymmetries. This work shows that the first-order toroidal force originating from the same anisotropic pressure tensor, despite having no flux surface average, can significantly modify the local perturbed force balance and thus must be included in perturbed equilibrium self-consistent with NTV. The force operator with an anisotropic pressure tensor is not self-adjoint when the NTV torque is finite and thus is solved directly for each component. This approach yields a modified, non-self-adjoint Euler-Lagrange equation that can be solved using a variety of common drift-kinetic models in generalized tokamak geometry. The resulting energy and torque integral provides a unique way to construct a torque response matrix, which contains all the information of self-consistent NTV torque profiles obtainable by applying non-axisymmetric fields to the plasma. This torque response matrix can then be used to systematically optimize non-axisymmetric field distributions for desired NTV profiles. Published by AIP Publishing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Validation of the kinetic-turbulent-neoclassical theory for edge intrinsic rotation in DIII-D

In a recent kinetic model of edge main-ion (deuterium) toroidal velocity, intrinsic rotation results from neoclassical orbits in an inhomogeneous turbulent field [T. Stoltzfus-Dueck, Phys. Rev. Lett. 108, 065002 (2012)]. This model predicts a value for the toroidal velocity that is co-current for a typical inboard X-point plasma at the core-edge boundary (ρ ~ 0.9). Using this model, the velocity prediction is tested on the DIII-D tokamak for a database of L-mode and H-mode plasmas with nominally low neutral beam torque, including both signs of plasma current. Values for the flux-surface-averaged main-ion rotation velocity in the database are obtained from the impurity carbon rotation by analytically calculating the main-ion—impurity neoclassical offset. Here, the deuterium rotation obtained in this manner has been validated by direct main-ion measurements for a limited number of cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

High-level theoretical characterization of the vinoxy radical ( • CH 2 CHO) + O 2 reaction

Numerous processes in atmospheric and combustion chemistry produce the vinoxy radical ( • CH 2 CHO). To understand the fate of this radical and to provide reliable energies needed for kinetic modeling of such processes, we have examined its reaction with O 2 using highly reliable theoretical methods. Utilizing the focal point approach, the energetics of this reaction and subsequent reactions were obtained using coupled-cluster theory with single, double, and perturbative triple excitations [CCSD(T)] extrapolated to the complete basis set limit. These extrapolated energies were appended with several corrections including a treatment of full triples and connected quadruple excitations, i.e., CCSDT(Q). In addition, this study models the initial vinoxy radical + O 2 reaction for the first time with multireference methods. We predict a barrier for this reaction of approximately 0.4 kcal mol –1 . This result agrees with experimental findings but is in disagreement with previous theoretical studies. The vinoxy radical + O 2 reaction produces a 2-oxoethylperoxy radical which can undergo a number of unimolecular reactions. Abstraction of a β-hydrogen (a 1,4-hydrogen shift) and dissociation back to reactants are predicted to be competitive to each other due to their similar barriers of 21.2 and 22.3 kcal mol –1 , respectively. The minimum-energy β-hydrogen abstraction pathway produces a hydroperoxy radical (QOOH) that eventually decomposes to formaldehyde, CO, and • OH. Two other unimolecular reactions of the peroxy radical are α-hydrogen abstraction (38.7 kcal mol –1 barrier) and HO 2 • elimination (43.5 kcal mol –1 barrier). Furthermore, these pathways lead to glyoxal + • OH and ketene + HO 2 • formation, respectively, but they are expected to be uncompetitive due to their high barriers.

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

Role of electron inertia and electron/ion finite Larmor radius effects in low-beta, magneto-Rayleigh-Taylor instability

The magneto-Rayleigh-Taylor (MRT) instability has been explored in great detail in previous work using magnetohydrodynamic and kinetic models for low-beta plasmas. The work presented here extends previous studies of this instability to regimes where finite-Larmor-Radius (FLR) effects may be important. Comparisons of the MRT instability are made using a 5-moment and a 10-moment two-fluid model, the two fluids being ions and electrons. The 5-moment model includes Hall stabilization, whereas the 10-moment model includes Hall and FLR stabilization. Results are presented for these two models using different electron mass to understand the role of electron iner- tia in the late-time nonlinear evolution of the MRT instability. For the 5-moment model, the late- time nonlinear MRT evolution does not significantly depend on the electron inertia. Yet, when FLR stabilization is important, the 10-moment results show that a lower ion-to-electron mass ratio (i.e., larger electron inertia) under-predicts the energy in high-wavenumber modes due to larger FLR stabilization.

nite Larmor radius↗
DOE OSTI2019

Collisional enhancement of energetic particle Alfvénic resonance width in tokamaks

The phase-space structure of resonances between fast ions and an Alfvénic mode and the associated modification of density profiles in tokamaks are studied as a function of particle collisions. Guiding-center simulations in a realistic tokamak equilibrium are employed to address the resonance broadening parametric dependencies with respect to changes in the pitch-angle scattering rate. The rate of collisional replenishment, along with resonance strength, given by the combination of eigenmode and resonance structures and equilibrium parameters, determines saturation amplitudes for a given damping rate. As seen from the distribution function flattening, collisions have an effect of broadening the resonances, while the absolute value of δf decreases with increasing collisionality. It is observed that the collisional broadening can be comparable to the collisionless resonance width due to the mode amplitude alone. The resonance broadening coefficients are compared with the existing theory based on analytically expected saturation levels, showing fair agreement. Furthermore the results can be useful in assisting reduced kinetic models, such as quasilinear models, when prescribing the effective resonance phase-space width, i.e., the mode-particle interaction platform, due to collisional or turbulent processes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Energetics and mechanisms for the acetonyl radical + O 2 reaction: An important system for atmospheric and combustion chemistry

The acetonyl radical (•CH 2 COCH 3 ) is relevant to atmospheric and combustion chemistry due to its prevalence in many important reaction mechanisms. One such reaction mechanism is the decomposition of Criegee intermediates in the atmosphere that can produce acetonyl radical and OH. In order to understand the fate of the acetonyl radical in these environments and to create more accurate kinetics models, we have examined the reaction system of the acetonyl radical with O 2 using highly reliable theoretical methods. Structures were optimized using coupled cluster theory with singles, doubles, and perturbative triples [CCSD(T)] with an atomic natural orbital (ANO0) basis set. Energetics were computed to chemical accuracy using the focal point approach involving perturbative treatment of quadruple excitations [CCSDT(Q)] and basis sets as large as cc-pV5Z. The addition of O 2 to the acetonyl radical produces the acetonylperoxy radical, and multireference computations on this reaction suggest it to be barrierless. Additionally, no submerged pathways were found for the unimolecular isomerization of the acetonylperoxy radical. Besides dissociation to reactants, the lowest energy pathway available for the acetonylperoxy radical is a 1-5 H shift from the methyl group to the peroxy group through a transition state that is 3.3 kcal mol -1 higher in energy than acetonyl radical + O 2 . The ultimate products from this pathway are the enol tautomer of the acetonyl radical along with O 2 . Multiple pathways that lead to OH formation are considered; however, all of these pathways are predicted to be energetically inaccessible, except at high temperatures.

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

Mean force kinetic theory applied to self-diffusion in supercritical Lennard-Jones fluids

A new kinetic-theory-based calculation of the self-diffusion coefficient for dense supercritical Lennard-Jones fluids is presented. The mean force kinetic theory, which was recently developed for transport in dense plasmas, is applied for the calculation of diffusion in dense neutral fluids. The calculation only requires the pair distribution function, a quantity that is readily calculable from equilibrium statistical mechanics for many systems, including the Lennard-Jones fluid. The self-diffusion coefficients are compared with calculations from molecular dynamics simulations, and good agreement at high density is demonstrated, even in the vicinity of the solid–fluid coexistence line. A comparison of different kinetic models with molecular dynamics simulations demonstrates that the transport coefficients have important contributions due to particle interaction via a potential of mean force and local correlations, which increase the collision rate. The new calculations compare well to those from free-volume theory and overcome a limitation of this theory that prevents its use in systems that interact via long range monotonic potentials. It is expected that this approach will also apply to other systems, including neutral-plasma and neutral-electrolyte mixtures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Study of vacuum ultraviolet emission in helium and helium/nitrogen mixtures

Helium is frequently used as a working medium for the generation of plasmas and is capable of energetic photon emissions. These energetic photon emissions are often attributed to the formation of helium excimer and subsequent photon emission. When the plasma device is exposed to another gas, such as nitrogen, this energetic photon emission can cause photoionization and further ionization wave penetration into the additional gas. Often ignored are the helium resonance emissions that are assumed to be radiation trapped and therefore not pertinent to photoionization. Here, experimental evidence for the presence of helium atomic emission in a pulsed discharge at ten’s of Torr is shown. Simulations of a discharge in similar conditions agree with the experimental measurements. In this context, the role of atomic and molecular helium light emission on photoionization of molecular nitrogen in an ionization wave is studied using a kinetic modeling approach that accounts for radiation dynamics in a developing low-temperature plasma. Three different mixtures of helium at a total pressure of 250 Torr are studied in simulation. Photoionization of the nitrogen molecule by vacuum ultraviolet helium emission is used as the only seed source ahead of the ionization front. It is found that even though radiation trapped, the atomic helium emission lines are the significant source of photoionization of nitrogen. The significant effect of radiation trapped photon emission on ionization wave dynamics demonstrates the need to consider these radiation dynamics in plasma reactors where self-absorbed radiation is ignored.

42 ENGINEERING↗
DOE OSTI2023

Ionization disequilibrium in K- and L-shell ions

Time-gated Sc K-shell and Ge L-shell spectra are presented from a range of characterized thermodynamic states spanning ion densities of 1019–1020cm−3 and plasma temperatures around 2000 eV. For the higher densities studied and temperatures from 1000 to 3000 eV, the Sc and Ge x-ray emission spectra are consistent with steady-state calculations from the modern atomic kinetics model SCRAM. At the lower ion densities achieved through plasma expansion, however, the model calculations require a higher plasma temperature to reproduce the observed Ge spectrum. We attribute this to ionization disequilibrium of the Sc because the ionization time scales exceed the hydrodynamic timescale when the inferred temperatures diverge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Effect of neutral interactions on parallel transport and blob dynamics in gyrokinetic scrape-off layer simulations

The effect of neutral interactions on scrape-off layer (SOL) turbulence is investigated in a continuum gyrokinetic code that has been coupled to a continuum kinetic model of neutral transport. This extends the work of a previous paper, which compared two NSTX SOL simulations in simple helical geometry, one with neutrals and one without. The former included electron-impact ionization, charge exchange, and wall recycling. Here, the case with neutrals is compared to a gyrokinetic-only simulation that includes an effective ionization source to separate the effect of sourcing from charge exchange collisions. It is observed that sourcing accounts for many features of the simulated SOL with neutrals, including density and temperature magnitudes and reduced normalized density fluctuations, but differences persist. In particular, a flatter density profile results due to changes in parallel transport when neutral collisions are included, illustrating the importance of neutral drag on global plasma properties. An analysis of coherent turbulent structures, or blobs, in these simulations demonstrates the case with neutrals has slower and larger blobs. Here, a series of seeded blob simulations corroborates the blob velocity observation. In general, the blob motion does not contribute significantly to radial transport in these simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Loss cone effects and monotonic sheath conditions of a partially magnetized plasma sheath

In this Letter, we propose the conditions for monotonic plasma sheaths adjacent to a floating wall in the presence of an applied, oblique magnetic field. The electron velocity distribution function (VDF) at the sheath edge obtained from a kinetic model exhibits a loss cone shaped truncation. Using an approximation of the truncated VDF, we derive an analytical framework of the sheath edge condition (i.e., Bohm condition), namely, the relation of ion injection velocity and sheath potential drop as a function of the magnetic field angle. The results show that the sheath edge velocity and total potential drop decrease for a steady-state sheath, which eventually collapses when the magnetic field lines become parallel to the wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗