Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase space 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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Extracting dynamical frequencies from invariants of motion in finite-dimensional nonlinear integrable systems

Integrable dynamical systems play an important role in many areas of science, including accelerator and plasma physics. An integrable dynamical system with n degrees of freedom possesses n nontrivial integrals of motion, and can be solved, in principle, by covering the phase space with one or more charts in which the dynamics can be described using action-angle coordinates. To obtain the frequencies of motion, both the transformation to action-angle coordinates and its inverse must be known in explicit form. However, no general algorithm exists for constructing this transformation explicitly from a set of n known (and generally coupled) integrals of motion. In this study we describe how one can determine the dynamical frequencies of the motion as functions of these n integrals in the absence of explicitly known action-angle variables, and we provide several examples.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Cluster Dynamics Simulations of Intra-Granular Fission Gas Bubble Size and Pressure Evolution in UO 2

Fission gases such as xenon (Xe) play a critical role in determining the behavior and response of nuclear fuel. Given that Xe has little solubility in UO 2 , it accumulates and forms bubbles, which significantly impact fuel performance. Intra- and inter-granular bubble nucleation and growth can lead to fuel swelling, and once bubbles interconnect at grain boundaries, fission gas can be released into the plenum. At low temperatures, limited uranium vacancy mobility can restrict swelling, therefore causing the bubbles to become highly pressurized. Consequently, this can induce micro-cracking, promote fission gas release (increasing the likelihood of cladding failure), and even lead to fuel pulverization under accident conditions such as a loss of coolant accident. As bubble evolution is strongly influenced by local temperature and fission rate, markedly different behavior occurs across the radial profile of the fuel pellet. Capturing the mechanisms that underpin bubble evolution is therefore important to predict these behaviors in the fuel. Previous models describing important mechanisms informed by lower length scale simulations have been developed under the NEAMS program. These can describe the evolution of a single bubble type (i.e., single value for radius and pressure) at each position in the pellet, for instance using the Centipede cluster dynamic code. However, in reality, a full distribution in bubble sizes and pressures exists within the microstructure at a given position in the pellet. To address this the cluster dynamics code Xolotl, which can predict Xe and vacancy phase space (i.e., bubble distributions) for intra-granular bubbles, has been used before. Prior work benchmarked the Xolotl code against the Centipede cluster dynamics code to ensure compatibility and to verify that mobile defect properties are adequately transferred between the two codes, along with some physics improvements. In this work, we go further by introducing a physics-based set of improvements that will allow us to accurately predict bubble size distributions and internal bubble pressures under representative UO 2 irradiation conditions. The improvements include (i) coupling bubble-defect reaction energies to a virial equation of state (EOS), (ii) including a bubble surface tension contribution, (iii) incorporating radiation-induced re-solution of Xe and vacancies, (iv) enabling pressure-driven dislocation loop punching through an effective emission of interstitial clusters informed by interstitial loop energetics, (v) accounting for radiation induced athermal diffusion of Xe, and (vi) implementing a Booth-type grain boundary sink representation for all mobile defects and defect clusters. After these modifications, we observe good agreement of Xolotl fission gas bubble size and concentration predictions with legacy experimental measurements. Additionally, it allows the distribution of Xe bubble pressures and radius to also be predicted and compared to data produced through the Advanced Fuels Campaign (AFC) program. Here, we have done this by running simulations under conditions similar to the AFC post-irradiation examination (PIE) samples irradiated at North Anna 2 light water reactor (LWR). Our results shows excellent agreement with these experimental measurements.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

A collision operator for describing dissipation in noncanonical phase space

