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

Variational design principles for nonequilibrium colloidal assembly

Using large deviation theory and principles of stochastic optimal control, we show that rare molecular dynamics trajectories conditioned on assembling a specific target structure encode a set of interactions and external forces that lead to enhanced stability of that structure. Such a relationship can be formulated into a variational principle, for which we have developed an associated optimization algorithm and have used it to determine optimal forces for targeted self-assembly within nonequilibrium steady-states. We illustrate this perspective on inverse design in a model of colloidal cluster assembly within linear shear flow. We find that colloidal clusters can be assembled with high yield using specific short-range interactions of tunable complexity. Shear decreases the yields of rigid clusters, while small values of shear increase the yields of nonrigid clusters. The enhancement or suppression of the yield due to shear is rationalized with a generalized linear response theory. Furthermore, by studying 21 unique clusters made of six, seven, or eight particles, we uncover basic design principles for targeted assembly out of equilibrium.

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

Potential energy surfaces for high-energy N + O 2 collisions

Potential energy surfaces for high-energy collisions between an oxygen molecule and a nitrogen atom are useful for modeling chemical dynamics in shock waves. In the present work, we present doublet, quartet, and sextet potential energy surfaces that are suitable for studying collisions of O 2 ( 3 Σ$^{–}_{g}$) with N( 4 S) in the electronically adiabatic approximation. Two sets of surfaces are developed, one using neural networks (NNs) with permutationally invariant polynomials (PIPs) and one with the least-squares many-body (MB) method, where a two-body part is an accurate diatomic potential and the three-body part is expressed with connected PIPs in mixed-exponential-Gaussian bond order variables (MEGs). We find, using the same dataset for both fits, that the fitting performance of the PIP-NN method is significantly better than that of the MB-PIP-MEG method, even though the MB-PIP-MEG fit uses a higher-order PIP than those used in previous MB-PIP-MEG fits of related systems (such as N 4 and N 2 O 2 ). However, the evaluation of the PIP-NN fit in trajectory calculations requires about 5 times more computer time than is required for the MB-PIP-MEG fit.

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

Particle resonances in toroidal fusion devices

Resonances of high energy particles in magnetic confinement devices due to electromagnetic instabilities can strongly modify the distribution, leading to a reduction in fusion power and even discharge termination and particle loss to the device walls through avalanche. The existence of a mode particle resonance depends on properties of the equilibrium, particle trajectories, and perturbation mode harmonic content. Resonance location is a function of particle energy and equilibrium field line helicity. Different methods for finding resonance location and energy dependence are developed. The properties of resonances are discussed using examples from magnetic fusion devices. Here, we show that if mode resonances exist at low particle energy, they very likely also exist at high energy, thus modifying high energy beam particles and fusion products. It is possible for a resonance to appear due to mode induced orbit helicity modification when it is forbidden at low mode amplitude.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Operator learning for predicting multiscale bubble growth dynamics

We report simulating and predicting multiscale problems that couple multiple physics and dynamics across many orders of spatiotemporal scales is a great challenge that has not been investigated systematically by deep neural networks (DNNs). Herein, we develop a framework based on operator regression, the so-called deep operator network (DeepONet), with the long-term objective to simplify multiscale modeling by avoiding the fragile and time-consuming “hand-shaking” interface algorithms for stitching together heterogeneous descriptions of multiscale phenomena. To this end, as a first step, we investigate if a DeepONet can learn the dynamics of different scale regimes, one at the deterministic macroscale and the other at the stochastic microscale regime with inherent thermal fluctuations. Specifically, we test the effectiveness and accuracy of the DeepONet in predicting multirate bubble growth dynamics, which is described by a Rayleigh–Plesset (R–P) equation at the macroscale and modeled as a stochastic nucleation and cavitation process at the microscale by dissipative particle dynamics (DPD). First, we generate data using the R–P equation for multirate bubble growth dynamics caused by randomly time-varying liquid pressures drawn from Gaussian random fields (GRFs). Our results show that properly trained DeepONets can accurately predict the macroscale bubble growth dynamics and can outperform long short-term memory networks. We also demonstrate that the DeepONet can extrapolate accurately outside the input distribution using only very few new measurements. Subsequently, we train the DeepONet with DPD data corresponding to stochastic bubble growth dynamics. Although the DPD data are noisy and we only collect sparse data points on the trajectories, the trained DeepONet model is able to predict accurately the mean bubble dynamics for time-varying GRF pressures. Taken together, our findings demonstrate that DeepONets can be employed to unify the macroscale and microscale models of the multirate bubble growth problem, hence providing new insight into the role of operator regression via DNNs in tackling realistic multiscale problems and in simplifying modeling with heterogeneous descriptions.

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

