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 325 records · Page 18

Edge localized mode characteristics and divertor heat flux during stationary and transient phase for CFETR hybrid scenario

The study of edge localized mode (ELM) behavior, stationary heat flux and transient heat flux during the ELM crash phase is performed for a China Fusion Engineering Test Reactor (CFETR) 1 GW hybrid mode operation scenario ($R$ = 7.2 m, ${B_{\text{T}}}$ = 6.5 T, ${I_{\text{p}}}$ = 13.78 MA). Modeling and simulation start with a scenario obtained by multi-code integrated modeling on the one modeling framework for integrated tasks framework. Linear stability and nonlinear simulations of ELM dynamics are carried out using the BOUT++ six-field reduced magnetohydrodynamic module, which show a much smaller ELM energy loss (Δ ELM ~ 0.13%) compared to that of a Type-I ELM. Parametric analysis of the weak linear growth rate and small ELM energy loss characteristics shed light on physics corresponding to a grassy ELM regime for CFETR 1 GW hybrid scenario. The transient heat flux on the divertor target during this small ELM phase is investigated using BOUT++. Here we found that upstream radial transport in the scrape-off-layer (SOL) induced by small ELMs is weak, which keeps it in the drift-dominated region. However, the heat flux width is still broadened to ${\lambda _{\text{q}}}$ = 4.64 mm by the increase of separatrix temperature during ELM nonlinear evolution. The impact of transient peak heat load and ELM energy fluence on tungsten melting and net erosion rate of divertor target is evaluated for the first-time using physics-based transport that connects the pedestal with the SOL. Energy fluence caused by a single ELM pulse is below the tungsten melting limit, while tungsten erosion would exceed the material requirements. We conclude that external mitigation methods, such as divertor detachment and advanced divertor geometry are likely needed for the steady state operation of CFETR.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗

Dynamic Catalysis Fundamentals: I. Fast calculation of limit cycles in dynamic catalysis

Dynamic catalysis—the forced oscillation of catalytic reaction coordinate potential energy surfaces (PES)—has recently emerged as a promising method for the acceleration of heterogeneously-catalyzed reactions. Theoretical study of enhancement of rates and supra-equilibrium product yield via dynamic catalysis has, to-date, been severely limited by onerous computational demands of forward integration of stiff, coupled ordinary differential equations (ODEs) that are necessary to quantitatively describe periodic cycling between PESs. Here, we establish a new approach that reduces, by ≳108×, the computational cost of finding the time-averaged rate at dynamic steady state (i.e. the limit cycle for linear and nonlinear systems of kinetic equations). Our developments are motivated by and conceived from physical and mathematical insight drawn from examination of a simple, didactic case study for which closed-form solutions of rate enhancement are derived in explicit terms of periods of oscillation and elementary step rate constants. Generalization of such closed-form solutions to more complex catalytic systems is achieved by introducing a periodic boundary condition requiring the dynamic steady state solution to have the same periodicity as the kinetic oscillations and solving the corresponding differential equations by linear algebra or Newton-Raphson-based approaches. The methodology is well-suited to extension to non-linear systems for which we detail the potential for multiple solutions or solutions with different periodicities. For linear and non-linear systems alike, the acute decrement in computational expense enables rapid optimization of oscillation waveforms and, consequently, accelerates understanding of the key catalyst properties that enable maximization of reaction rates, conversions, and selectivities during dynamic catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Domain wall motion across microstructural features in polycrystalline ferroelectric films

