Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “chemical dynamics”

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 235 records · Page 13

PolyODENet: Deriving mass-action rate equations from incomplete transient kinetics data

Kinetics of a reaction network that follows mass-action rate laws can be described with a system of ordinary differential equations (ODEs) with polynomial right-hand side. However, it is challenging to derive such kinetic differential equations from transient kinetic data without knowing the reaction network, especially when the data are incomplete due to experimental limitations. We introduce a program, PolyODENet, toward this goal. Based on the machine-learning method Neural ODE, PolyODENet defines a generative model and predicts concentrations at arbitrary time. As such, it is possible to include unmeasurable intermediate species in the kinetic equations. Importantly, we have implemented various measures to apply physical constraints and chemical knowledge in the training to regularize the solution space. Here, using simple catalytic reaction models, we demonstrate that PolyODENet can predict reaction profiles of unknown species and doing so even reveal hidden parts of reaction mechanisms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Theoretical and experimental studies of atmospheric structure and dynamics, using high altitude chemical release, Radio meteor, and meteorological rocket network and other data

Data collected by the Georgia Tech Radio Meteor Wind Facility during the fall and winter of 1975 are analyzed indicating a relationship between lower thermospheric circulation at mid latitudes and polar stratospheric dynamics. Techniques of measurement of mixing processes in the upper atmosphere and the interpretation of those measurements are described along with a diffusion simulation program based on the Global Reference Atmosphere program.

Edwards, H. D.↗

In Situ Probe Science at Saturn

A fundamental goal of solar system exploration is to understand the origin of the solar sys-tem, the initial stages, conditions, and processes by which the solar system formed, how the formation pro-cess was initiated, and the nature of the interstellar seed material from which the solar system was born. Key to understanding solar system formation and subsequent dynamical and chemical evolution is the origin and evolution of the giant planets and their atmospheres. Several theories have been put forward to explain the process of solar system formation, and the origin and evolution of the giant planets and their atmospheres. Each theory offers quantifiable predictions of the abundances of noble gases He, Ne, Ar, Kr, and Xe, and abundances of key isotopic ratios 4He3He, DH, 15N14N, 18O16O, and 13C12C. Detection of certain dis-equilibrium species, diagnostic of deeper internal pro-cesses and dynamics of the atmosphere, would also help discriminate between competing theories. Measurements of the critical abundance profiles of these key constituents into the deeper well-mixed at-mosphere must be complemented by measurements of the profiles of atmospheric structure and dynamics at high vertical resolution and also require in situ explora-tion. The atmospheres of the giant planets can also serve as laboratories to better understand the atmospheric chem-istries, dynamics, processes, and climates on all planets including Earth, and offer a context and provide a ground truth for exoplanets and exoplanetary systems. Additionally, Giant planets have long been thought to play a critical role in the development of potentially habitable planetary systems. In the context of giant planet science provided by the Galileo, Juno, and Cassini missions to Jupiter and Sat-urn, a small, relatively shallow Saturn probe capable of measuring abundances and isotopic ratios of key at-mospheric constituents, and atmospheric structure in-cluding pressures, temperatures, dynamics, and cloud locations and properties not accessible by remote sens-ing can serve to test competing theories of solar system and giant planet origin, chemical, and dynamical evolution.

Interstellar↗

Modeling atmospheric data and identifying dynamics Temporal data-driven modeling of air pollutants

Atmospheric modelling has recently experienced a surge with the advent of deep learning. Most of these models, however, predict concentrations of pollutants following a data-driven approach in which the physical laws that govern their behaviors and relationships remain hidden. Furthermore, with the aid of real-world air quality data collected hourly in different stations throughout Madrid, we present a case study using a series of data-driven techniques with the following goals: (1) Find systems of ordinary differential equations that model the concentration of pollutants and their changes over time; (2) assess the performance and limitations of our model using stability analysis; (3) reconstruct the time series of chemical pollutants not measured in certain stations using delay coordinate embedding results.

54 ENVIRONMENTAL SCIENCES↗

