Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Linear 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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Control simulations of many-body quantum systems by a synergism of discrete real-time learning and optimal control theory

We present a self-consistent algorithm for optimal control simulations of many-body quantum systems. The algorithm features a two-step synergism that combines discrete real-time machine learning (DRTL) with Quantum Optimal Control Theory (QOCT) using the time-dependent Schrödinger equation. Specifically, in step (1), DRTL is employed to identify a compact working space (i.e., the important portion of the Hilbert space) for the time evolution of the many-body quantum system in the presence of a control field (i.e., the initial or previously updated field), and in step (2), QOCT utilizes the DRTL-determined working space to find a newly updated control field for a chosen objective. Steps 1 and 2 are iterated until a self-consistent control objective value is reached such that the resulting optimal control field yields the same targeted objective value when the corresponding working space is systematically enlarged. Furthermore, to demonstrate this two-step self-consistent DRTL-QOCT synergistic algorithm, we perform optimal control simulations of strongly interacting 1D as well as 2D Heisenberg spin systems. In both scenarios, only a single spin (at the left end site for 1D and the upper left corner site for 2D) is driven by the time-dependent control fields to create an excitation at the opposite site as the target. It is found that, starting from all spin-down zero excitation states, the synergistic method is able to identify working spaces and convergence of the desired controlled dynamics with just a few iterations of the overall algorithm. In the cases studied, the dimensionality of the working space scales only quasi-linearly with the number of spins.

Artificial neural networks↗
DOE OSTI2021

Edge-localized-mode simulation in CFETR steady-state scenario

We report the EPED1 model and self-consistent core-pedestal coupling in integrated modeling are used to design the pedestal structure of the China Fusion Engineering Testing Reactor (CFETR) steady-state scenario. The key parameters, such as β p and q 95 , are based on the grassy edge-localized-mode (ELM) experimental database. In this work, we use the BOUT++ six-field two-fluid code to simulate the onset of the ELM in the CFETR steady-state scenario. The ELM size is around 0.2% in nonlinear simulations, which is in the experimental range of the grassy ELM discharges, 0.1%–1% observed in multiple tokamak devices. Linear and nonlinear simulations show that the dominant high-n ballooning modes peak around n = 40. Compared to type-I ELM crashing dynamics, grassy ELM crashing has a smaller initial crash and is then followed by three phases of turbulence spreading, which are dominated by multi-modes, a high-n mode of n = 45 and low-n mode of n = 5, respectively. In contras to type-I ELM, the perturbation of the high-n mode has a narrow width around ψ = 0.95, and magnetic island formation and reconnection occur only beyond ψ = 0.95, leading to a small initial crash. Mode–mode interaction in the multi-mode coexistence stage stops the growth of individual modes and reduces the transport of particles and heat, and these are the two reasons why the ELM size is small. In–out asymmetry of transient heat flux with a ratio of E out /E in = 3.5 is found during grassy ELM crash. The rise and delay times of the heat flux match the calculation from the free-streaming model. To evaluate the erosion of the divertor target, the energy fluence at the outer divertor target is calculated, which is 0.029 MJ m -2 , 5.5 times smaller than the tungsten melting limit 0.16 MJ m -2 . The calculated energy fluency still follows the experimental scaling law from type-I ELM experiments. The fluctuation eddies in the toroidal direction show a filament structure at the outer mid-plane. Parallel heat flux patterns with a toroidal mode number n = 10 are found at the outer divertor with an amplitude of 680 MW m -2 .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Photoabsorption, photoionization, and Auger processes at the carbon K edge in C H 3 I