The phase space of a noncanonical Hamiltonian system is partially inaccessible due to dynamical constraints (Casimir invariants) arising from the kernel of the Poisson tensor. When an ensemble of noncanonical Hamiltonian systems is allowed to interact, dissipative processes eventually break the phase space constraints, resulting in a thermodynamic equilibrium described by a Maxwell–Boltzmann distribution. However, the time scale required to reach Maxwell–Boltzmann statistics is often much longer than the time scale over which a given system achieves a state of thermal equilibrium. Examples include diffusion in rigid mechanical systems, as well as collisionless relaxation in magnetized plasmas and stellar systems, where the interval between binary Coulomb or gravitational collisions can be longer than the time scale over which stable structures are self-organized. Here, we focus on self-organizing phenomena over spacetime scales such that particle interactions respect the noncanonical Hamiltonian structure, but yet act to create a state of thermodynamic equilibrium. We derive a collision operator for general noncanonical Hamiltonian systems, applicable to fast, localized interactions. This collision operator depends on the interaction exchanged by colliding particles and on the Poisson tensor encoding the noncanonical phase space structure, is consistent with entropy growth and conservation of particle number and energy, preserves the interior Casimir invariants, reduces to the Landau collision operator in the limit of grazing binary Coulomb collisions in canonical phase space, and exhibits a metriplectic structure. We further show how thermodynamic equilibria depart from Maxwell–Boltzmann statistics due to the noncanonical phase space structure, and how self-organization and collisionless relaxation in magnetized plasmas and stellar systems can be described through the derived collision operator.

Boltzmann equation↗
DOE OSTI2025

A Perspective on Quantum Computing Applications in Quantum Chemistry Using 25-100 Logical Qubits

The intersection of quantum computing and quantum chemistry represents a promising frontier for achieving quantum utility in domains of both scientific and societal relevance. Owing to the exponential growth of classical resource requirements for simulating quantum systems, quantum chemistry has long been recognized as a natural candidate for quantum computation. This perspective focuses on identifying scientifically meaningful use cases where early fault-tolerant quantum computers, which are considered to be equipped with approximately 25-100 logical qubits, could deliver tangible impact. While recent advances in classical computing have pushed the boundaries of tractable simulations to unprecedented scales, this logical-qubit regime represents the first window where quantum devices can pursue qualitatively distinct strategies, such as polynomial-scaling phase estimation, direct simulation of quantum dynamics, and active-space embedding, that remain challenging for classical solvers, such as multireference charge-transfer and conical-intersection states central to photochemistry and materials design. We highlight near-term opportunities in algorithm and software design, discuss representative chemical problems suited for quantum acceleration, and propose strategic roadmaps and collaborative pathways for advancing practical quantum utility in quantum chemistry.

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

Exact Coherent Structures and Phase Space Geometry of Preturbulent 2D Active Nematic Channel Flow

Confined active nematics exhibit rich dynamical behavior, including spontaneous flows, periodic defect dynamics, and chaotic “active turbulence.” Here, we study these phenomena using the framework of exact coherent structures, which has been successful in characterizing the routes to high Reynolds number turbulence of passive fluids. Exact coherent structures are stationary, periodic, quasiperiodic, or traveling wave solutions of the hydrodynamic equations that, together with their invariant manifolds, serve as an organizing template of the dynamics. We compute the dominant exact coherent structures and connecting orbits in a preturbulent active nematic channel flow, which enables a fully nonlinear but highly reduced-order description in terms of a directed graph. Using this reduced representation, we compute instantaneous perturbations that switch the system between disparate spatiotemporal states occupying distant regions of the infinite-dimensional phase space. Overall, our results lay the groundwork for a systematic means of understanding and controlling active nematic flows in the moderate- to high-activity regime.

36 MATERIALS SCIENCE↗
DOE OSTI2021

On the Formation of Phantom Electron Phase Space Density Peaks in Single Spacecraft Radiation Belt Data

This study examines the rapid losses and acceleration of trapped relativistic and ultrarelativistic electron populations in the Van Allen radiation belt during the September 7–9, 2017 geomagnetic storm. By analyzing the dynamics of the last closed drift shell (LCDS), and the electron flux and phase space density (PSD), we show that the electron dropouts are consistent with magnetopause shadowing and outward radial diffusion to the compressed LCDS. During the recovery phase, an in-bound pass of Van Allen Probe A shows an apparent local peak in PSD, which does not exist in reality. Here, a careful analysis of the multipoint measurements by the Van Allen Probes reveals instead how the apparent PSD peak arises from aliasing monotonic PSD profiles which are rapidly increasing due to acceleration from very fast inwards radial diffusion. In the absence of such multisatellite conjunctions during fast acceleration events, such peaks might otherwise be associated with local acceleration processes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Three-Dimensional Dynamics of a Magnetic Hopfion Driven by Spin Transfer Torque

