Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Equations of motion”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Spontaneous Exciton Collapse in a Strongly Flattened Ellipsoidal InSb Quantum Dot

Abstract Electronic and excitonic states in an InSb strongly flattened ellipsoidal quantum dot (QD) with complicated dispersion law are theoretically investigated within the framework of the geometric adiabatic approximation in the strong, intermediate, and weak quantum confinement regimes. For the lower levels of the spectrum, the square root dependence of energy on QD sizes is revealed in the case of Kane’s dispersion law. The obtained results are compared to the case of a parabolic (standard) dispersion law of charge carriers. The possibility of the accidental exciton instability is revealed for the intermediate quantum confinement regime. For the weak quantum confinement regime, the motion of the exciton's center-of-gravity is quantized, which leads to the appearance of additional Coulomb-like sub-levels. It is revealed that in the case of the Kane dispersion law, the Coulomb levels shift into the depth of the forbidden band gap, moving away from the quantum confined level, whereas in the case of the parabolic dispersion law, the opposite picture is observed. The corresponding selection rules of quantum transitions for the interband absorption of light are obtained. New selection rules of quantum transitions between levels conditioned by 2D exciton center of mass vertical motion quantization in a QD are revealed. The absorption threshold behavior characteristics depending on the QDs geometrical sizes are also revealed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A Quantum-mechanical Treatment of Electron Broadening in Strong Magnetic Fields

Abstract Spectral line-shape fitting is an extremely useful tool in determining the gravity of white dwarf stars. This method is so far limited to nonmagnetic white dwarfs largely because the theory of line broadening in high magnetic fields is not as complete as in the nonmagnetic case. Current Stark+Zeeman models treat plasma particles classically and ignore the motion of the nucleus. We develop the formalism for a quantum-mechanical treatment of the perturbing electrons and include the nuclear motion as part of the broadening and explore their relative importance. The conditions we explore are those found in white dwarf and neutron star atmospheres. We find that, contrary to previous studies, the quantized perturbing electrons create more broadening than perturbers on a straight-path trajectory. Additionally, the quantization of the plasma electrons gives rise to resonances away from the line center. The nuclear motion creates an additional electric field, which also leads to an increase in line broadening; however, this effect in neutron star atmospheres is not as large as previously estimated. This suggests that neutron star spectral lines are sensitive to density and that their mass and radius can be obtained from spectral line fitting, which would help constrain the neutron star equation of state.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Price equation captures the role of drug interactions and collateral effects in the evolution of multidrug resistance

Bacterial adaptation to antibiotic combinations depends on the joint inhibitory effects of the two drugs (drug interaction [DI]) and how resistance to one drug impacts resistance to the other (collateral effects [CE]). Here we model these evolutionary dynamics on two-dimensional phenotype spaces that leverage scaling relations between the drug-response surfaces of drug-sensitive (ancestral) and drug-resistant (mutant) populations. We show that evolved resistance to the component drugs – and in turn, the adaptation of growth rate – is governed by a Price equation whose covariance terms encode geometric features of both the two-drug-response surface (DI) in ancestral cells and the correlations between resistance levels to those drugs (CE). Within this framework, mean evolutionary trajectories reduce to a type of weighted gradient dynamics, with the drug interaction dictating the shape of the underlying landscape and the collateral effects constraining the motion on those landscapes. We also demonstrate how constraints on available mutational pathways can be incorporated into the framework, adding a third key driver of evolution. Our results clarify the complex relationship between drug interactions and collateral effects in multidrug environments and illustrate how specific dosage combinations can shift the weighting of these two effects, leading to different and temporally explicit selective outcomes.

Gjini, Erida↗
DOE OSTI2020

Long range frequency chirping of Alfvén eigenmodes