The effect of breaking gravity waves on the dynamics and chemistry of the mesosphere and lower thermosphere (invited review)

The influence of breaking gravity waves on the dynamics and chemical composition of the 60 to 110 km region is investigated with a two dimensional model that includes a parameterization of gravity wave momentum deposition and diffusion. The dynamical model is described by Garcia and Solomon (1983) and Solomon and Garcia (1983) and includes a complete chemical scheme for the mesosphere and lower thermosphere. The parameterization of Lindzen (1981) is used to calculate the momentum deposited and the turbulent diffusion produced by the gravity waves. It is found that wave momentum deposition drives a very vigorous mean meridional circulation, produces a very cold summer mesopause and reverse the zonal wind jets above about 85 km. The seasonal variation of the turbulent diffusion coefficient is consistent with the behavior of mesospheric turbulences inferred from MST radar echoes. The large degree of consistency between model results and various types of dynamical and chemical data supports very strongly the hypothesis that breaking gravity waves play a major role in determining the zonally-averaged dynamical and chemical structure of the 60 to 110 km region of the atmosphere.

Garcia, R. R.↗

On the Dynamics of Computing a Chemically Relaxed Nonequilibrium Flow

In order to gain insights into the strong dependence of numerical solutions on initial data for finite time steps, a set of nonlinear test problems rich enough to capture the behavior of difference schemes were recently identified and the numerical basins of attraction for these problems were computed using commonly used time discretizations in CFD. Our study revealed a wealth of surprisingly nonlinear behavior of numerical schemes that were not observed before, in particular for the implicit time discretizations that are commonly used in CFD. The goal of this work is to apply these tools to study a practical model from non-equilibrium flowfield relaxation. This type of problem arises in chemically nonequilibrium hypersonic flows such as in a shock tube experiment or an expanding nozzle. Here we consider a reacting mixture of (N2, N) for an inviscid one-dimensional steady model. Preliminary numerical results indicate that, aside from the possibility of spurious numerical solutions being introduced by the time discretizations, limitations on the model for physical or accurate solutions may also play a part in the dynamics observed.

Sweby, P. K.↗

Ultrafast time-resolved x-ray absorption spectroscopy of ionized urea and its dimer through ab initio nonadiabatic dynamics

Investigating the early dynamics of chemical systems following ionization is essential for our understanding of radiation damage. However, experimental as well as theoretical investigations are very challenging due to the complex nature of these processes. Time-resolved x-ray absorption spectroscopy on a femtosecond timescale, in combination with appropriate simulations, is able to provide crucial insights into the ultrafast processes that occur upon ionization due to its element-specific probing nature. In this theoretical study, we investigate the ultrafast dynamics of valence-ionized states of urea and its dimer employing Tully's fewest switches surface hopping approach using Koopmans' theorem to describe the ionized system. We demonstrate that following valence ionization through a pump pulse, the time-resolved x-ray absorption spectra at the carbon, nitrogen, and oxygen K-edges reveal rich insights into the dynamics. Excited states of the ionized system give rise to time-delayed blueshifts in the x-ray absorption spectra as a result of electronic relaxation dynamics through nonadiabatic transitions. Moreover, our statistical analysis reveals specific structural dynamics in the molecule that induce time-dependent changes in the spectra. For the urea monomer, we elucidate the possibility to trace effects of specific molecular vibrations in the time-resolved x-ray absorption spectra. For the urea dimer, where ionization triggers a proton transfer reaction, we show how the x-ray absorption spectra can reveal specific details on the progress of proton transfer.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Multichannel dynamics in the OH+ n -butane reaction revealed by crossed-beam slice imaging and quasiclassical trajectory calculations