This paper describes the effect of microstructural features such as grain boundaries and triple points on the pinning of domain wall motion in perovskite Pb(Zr 0.52 Ti 0.48 )O 3 (PZT) films on polycrystalline SrRuO 3 /SrTiO 3 substrates. Here, spatial variability in the collective domain wall dynamics was assessed using non-linearity mapping via Band Excitation Piezoresponse Force Microscopy (BE-PFM). Collocating the non-linearity maps with triple point locations (as visualized by EBSD) allowed for exploration of the effects that local microstructure (e.g., grain boundary) have on domain wall motion. It was found that the extrinsic behavior varied with both the misorientation angle and the proximity to the grain boundary. The width of influence of individual grain boundaries on the motion of domain walls was a function of the character of the grain boundary; random grain boundaries exhibit deeper minima in α d /d 33,initial and larger widths of influence (up to 905 nm) compared to coincident site lattice (CSL) boundaries (up to 572 nm). Additionally, triple points containing larger numbers of random boundaries exhibited non-Rayleigh behavior to greater distances, suggesting that the triple point provides either a deep potential minimum or a region where domain wall motion is unfavorable.

36 MATERIALS SCIENCE↗

Solvent-Treated Zirconium-Based Nanoporous UiO-66 Metal–Organic Frameworks for Enhanced CO 2 Capture

Metal–organic framework (MOF)-derived nanoscale porous materials are widely deployed in carbon capture, but the CO 2 capacity is varied by different post-treatments, and its mechanism is still unclear. In this work, we prepared UiO-66 and treated it with methanol solvent and thermal activation approaches, which showed ~3 times enhanced CO 2 capacities from 15.1 to 45.0 mg/g and excellent recyclability of 10 cycles. The methanol treatment efficiently removes the residual guest molecules, including N,N-dimethylformamide, dangling organic linkers, and their derivatives, in the micropores (~0.8 nm) of UiO-66 and improves the surface area, pore volume, and void fraction to enhance the CO 2 capacity close to the ideal UiO-66 materials. The molecular dynamics simulation also proved a good linear relationship between the surface areas, void fraction, and CO 2 capacity. This work provides a deep understanding of the MOF’s activation mechanism and its applications in CO 2 capture.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Pressure-tuned quantum criticality in the large-D antiferromagnet DTN

Strongly correlated spin systems can be driven to quantum critical points via various routes. In particular, gapped quantum antiferromagnets can undergo phase transitions into a magnetically ordered state with applied pressure or magnetic field, acting as tuning parameters. These transitions are characterized by z = 1 or z = 2 dynamical critical exponents, determined by the linear and quadratic low-energy dispersion of spin excitations, respectively. Employing high-frequency susceptibility and ultrasound techniques, we demonstrate that the tetragonal easy-plane quantum antiferromagnet NiCl 2 · 4SC(NH 2 ) 2 (aka DTN) undergoes a spin-gap closure transition at about 4.2 kbar, resulting in a pressure-induced magnetic ordering. The studies are complemented by high-pressure-electron spin-resonance measurements confirming the proposed scenario. Powder neutron diffraction measurements revealed that no lattice distortion occurs at this pressure and the high spin symmetry is preserved, establishing DTN as a perfect platform to investigate z = 1 quantum critical phenomena. The experimental observations are supported by DMRG calculations, allowing us to quantitatively describe the pressure-driven evolution of critical fields and spin-Hamiltonian parameters in DTN.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Complex motion of steerable vesicular robots filled with active colloidal rods

Abstract While the collective motion of active particles has been studied extensively, effective strategies to navigate particle swarms without external guidance remain elusive. We introduce a method to control the trajectories of two-dimensional swarms of active rod-like particles by confining the particles to rigid bounding membranes (vesicles) with non-uniform curvature. We show that the propelling agents spontaneously form clusters at the membrane wall and collectively propel the vesicle, turning it into an active superstructure. To further guide the motion of the superstructure, we add discontinuous features to the rigid membrane boundary in the form of a kinked tip, which acts as a steering component to direct the motion of the vesicle. We report that the system’s geometrical and material properties, such as the aspect ratio and Péclet number of the active rods as well as the kink angle and flexibility of the membrane, determine the stacking of active particles close to the kinked confinement and induce a diverse set of dynamical behaviors of the superstructure, including linear and circular motion both in the direction of, and opposite to, the kink. From a systematic study of these various behaviors, we design vesicles with switchable and reversible locomotions by tuning the confinement parameters. The observed phenomena suggest a promising mechanism for particle transportation and could be used as a basic element to navigate active matter through complex and tortuous environments.