Magnetic hopfion is a three-dimensional (3D) topological soliton with novel spin structure that would enable exotic dynamics. Herein, we observe and research the current-driven 3D dynamics of a magnetic hopfion with a unit Hopf index in a frustrated magnet. Attributed to the spin Berry phase and symmetry of the hopfion, the phase space entangles multiple collective coordinates, thus the hopfion exhibits rich dynamics including longitudinal motion along the current direction, transverse motion perpendicular to the current direction, rotational motion, and dilation. Furthermore, the characteristics of hopfion dynamics is determined by the ratio between the nonadiabatic spin transfer torque parameter and the damping parameter. Such peculiar 3D dynamics of magnetic hopfion could shed light on understanding the universal physics of hopfions in different systems and boost the prosperous development of 3D spintronics.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Ultrafast Dynamics of Excitonic Complexes in Electrolyte-Gated Monolayer MoS 2

Understanding and controlling ultrafast excitons and trion dynamics in monolayer transition metal dichalcogenides (TMDs) is critical for optoelectronic and photonic applications. These dynamics depend strongly on carrier density, but most studies use fluence-dependence to modulate electron-hole populations rather than direct electrostatic gating that is more relevant to optoelectronic devices. We utilize electrolyte gating to tune the ground-state carrier density of monolayer MoS 2 and probe excitonic dynamics using transient absorption spectroscopy, thereby modifying absorption through bandgap renormalization, screening, phase-space filling, and exciton-trion crossover. Spectral deconvolution reveals distinct but coupled exciton and trion dynamics. Excitons form independently of voltage but relax faster with increasing carrier density, consistent with Auger-like processes, while trions decay more rapidly through multiple channels. At higher voltages, trion relaxation shifts from a sub-picosecond cooling to slower trapping-assisted processes. Furthermore, our results provide missing insights into the interplay between excitons and trions in monolayer MoS 2 , relevant for TMD-based optoelectronics.

chalcogenides↗
DOE OSTI2023

Dynamic Aperture of the PIP-II Accumulator Ring (PAR)

The Proton Improvement Plan-II (PIP-II) will include a superconducting linear accelerator which will accelerate H- ions to 800MeV, at an intensity never seen before. A proposed PIP-II Accumulator Ring (PAR) can reduce the injection losses to the booster, enhance the ramp up of the Long Baseline Neutrino Facility (LBNF) and support experiments in Dark Sector Physics. When the beam is circulating in PAR at its full intensity, it will then be transferred to the Booster via single turn injection. Strong PAR focusing results in larger phase advances, shorter focal lengths and higher machine tunes, which creates robust beam envelope control. Particles will be able to survive modest excursions from the elliptical phase-space trajectories of linear optics, and any loss of particles would be due to reduced dynamic aperture. The values for tunes and chromaticities are flexible and will continue to evolve as the lattice matures. Since there will be non-linear terms impacting the particles’ trajectories, tracking will include realistic magnetic field errors and installation alignment errors to determine the dynamic aperture as a function of the machine tunes. This will allow us to determine the optimum operating point in phase-space.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Kinetically constrained freezing transition in a dipole-conserving system

Here, we study a stochastic lattice gas of particles in one dimension with strictly finite-range interactions that respect the fractonlike conservation laws of total charge and dipole moment. As the charge density is varied, the connectivity of the system's charge configurations under the dynamics changes qualitatively. We find two distinct phases: Near half filling the system thermalizes subdiffusively, with almost all configurations belonging to a single dynamically connected sector. As the charge density is tuned away from half filling there is a phase transition to a frozen phase, where locally active finite bubbles cannot exchange particles and the system fails to thermalize. The two phases exemplify what has recently been referred to as weak and strong Hilbert space fragmentation, respectively. We study the static and dynamic scaling properties of this weak-to-strong fragmentation phase transition in a kinetically constrained classical Markov circuit model, obtaining some conjectured exact critical exponents.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Asymptotic optimality of twist-untwist protocols for Heisenberg scaling in atom-based sensing