Multidimensional reactions present various channels that can exhibit very different dynamics and give products of varying subsequent reactivity. Here, we present a combination of experiment and theory to reveal the dynamics of hydrogen abstraction by OH radical at primary and secondary sites in n-butane at a collision energy of 8 kcal/mol. Crossed molecular beam slice imaging experiments unequivocally probe the secondary abstraction channel showing backward angular distributions with mild energy release to product translation, which are accurately captured by trajectory calculations using a specific-reaction-parameter Hamiltonian. Additionally, experiments containing both reaction channels indicate a less marked backward character in the angular distribution, whose origin is shown by trajectory calculations to appear as an evolution toward more sideways scattering from the secondary to primary channel. While the two channels have markedly different angular distributions, their energy release is largely comparable, showing ample energy release into the water product. The synergistic combination of crossed-beam imaging and trajectories opens the door to detailed reaction-dynamics studies of chemical reactions with ever-increasing complexity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The quantum dynamics of electronically nonadiabatic chemical reactions

Considerable progress was achieved on the quantum mechanical treatment of electronically nonadiabatic collisions involving energy transfer and chemical reaction in the collision of an electronically excited atom with a molecule. In the first step, a new diabatic representation for the coupled potential energy surfaces was created. A two-state diabatic representation was developed which was designed to realistically reproduce the two lowest adiabatic states of the valence bond model and also to have the following three desirable features: (1) it is more economical to evaluate; (2) it is more portable; and (3) all spline fits are replaced by analytic functions. The new representation consists of a set of two coupled diabatic potential energy surfaces plus a coupling surface. It is suitable for dynamics calculations on both the electronic quenching and reaction processes in collisions of Na(3p2p) with H2. The new two-state representation was obtained by a three-step process from a modified eight-state diatomics-in-molecules (DIM) representation of Blais. The second step required the development of new dynamical methods. A formalism was developed for treating reactions with very general basis functions including electronically excited states. Our formalism is based on the generalized Newton, scattered wave, and outgoing wave variational principles that were used previously for reactive collisions on a single potential energy surface, and it incorporates three new features: (1) the basis functions include electronic degrees of freedom, as required to treat reactions involving electronic excitation and two or more coupled potential energy surfaces; (2) the primitive electronic basis is assumed to be diabatic, and it is not assumed that it diagonalizes the electronic Hamiltonian even asymptotically; and (3) contracted basis functions for vibrational-rotational-orbital degrees of freedom are included in a very general way, similar to previous prescriptions for locally adiabatic functions in various quantum scattering algorithms.

Truhlar, Donald G.↗

A window on the composition of the early solar nebula: 2014MU69, Pluto, and Phoebe

The initial chemical composition of any solar nebula will depend upon the degree to which 1) organic and ice components form on dust grains, 2) organic and molecular species form in the gas phase, 3) organics and ices are exchanged between the gas and solid state, and 4) the precursor and newly formed (more complex) materials survive and are modified in the developing planetary system. Infrared and radio observations of star-forming regions reveal that complex chemistry occurs on icy grains, sometimes before stars even form. Additional processing, through the protosolar disk and within the solar nebula further modifies most, but probably not all, of the initial materials. In fact, the modern Solar System still carries a fraction of its interstellar inheritance [Alexander et al., 2017]. Here we focus on three examples of small bodies in our Solar System, each containing chemical and dynamical clues to its origin and evolution: the small, cold-classical Kuiper Belt object (KBO) 2014MU69 , Pluto and Saturn’s moon, Phoebe. The New Horizons flyby of 2014 MU69 has given the first view of an unaltered body composed of material originally in the solar nebula at ~45 AU. The spectrum reveals methanol ice (not commonly found), a possible detection of water ice (Stern et al. 2019). Pluto’s internal and surface inventory of volatiles and complex organics, together with active geological processes including cryo-volcanism, indicate a surprising level of activity on a body in the outermost region of the Solar System, and the fluid that emerges from subsurface reservoirs may contain material inherited from the solar nebula (Cruikshank et al. 2019). Meanwhile, Saturn's captured moon, Phoebe, carries high D/H in H2O [Clark et al. 2018)], and complex organics (Cruikshank et al. 2008), both consistent with its formation in, and inheritance from, the outer solar nebula. Together, these objects provide windows on the origin and evolution of our Solar System and constraints to be considered in future chemical and physical models of PPDs. The Spitzer Space telescope discovered the ring of dust around Saturn that emanated from the bombardment of Phoebe in the recent past, dislodging primitive Solar System material from deep within that satellite. In addition, Spitzer has provided detailed studies of the chemistry of star forming regions where methanol and other ices are readily observed. As a result of these and other observations, protoplanetary disk models are beginning to combine both dynamical and chemical complexities into their codes, resulting in new insights for exoplanetary systems.