42 ENGINEERING↗

Hybrid particle-spectral method for kinetic plasma simulations

A hybrid model for numerical solutions of the Vlasov–Poisson equations is presented, which blends spectral and particle approaches. The model splits the distribution function for plasma species into both spectral and particle representations in the velocity space to combine the advantages of each approach. The spectral representation leverages asymmetrically weighted Hermite basis, whereas the particle representation leverages the particle-in-cell method. Configuration phase space is decomposed with the Fourier method, which is well suited for periodic problems. We derive conservation equations for mass, momentum, and energy for the proposed combined method. It is shown that the coupling error between the two methods is absent in the semi-discrete setting (not taking into account time discretization). Finally, numerical test cases are presented simulating a weak electron beam interaction with plasma, leading to beam–plasma instability. The initially localized electron beam evolved into a highly non-equilibrium distribution function in the velocity space. A small growth rate and the resonance nature of instability make it difficult to obtain accurate solutions for purely particle methods due to noise, which falls as ∼1/Np with a number of particles. At the same time, purely spectral methods may require a large number of modes to capture the highly non-equilibrium state of the evolved beam. We show that the hybrid method is well suited for such problems: it reproduces the linear stage as well as nonlinear dynamics with sufficient accuracy using a highly non-equilibrium distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A novel method for thermal noise reduction, enabling measurements of broadband, low-amplitude electron temperature fluctuations using individual radiometer channels

A new analysis method has been developed for measurements of broadband, low-amplitude turbulent electron temperature fluctuations in fusion plasmas using individual radiometer channels of a correlation electron cyclotron emission diagnostic. This method takes advantage of differences in the correlation time of thermal noise compared to the correlation time of plasma fluctuations in fusion reactors. The validation of this single-channel method is demonstrated using comparisons with the standard dual-channel radiometer spectral decorrelation method for measurements of turbulent electron temperature fluctuations in the core and edge of low confinement (L), improved confinement (I), and high confinement (H)-mode plasmas at the ASDEX Upgrade tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NeuroFEM

SAND2025-00525O NeuroFEM is a software tool that demonstrates a neuromorphic algorithm for solving finite element problems. It sets up a 2D finite element problem for the Poisson equation on a disk, constructs synaptic matrices, and simulates neural dynamics to solve the resulting sparse linear system. The software illustrates how the algorithm converges to the solution and plots the results, showcasing a neuromorphic counterpart to traditional methods like Conjugate Gradient or GMRES. This tool is designed to highlight the potential of neuromorphic algorithms for solving sparse linear systems, which are prevalent in various computational applications. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

SciDAC↗

Strengthening the Security of Operational Technology: Understanding Contemporary Bill of Materials

