Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear dynamical system”

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 721 records · Page 40

Nonlinear dynamics of topological Dirac fermions in 2D spin-orbit coupled materials

Abstract The graphene family materials are two-dimensional staggered monolayers with a gapped energy band structure due to intrinsic spin-orbit coupling. The mass gaps in these materials can be manipulated on-demand via biasing with a static electric field, an off-resonance circularly polarized laser, or an exchange interaction field, allowing the monolayer to be driven through a multitude of topological phase transitions. We investigate the dynamics of spin-orbit coupled graphene family materials to unveil topological phase transition fingerprints embedded in the nonlinear regime and show how these signatures manifest in the nonlinear Kerr effect and in third-harmonic generation processes. We show that the resonant nonlinear spectral response of topological fermions can be traced to specific Dirac cones in these materials, enabling characterization of topological invariants in any phase by detecting the cross-polarized component of the electromagnetic field. By shedding light on the unique processes involved in harmonic generation via topological phenomena our findings open an encouraging path towards the development of novel nonlinear systems based on two-dimensional semiconductors of the graphene family.

36 MATERIALS SCIENCE↗

Ocean Energy Systems Wave Energy Modeling Task 10.4: Numerical Modeling of a Fixed Oscillating Water Column

This paper reports on an ongoing international effort to establish guidelines for numerical modeling of wave energy converters, initiated by the International Energy Agency Technology Collaboration Program for Ocean Energy Systems. Initial results for point absorbers were presented in previous work, and here we present results for a breakwater-mounted Oscillating Water Column (OWC) device. The experimental model is at scale 1:4 relative to a full-scale installation in a water depth of 12.8 m. The power-extracting air turbine is modeled by an orifice plate of 1–2% of the internal chamber surface area. Measurements of chamber surface elevation, air flow through the orifice, and pressure difference across the orifice are compared with numerical calculations using both weakly-nonlinear potential flow theory and computational fluid dynamics. Both compressible- and incompressible-flow models are considered, and the effects of air compressibility are found to have a significant influence on the motion of the internal chamber surface. Recommendations are made for reducing uncertainties in future experimental campaigns, which are critical to enable firm conclusions to be drawn about the relative accuracy of the numerical models. It is well-known that boundary element method solutions of the linear potential flow problem (e.g., WAMIT) are singular at infinite frequency when panels are placed directly on the free surface. This is problematic for time-domain solutions where the value of the added mass matrix at infinite frequency is critical, especially for OWC chambers, which are modeled by zero-mass elements on the free surface. A straightforward rational procedure is described to replace ad-hoc solutions to this problem that have been proposed in the literature.

16 TIDAL AND WAVE POWER↗

Fracton Hydrodynamics without Time-Reversal Symmetry

Here, we present an effective field theory for the nonlinear fluctuating hydrodynamics of a single conserved charge with or without time-reversal symmetry, based on the Martin-Siggia-Rose formalism. Applying this formalism to fluids with only charge and multipole conservation, and with broken time-reversal symmetry, we predict infinitely many new dynamical universality classes, including some with arbitrarily large upper critical dimensions. Using large scale simulations of classical Markov chains, we find numerical evidence for a breakdown of hydrodynamics in quadrupole-conserving models with broken time-reversal symmetry in one spatial dimension. Our framework can be applied to the hydrodynamics around stationary states of open systems, broadening the applicability of previously developed ideas and methods to a wide range of systems in driven and active matter.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Vortex dynamics induced by scanning SQUID susceptometry

In this work, we measured the local magnetic response of a niobium thin film by applying a millitesla-scale AC magnetic field using a micron-scale field coil and detecting the response with a micron-scale pickup loop in a scanning superconducting quantum interference device (SQUID) susceptometry measurement. Near the film's critical temperature, we observed a steplike nonlinear and dissipative magnetic response due to the dynamics of a small number of vortex-antivortex pairs induced in the film by the local applied AC field. We modeled the dynamics of the measurement using a combined two-dimensional London-Maxwell and time-dependent Ginzburg-Landau approach, allowing us to construct a detailed real-space picture of the vortex motion causing the observed dissipative response. This work pushes scanning SQUID susceptometry of two-dimensional superconductors beyond the regime of linear response and lays the foundation for microscopic studies of vortex dynamics and pinning in superconducting devices and more exotic materials systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Collisional simulations of the modulator section in coherent electron cooling