Disruption avoidance via radio frequency current condensation in magnetic islands produced by off-normal events

This article discusses the use of radio frequency (RF) current drive to stabilize large islands, focusing on nonlinear effects that appear when relatively high powers are used to stabilize large islands. We are interested in developing a capability to stabilize large islands via RF driven currents to avoid the need for mitigation to the extent possible. As tokamaks are designed and built with increasing levels of stored energy in the plasma, disruptions become increasingly dangerous. It has been reported that 95% of the disruptions in the Joint European Torus tokamak with the ITER-like wall are preceded by the growth of large locked islands. These large islands are mostly produced by off-normal events other than neoclassical tearing modes. This paper presents theory and modeling for a nonlinear “RF current condensation” effect that can concentrate the RF driven current near the center of a large island, thereby increasing the efficiency of the stabilization. A nonlinear shadowing effect can hinder the stabilization of islands if the aiming of the ray trajectories does not properly consider the nonlinear effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Ray-tracing analysis for cross-polarization scattering diagnostic on MAST-upgrade spherical tokamak

A combined Doppler backscattering/cross-polarization scattering (DBS/CPS) system is being deployed on MAST-U for simultaneous measurements of local density turbulence, turbulence flows, and magnetic turbulence. In this design, CPS shares the probing beam with the DBS and uses a separate parallel-viewing receiver system. In this study, we utilize a modified GENRAY 3D ray-tracing code to simulate the propagation of the probing and scattered beams. The contributions of different scattering locations along the entire beam trajectories are considered, and the corresponding local B̃ wavenumbers are estimated using the wavevector matching criterion. The wavenumber ranges of the local B̃ that are detectable to the CPS system are explored for simulated L- and H-mode plasmas.

Hong, R. (ORCID:0000000347508015)↗
DOE OSTI2021

Observation of a highly conductive warm dense state of water with ultrafast pump–probe free-electron-laser measurements

The electrical conductivity of water under extreme temperatures and densities plays a central role in modeling planetary magnetic fields. Experimental data are vital to test theories of high-energy-density water and assess the possible development and presence of extraterrestrial life. These states are also important in biology and chemistry studies when specimens in water are confined and excited using ultrafast optical or free-electron lasers (FELs). Here we utilize femtosecond optical lasers to measure the transient reflection and transmission of ultrathin water sheet samples uniformly heated by a 13.6 nm FEL approaching a highly conducting state at electron temperatures exceeding 20 000 K. The experiment probes the trajectory of water through the high-energy-density phase space and provides insights into changes in the index of refraction, charge carrier densities, and AC electrical conductivity at optical frequencies. At excitation energy densities exceeding 10 MJ/kg, the index of refraction falls to n = 0.7, and the thermally excited free-carrier density reaches ne = 5 × 1027 m−3, which is over an order of magnitude higher than that of the electron carriers produced by direct photoionization. Significant specular reflection is observed owing to critical electron density shielding of electromagnetic waves. The measured optical conductivity reaches 2 × 104 S/m, a value that is one to two orders of magnitude lower than those of simple metals in a liquid state. At electron temperatures below 15 000 K, the experimental results agree well with the theoretical calculations using density-functional theory/molecular-dynamics simulations. With increasing temperature, the electron density increases and the system approaches a Fermi distribution. In this regime, the conductivities agree better with predictions from the Ziman theory of liquid metals.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Cylinders and jets in crossflow: Wake formations as a result of varying geometric proximities

