Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TRAJECTORIES”

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,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Hotspot formation due to shock-induced pore collapse in 1,3,5,7-tetranitro-1,3,5,7-tetrazoctane (HMX): Role of pore shape and shock strength in collapse mechanism and temperature

The shock to detonation transition in heterogeneous high energy density solids starts with the spatial localization of mechanical energy into so-called hotspots that form due to the interaction between the leading wave and microstructural features and defects. We used large-scale molecular dynamics to characterize the hotspots resulting from the shock-induced collapse of cylindrical voids and elongated cracks focusing on the effect of shock strength, defect shape, and size. The temperature fields resulting from the collapse of cracks elongated along the shock direction show significantly higher sensitivity to both shock strength and size than cylindrical voids. Cracks 80 nm in length result in temperatures almost three times higher than voids 80 nm in diameter, reaching values corresponding to the ideal case of isentropic recompression of a gas. The molecular dynamics trajectories reveal the atomic origin of this contrasting behavior. While circular voids undergo a transition from viscoelastic pore collapse to a hydrodynamic regime with increasing shock strength, shock focusing in elongated cracks results in jetting and vaporization which, upon recompression, leads to increased heating.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Excited state dynamics of cis , cis -1,3-cyclooctadiene: UV pump VUV probe time-resolved photoelectron spectroscopy

We present UV pump, vacuum ultraviolet probe time-resolved photoelectron spectroscopy measurements of the excited state dynamics of cis,cis-1,3-cyclooctadiene. A 4.75 eV deep UV pump pulse launches a vibrational wave packet on the first electronically excited state, and the ensuing dynamics are probed via ionization using a 7.92 eV probe pulse. The experimental results indicate that the wave packet undergoes rapid internal conversion to the ground state in under 100 fs. Comparing the measurements with electronic structure and trajectory surface hopping calculations, we are able to interpret the features in the measured photoelectron spectra in terms of ionization to several states of the molecular cation.

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

Molecular simulation of liquid–vapor coexistence for NaCl: Full-charge vs scaled-charge interaction models

Scaled-charge models have been recently introduced for molecular simulations of electrolyte solutions and molten salts to attempt to implicitly represent polarizability. Although these models have been found to accurately predict electrolyte solution dynamic properties, they have not been tested for coexistence properties, such as the vapor pressure of the melt. In this work, we evaluate the vapor pressure of a scaled-charge sodium chloride (NaCl) force field and compare the results against experiments and a non-polarizable full-charge force field. The scaled-charge force field predicts a higher vapor pressure than found in experiments, due to its overprediction of the liquid-phase chemical potential. Reanalyzing the trajectories generated from the scaled-charge model with full charges improves the estimation of the liquid-phase chemical potential but not the vapor pressure.

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

Effect of Many Modes on Self-Polarization and Photochemical Suppression in Cavities

The standard description of cavity-modified molecular reactions typically involves a single (resonant) mode, while in reality the quantum cavity supports a range of photon modes. In this work we demonstrate that as more photon modes are accounted for, physico-chemical phenomena can dramatically change, as illustrated by the cavity-induced suppression of the important and ubiquitous process of proton-coupled electron-transfer. Using a multi-trajectory Ehrenfest treatment for the photonmodes, we find that self-polarization effects become essential, and we introduce the concept of selfpolarization- modified Born-Oppenheimer surfaces as a new construct to analyze dynamics. As the number of cavity photon modes increases, the increasing deviation of these surfaces from the cavityfree Born-Oppenheimer surfaces, together with the interplay between photon emission and absorption inside the widening bands of these surfaces, leads to enhanced suppression. The present findings are general and will have implications for the description and control of cavity-driven physical processes of molecules, nanostructures and solids embedded in cavities.

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

Long short-term memory embedded nudging schemes for nonlinear data assimilation of geophysical flows

Reduced rank nonlinear filters are increasingly utilized in data assimilation of geophysical flows, but often require a set of ensemble forward simulations to estimate forecast covariance. On the other hand, predictor-corrector type nudging approaches are still attractive due to their simplicity of implementation when more complex methods need to be avoided. However, optimal estimate of nudging gain matrix might be cumbersome. In this paper, we put forth a fully nonintrusive recurrent neural network approach based on a long short-term memory (LSTM) embedding architecture to estimate the nudging term, which plays a role not only to force the state trajectories to the observations but also acts as a stabilizer. Furthermore, our approach relies on the power of archival data and the trained model can be retrained effectively due to power of transfer learning in any neural network applications. In order to verify the feasibility of the proposed approach, we perform twin experiments using Lorenz 96 system. Our results demonstrate that the proposed LSTM nudging approach yields more accurate estimates than both extended Kalman filter (EKF) and ensemble Kalman filter (EnKF) when only sparse observations are available. With the availability of emerging AI-friendly and modular hardware technologies and heterogeneous computing platforms, we articulate that our simplistic nudging framework turns out to be computationally more efficient than either the EKF or EnKF approaches.