The evolution of cyber-physical infrastructure has made its security more challenging. The last few years have witnessed a convergence of hardware and software segments in various domains, including operational technology (OT) which is responsible for carrying out critical tasks such as monitoring and controlling power grids, nuclear plants, transportation, and emergency services. Both hardware and software encapsulate numerous open source and proprietary subcomponents, making it crucial for end-users to understand the composition of the products they are using. For example, wind turbines incorporate thousands of lines of code (software) used for the turbine's design, planning, operation, and analytics in addition to the numerous hardware subcomponents that construct it. Due to the highly complex nature of software and hardware, knowledge of the components and subcomponents is required to mitigate cyber vulnerabilities and defend against cyberattacks. There has also been a transformation from a traditional linear supply chain into a global, dynamic, diverse, and interconnected system. The digitization of the supply chain makes it easier to find and exploit vulnerabilities. Critical infrastructures (e.g., power grids, oil, natural gas, water, and wastewater) rely on OT to function, and if the OT is compromised, equipment damage and potential interruption of services could result. A significant security measure to protect OT systems from disruption is to develop a supply chain bill of materials (BoM) corresponding to the software and hardware used in OT, along with attestations amongst vendors and asset owners. A supply chain BoM is a proactive way to understand the inherent vulnerabilities in the system and mitigate them in advance of being exploited. BoMs bolster the trust placed in the digital infrastructure and enhance software supply chain security by sustaining the management of component obsolescence and compliance, along with the seclusion of unsafe segments of a specific product. Adopting BoM tools is becoming increasingly important across various government sectors, as evidenced by the recent U.S. executive order on cybersecurity (NIST 2021). This paper aims to classify BoMs based on structure, functionality, component type, and architecture. The work also discusses case studies to further highlight the benefits of BoMs. In addition, it identifies missing pieces in existing BoM implementations so that future research may identify bounds on where it could expect to make improvements and directly enable researchers to identify promising areas for exploration. Further, the authors provide valuable recommendations to tool developers, researchers, and standardizing organizations (policymakers), additionally benefitting critical infrastructure owners and government executives. This aids in paving a path for future work, thereby, providing suggestions to determine a tool for consumers that best suit their needs.

97 MATHEMATICS AND COMPUTING↗

Exploring the Dynamics of Transverse Inter-planar Coupling in the SC Section of the PIP-II Linac

This study investigates the crucial role that an accurate understanding of inter-planar coupling in the transverse plane plays in regulating charged particle dynamics in a high-intensity linear accelerator and minimizing foil/septum impacts during injection from the linac to a ring. We in-depth analyze the emergence and evolution of transverse inter-planar coupling through multiple active lattice elements, taking into account space charge and field nonlinearities in the superconducting section of the PIP-II linac. The article compares various analytical, numerical, and experimental techniques for measuring transverse coupling using beam and lattice matrices and provides insight into effective strategies for its mitigation prior to ring injection.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Hybrid particle-spectral method for kinetic plasma simulations

A hybrid model for numerical solutions of the Vlasov–Poisson equations is presented, which blends spectral and particle approaches. The model splits the distribution function for plasma species into both spectral and particle representations in the velocity space to combine the advantages of each approach. The spectral representation leverages asymmetrically weighted Hermite basis, whereas the particle representation leverages the particle-in-cell method. Configuration phase space is decomposed with the Fourier method, which is well suited for periodic problems. We derive conservation equations for mass, momentum, and energy for the proposed combined method. It is shown that the coupling error between the two methods is absent in the semi-discrete setting (not taking into account time discretization). Finally, numerical test cases are presented simulating a weak electron beam interaction with plasma, leading to beam–plasma instability. The initially localized electron beam evolved into a highly non-equilibrium distribution function in the velocity space. A small growth rate and the resonance nature of instability make it difficult to obtain accurate solutions for purely particle methods due to noise, which falls as ~1/$\sqrt{N_p}$ with a number of particles. At the same time, purely spectral methods may require a large number of modes to capture the highly nonequilibrium state of the evolved beam. We show that the hybrid method is well suited for such problems: it reproduces the linear stage as well as nonlinear dynamics with sufficient accuracy using a highly non-equilibrium distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Cross-beam energy transfer in direct-drive ICF. I. Nonlinear and kinetic effects