The dynamics of photoabsorption, photoionization, and the associated Auger decay have been investigated at the carbon K edge in methyl iodide (CH 3 I) using linearly polarized synchrotron radiation. Ion yield measurements were used to investigate transitions in the pre-edge region due to excitations into either unoccupied valence or Rydberg states. The assignment of these transitions was achieved through comparison with theoretical x-ray absorption spectra calculated using time-dependent density functional theory, within the Tamm-Dancoff approximation. Several of the Rydberg states belonging to series converging onto the C1s ionization limit exhibit significant vibrational structure that is also interpreted using theoretical calculations. The C1s in CH 3 I photoelectron spectrum was measured, and the observed vibrational structure was assigned with the aid of theoretical predictions. Polarization dependent, resonantly excited, valence shell photoelectron spectra were recorded at photon energies coinciding with the C1s → σ*, C1s → 6sa1 and C1s → 6pe transitions in CH 3 I , thereby allowing photoelectron angular distributions to be determined. The nonresonantly excited C(KVV) Auger electron spectrum was measured and some of the features observed at high kinetic energies were attributed to transitions into valence orbitals possessing significant iodine character. The contributions of participator and spectator Auger decay to the resonantly excited photoelectron spectra have been assessed. The influence of participator decay appears minor whereas spectator decay results in the enhanced population of satellite states.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Quantum stabilization of unexpected ordered phases on the honeycomb lattice

In this work, I discuss and showcase the utility of the method called minimally-augmented spin-wave theory (MAGSWT) as a relatively simple semi-analytical approach to the phase diagrams of quantum magnets. It complements numerical methods by providing physical insight into which states are competitive and by yielding approximate phase boundaries that agree well with much more numerically intensive calculations. This approach enabled me to construct the faithful phase diagrams of two paradigmatic honeycomb-lattice models in the quantum S=1/2 limit, the J 1 –J 3 FM-AF and J 1 –J 2 AF models, for the collinear quantum phases that replace or extend the classical ones. The results are in good qualitative and semi-quantitative agreement with state-of-the-art numerical studies for these models, correctly capturing the emergence of unexpected quantum phases and the suppression of classically favored spiral orders by quantum fluctuations. This study provides a much-needed important guidance to the ongoing theoretical and experimental searches of the unconventional quantum states. This work will be of significant and timely interest to both theorists and experimentalists in the field of quantum magnetism, broadly defined.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

A Generalized Framework for Service Restoration in a Resilient Power Distribution System

An electric power grid is one of the complex infrastructures, and because of its complex nature, there is simply no way that outages can be completely avoided. Thus, a modern society that depends on reliable electric supply requires a resilient electric system that can recover from disruptions while integrating emerging smart grid technologies. Here, this article presents a novel approach for service restoration in a modern power distribution system with controllable switches and distributed generation (DG) resources for any kind of outage. The proposed framework supports both the traditional service restoration using feeder reconfiguration and the grid-forming DG-assisted intentional islands that are dynamically sized using algorithms based on the fault scenario, available resources, and priority of loads. The problem is formulated as a mixed-integer linear program that incorporates critical system connectivity and operating constraints. Simulations are performed to demonstrate the effectiveness of the proposed approach using a large-scale four-feeder 1069-bus three-phase unbalanced distribution test system. It is demonstrated that the framework is effective in utilizing all available resources in quickly restoring the power supply to improve resiliency during extreme events and is scalable for a large-scale unbalanced power distribution system.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Distributed State Estimation Over Time-Varying Graphs: Exploiting the Age-of-Information

Here, we study the problem of designing a distributed observer for an LTI system over a time-varying communication graph. The limited existing work on this topic imposes various restrictions either on the observation model or on the sequence of communication graphs. In contrast, we propose a single-time-scale distributed observer that works under mild assumptions. Specifically, our communication model only requires strong-connectivity to be preserved over non-overlapping, contiguous intervals that are even allowed to grow unbounded over time. We show that under suitable conditions that bound the growth of such intervals, joint observability is sufficient to track the state of any discrete-time LTI system exponentially fast, at any desired rate. We also develop a variant of our algorithm that is provably robust to worst-case adversarial attacks, provided the sequence of graphs is sufficiently connected over time. The key to our approach is the notion of a "freshness-index" that keeps track of the age-of-information being diffused across the network. Such indices enable nodes to reject stale estimates of the state, and, in turn, contribute to stability of the error dynamics.