In this work, a theoretical framework has been developed for an NBI scenario to model the hard non-linear evolution of global Alfvén eigenmodes (GAEs) where the adiabatic motion of phase-space structures (holes and clumps), associated with frequency chirping, occurs in generalized phase-space of slowing down energetic particles. The radial profile of the GAE is expanded using finite elements, which allows an update of the mode structure as the mode frequency chirps. Constants of motion are introduced to track the dynamics of energetic particles during frequency chirping by implementing proper action-angle variables and canonical transformations which reduce the dynamics essentially to 1D. Consequently, we specify whether the particles are drifting inward/outward as the frequency deviates from the initial MHD eigenfrequency. Using the principle of least action, we have derived the non-linear equation describing the evolution of the radial profile by varying the total Lagrangian of the system with respect to the weights of the finite elements. For the choice of parameters in this work, it is shown that the peak of the radial profile is shifted and also broadens due to frequency chirping. The time rate of frequency change is also calculated using the energy balance and we show that the adiabatic condition remains valid once it is satisfied. This model clearly illustrates the theoretical treatment to study the long-range adiabatic frequency sweeping events observed for Alfvén gap modes in real experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Examining astrophysical gas cloud collapse using an optical depth-scaled, x-ray-irradiated, carbon-foam sphere

When stellar radiation interacts with a molecular cloud, the cloud's fate depends on the strength of the incident radiation and the radiation's mean-free-path within the cloud [F. Bertoldi, Astrophys. J. 346, 735–755 (1989)]. Under the right conditions, the radiation compresses the cloud and a star formation may occur. Where and when the stellar formation occurs in the cloud's collapse are open questions. Direct observation of the complete star–cloud lifecycle is nearly impossible due to the immense timescales and distances over which the interaction occurs. Laboratory astrophysics offers a way to investigate such a system by scaling the important astrophysical parameters to the laboratory. This work describes laboratory experiments to study the radiation-driven implosion of clouds, using x rays from a laser-irradiated, thin, gold foil as a surrogate star and a carbon-foam sphere as a surrogate cloud. An optically thick system, theoretically corresponding to a star-forming regime, was selected by choice of the foam density. Gold foil and sphere motions were imaged by x-ray radiography. Radiographic images show the formation of an interface between rarefied gold and carbon plasmas, a shock moving into the sphere, and a blunting of the initial sphere's shape. Measurements show that the shock moved linearly around 64 μm/ns into the sphere, and the gold–carbon interface formed by 2 ns at the sphere edge remained stationary. The deformation of the sphere was driven by the incident radiation and not by mechanical pressures applied by gold plasma. The blunting of the sphere was likely due to the geometric reduction of flux near the sphere's poles. Higher x-ray flux near the sphere's equator caused high compression and a faster shock, which flattened the sphere. We will discuss the results and implications of our observations.