The first section of any coherent electron cooling (CeC) system is the modulator, where the density of the electron beam is modulated by the copropagating ion beam. This density modulation is a result of Coulomb collisions between the individual particles of the two beams. The pairwise, stochastic part of the interactions impacts the overall performance of the CeC process. We present the first simulations of the density modulations of the electron beams from a collisional picture of the dynamics, considering the proof-of-principle CeC experiments at Brookhaven National Laboratory. These simulations were performed using PHAD, which is the first efficient, large-scale collisional numerical method in beam physics that we have previously developed and benchmarked. Realistic beam distributions and external fields have been optimized to provide strong modulation signals necessary for variations of coherent electron cooling systems. Cooling performance limits and potential collisionless simulation pitfalls are pointed out. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗

Time-dependent saturation and physics-based nonlinear model of cross-beam energy transfer

The nonlinear physics of cross-beam energy transfer (CBET) for multi-speckled laser beams is examined using large-scale particle-in-cell simulations for a range of laser and plasma conditions relevant to indirect-drive inertial confinement fusion (ICF) experiments. The time-dependent growth and saturation of CBET involve complex, nonlinear ion and electron dynamics, including ion trapping-induced enhancement and detuning, ion acoustic wave (IAW) nonlinearity, oblique forward stimulated Raman scattering (FSRS), and backward stimulated Brillouin scattering (BSBS) in a CBET-amplified seed beam. Ion-trapping-induced detuning of CBET is captured in the kinetic linear response by a new δf-Gaussian-mixture algorithm, enabling an accurate characterization of trapping-induced non-Maxwellian distributions. Ion trapping induces nonlinear processes, such as changes to the IAW dispersion and nonlinearities (e.g., bowing and self-focusing), which, together with pump depletion, FSRS, and BSBS, determine the time-dependent nature and level of CBET gain as the system approaches a steady state. Using VPIC simulations at intensities at and above the onset threshold for ion trapping and the insight from the time-dependent saturation analyses, we construct a nonlinear CBET model from local laser and plasma conditions that predicts the CBET gain and the energy deposition into the plasma. This model is intended to provide a more accurate, physics-based description of CBET saturation over a wide range of conditions encountered in ICF hohlraums compared with linear CBET gain models with ad hoc saturation clamps often used in laser ray-based methods in multi-physics codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nonlinear burn control in ITER using adaptive allocation of actuators with uncertain dynamics

Abstract ITER will be the first tokamak to sustain a fusion-producing, or burning, plasma. If the plasma temperature were to inadvertently rise in this burning regime, the positive correlation between temperature and the fusion reaction rate would establish a destabilizing positive feedback loop. Careful regulation of the plasma’s temperature and density, or burn control, is required to prevent these potentially reactor-damaging thermal excursions, neutralize disturbances and improve performance. In this work, a Lyapunov-based burn controller is designed using a full zero-dimensional nonlinear model. An adaptive estimator manages destabilizing uncertainties in the plasma confinement properties and the particle recycling conditions (caused by plasma–wall interactions). The controller regulates the plasma density with requests for deuterium and tritium particle injections. In ITER-like plasmas, the fusion-born alpha particles will primarily heat the plasma electrons, resulting in different electron and ion temperatures in the core. By considering separate response models for the electron and ion energies, the proposed controller can independently regulate the electron and ion temperatures by requesting that different amounts of auxiliary power be delivered to the electrons and ions. These two commands for a specific control effort (electron and ion heating) are sent to an actuator allocation module that optimally maps them to the heating actuators available to ITER: an electron cyclotron heating system (20 MW), an ion cyclotron heating system (20 MW), and two neutral beam injectors (16.5 MW each). Two different actuator allocators are presented in this work. The first actuator allocator finds the optimal mapping by solving a convex quadratic program that includes actuator saturation and rate limits. It is nonadaptive and assumes that the mapping between the commanded control efforts and the allocated actuators (i.e. the effector model) contains no uncertainties. The second actuator allocation module has an adaptive estimator to handle uncertainties in the effector model. This uncertainty includes actuator efficiencies, the fractions of neutral beam heating that are deposited into the plasma electrons and ions, and the tritium concentration of the fueling pellets. Furthermore, the adaptive allocator considers actuator dynamics (actuation lag) that contain uncertainty. This adaptive allocation algorithm is more computationally efficient than the aforementioned nonadaptive allocator because it is computed using dynamic update laws so that finding the solution to a static optimization problem is not required at every time step. A simulation study assesses the performance of the proposed adaptive burn controller augmented with each of the actuator allocation modules.