42 ENGINEERING↗
DOE OSTI2020

A novel, magnetically driven convergent Richtmyer–Meshkov platform

Here, we introduce a novel experimental platform for the study of the Richtmyer–Meshkov instability in a cylindrically converging geometry using a magnetically driven cylindrical piston. Magnetically driven solid liner implosions are used to launch a shock into a liquid deuterium working fluid and, ultimately, into an on-axis rod with a pre-imposed perturbation. The shock front trajectory is tracked through the working fluid and up to the point of impacting the rod through the use of on axis photonic Doppler velocimetry. This configuration allows for precise characterization of the shock state as it impacts the perturbed rod interface. Monochromatic x-ray radiography is used to measure the post-shock interface evolution and rod density profile. The ALEGRA MHD model is used to simulate the dynamics of the experiment in one dimension. We show that late in time the perturbation growth becomes non-linear as evidenced by the observation of high-order harmonics, up to n = 5. Two dimensional simulations performed using a combination of the GORGON MHD code and the xRAGE radiation hydrodynamics code suggest that the late time non-linear growth is modified by convergence effects as the bubbles and spikes experience differences in the pressure of the background flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Mixed quantum–classical approach to model non-adiabatic electron–nuclear dynamics: Detailed balance and improved surface hopping method

We develop a density matrix formalism to describe coupled electron–nuclear dynamics. To this end, we introduce an effective Hamiltonian formalism that describes electronic transitions and small (quantum) nuclear fluctuations along a classical trajectory of the nuclei. Using this Hamiltonian, we derive equations of motion for the electronic occupation numbers and for the nuclear coordinates and momenta. We show that, in the limit, when the number of nuclear degrees of freedom coupled to a given electronic transition is sufficiently high (i.e., the strong decoherence limit), the equations of motion for the electronic occupation numbers become Markovian. Furthermore, the transition rates in these (rate) equations are asymmetric with respect to the lower-to-higher energy transitions and vice versa. In thermal equilibrium, such asymmetry corresponds to the detailed balance condition. We also study the equations for the electronic occupations in the non-Markovian regime and develop a surface hopping algorithm based on our formalism. To treat the decoherence effects, we introduce additional “virtual” nuclear wave packets whose interference with the “real” (physical) wave packets leads to the reduction in coupling between the electronic states (i.e., decoherence) as well as to the phase shifts that improve the accuracy of the numerical approach. Remarkably, the same phase shifts lead to the detailed balance condition in the strong decoherence limit.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Evolutionary reinforcement learning of dynamical large deviations

In this work, we show how to bound and calculate the likelihood of dynamical large deviations using evolutionary reinforcement learning. An agent, a stochastic model, propagates a continuous-time Monte Carlo trajectory and receives a reward conditioned upon the values of certain path-extensive quantities. Evolution produces progressively fitter agents, potentially allowing the calculation of a piece of a large-deviation rate function for a particular model and path-extensive quantity. For models with small state spaces, the evolutionary process acts directly on rates, and for models with large state spaces, the process acts on the weights of a neural network that parameterizes the model's rates. This approach shows how path-extensive physics problems can be considered within a framework widely used in machine learning.

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

Using solenoid as multipurpose tool for measuring beam parameters

Solenoids are frequently used for focusing low-energy beams. In this paper, we show how they can serve as multipurpose diagnostics tools to measure various beam parameters, including energy, emittance, the second moments of the transverse distribution, and the beam position and angle with respect to the solenoid’s axis. The energy measurement is based on rotation of the plane of the transverse motion, as opposed to generating dispersion with a dipole. Measurement of the beam trajectory with respect to the solenoid axis is done by analyzing the beam orbit downstream of the solenoid while varying its current. The second moments are calculated by analyzing the beam image on a profile monitor while accounting for the beam rotation caused by the solenoid. We describe in detail the corresponding procedures and the experimental results of these measurements.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Raman studies of hydrogen trapped in As 4 O 6 ·2H 2 at high pressure and low temperature