The combined flow physics of several canonical flow configurations is experimentally studied. Here, we analyze an array of jets issuing into a crossflow, then immediately navigating past an array of cylinders. This is achieved with a 2 × 3 triangular pattern of jets and symmetric cylinders at three jets to crossflow velocity ratios, enabling near-complete optical access of each jet, with velocities measured by time-resolved particle image velocimetry. Jet trajectories reveal that each configuration adheres to a power-law trend and that greater penetration is achieved by the downstream and confined jets compared to the more conventional upstream one. Recirculation regions of the upstream and downstream jets are nearly independent, with the confined jet encompassing regions of overlap with both. Turbulent statistics reveal the influence of geometric placement and velocity ratio on the time-averaged velocity, anisotropy, and Reynolds stresses incurred by each jet. Galilean decomposition utilizes a supplemental crossflow-only velocity field to delineate the influence of each jet's low- and high-pressure regions on the otherwise uniform stream. Proper orthogonal decomposition suggests that increased jet penetration decreases the number of modes required for truncation in the investigated spanwise plane. Vortex identification algorithms are applied to the reconstructed flow fields, reaffirming that with increasing velocity ratio, the jets generate vortices of their own in similar statistical formations as the cylinders. This investigation provides a foundation to aid future modeling efforts toward characterizing flow physics of importance in designing and passively controlling transverse jets with varying blockage proximities in a crossflow.

Kristo, Paul J. (ORCID:0000000225605175)↗
DOE OSTI2021

Dynamics and energy coupling of gas puff Z-pinches on a fast linear transformer driver

Gas puff Z-pinch experiments with annular Ar and Ne gas shells have been conducted on the Compact Experimental System for Z-pinch and Ablation Research (CESZAR) linear transformer driver (LTD) with 500 kA current and 160 ns rise time. Here, we present results from the first systematic gas puff Z-pinch experiments using a fast (≤ 200 ns) LTD as a driver, in which we show that 7% of the stored energy in the capacitors is coupled to plasma kinetic energy as estimated via self-emission and laser schlieren images. 0D and 1D simulations—which do not allow instability growth and thus reach greater maximum average velocities—using initial conditions inferred from experimental implosion trajectories predict coupling in excess of 10% of the stored energy. Additionally, the Ar and Ne implosions were comparably massed and thus achieved similar maximum kinetic energies, though the Ne pinches were more stable and the x-ray pulses were longer and produced higher yield: 2–5 ns and 0.21–0.52 J (0.15–0.37 J/cm) of Ar K-shell and 12–25 ns and 2.2–3.9 J (1.6–2.6 J/cm) of Ne K-shell, respectively. The difference in stability is most likely attributed to variations in initial conditions such as density distribution and gas breakdown initiation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

High-resolution planar electron beam induced current in bulk diodes using high-energy electrons

Understanding the impact of high-energy electron radiation on device characteristics remains critical for the expanding use of semiconductor electronics in space-borne applications and other radiation harsh environments. Here in this paper, we report on in situ measurements of high-energy electron radiation effects on the hole diffusion length in low threading dislocation density homoepitaxial bulk n-GaN Schottky diodes using electron beam induced current (EBIC) in high-voltage scanning electron microscopy mode. Despite the large interaction volume in this system, quantitative EBIC imaging is possible due to the sustained collimation of the incident electron beam. This approach enables direct measurement of electron radiation effects without having to thin the specimen. Using a combination of experimental EBIC measurements and Monte Carlo simulations of electron trajectories, we determine a hole diffusion length of 264 ± 11 nm for n-GaN. Irradiation with 200 kV electron beam with an accumulated dose of 24 × 10 16 electrons cm -2 led to an approximate 35% decrease in the minority carrier diffusion length.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Expanding the momentum field of view in angle-resolved photoemission systems with hemispherical analyzers

In photoelectron spectroscopy, the measured electron momentum range is intrinsically related to the excitation photon energy. Low photon energies <10 eV are commonly encountered in laser-based photoemission and lead to a momentum range that is smaller than the Brillouin zones of most materials. This can become a limiting factor when studying condensed matter with laser-based photoemission. An additional restriction is introduced by widely used hemispherical analyzers that record only electrons photoemitted in a solid angle set by the aperture size at the analyzer entrance. Here, we present an upgrade to increase the effective solid angle that is measured with a hemispherical analyzer. We achieve this by accelerating the photoelectrons toward the analyzer with an electric field that is generated by a bias voltage on the sample. Our experimental geometry is comparable to a parallel plate capacitor, and therefore, we approximate the electric field to be uniform along the photoelectron trajectory. With this assumption, we developed an analytic, parameter-free model that relates the measured angles to the electron momenta in the solid and verify its validity by comparing with experimental results on the charge density wave material TbTe3. By providing a larger field of view in momentum space, our approach using a bias potential considerably expands the flexibility of laser-based photoemission setups.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A kinetic model of friction in strongly coupled strongly magnetized plasmas