42 ENGINEERING↗
DOE OSTI2023

The Optical Signatures of Stochastic Processes in Many-Body Exciton Scattering

We review our recent quantum stochastic model for spectroscopic lineshapes in the presence of a coevolving and nonstationary background population of excitations. Starting from a field theory description for interacting bosonic excitons, we derive a reduced model whereby optical excitons are coupled to an incoherent background via scattering as mediated by their screened Coulomb coupling. The Heisenberg equations of motion for the optical excitons are then driven by an auxiliary stochastic population variable, which we take to be the solution of an Ornstein–Uhlenbeck process. Here, we present an overview of the theoretical techniques we have developed as applied to predicting coherent nonlinear spectroscopic signals. We show how direct (Coulomb) and exchange coupling to the bath give rise to distinct spectral signatures and discuss mathematical limits on inverting spectral signatures to extract the background density of states.

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

The Radiative Magneto-Hydrodynamics of Colliding Flow: Instabilities, Reconnection and Exoplanet Atmosphere Connections

Hypersonic flows of magnetized plasmas naturally produce strong shocks when obstructions or self‐interactions (flow collisions) occur. Such complex flows occur in a wide variety of HEDP settings such as Z‐pinches and laser driven implosion experiments that include magnetic fields via explosive Helmholtz coils (MIFEDS) When thermal energy behind the shocks can be lost to optically thin radiation, the flow dynamics becomes even more complex, admitting a range of new instabilities such as the Nonlinear Thin Shell Instability (NTSI). Because of their inherent non‐ linearity, multi‐dimensionality and time‐dependence (i.e. turbulence), the ability to articulate the underlying physics of colliding hypersonic magnetized flows represents a significant frontier in HEDP studies. The addition of radiative cooling to these flows propels the study even further into the domain of extreme plasma physics. Exploration of these forms of HED plasmas is of interest to basic plasma physics as well as of interest for their application to Inertial Confinement Fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Simulation, Structural, Thermal and Mechanical Properties of the FeTiTaVW High Entropy Alloy

Developing new materials to be applied in extreme environments is an opportunity and a challenge for the future. High entropy alloys are new materials that seem promising approaches to work in nuclear fusion reactors. In this work, FeTaTiVW high entropy alloys were developed and characterized with Molecular Dynamic and Hybrid Molecular Dynamic Monte Carlo simulations. The simulation results show that phase separation originates a lower potential energy per atom and a high level of segregation compared to those of a uniform solid solution. Moreover, the experimental diffractogram of the milled powder shows the formation of a body-centred cubic-type structure and the presence of TiO2. In addition, the microstructure of the consolidated material evidenced three phases: W-rich, Ti-rich, and a phase with all the elements. This phase separation observed in the microstructure agrees with the Hybrid Molecular Dynamic Monte Carlo simulation. Moreover, the consolidated material’s thermal conductivity and specific heat are almost constant from 25 °C to 1000 °C, and linear expansion increases with increasing temperature. On the other hand, specific heat and thermal expansion values are in between CuCrZr and W values (materials chosen for the reactor walls). The FeTaTiVW high entropy alloy evidences a ductile behaviour at 1000 °C. Therefore, the promising thermal properties of this system can be attributed to the multiple phases and systems with different compositions of the same elements, which is exciting for future developments.

Martins, Ricardo (ORCID:0000000177685704)↗
DOE OSTI2022

Perturbative theoretical model of electronic transient circular dichroism spectroscopy of molecular aggregates