Twist-untwist protocols for quantum metrology consist of a serial application of (1) unitary nonlinear dynamics (e.g., spin squeezing or Kerr nonlinearity), (2) parameterized dynamics U ( φ ) (e.g., a collective rotation or phase space displacement), and (3) time reversed application of step 1. Such protocols are known to produce states that allow Heisenberg scaling for experimentally accessible estimators of φ even when the nonlinearities are applied for times much shorter than required to produce Schrödinger cat states. In this work, we prove that, asymptotically in the number of particles, twist-untwist protocols provide the lowest estimation error among quantum metrology protocols that utilize two calls to a weakly nonlinear evolution and a readout involving only first and second moments of a total spin operator n → · J → . We consider the following physical settings: all-to-all interactions generated by one-axis twisting J z 2 (e.g., interacting Bose gases), constant finite range spin-spin interactions of distinguishable or bosonic atoms (e.g., trapped ions or Rydberg atoms, or lattice bosons). In these settings, we further show that the optimal twist-untwist protocols asymptotically achieve 85% and 92% of the respective quantum Cramér-Rao bounds. We show that the error of a twist-untwist protocol can be decreased by a factor of L without an increase in the noise of the spin measurement if the twist-untwist protocol can be noiselessly iterated as an L layer quantum alternating operator ansatz.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Atoms to fibers: Identifying novel processing methods in the synthesis of pitch-based carbon fibers

Understanding and optimizing the key mechanisms used in the synthesis of pitch-based carbon fibers (CFs) are challenging, because unlike polyacrylonitrile-based CFs, the feedstock for pitch-based CFs is chemically hetero- geneous, resulting in complex fabrication leading to inconsistency in the final properties. In this work, we use molecular dynamics simulations to explore the processing and chemical phase space through a framework of CF models to identify their effects on elastic performance. The results are in excellent agreement with experiments. We find that density, followed by alignment, and functionality of the molecular constituents dictate the CF mechanical properties more strongly than their size and shape. Last, we propose a previously unexplored fabrication route for high-modulus CFs. Unlike graphitization, this results in increased sp 3 fraction, achieved via generating high-density CFs. In addition, the high sp 3 fraction leads to the fabrication of CFs with isometric compressive and tensile moduli, enabling their potential applications for compressive loading.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

A Sparse-Grid Probabilistic Scheme for Approximation of the Runaway Probability of Electrons in Fusion Tokamak Simulation

Runaway electrons (RE) generated during magnetic disruptions present a major threat to the safe operation of plasma nuclear fusion reactors. A critical aspect of understanding RE dynamics is to calculate the runaway probability, i.e., the probability that an electron in the phase space will runaway on, or before, a prescribed time. Such probability can be obtained by solving the adjoint equation of the underlying Fokker-Planck equation that controls the electron dynamics. In this effort, we present a sparse-grid probabilistic scheme for computing the runaway probability. The key ingredient of our approach is to represent the solution of the adjoint equation as a conditional expectation, such that discretizing the differential operator reduces to the approximation of a set of integrals. Adaptive sparse grid interpolation is utilized to approximate the map from the phase space to the runaway probability. The main novelties of this effort are the integration of the sparse-grid method into the probabilistic numerical scheme for computing escape probability, and the application of the proposed method in computing RE probabilities. Two numerical examples are given to illustrate that the proposed method can achieve O(Δt) convergence, and that the local anisotropic adaptive refinement strategy (M. Stoyanov, Adaptive sparse grid construction in a context of local anisotropy and multiple hierarchical parents. In: Sparse Grids and Applications-Miami 2016, Springer, Berlin, 2018, pp. 175–199) can effectively handle the sharp transition layer between the runaway and non-runaway regions.

Yang, Minglei↗
DOE OSTI2022

Unsupervised discovery of nonlinear plasma physics using differentiable kinetic simulations

Plasma supports collective modes and particle–wave interactions that lead to complex behaviour in, for example, inertial fusion energy applications. While plasma can sometimes be modelled as a charged fluid, a kinetic description is often crucial for studying nonlinear effects in the higher-dimensional momentum–position phase space that describes the full complexity of the plasma dynamics. We create a differentiable solver for the three-dimensional partial-differential equation describing the plasma kinetics and introduce a domain-specific objective function. Using this framework, we perform gradient-based optimization of neural networks that provide forcing function parameters to the differentiable solver given a set of initial conditions. We apply this to an inertial-fusion-relevant configuration and find that the optimization process exploits a novel physical effect.