Plasmas that are strongly magnetized in the sense that the gyrofrequency exceeds the plasma frequency exhibit novel transport properties that are not well understood. As a representative example, we compute the friction force acting on a massive test charge moving through a strongly coupled and strongly magnetized one-component plasma using a generalized Boltzmann kinetic theory. Recent works studying the weakly coupled regime have shown that strong magnetization leads to a transverse component of the friction force that is perpendicular to both the Lorentz force and velocity of the test charge, in addition to the stopping power component aligned antiparallel to the velocity. Recent molecular dynamics simulations have also shown that strong Coulomb coupling in addition to strong magnetization gives rise to a third component of the friction force in the direction of the Lorentz force. Here, we show that the generalized Boltzmann kinetic theory captures these effects and generally agrees well with the molecular dynamics simulations over a broad range of Coulomb coupling and magnetization strength regimes. The theory is also used to show that the “gyro” component of the friction in the direction of the Lorentz force arises due to asymmetries associated with gyromotion during short-range collisions. Computing the average motion of the test charge through the background plasma, the transverse force is found to strongly influence the trajectory by changing the gyroradius and the gyrofriction force is found to slightly change the gyrofrequency of the test charge resulting in a phase shift.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Quasioptical modeling of wave beams with and without mode conversion. IV. Numerical simulations of waves in dissipative media

We report the first quasi-optical simulations of wave beams in a hot plasma using the quasi-optical code PARADE (PAraxial RAy DEscription) [K. Yanagihara, I. Y. Dodin, and S. Kubo, Phys. Plasmas 26, 072112 (2019)]. This code is unique in that it accounts for inhomogeneity of the dissipation-rate across the beam and mode conversion simultaneously. We show that the dissipation-rate inhomogeneity shifts beams relative to their trajectories in cold plasma and that the two electromagnetic modes are coupled via this process, an effect that was ignored in the past. Additionally, we propose a simplified approach to account for the dissipation-rate inhomogeneity. This approach is computationally inexpensive and simplifies the analysis of actual experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Particle emission with identification from an atmospheric pressure plasma liquid interface

A conspicuous emission phenomenon of luminous particles from the liquid anode surface of an atmospheric pressure DC glow discharge is reported. The emission has been shown to occur when the plasma forms a pattern on the surface of the liquid electrolyte. Here, the spatial-temporal evolutions of the trajectories of emitted particles were studied using a high-speed imaging system. Particles were sampled in flight using a witness plate for analysis. Critical particle characteristics were examined using electron microscope technologies. The morphology of the resulting splats and compositions suggest that they are molten droplets with a great deal of structure, including evidence of nanoprecipitation. A theoretical model was employed to estimate the particle size from the measured splats, which provides a radius range within a hundred micrometers. This experiment leads to the postulation of a mechanism in which particle emission at the plasma liquid interface is likely driven by the Taylor cone effect.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Scaling of silicon nanoparticle growth in low temperature flowing plasmas

Low temperature plasmas are an emerging method to synthesize high quality nanoparticles (NPs). An established and successful technique to produce NPs is using a capacitively coupled plasma (CCP) in cylindrical geometry. Although a robust synthesis technique, optimizing or specifying NP properties using CCPs, is challenging. In this paper, results from a computational investigation for the growth of silicon NPs in flowing inductively coupled plasmas (ICPs) using Ar/SiH 4 gas mixtures of up to a few Torr are discussed. ICPs produce more locally constrained and quiescent plasma potentials. These positive plasma potentials produce an electrostatic trap for negatively charged NPs, which can significantly extend the residence time of NPs in the plasma, which in turn provides a controllable period for particle growth. The computational platforms used in this study consist of a two-dimensional plasma hydrodynamics model, a three-dimensional nanoparticle growth and trajectory tracking model, and a molecular dynamics simulation for deriving reactive sticking coefficients of silane radicals on Si NPs. Trends for the nanoparticle growth as a function of SiH 4 inlet fraction, gas residence time, energy deposition per particle, pressure, and reactor diameter are discussed. The general path for particle synthesis is the trapping of small NPs in the positive electrostatic potential, followed by entrainment in the gas flow upon reaching a critical particle size. Optimizing or controlling NP synthesis then depends on the spatial distribution of plasma potential, the density of growth species, and the relative time that particles spend in the electrostatic trap and flowing through higher densities of growth species upon leaving the trap.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Influence of non-adiabatic effects on linear absorption spectra in the condensed phase: Methylene blue