In this work, Raman spectroscopic measurements of the arsenolite–hydrogen inclusion compound As 4 O 6 ·2H 2 were performed in diamond anvil cells at high pressure and variable temperature down to 80 K. The experimental results were complemented by ab initio molecular dynamics simulations and phonon calculations. Observation of three hydrogen vibrons in As 4 O 6 ·2H 2 is reported in the entire temperature and pressure range studied (up to 24 GPa). While the experiments performed with protium and deuterium at variable temperatures allowed for the assignment of two vibrons as Q 1 (1) and Q 1 (0) transitions of ortho and para spin isomers of hydrogen trapped in the inclusion compound, the origin of the third vibron could not be unequivocally established. Low-temperature spectra revealed that the lowest-frequency vibron is actually composed of two overlapping bands of A g and T 2g symmetries dominated by H 2 stretching modes as predicted by our previous density functional theory calculations. We observed low-frequency modes of As 4 O 6 ·2H 2 vibrations dominated by H 2 “librations,” which were missed in a previous study. A low-temperature fine structure was observed for the J = 0 → 2 and J = 1 → 3 manifolds of hydrogen trapped in As 4 O 6 ·2H 2 , indicating the lifting of degeneracy due to an anisotropic environment. A non-spherical distribution was captured by molecular dynamics simulations, which revealed that the trajectory of H 2 molecules is skewed along the crystallographic $\langle$111$\rangle$ direction. Last but not least, low-temperature synchrotron powder x-ray diffraction measurements on As 4 O 6 ·2H 2 revealed that the bulk structure of the compound is preserved down to 5 K at 1.6 GPa.

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

Nonlinear fluorescence spectroscopy of layered perovskite quantum wells

Interest in layered organohalide perovskites is motivated by their potential for use in optoelectronic devices. In these systems, the smallest and largest quantum wells are primarily concentrated near the glass and air interfaces of a film, thereby establishing a gradient in the average values of the bandgaps. It has been suggested that this layered architecture promotes the funneling of electronic excitations through space in a manner similar to light-harvesting processes in photosynthetic antennae. Whereas energy and charge transfer are difficult to distinguish by conventional transient absorption techniques, it has recently been shown that these competing relaxation mechanisms may be separately targeted with nonlinear fluorescence (NLFL) and photocurrent “action spectroscopies.” Here, we present perturbative rate functions to describe NLFL experiments conducted on layered perovskite systems. The formulas reproduce the patterns of resonances observed in experimental measurements and show how signatures of energy transfer manifest in two-dimensional spectra. Overall, this work suggests that NLFL spectroscopy may be used to fully reveal the trajectories of electronic excitations by correlating ultrafast energy transfer pathways to fluorescence emission from the thickest quantum wells.

Ouyang, Zhenyu↗
DOE OSTI2020

A generalized Boltzmann kinetic theory for strongly magnetized plasmas with application to friction

Coulomb collisions in plasmas are typically modeled using the Boltzmann collision operator, or its variants, which apply to weakly magnetized plasmas in which the typical gyroradius of particles significantly exceeds the Debye length. Conversely, O’Neil has developed a kinetic theory to treat plasmas that are so strongly magnetized that the typical gyroradius of particles is much smaller than the distance of the closest approach in a binary collision. Here, we develop a generalized collision operator that applies across the full range of magnetization strength. Since there is no closed-form solution for the scattering cross section when plasma is strongly magnetized, the input to the collision operator is obtained by numerically calculating particle trajectories. To demonstrate novel physics associated with strong magnetization, it is used to compute the friction force on a massive test charge. In addition to the traditional stopping power component, this is found to exhibit a transverse component that is perpendicular to both the velocity and Lorentz force vectors in the strongly magnetized regime, as was predicted recently using linear response theory. As a result, good agreement is found between the collision theory and linear response theory in the regime in which both apply, but the new collision theory also applies to stronger magnetization strength regimes than the linear response theory is expected to apply in.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Excited-state dynamics of CH 2 I 2 and CH 2 IBr studied with UV-pump VUV-probe momentum-resolved photoion spectroscopy

We perform time-resolved ionization spectroscopy measurements of the excited state dynamics of CH 2 I 2 and CH 2 IBr following photoexcitation in the deep UV. The fragment ions produced by ionization with a vacuum-ultraviolet probe pulse are measured with velocity map imaging, and the momentum resolved yields are compared with trajectory surface hopping calculations of the measurement observable. Together with recent time-resolved photoelectron spectroscopy measurements of the same dynamics, these results provide a detailed picture of the coupled electronic and nuclear dynamics involved. Overall, our measurements highlight the non-adiabatic coupling between electronic states, which leads to notable differences in the dissociation dynamics for the two molecules.

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

Implosion dynamics of triple-nozzle gas-puff z pinches on COBRA

Experiments on the 1-MA, 220-ns COBRA generator at Cornell University were conducted to provide detailed measurements of structured cylindrical gas-puff z pinches. In the experiments, a 7 cm diameter triple-nozzle gas valve assembly with concentric outer and inner annular nozzles and a central gas jet initialize the z-pinch load with various working gases, radial density profiles, and externally applied axial magnetic fields. Planar laser-induced fluorescence provides a measure of the initial neutral gas density of the load, while three-frame laser shearing interferometry and multi-frame extreme ultraviolet (XUV) cameras reveal the formation and propagation of a magneto-Rayleigh–Taylor (MRT) unstable shock layer. Implosion trajectories are compared to simple, experimentally informed models and found to be in good agreement. Differences in the structure of the accelerating plasma sheath and evolution of the MRT instability are observed for different gas species and axial magnetic field strengths, correlating with differences in pinch uniformity and x-ray emission. Here, the average instability growth is compared to linear MRT theory predictions using the instantaneous acceleration of the best-fit implosion models and characteristic instability wavelength, with the effective Atwood number and seed perturbation size as fit parameters. For high density argon center jets, ionization prior to the arrival of the imploding plasma sheath suggests a heating mechanism consistent with photoionization by XUV self-emission.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Case studies of the time-dependent potential energy surface for dynamics in cavities