Yvonne J. Pendleton↗

Agent-Based Chemical Plume Tracing Using Fluid Dynamics

This paper presents a rigorous evaluation of a novel, distributed chemical plume tracing algorithm. The algorithm is a combination of the best aspects of the two most popular predecessors for this task. Furthermore, it is based on solid, formal principles from the field of fluid mechanics. The algorithm is applied by a network of mobile sensing agents (e.g., robots or micro-air vehicles) that sense the ambient fluid velocity and chemical concentration, and calculate derivatives. The algorithm drives the robotic network to the source of the toxic plume, where measures can be taken to disable the source emitter. This work is part of a much larger effort in research and development of a physics-based approach to developing networks of mobile sensing agents for monitoring, tracking, reporting and responding to hazardous conditions.

Zarzhitsky, Dimitri↗

Expansion dynamics and chemistry evolution in ultrafast laser filament produced plasmas

Laser ablation in conjunction with optical emission spectroscopy is a potential non-contact, standoff detection method for all elements in the periodic table and certain isotopes such as radionuclides. Currently, significant development efforts are on-going to use ultrafast laser filaments for remote sensing. The application of filaments is of particular interest in extending the range of stand-off capability associated with elemental and isotopic detection via laser-induced breakdown spectroscopy. In this study, we characterize the expansion dynamics and chemical evolution of filament-produced uranium (U) plasmas. Laser filaments are generated in the laboratory by loosely focusing 35 femtosecond (fs), 6 milli Joule (mJ) pulses in air. Two-dimensional spectrally-integrated and time-resolved imaging was performed to study hydrodynamics and evolution of U atomic and UO molecular emission in filament-produced U plasmas. Our results highlight that filament ablation of U plasmas gives a cylindrical plasma plume morphology with an appearance of plume splitting into slow and fast moving components at later times of its evolution. Emission from the slow-moving region shows no distinct spectral features (i.e. broadband-like) and is contributed in part by nanoparticles generated during ultrafast laser ablation. Additionally, we find U atoms and U oxide molecules (i.e. UO, UxOy) co-exist in the filament produced plasma, which can be attributed to the generation of low-temperature plasma conditions during filament ablation

(140.3440) Laser-induced breakdown, (280.5395) Pla↗

ALMA Spectral Line and Imaging Survey of a Low and a High Mass-Loss Rate AGB Star Between 335 and 362 GHz