Modeling linear absorption spectra of solvated chromophores is highly challenging as contributions are present both from coupling of the electronic states to nuclear vibrations and from solute–solvent interactions. In systems where excited states intersect in the Condon region, significant non-adiabatic contributions to absorption line shapes can also be observed. Here, we introduce a robust approach to model linear absorption spectra accounting for both environmental and non-adiabatic effects from first principles. This model parameterizes a linear vibronic coupling (LVC) Hamiltonian directly from energy gap fluctuations calculated along molecular dynamics (MD) trajectories of the chromophore in solution, accounting for both anharmonicity in the potential and direct solute–solvent interactions. The resulting system dynamics described by the LVC Hamiltonian are solved exactly using the thermalized time-evolving density operator with orthogonal polynomials algorithm (T-TEDOPA). The approach is applied to the linear absorption spectrum of methylene blue in water. We show that the strong shoulder in the experimental spectrum is caused by vibrationally driven population transfer between the bright S1 and the dark S2 states. The treatment of the solvent environment is one of many factors that strongly influence the population transfer and line shape; accurate modeling can only be achieved through the use of explicit quantum mechanical solvation. The efficiency of T-TEDOPA, combined with LVC Hamiltonian parameterizations from MD, leads to an attractive method for describing a large variety of systems in complex environments from first principles.

Dunnett, Angus J. (ORCID:0000000285797826)↗
DOE OSTI2021

The microscopic mechanism of bulk melting of ice

We study the initial stages of homogeneous melting of a hexagonal ice crystal at coexistence and at moderate superheating. Our trajectory-based computer simulation approach provides a comprehensive picture of the events that lead to melting, from the initial accumulation of 5+7 defects, via the formation of L-D and interstitial-vacancy pairs, to the formation of a liquid nucleus. Of the different types of defects that we observe to be involved in melting, a particular kind of 5+7 type defect (type 5) plays a prominent role as it often forms prior to the formation of the initial liquid nucleus and close to the site where the nucleus forms. Hence, like other solids, ice homogeneously melts via the prior accumulation of defects.

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

High-speed impact of micron-sized diesel drop trains—Splashing dynamics, secondary droplet formation, and effects of pre-existing film thickness

This work investigates surface impingement of monodispersed diesel drop trains using computational simulations. The three-dimensional simulations include high-speed impact of micrometer-sized fuel drops onto initially dry and wetted stainless steel substrates. The drop size and impact velocity are representative of fuel injection conditions in internal combustion engines (ICEs). The drop trains serve as a simplified representation of fuel spray. To accurately capture the film that naturally develops on initially dry substrates, a dynamic diesel-stainless steel contact angle model was developed using experiments of single drop impact under ICE representative temperature. Simulations are presented for a highly splashing case and the effects of a pre-existing film on splashing dynamics were investigated, including the temporal evolution of splashed mass and film thickness. It was concluded that for thin films, the effect of pre-existing film thickness is less significant than in single drop impingement. Secondary droplet characterization was performed on simulation results, leading to instantaneous and time-averaged distributions of secondary droplet size, velocity magnitude, and trajectory angle. It was found that for each drop impingement, approximately 58% of the splashed mass is from the impinging drop itself, while the remainder of the splashed mass is composed of film liquid. From a high resolution case, details of secondary droplet formation are observed and three distinct phases of secondary droplet formation are identified. The detailed analysis of drop train impingement under engine-relevant conditions serves as a first step toward a robust understanding of fuel impingement and the development of cleaner and more efficient ICEs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