A chiral analog of transient absorption spectroscopy, transient circular dichroism (TCD) spectroscopy is an emerging time-resolved method. Both spectroscopic methods can probe the electronic transitions of a sample, and TCD is additionally sensitive to the dynamic aspects of chirality, such as those induced by molecular excitons. Here, we develop a theoretical description of TCD for electronic multi-level models in which the pump pulse is linearly polarized and probe pulse is alternately left- and right-circularly polarized. We derive effective response functions analogous to those often used to describe other four-wave mixing methods and then simulate and analyze TCD spectra for three representative multi-level electronic model systems. We elaborate on the presence and detection of the spectral signatures of electronic coherences.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Coupling flux balance analysis with reactive transport modeling through machine learning for rapid and stable simulation of microbial metabolic switching

Integrating genome-scale metabolic networks with reactive transport models (RTMs) provides a detailed description of the dynamic changes in microbial growth and metabolism. Despite promising demonstrations in the past, computational inefficiency has been pointed out as a critical issue to overcome because it requires repeated application of linear programming (LP) to obtain flux balance analysis (FBA) solutions in every time step and spatial grid. To address this challenge, we propose a new simulation method where we train and validate artificial neural networks (ANNs) using randomly sampled FBA solutions and incorporate the resulting surrogate FBA model (represented as algebraic equations) into RTMs as source/sink terms. We demonstrate the efficiency of our method via a case study of Shewanella oneidensis MR-1. During aerobic growth on lactate, S. oneidensis produces metabolic byproducts (such as pyruvate and acetate), which are subsequently consumed as alternative carbon sources when the preferred nutrients are depleted. To effectively simulate these complex dynamics, we used a cybernetic approach that models metabolic switches as the outcome of dynamic competition among multiple growth options. In both zero-dimensional batch and one-dimensional column configurations, the ANN-based surrogate models achieved substantial reduction of computational time by several orders of magnitude compared to the original LP-based FBA models. Moreover, the ANN models produced robust solutions without any special measures to prevent numerical instability. These developments significantly promote our ability to utilize genome-scale networks in complex, multi-physics, and multi-dimensional ecosystem modeling.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

A smeared crack modeling framework accommodating multi-directional fracture at finite strains

A generic smeared crack modeling framework predicated on the deformation gradient decomposition (DGD) approach is proposed for use in dynamic fracture problems at finite strains, accommodating failure along multiple mutually orthogonal fracture planes embedded within an independently defined bulk material model. Within this constitutive framework, the traction equilibrium conditions imposed at each failure surface are used to determine the associated crack displacements stored as internal state variables. In general, the enforcement of interfacial equilibrium entails the implicit solution of a non-linear system of equations within the constitutive update procedure. However, if inertial effects arising due to the relative motion of the fractured material are incorporated within the model, the traction equilibrium conditions are shown to give rise to corresponding dynamic equations of motion governing the time-evolution of the crack opening displacements. For dynamic problems, an explicit time-integration procedure is devised to efficiently update the material state, subject to a set of internal frictionless contact constraints to prevent material inversion. Finally, the efficacy of the proposed modeling framework is investigated through several benchmark dynamic fracture problems run within the explicit finite element code DYNA3D.

42 ENGINEERING↗
DOE OSTI2023

Softening by charging: how collective modes of ionic association in concentrated redoxmer/electrolyte solutions define the structural and dynamic properties in different states of charge

Understanding the physical and chemical processes occurring in concentrated electrolyte solutions is required to achieve redox flow batteries with high energy density. Highly concentrated electrolyte solutions are often studied in which collective crowded interactions between molecules and ions become predominant. In this work, experimental and computational methods were used to examine non-aqueous electrolyte solutions in two different states of charge as a function of redoxmer concentration. As the latter increases and the ionic association strengthens, the electric conductivity passes through a maximum and the solution increasingly gels, which is seen through a rapid non-linear increase in viscosity. We establish that the structural rigidity of ionic networks is closely connected with this loss of fluidity and show that charging generally yields softer ionic assemblies with weaker attractive forces and improved dynamical properties.

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

Dynamic wavelength control of laser pulse profiles at picosecond to nanosecond timescales