Physics↗

Identification of the weak-to-strong transition in Alfvénic turbulence from space plasma

Abstract Plasma turbulence is a ubiquitous dynamical process that transfers energy across many spatial and temporal scales in astrophysical and space plasma systems. Although the theory of anisotropic magnetohydrodynamic (MHD) turbulence has successfully described natural phenomena, its core prediction of an Alfvénic transition from weak to strong MHD turbulence when energy cascades from large to small scales has not been observationally confirmed. Here we report evidence for the Alfvénic weak-to-strong transition in small-amplitude, turbulent MHD fluctuations in Earth’s magnetosheath using data from four Cluster spacecraft. Our observations demonstrate the universal existence of strong turbulence accompanied by weak turbulent fluctuations on large scales. Moreover, we find that the nonlinear interactions of MHD turbulence are crucial to the energy cascade, as they broaden the cascade directions and fluctuating frequencies. The observed connection between weak and strong MHD turbulence systems may be present in many astrophysical environments, such as star formation, energetic particle transport, turbulent dynamos, and solar corona or solar wind heating.

79 ASTRONOMY AND ASTROPHYSICS↗

Table-top extreme ultraviolet second harmonic generation

The lack of available table-top extreme ultraviolet (XUV) sources with high enough fluxes and coherence properties has limited the availability of nonlinear XUV and x-ray spectroscopies to free-electron lasers (FELs). Here, we demonstrate second harmonic generation (SHG) on a table-top XUV source by observing SHG near the Ti M 2,3 edge with a high-harmonic seeded soft x-ray laser. Furthermore, this experiment represents the first SHG experiment in the XUV. First-principles electronic structure calculations suggest the surface specificity and separate the observed signal into its resonant and nonresonant contributions. The realization of XUV-SHG on a table-top source opens up more accessible opportunities for the study of element-specific dynamics in multicomponent systems where surface, interfacial, and bulk-phase asymmetries play a driving role.

47 OTHER INSTRUMENTATION↗

Visualizing energy transfer between redox-active colloids

Redox-active colloids (RACs) represent a novel class of energy carriers that exchange electrical energy upon contact. Understanding contact-mediated electron transfer dynamics in RACs offers insights into physical contact events in colloidal suspensions and enables quantification of electrical energy transport in nonconjugated polymers. Redox-based electron transport was directly observed in monolayers of micron-sized RACs containing ethyl-viologen side groups via fluorescence microscopy through an unexpected nonlinear electrofluorochromism that is quantitatively coupled to the redox state of the colloid. Via imaging studies, using this electrofluorochromism, the apparent charge transfer diffusion coefficient D CT of the RAC was easily determined. The visualization of energy transport within suspensions of redox-active colloids was also demonstrated. Our work elucidates fundamental mechanisms of energy transport in colloidal systems, informs the development of next-generation redox flow batteries, and may inspire new designs of smart active soft matter including conductive polymers for applications ranging from electrochemical sensors and organic electronics to colloidal robotics.