Here, the exact time-dependent potential energy surface driving the nuclear dynamics was recently shown to be a useful tool to understand and interpret the coupling of nuclei, electrons, and photons in cavity settings. Here, we provide a detailed analysis of its structure for exactly solvable systems that model two phenomena: cavity-induced suppression of proton-coupled electron-transfer and its dependence on the initial state, and cavity-induced electronic excitation. We demonstrate the inadequacy of simply using a weighted average of polaritonic surfaces to determine the dynamics. Such a weighted average misses a crucial term that redistributes energy between the nuclear and the polaritonic systems, and this term can in fact become a predominant term in determining the nuclear dynamics when several polaritonic surfaces are involved. Evolving an ensemble of classical trajectories on the exact potential energy surface reproduces the nuclear wavepacket quite accurately, while evolving on the weighted polaritonic surface fails after a short period of time. The implications and prospects for application of mixed quantum-classical methods based on this surface are discussed.

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

Suppressing the enhancement of stimulated Raman scattering in inhomogeneous plasmas by tuning the modulation frequency of a broadband laser

The stimulated Raman scattering (SRS) instability can inhibit the performance of laser driven inertial confinement fusion (ICF) implosions by scattering light into unwanted directions or by generating hot electrons that preheat the target fuel. In principle, ICF target designs can avoid parameter regimes conducive to large, linear SRS gains. In practice, kinetic inflation—the nonlinear enhancement of SRS due to electron trapping in the excited plasma wave—makes this difficult. Here we show that laser bandwidth in the form of frequency modulation can either decrease or increase the inflationary SRS (iSRS) threshold in inhomogeneous plasmas depending on the maximum chirp of the laser pulse. The threshold, mapped out by a series of particle-in-cell simulations, exhibits a minimum when the frequency change within the pulse cancels the spatial detuning due to density inhomogeneities along the trajectory of the scattered light. In conclusion, by tuning the pump laser parameters away from this minimum, the iSRS threshold can be larger than at zero bandwidth, providing a path to mitigating kinetic inflation in ignition relevant plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Guiding center and gyrokinetic orbit theory for large electric field gradients and strong shear flows

The guiding center and gyrokinetic theory of magnetized particle motion is extended to the regime of large electric field gradients perpendicular to the magnetic field. A gradient in the electric field directly modifies the oscillation frequency and causes the Larmor orbits to deform from circular to elliptical trajectories. In order to retain a good adiabatic invariant, there can only be strong dependence on a single coordinate at lowest order, so that resonances do not generate chaotic motion that destroys the invariant. When the gradient across magnetic flux surfaces is dominant, the guiding center drift velocity becomes anisotropic in response to external forces and additional curvature drifts must be included. The electric polarization density remains gyrotropic, but both the polarization and magnetization are modified by the change in gyrofrequency. The theory can be applied to shear flows that are even stronger than those observed in the edge transport barrier of a high-performance tokamak (H-mode) pedestal, even if the toroidal field is as small as or even smaller than the poloidal field. Yet, the theory retains a mathematical form that is similar to the standard case and can readily be implemented within existing simulation tools.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

State predictive information bottleneck

We report the ability to make sense of the massive amounts of high-dimensional data generated from molecular dynamics simulations is heavily dependent on the knowledge of a low-dimensional manifold (parameterized by a reaction coordinate or RC) that typically distinguishes between relevant metastable states, and which captures the relevant slow dynamics of interest. Methods based on machine learning and artificial intelligence have been proposed over the years to deal with learning such low-dimensional manifolds, but they are often criticized for a disconnect from more traditional and physically interpretable approaches. To deal with such concerns, in this work we propose a deep learning based state predictive information bottleneck approach to learn the RC from high-dimensional molecular simulation trajectories. We demonstrate analytically and numerically how the RC learnt in this approach is connected to the committor in chemical physics and can be used to accurately identify transition states. A crucial hyperparameter in this approach is the time delay or how far into the future the algorithm should make predictions about. Through careful comparisons for benchmark systems, we demonstrate that this hyperparameter choice gives useful control over how coarse-grained we want the metastable state classification of the system to be. We thus believe that this work represents a step forward in systematic application of deep learning based ideas to molecular simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