We report on a novel combined laser pulse shaping and dynamic wavelength encoding capability based on a simple architecture implementing direct space to time mapping. There are several potential applications that can be enabled by the ability to control the instantaneous intensity or wavelength of an optical waveform on a picosecond-to-nanosecond timescale. To our knowledge, no known methods can access this temporal regime with a practical architecture. Here, we demonstrate an extension of the Space–Time Induced Linearly Encoded Transcription for Temporal Optimization (STILETTO) technique that can generate optical waveforms with a programmable instantaneous wavelength vs. time. We experimentally demonstrate the technique by generating self-gated spectrograms and show that it can encode dynamic wavelength vs time profiles at timescales not achievable by any other known method.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

Dynamics of McMillan mappings III. Symmetric map with mixed nonlinearity

This article extends the study of the dynamical properties of the symmetric McMillan map, emphasizing its utility in understanding and modeling complex nonlinear systems. Although the map features six parameters, we demonstrate that only two are irreducible: the linearized rotation number at the fixed point and a nonlinear parameter representing the ratio of terms in the biquadratic invariant. Through a detailed analysis, we classify regimes of stable motion, provide exact solutions to the mapping equations, and derive a canonical set of action-angle variables, offering analytical expressions for the rotation number and nonlinear tune shift. We further establish connections between general standard-form mappings and the symmetric McMillan map, using the area-preserving Hénon map and accelerator lattices with thin sextupole magnet as representative case studies. Our results show that, despite being a second-order approximation, the symmetric McMillan map provides a highly accurate depiction of dynamics across a wide range of system parameters, demonstrating its practical relevance in both theoretical and applied contexts.

43 PARTICLE ACCELERATORS↗
DOE OSTI2019

First principles simulation of ultracold ion crystals in a Penning trap with Doppler cooling and a rotating wall potential

A direct numerical simulation of many interacting ions in a Penning trap with a rotating wall is presented. The ion dynamics is modeled classically. Both axial and planar Doppler laser cooling processes are modeled using stochastic momentum impulses based on two-level atomic scattering rates. The plasmas being modeled are ultracold two-dimensional crystals made up of hundreds of ions. We compare Doppler cooled results directly to a previous linear eigenmodes analysis. Agreement in both frequency and mode structure is obtained. Additionally, when Doppler laser cooling is applied, the laser cooled steady state plasma axial temperature agrees with the Doppler cooling limit. Numerical simulations using the approach described and benchmarked here will provide insights into the dynamics of large trapped-ion crystals, improving their performance as a platform for quantum simulation and sensing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Impact of pedestal density gradient and collisionality on ELM dynamics

BOUT++ turbulence simulations are conducted to capture the underlying physics of small ELM characteristics achieved by increasing separatrix density via controlling strike points from vertical to horizontal divertor plates for three EAST discharges. BOUT++ linear simulations show that the most unstable modes change from high-n ideal ballooning modes to intermediate-n peeling–ballooning modes and eventually to peeling–ballooning stable plasmas in the pedestal. Nonlinear simulations show that the fluctuation is saturated at a high level for the lowest separatrix density. The ELM size decreases with increasing separatrix density, until the fraction of this energy lost during the ELM crash becomes less than 1% of the pedestal stored energy, leading to small ELMs. Simulations indicate that small ELMs can be triggered either by the marginally peeling–ballooning instability near the peak pressure gradient position inside the pedestal or by a local instability in the pedestal foot with a larger separatrix density gradient. The pedestal collisionality scan for type-I ELMs with steep pedestal density gradient shows that both linear growth rate and ELM size decrease with increasing collisionality. On the contrary, the pedestal collisionality and pedestal density width scan with a weak pedestal density gradient indicate small ELMs can either be triggered by a high-n ballooning mode or by a low-n peeling mode in a low collisionality region 0.04–0.1. We report the simulations indicate the weaker the linear unstable modes near marginal stability with small linear growth rate, the lower nonlinearly saturated fluctuation intensity and the smaller turbulence spreading from the linear unstable zone to stable zone in the nonlinear saturation phase, leading to small ELMs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