plasma nonlinear phenomena↗
DOE OSTI2021

Accurate and confident prediction of electron beam longitudinal properties using spectral virtual diagnostics

Abstract Longitudinal phase space (LPS) provides a critical information about electron beam dynamics for various scientific applications. For example, it can give insight into the high-brightness X-ray radiation from a free electron laser. Existing diagnostics are invasive, and often times cannot operate at the required resolution. In this work we present a machine learning-based Virtual Diagnostic (VD) tool to accurately predict the LPS for every shot using spectral information collected non-destructively from the radiation of relativistic electron beam. We demonstrate the tool’s accuracy for three different case studies with experimental or simulated data. For each case, we introduce a method to increase the confidence in the VD tool. We anticipate that spectral VD would improve the setup and understanding of experimental configurations at DOE’s user facilities as well as data sorting and analysis. The spectral VD can provide confident knowledge of the longitudinal bunch properties at the next generation of high-repetition rate linear accelerators while reducing the load on data storage, readout and streaming requirements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Dynamically generated momentum space shell structure of quarkyonic matter via an excluded volume model

The phase-space structure of zero-temperature quarkyonic matter is a Fermi sphere of quark matter surrounded by a shell of nucleonic matter. Here, we construct a quasiparticle model of quarkyonic matter based on the constituent quark model, where the quark and nucleon masses are related by $m_Q = m_N/N_c$, and $N_c$ is the number of quark colors. The region of occupied states is for quarks $k_Q < k_F/N_c$ and for nucleons $k_F < k_N < k_F + Δ$. We first consider the general problem of quarkyonic matter with hard-core nucleon interactions. We then specialize to a quasiparticle model where the hard-core nucleon interactions are accounted for by an excluded volume. In this model, we show that the nucleonic shell forms past some critical density related to the hard-core size and for large densities becomes a thin shell. We explore the basic features of such a model and argue this model has the semiquantitative behavior needed to describe neutron stars.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Machine learning for collective variable discovery and enhanced sampling in biomolecular simulation

Classical molecular dynamics simulates the time evolution of molecular systems through the phase space spanned by the positions and velocities of the constituent atoms. Molecular-level thermodynamic, kinetic, and structural data extracted from the resulting trajectories provide valuable information for the understanding, engineering, and design of biological and molecular materials. The cost of simulating many-body atomic systems makes simulations of large molecules prohibitively expensive, and the high-dimensionality of the resulting trajectories presents a challenge for analysis. Driven by advances in algorithms, hardware, and data availability, there has been a flare of interest in recent years in the applications of machine learning – especially deep learning – to molecular simulation. Furthermore, these techniques have demonstrated great power and flexibility in both extracting mechanistic understanding of the important nonlinear collective variables governing the dynamics of a molecular system, and in furnishing good low-dimensional system representations with which to perform enhanced sampling or develop long-timescale dynamical models. It is the purpose of this article to introduce the key machine learning approaches, describe how they are married with statistical mechanical theory into domain-specific tools, and detail applications of these approaches in understanding and accelerating biomolecular simulation.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Random Close Packing as a Dynamical Phase Transition

Sphere packing is an ancient problem. The densest packing is known to be a face-centered cubic (FCC) crystal, with space-filling fraction Φ FCC = π / √ 18 ≈ 0.74 . The densest “random packing,” random close packing (RCP), is yet ill defined, although many experiments and simulations agree on a value Φ RCP ≈ 0.64 . We introduce a simple absorbing-state model, biased random organization (BRO), which exhibits a Manna class dynamical phase transition between absorbing and active states that has as its densest critical point Φ cmax ≈ 0.64 ≈ Φ RCP and, like other Manna class models, is hyperuniform at criticality. The configurations we obtain from BRO appear to be structurally identical to RCP configurations from other protocols. This leads us to conjecture that the highest-density absorbing state for an isotropic biased random organization model produces an ensemble of configurations that characterizes the state conventionally known as RCP.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