Context. Low and intermediate mass stars are known to power strong stellar winds when evolving through the asymptotic giant branch (AGB) phase. Initial mass, luminosity, temperature, and composition determine the pulsation characteristics of the star and the dust species formed in the pulsating photospheric layers. Radiation pressure on these grains triggers the onset of a stellar wind. However, as of today, we still cannot predict the wind mass-loss rates and wind velocities from first principles neither do we know which species are the first to condense in the upper atmospheric regions. Aims. We aim to characterise the dominant physical, dynamical, and chemical processes in the inner wind region of two archetypical oxygen-rich (C/O < 1) AGB stars, that is, the low mass-loss rate AGB star R Dor (M’ ∼ 1 × 10(exp −7) M/yr) and the high mass-loss rate AGB star IK Tau (M’ ∼ 5 × 10(exp −6) M/yr). The purpose of this study is to observe the key molecular species contributing to the formation of dust grains and to cross-link the observed line brightnesses of several species to the global and local properties of the star and its wind. Methods. A spectral line and imaging survey of IK Tau and R Dor was made with ALMA between 335 and 362 GHz (band 7) at a spatial resolution of ∼150 mas, which corresponds to the locus of the main dust formation region of both targets. Results. Some two hundred spectral features from 15 molecules (and their isotopologues) were observed, including rotational lines in both the ground and vibrationally excited states (up to v = 5 for SiO). Detected species include the gaseous precursors of dust grains such as SiO, AlO, AlOH, TiO, and TiO2. We present a spectral atlas for both stars and the parameters of all detected spectral features. A clear dichotomy for the sulphur chemistry is seen: while CS, SiS, SO, and SO2 are abundantly present in IK Tau, only SO and SO2 are detected in R Dor. Also other species such as NaCl, NS, AlO, and AlOH display a completely different behaviour. From some selected species, the minor isotopologues can be used to assess the isotopic ratios. The channel maps of many species prove that both large and small-scale inhomogeneities persist in the inner wind of both stars in the form of blobs, arcs, and/or a disk. The high sensitivity of ALMA allows us to spot the impact of these correlated density structures in the spectral line profiles. The spectral lines often display a half width at zero intensity much larger than expected from the terminal velocity, v∞, previously derived for both objects (36 km/s versus v(∞) ∼ 17.7 km/s for IK Tau and 23 km/s versus v(∞) ∼ 5.5 km/s for R Dor). Both a more complex 3D morphology and a more forceful wind acceleration of the (underlying) isotropic wind can explain this trend. The formation of fractal grains in the region beyond ∼400 mas can potentially account for the latter scenario. From the continuum map, we deduce a dust mass of ∼3.7 × 10(exp −7) Mand ∼2 × 10(exp −8) M for IK Tau and R Dor, respectively. Conclusions. The observations presented here provide important constraints on the properties of these two oxygen-dominated AGB stellar winds. In particular, the ALMA data prove that both the dynamical and chemical properties are vastly different for this high mass-loss rate (IK Tau) and low mass-loss rate (R Dor) star.

L. Decin↗

ALMA Spectral Line and Imaging Survey of a Low and a High Mass-Loss Rate AGB Star Between 335 and 362 GHz

Context. Low and intermediate mass stars are known to power strong stellar winds when evolving through the asymptotic giant branch (AGB) phase. Initial mass, luminosity, temperature, and composition determine the pulsation characteristics of the star and the dust species formed in the pulsating photospheric layers. Radiation pressure on these grains triggers the onset of a stellar wind. However, as of today, we still cannot predict the wind mass-loss rates and wind velocities from first principles neither do we know which species are the first to condense in the upper atmospheric regions. Aims. We aim to characterise the dominant physical, dynamical, and chemical processes in the inner wind region of two archetypical oxygen-rich (C/O<1) AGB stars, that is, the low mass-loss rate AGB star R Dor (M~ 1x10(exp 7) M☉/yr) and the high mass-loss rate AGB star IK Tau (M~5x10(exp 6) M☉/yr). The purpose of this study is to observe the key molecular species contributing to the formation of dust grains and to cross-link the observed line brightnesses of several species to the global and local properties of the star and its wind. Methods. A spectral line and imaging survey of IK Tau and R Dor was made with ALMA between 335 and 362 GHz (band 7) at a spatial resolution of ~150 mas, which corresponds to the locus of the main dust formation region of both targets. Results. Some two hundred spectral features from 15 molecules (and their isotopologues) were observed, including rotational lines in both the ground and vibrationally excited states (up to v=5 for SiO). Detected species include the gaseous precursors of dust grains such as SiO, AlO, AlOH, TiO, and TiO2.We present a spectral atlas for both stars and the parameters of all detected spectral features. A clear dichotomy for the sulphur chemistry is seen: while CS, SiS, SO, and SO2 are abundantly present in IK Tau, only SO and SO2 are detected in R Dor. Also other species such as NaCl, NS, AlO, and AlOH display a completely different behaviour. From some selected species, the minor isotopologues can be used to assess the isotopic ratios. The channel maps of many species prove that both large and small-scale inhomogeneities persist in the inner wind of both stars in the form of blobs, arcs, and/or a disk. The high sensitivity of ALMA allows us to spot the impact of these correlated density structures in the spectral line profiles. The spectral lines often display a half width at zero intensity much larger than expected from the terminal velocity, v1, previously derived for both objects (36 km/s versus v(∞)~17.7 km/s for IK Tau and 23 km/s versus v(∞)~5.5 km/s for R Dor). Both a more complex 3D morphology and a more forceful wind acceleration of the (underlying) isotropic wind can explain this trend. The formation of fractal grains in the region beyond _400 mas can potentially account for the latter scenario. From the continuum map, we deduce a dust mass of ~3.7 x 10(exp 7) M☉ and ~2 x 10(ex0 8) M☉ for IK Tau and R Dor, respectively. Conclusions. The observations presented here provide important constraints on the properties of these two oxygen-dominated AGB stellar winds. In particular, the ALMA data prove that both the dynamical and chemical properties are vastly different for this high mass-loss rate (IK Tau) and low mass-loss rate (R Dor) star.