VanDervort, R. W. [University of Michigan 1 , Ann ↗
DOE OSTI2022

Effect of Heave Plate Hydrodynamic Force Parameterization on a Two-Body Wave Energy Converter

Heave plates are one approach to generating the reaction force necessary to harvest energy from ocean waves. In a Morison equation description of the hydrodynamic force, the components of drag and added mass depend primarily on the heave plate oscillation. These terms may be parameterized in three ways: (1) as a single coefficient invariant across sea state, most accurate at the reference sea state, (2) coefficients dependent on the oscillation amplitude, but invariant in phase, that are most accurate for relatively small amplitude motions, and (3) coefficients dependent on both oscillation amplitude and phase, which are accurate for all oscillation amplitudes. We validate a MATLAB model for a two-body point absorber wave energy converter against field data and a dynamical model constructed in ProteusDS. We then use the MATLAB model to evaluate the effect of these parameterizations on estimates of heave plate motion, tension between the float and heave plate, and wave energy converter electrical power output. We find that power predictions using amplitude-dependent coefficients differ by up to 30% from models using invariant coefficients for regular waves ranging in height from 0.5 to 1.9 m. Amplitude- and phase-dependent coefficients, however, yield less than a 5% change when compared with coefficients dependent on amplitude only. This suggests that amplitude-dependent coefficients can be important for accurate wave energy converter modeling, but the added complexity of phase-dependent coefficients yields little further benefit. We show similar, though less pronounced, trends in maximum tether tension, but note that heave plate motion has only a weak dependence on coefficient fidelity. Finally, we emphasize the importance of using experimentally derived added mass over that calculated from boundary element methods, which can lead to substantial under-prediction of power output and peak tether tension.

dynamical model↗
DOE OSTI2020

Effect of finite ion mass on relativistic self-induced transparency of plasma layers with a sharp boundary

We consider the effect of ion mobility on the threshold of relativistic self-induced transparency under irradiation of thin plasma layers by circularly polarised laser radiation. An analytical model of the motion of ions during the removal of electrons from the layer surface by the ponderomotive force of laser radiation is constructed. The model is used to analyse the motion of probe electrons in the resulting electromagnetic field. It is shown that the higher the ion mobility and the longer the laser pulse, the more stable the plasma layer to longitudinal perturbations of a finite value and the higher the threshold of self-induced transparency. These conclusions are verified by one-dimensional numerical simulation of a complete system of kinetic equations for plasma and Maxwell’s equations. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Dynamics of metastable contact soliton dissipative exchange flows in one-dimensional ferromagnetic channels

Dissipative exchange flows (DEFs) are large-amplitude boundary value solutions of ferromagnetic channels. In their low-injection limit, DEFs reduce to spin superfluids. However, in the strong injection limit, nonlinearities dominate close to the injection site and a soliton is formed; this solution has been termed a contact soliton dissipative exchange flow (CS-DEF). Here, in this work, we numerically investigate CS-DEF solutions in a moderate injection regime and a finite injection width. We find a solution where two metastable solitons coexist in the injection region. This solution is metastable in the sense that any perturbation to the system will eject one of the solitons out of the injection region. Moreover, soliton dynamics can be excited when two injection regions are separated by a certain distance. We find that the ensuing DEF between the solitons induces a steady-state dynamics in which metastable solitons are continually ejected and nucleated. Furthermore, and depending on the relative signs of the spin injections, the soliton dynamics possess a particular handedness and frequency related to the spin transfer torque delivered by the DEF. Our results provide insights into the transport of spin current by DEFs - where the interaction between DEFs and solitons suggests a mechanism for detaching contact-solitons from the injection boundary. Although this study focuses on the "nonlocal" interaction between solitons, it may lead to the investigation of new mechanisms for inserting solitons in a DEF, e.g., for discrete motion and transport of information over long distances.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Ab initio generalized Langevin equation

We introduce a machine learning–based approach called ab initio generalized Langevin equation (AIGLE) to model the dynamics of slow collective variables (CVs) in materials and molecules. In this scheme, the parameters are learned from atomistic simulations based on ab initio quantum mechanical models. Force field, memory kernel, and noise generator are constructed in the context of the Mori–Zwanzig formalism, under the constraint of the fluctuation–dissipation theorem. Combined with deep potential molecular dynamics and electronic density functional theory, this approach opens the way to multiscale modeling in a variety of situations. Here, we demonstrate this capability with a study of two mesoscale processes in crystalline lead titanate, namely the field-driven dynamics of a planar ferroelectric domain wall, and the dynamics of an extensive lattice of coarse-grained electric dipoles. In the first case, AIGLE extends the reach of ab initio simulations to a regime of noise-driven motions not accessible to molecular dynamics. In the second case, AIGLE deals with an extensive set of CVs by adopting a local approximation for the memory kernel and retaining only short-range noise correlations. The scheme is computationally more efficient than molecular dynamics by several orders of magnitude and mimics the microscopic dynamics at low frequencies where it reproduces accurately the dominant far-infrared absorption frequency.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Holographic Space-Time and Quantum Information

The formalism of Holographic Space-time (HST) is a translation of the principles of Lorentzian geometry into the language of quantum information. Intervals along time-like trajectories, and their associated causal diamonds, completely characterize a Lorentzian geometry. The Bekenstein-Hawking-Gibbons-’t Hooft-Jacobson-Fischler-Susskind-Bousso Covariant Entropy Principle, equates the logarithm of the dimension of the Hilbert space associated with a diamond to one quarter of the area of the diamond’s holographic screen, measured in Planck units. The most convincing argument for this principle is Jacobson’s derivation of Einstein’s equations as the hydrodynamic expression of this entropy law. In that context, the null energy condition (NEC) is seen to be the analog of the local law of entropy increase. The quantum version of Einstein’s relativity principle is a set of constraints on the mutual quantum information shared by causal diamonds along different time-like trajectories. The implementation of this constraint for trajectories in relative motion is the greatest unsolved problem in HST. The other key feature of HST is its claim that, for non-negative cosmological constant or causal diamonds much smaller than the asymptotic radius of curvature for negative c.c., the degrees of freedom localized in the bulk of a diamond are constrained states of variables defined on the holographic screen. This principle gives a simple explanation of otherwise puzzling features of BH entropy formulae, and resolves the firewall problem for black holes in Minkowski space. It motivates a covariant version of the CKN [1] bound on the regime of validity of quantum field theory (QFT) and a detailed picture of the way in which QFT emerges as an approximation to the exact theory

tensor network↗
DOE OSTI2024

Elucidating the Mechanism of Helium Evaporation from Liquid Water

Here, we investigate the evaporation of trace amounts of helium solvated in liquid water using molecular dynamics simulations and theory. Consistent with experimental observations, we find a super-Maxwellian distribution of kinetic energies of evaporated helium. This excess of kinetic energy over typical thermal expectations is explained by an effective continuum theory of evaporation based on a Fokker–Planck equation, parametrized molecularly by a potential of mean force and position-dependent friction. Using this description, we find that helium evaporation is strongly influenced by the friction near the interface, which is anomalously small near the Gibbs dividing surface due to the ability of the liquid–vapor interface to deform around the gas particle. Our reduced description provides a mechanistic interpretation of trace gas evaporation as the motion of an underdamped particle in a potential subject to a viscous environment that varies rapidly across the air–water interface. From it we predict the temperature dependence of the excess kinetic energy of evaporation, which is yet to be measured.

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

Fluid-like shock compression of dilute polymer nanocomposites

The shock Hugoniot of heterogenous mixtures of discrete particles has been experimentally investigated for porous agglomerates and fluids, while the study of full-density solids has been primarily limited to compressed powders and high-dimensional composites. By dispersing ceria nanoparticles in a polyethylene matrix, we are able to examine the hydrodynamic behavior of a nonporous, heterogenous solid in thermal equilibrium during weak shock compression. Phase-driven discontinuities in the Hugoniot particle velocity–shock velocity (u−D) relationship of pure polyethylene are replicated in the nanocomposites but are shifted to lower velocity and to higher pressure with higher particle concentration. The results are explained using an isothermal, two-velocity fluid model under the hydrodynamic approximation. The model, which assumes a theoretical equation-of-state for ceria and either a low-order or high-order fit to the measured polyethylene Hugoniot, reasonably predicts the Hugoniot for two different polyethylene/ceria nanocomposites. Using the model, the mixture Hugoniot is shown to be insensitive to the Hugoniot of the stiffer constituent when the moduli are sufficiently disparate, while dependence on particle density and volume fraction is preserved through fluid-like motion.

Moore, Nathan W. [Sandia National Laboratories (SN↗
DOE OSTI2025

Mapping Spatiotemporal Solvent Velocity from Measured Concentration Gradients in a Polarized Electrolyte

The electric-field induced motion of neutral species impedes the efficacy of electrochemical devices. By combining operando X-ray transmission measurements with continuum mechanics, we have developed a methodology for determining the velocity of neutral solvent molecules under an applied field. The X-ray transmission experiments were used to determine ion concentration profiles as a function of space and time in a polymer electrolyte. The unsteady state solvent mass balance equation was solved numerically with experimental concentration profiles to map spatiotemporal solvent velocities. We compare our experimentally derived results with predictions made with concentrated solution theory. We use the cation transference number as the only adjustable parameter to match experimental measurements of both concentration and solvent velocity. Our approach may be used to determine solvent velocity with any operando technique used to measure time-dependent ion concentration profiles.

Abdo, Emily E↗
DOE OSTI2023

Forecasts on the Dark Matter Density Profiles of Dwarf Spheroidal Galaxies with Current and Future Kinematic Observations

We forecast parameter uncertainties on the mass profile of a typical Milky Way dwarf spheroidal galaxy (dSph) using the spherical Jeans equation and Fisher matrix formalism. For a Draco-like system we show that radial velocity measurements for 1000 individual stars can constrain the mass contained within the effective radius of a dSph to within 5%. This is consistent with constraints extracted from current observational data. We compare two systems, a cusp and core, and demonstrate that a minimum sample of 100,000 (10,000) stars with both radial and proper motions measurements is required to disentangle their inner slopes at the 2σ (1σ) level. If using the log-slope measured at the half-light radius as a proxy for differentiating between a core or cusp slope, only 1000 line-of-sight and proper motions measurements are required; however, we show this choice of radius does not always unambiguously differentiate between core and cusped profiles. Once observational errors are below half the value of the intrinsic dispersion, improving the observational precision yields little change in the density profile uncertainties. The choice of priors in our profile shape analysis plays a crucial role when the number of stars in a system is less than 100 but does not affect the resulting uncertainties for larger kinematic samples. Our predicted 2D confidence regions agree well with those from a full likelihood analysis run on a mock kinematic data set taken from the Gaia Challenge, validating our Fisher predictions. Our methodology is flexible, allowing us to predict density profile uncertainties for a wide range of current and future kinematic data sets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

View factor estimation of hot spot velocities in inertial confinement fusion implosions at the National Ignition Facility

nertial confinement fusion (ICF) experiments at the National Ignition Facility (NIF) suffer from asymmetries in thex-ray drive that degrade capsule performance compared to expectations for a symmetric one-dimensional implosion.Mode 1, or pole-to-pole, drive asymmetry can reduce confinement and implosion efficiency, driving a bulk motion ofthe hot spot that is detectable with neutron diagnostics. Understanding and removing sources of mode 1 asymmetry inICF implosions is important for improving performance, and the three-dimensional nature of the problem makes high-resolution radiation-hydrodynamic modeling extremely computationally expensive. This work describes a reducedorder view factor model that calculates the drive asymmetry induced by beam-to-beam variations in laser delivery andhohlraum diagnostic windows along the equator. Capsule response is estimated by coupling to a Green’s function thatrelates final hot spot velocity to the applied time-varying mode 1 asymmetry. The model makes several predictionsabout the impact of mode 1 drivers such as laser delivery and target misalignment, and achieves good agreementin both magnitude and vector direction for several shots in three families of high-performance platforms. However,notable discrepancies suggest other potential sources of mode 1 asymmetry not captured by the model are also at play

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Ultrafast x-ray imaging of coherently controlled molecular dynamics in real space and time

Coherent control aims to manipulate chemical processes on the latent length and timescales of atoms and bonds, scales that are intrinsic to molecular dynamics but not directly resolved by most experimental probes. As a result, existing coherent-control experiments, which overwhelmingly rely on spectroscopic observables, leave a crucial blind spot: the direct, real-space recovery of all atomic and molecular rearrangements in response to coherently controlled excitations. Here, we overcome this limitation by integrating a Tannor-Kosloff-Rice pump-control-probe scheme with ultrafast X-ray scattering to capture snapshots of the wavepacket motion in a benchmark molecular system. We demonstrate this by photoexciting diatomic iodine vapor with a visible pump pulse, selectively steering the wavepacket toward ground-state recombination or dissociative pathways with a time-delayed near-IR control pulse, and recording the dynamics with angstrom and femtosecond precision with an ultrashort hard X-ray probe pulse. By comparing these structural observations with numerical solutions of the time-dependent Schrödinger equation, we reveal how coherent control actively reshapes the molecular charge density distribution. Our results pave the way for leveraging structural feedback as a control handle and provide a fundamental microscopic visualization of quantum decoherence and energy redistribution at the atomic level.

Hopper, Thomas R. [SLAC National Accelerator Labor↗
DOE OSTI2021

In-Home Mobility Frequency and Stability in Older Adults Living Alone With or Without MCI: Introduction of New Metrics

Background: Older adults spend a considerable amount of time inside their residences; however, most research investigates out-of-home mobility and its health correlates. We measured indoor mobility using room-to-room transitions, tested their psychometric properties, and correlated indoor mobility with cognitive and functional status. Materials and Methods: Community-dwelling older adults living alone ( n = 139; age = 78.1 ± 8.6 years) from the Oregon Center for Aging & Technology (ORCATECH) and Minority Aging Research Study (MARS) were included in the study. Two indoor mobility features were developed using non-parametric parameters (frequency; stability): Indoor mobility frequency (room-to-room transitions/day) was detected using passive infrared (PIR) motion sensors fixed on the walls in four geographic locations (bathroom; bedroom; kitchen; living room) and using door contact sensors attached to the egress door in the entrance. Indoor mobility stability was estimated by variances of number of room-to-room transitions over a week. Test-retest reliability (Intra-class coefficient, ICC) and the minimal clinically important difference (MCID) defined as the standard error of measurement (SEM) were generated. Generalized estimating equations models related mobility features with mild cognitive impairment (MCI) and functional status (gait speed). Results: An average of 206 days (±127) of sensor data were analyzed per individual. Indoor mobility frequency and stability showed good to excellent test-retest reliability (ICCs = 0.91[0.88–0.94]; 0.59[0.48–0.70]). The MCIDs of mobility frequency and mobility stability were 18 and 0.09, respectively. On average, a higher indoor mobility frequency was associated with faster gait speed (β = 0.53, p = 0.04), suggesting an increase of 5.3 room-to-room transitions per day was associated with an increase of 10 cm/s gait speed. A decrease in mobility stability was associated with MCI (β = −0.04, p = 0.03). Discussion: Mobility frequency and stability in the home are clinically meaningful and reliable features. Pervasive-sensing systems deployed in homes can objectively reveal cognitive and functional status in older adults who live alone.

Wu, Chao-Yi↗
DOE OSTI2022

Simulation of dislocation evolution in microparticle impacts over a wide range of impact velocities

We report the dynamic loading environment created by high velocity microparticle impacts is very difficult to model. The extremely high strain rates, large deformations, and high temperatures affect the flow stress of the material in ways that standard flow stress models cannot capture. Plastic deformation and strain hardening in metals is controlled by the motion of dislocations. Dislocations can be nucleated, stored as forest dislocations, or be annihilated as loading progresses. A comprehensive accounting of dislocation density change is needed to accurately describe rate and temperature dependent dislocation glide and evolution across the wide rage of loading conditions present in microparticle impact problems. Therefore, we implement and apply a newly proposed flow stress model (Hunter and Preston, 2015, 2022) that is then coupled with mobile and immobile dislocation density evolution equations. This model is implemented in Los Alamos National Laboratory’s hydrodynamics code, FLAG, to model copper-on-copper microparticle impacts. This new strength model allows for accurate simulation of particle rebound and flattening across a wide range of impact velocities.

36 MATERIALS SCIENCE↗