Qu, Alan Subing [Univ. of Illinois at Urbana-Champ↗

Identification of the weak-to-strong transition in Alfvénic turbulence from space plasma

Plasma turbulence is a ubiquitous dynamical process that transfers energy across many spatial and temporal scales in astrophysical and space plasma systems. Although the theory of anisotropic magnetohydrodynamic (MHD) turbulence has successfully described natural phenomena, its core prediction of an Alfvénic transition from weak to strong MHD turbulence when energy cascades from large to small scales has not been observationally confirmed. Here we report evidence for the Alfvénic weak-to-strong transition in small-amplitude, turbulent MHD fluctuations in Earth’s magnetosheath using data from four Cluster spacecraft. Our observations demonstrate the universal existence of strong turbulence accompanied by weak turbulent fluctuations on large scales. Moreover, we find that the nonlinear interactions of MHD turbulence are crucial to the energy cascade, as they broaden the cascade directions and fluctuating frequencies. The observed connection between weak and strong MHD turbulence systems may be present in many astrophysical environments, such as star formation, energetic particle transport, turbulent dynamos, and solar corona or solar wind heating.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Shadow Lagrangian dynamics for superfluidity

Motivated by a similar approach for Born-Oppenheimer molecular dynamics, this paper proposes an extended "shadow" Lagrangian density for quantum states of superfluids. The extended Lagrangian contains an additional field variable that is forced to follow the wave function of the quantum state through a rapidly oscillating extended harmonic oscillator. By considering the adiabatic limit for large frequencies of the harmonic oscillator, we can derive the two equations of motions, a Schrödinger-type equation for the quantum state and a wave equation for the extended field variable. The equations are coupled in a nonlinear way, but each equation individually is linear with respect to the variable that it defines. The computational advantage of this new system is that it can be easily discretized using linear time stepping methods, where we propose to use a Crank-Nicolson-type approach for the Schrödinger equation and an extended leapfrog scheme for the wave equation. Furthermore, the difference between the quantum state and the extended field variable defines a consistency error that should go to zero if the frequency tends to infinity. By coupling the time-step size in our discretization to the frequency of the harmonic oscillator we can extract an easily computable consistency error indicator that can be used to estimate the numerical error without additional costs. The findings are illustrated in numerical experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

HYDRODYNAMIC ANALYSIS OF MACROALGAE LOCAL MODEL USING COMPUTATIONAL FLUID DYNAMICS

In this article, a local scale, fully nonlinear coupled fluid-structural interaction (FSI) sugar kelp model has been developed using a computational fluid dynamics (CFD) method. In this model, to be consistent with available experimental data, the sugar kelp is approximated as elongated rectangles with smoothed isosceles triangles at the ends and a single kelp model with one end fixed in a channel with constant current model is developed. Several different current speeds are simulated, and the resulting drag forces and calculated drag coefficients are validated by comparison with experimental data from the literature. In a previous study, a global scale model was developed using a computational structural dynamics (CSD) method to simulate macroalgae farming system and guide the system configuration design. In the global scale model, the hydrodynamic forces are calculated using Morison’s equation and the kinematics and dynamics of the sugar kelp are simplified and the group of kelps attached to the long line is modeled as a slender structure with the same length and an effective diameter such that the volumes are consistent with the real physical system. This simplified model matches the weight and buoyancy but adjusting the hydrodynamic properties when the general hydrodynamic coefficients are employed. Therefore, optimal hydrodynamic coefficients used in global scale model were determined to obtain the hydrodynamic force more accurately. The validated local scale model is then be applied to determine the hydrodynamic coefficients of the simplified sugar kelp model for global dynamic analysis.

hydrodynamic modeling, macroalgae, CFD model↗

Solving the Dynamics-Aware Economic Dispatch Problem with the Koopman Operator

The dynamics-aware economic dispatch (DED) problem embeds low-level generator dynamics and operational constraints to enable near real-time scheduling of generation units in a power network. DED produces a more dynamic supervisory control policy than traditional economic dispatch (T-ED) that reduces overall generation costs. However, in contrast to T-ED, DED is a nonlinear, non-convex optimization problem that is computationally prohibitive to solve. We introduce a machine learning-based operator-theoretic approach for solving the DED problem efficiently. Specifically, we develop a novel discrete-time Koopman Operator (KO) formulation that embeds domain information into the structure of the KO to learn high-fidelity approximations of the generator dynamics. Using the KO approximation, the DED problem can be reformulated as a computationally tractable linear program (abbreviated DED-KO). We demonstrate the high solution quality and computational-time savings of the DED-KO model over the original DED formulation on a 9-bus test system.

King, Ethan↗

Hybrid model predictive control techniques for safety factor profile and stored energy regulation while incorporating NBI constraints

Abstract A novel hybrid Model Predictive Control (MPC) algorithm has been designed for simultaneous safety factor ( q ) profile and stored energy ( w ) control while incorporating the pulse-width-modulation constraints associated with the neutral beam injection (NBI) system. Regulation of the q -profile has been extensively shown to be a key factor for improved confinement as well as non-inductive sustainment of the plasma current. Simultaneous control of w is necessary to prevent the triggering of pressure-driven magnetohydrodynamic instabilities as the controller shapes the q profile. Conventional MPC schemes proposed for q -profile control have considered the NBI powers as continuous-time signals, ignoring the discrete-time nature of these actuators and leading in some cases to performance loss. The hybrid MPC scheme in this work has the capability of incorporating the discrete-time actuator dynamics as additional constraints. In nonlinear simulations, the proposed hybrid MPC scheme demonstrates improved q -profile+ w control performance for NSTX-U operating scenarios.

Physics↗

Oscillation Risks of Grid-Following and Grid-Forming Inverter-Based Resources in Series-Compensated Networks

Here, this paper investigates the dynamic behavior of a grid-connected inverter-based resource (IBR) when connected radially to a series compensated line. Potential interactions between the series compensation and the IBR have been identified for both types: grid-following (GFL) or grid-forming (GFM). The study begins with electromagnetic transient (EMT) simulations to demonstrate stability issues. Subsequently, nonlinear analytical models are formulated in the dq frame, validated against the EMT simulation, and leveraged to assess eigenvalues and participation factors. Influencing factors of the dominant oscillation modes have been identified. The analysis results show that series compensation may make a mode associated with the synchronization unit unstable. Furthermore, customized feedback systems are built for the synchronizing loop. Series compensation can increase the sensitivity of the voltage phase angle towards the synchronizing angle, and introduce phase lag in the real power response towards the synchronizing angle. These factors may cause interactions with the phase-locked loop in GFL-IBR systems and with power-based synchronization in GFM-IBR systems, potentially leading to instability.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Modeling atmospheric data and identifying dynamics Temporal data-driven modeling of air pollutants

Atmospheric modelling has recently experienced a surge with the advent of deep learning. Most of these models, however, predict concentrations of pollutants following a data-driven approach in which the physical laws that govern their behaviors and relationships remain hidden. Furthermore, with the aid of real-world air quality data collected hourly in different stations throughout Madrid, we present a case study using a series of data-driven techniques with the following goals: (1) Find systems of ordinary differential equations that model the concentration of pollutants and their changes over time; (2) assess the performance and limitations of our model using stability analysis; (3) reconstruct the time series of chemical pollutants not measured in certain stations using delay coordinate embedding results.

54 ENVIRONMENTAL SCIENCES↗

CMLM (Co-Optimized Machine-Learned Manifolds) [SWR-23-41]

Co-optimized Machine-Learned Manifolds (CMLM) is a data-driven approach for developing reduced-order manifold models for high-dimensional chemically reacting systems. It involves a specially designed neural network, the training of which simultaneously optimizes linear combinations of species that define the manifold, nonlinear mapping to outputs of interest such as reaction rates, and (optionally) subfilter closure for large eddy simulation. This software package provides an implementation of the CMLM approach in Python using the PyTorch machine learning library. A few example cases are included, showing how the tool can be applied to different types of data from 0D and 1D reacting simulations performed using Cantera. The neural networks can be saved in a format that is readable by the Pele suite of combustion solvers for use in reacting computational fluid dynamics simulations. This software repository contains several python scripts to perform various tasks associated with the Co-optimized Machine Learned Manifolds (CMLM) model, which is described in Perry, Henry de Frahan, and Yellapantula, CNF, 2022 (https://doi.org/10.1016/j.combustflame.2022.112286). This includes not only the code that defines the CMLM model, but also scripts to generate suitable training data, scripts to pre-process the data, scripts to train the CMLM model, and scripts to plot the output, as well as various other helper files. The scripts depend on several commonly used python libraries for data analysis and chemical reaction computations. The trained models that result from this tool are designed to work with the an interface being implemented in the Pele suite of reacting flow solvers (https://github.com/AMReX-Combustion).

Perry, Bruce↗