Stars: AGB and post-AGB, Stars: mass loss, Stars: ↗

Chemical structure and curing dynamics of bisphenol S, PEEK™-like, and resveratrol phthalonitrile thermoset resins

This effort provides a multifaceted analysis of the structural changes and material dynamics of thermally driven softening and curing of three distinct phthalonitrile (PN) resins that cross-link into thermally stable and oxidation-resistant thermosets. Although this material system has yielded a large subset of fire-retardant composites that require facile processing and low-temperature curing, to date, insufficient information had been available on the fundamental processes that drive their softening and curing stages. Our approach conducted a complementary analysis the chemistry, monomer mobility, and rheology of three PN polymers in order to correlate the curing processes with corresponding structural and behavioral transformations of thermosets. Here, we focused on PNs with a bisphenol S backbone, a bisphenol A (PEEK™-like) backbone, and a resveratrol backbone. We relied on quasi-elastic neutron scattering (QENS) in order to analyze the in situ dynamics and self-diffusion properties of PN monomers, and to track changes in their mobilities during cross-linking and staging. Our analysis facilitates proper control over the staging and final curing of these resins and enables more efficient processing of these thermosets.

36 MATERIALS SCIENCE↗

Temperature and density dependent pair potential for deuterium under shock

Large-scale classical molecular dynamics (CMD) simulations naturally include the microscopic physics necessary for atomistic modeling of shock release at the ablator-fuel interface in an inertial confinement fusion (ICF) capsule. Here, the multi-megabar shocks utilized in ICF experiments can drive the deuterium fuel from ambient to electron volt temperatures (T) and multi-fold compression. Modeling interatomic interactions over such an extreme range of conditions is challenging for empirical bond order potentials. We generate a pair potential for deuterium with explicit temperature and mass density dependence from ab initio density functional theory molecular dynamics using the iterative Boltzmann inversion method. This potential accurately reproduces the radial distribution functions and pressures from DFT in CMD equilibrium simulations across a wide range of thermodynamic conditions, yet fails to return the expected Hugoniot relations when used in direct CMD shock simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Energy conservation in real-time nuclear–electronic orbital Ehrenfest dynamics

Real-time nuclear–electronic orbital Ehrenfest (RT-NEO-Ehrenfest) dynamics methods provide a first-principles approach for describing nonadiabatic molecular processes with nuclear quantum effects. For an efficient description of proton transfer within RT-NEO-Ehrenfest dynamics, the basis function center associated with the quantum proton can be allowed to move classically. Here, this traveling proton basis (TPB) approach effectively captures proton quantum dynamics, although its energy conservation behavior is not yet fully satisfactory. Two recently proposed TPB approaches, in principle, conserve the extended energy, which includes both the system energy and the kinetic energy associated with the proton basis function center. Herein, a thermostatted TPB approach is proposed to improve the conservation of the system energy, excluding the kinetic energy associated with the proton basis function center. In this approach, the quantum proton dynamics are modulated by dynamically rescaling the proton momentum operator to maintain the system energy conservation. With the excited-state intramolecular proton transfer of o-hydroxybenzaldehyde as an example, this approach is shown to significantly improve the system energy conservation while preserving the accuracy of the quantum proton dynamics as achieved in the original TPB approach.

Ehrenfest dynamics↗