Results are presented from a series of simulations examining the susceptibility of the cross-beam energy transfer (CBET) instability to nonlinear processes in the context of direct-drive inertial confinement fusion experiments on the OMEGA laser facility. These form the basis for the second paper of this series [A. G. Seaton, L. Yin, R. Follett, B. J. Albright, and A. Le, “Cross-beam energy transfer in direct-drive ICF. II. Theory and simulation of mitigation through increased laser bandwidth,” Phys. Plasmas 29, 042707 (2022)], where we examine the efficacy of increases in laser bandwidth at suppressing CBET. We choose laser and plasma conditions for the simulations that are favorable to CBET and promote nonlinearity. Through a comparison of outputs from the particle-in-cell code vector particle in cell (VPIC) and the linearized fluid code laser-plasma simulation environment (LPSE), a series of nonlinear effects have been identified in the kinetic simulations that include particle trapping, the two-ion wave decay, and ion-acoustic wave self-focusing. These effects produce time-dependent energy transfer, in contrast to the linearized fluid simulations in which a steady state is reached after an initial transient. Ion trapping is shown to allow for increased energy transfer relative to fluid simulations, with the remaining nonlinear processes acting to reduce the energy transfer. Nonlinear dynamics is contrasted for low- and high-intensity beams as well as between speckled and planar beams. For the parameters under consideration, beam profile has a significant effect on nonlinear dynamics, though the greatest sensitivity is to beam intensity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Data-Driven Modeling and Correction of Vehicle Dynamics

We develop a data-driven framework for learning and correcting nonautonomous vehicle dynamics. Physics-based vehicle models are often simplified for tractability and therefore exhibit inherent model-form uncertainty, motivating the need for data-driven correction. Moreover, nonautonomous dynamics are governed by time-dependent control inputs, which pose challenges in learning predictive models directly from temporal snapshot data. To address these, we reformulate the vehicle dynamics via a local parameterization of the time-dependent inputs, yielding a modified system composed ofa sequence of local parametric dynamical systems. Here, we approximate these parametric systems using two complementary approaches. First, we employ the dimension reduction and interpolation in parameter space (DRIPS) methodology to construct efficient linear surrogate models, equipped with lifted observable spaces and manifold-based operator interpolation. This enables data-efficient learning of vehicle models whose dynamics admit accurate linear representations in the lifted spaces. Second, for more strongly nonlinear systems, we employ flow map learning (FML), a deep neural network (DNN) approach that approximates the parametric evolution map without requiring special treatment of nonlinearities. We further extend FML with a transfer-learning-based model correction procedure, enabling the correction of misspecified prior models using only a sparse set of high-fidelity or experimental measurements, without assuming a prescribed form for the correction term. Through a suite of numerical experiments on unicycle, simplified bicycle, and slip-based bicycle models, we demonstrate that DRIPS offers robust and highly data-efficient learning of nonautonomous vehicle dynamics, while FML provides expressive nonlinear modeling and effective correction of model-form errors under severe data scarcity.

data-driven modeling↗

Numerical modeling of a proton spin-flipping system in the spin transparency mode at an integer spin resonance in JINR's Nuclotron

In this paper we propose a lattice insertion for the Nuclotron ring called a “spin navigator” that can adjust any direction of the proton polarization in the orbital plane using weak solenoids. The polarization control is realized in the spin transparency mode at the energy of 108 MeV, which corresponds to the integer spin resonance γ G = 2. The requirements on the navigator solenoid fields are specified considering the criteria for stability of the spin motion during any manipulation of the polarization direction in an experiment. Additionally, this paper presents the results of numerical modeling of the proton spin dynamics in the Nuclotron ring operated in the spin transparency mode. The verified spin navigator is aimed at an experimental study of a spin-flipping system using the Nuclotron ring. The results are relevant to the NICA (JINR), EIC (BNL) and COSY (FZJ) facilities where the spin transparency mode can be applied for polarization control.

47 OTHER INSTRUMENTATION↗

Real-Time Modal Analysis of Electric Power Grids - The Need for Dynamic State Estimation

We articulate the reason why dynamic state estimation is needed to push the boundaries of the modal analysis of electric power grids in real-time operation. Then, we demonstrate how to unravel linear and nonlinear modes by using the extended dynamic mode decomposition along with estimates of the synchronous generators’ rotor angles and rotor speed deviations from nominal speed. The estimated modes are associated with electromechanical oscillations that take place continuously in electric power grids because of imbalances between power generation and demand. The numerical simulations are performed on a synthetic, albeit realistic, 2,000-bus network that was designed to resemble the electric power grid of Texas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