Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear dynamical systems”

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 55 records · Page 3

Thermal Reservoir Networks for Modularly Expandable Thermal Microgrids

The Department of Defense (DoD) faces the substantial challenge of cost-effectively retrofitting one to two installations per month, each comprising approximately 1,000 buildings, to improve resilience, reduce energy consumption, and enhance energy supply security. Achieving these objectives requires optimal system selection and effective risk mitigation during system integration. To address this need, we introduce Platform-Based Design (PBD), a structured, hierarchical methodology adapted from other industrial sectors to the domain of energy system retrofits. We demonstrate the effectiveness of PBD through a techno-economic feasibility study comparing geothermal-coupled thermal energy networks (TENs) with conventional energy systems for heating, cooling, and powering 17 buildings at Joint Base Andrews (JBA) in Maryland. Our analysis illustrates that the PBD approach enables rigorous, data-driven, sequential decision making, resulting in a family of Pareto-optimal systems, among which the TEN emerged as the most promising solution. The selected TEN design integrates geothermal borefields, heat recovery heat pumps, photovoltaic (PV) arrays, and battery storage. Compared to the baseline system – gas heating combined with air-source chillers – the proposed TEN reduces annual imported energy by 74% and peak electricity demand by 45%, achieves a levelized cost of energy of $\$0.210$/kWh, and substantially enhances resilience. Life-cycle costs increase by approximately 6%, and initial investment costs are about 2.5 times higher than the baseline. However, if central plant infrastructure, district loops, and utility-scale PV and battery systems are privately funded and operated, the initial investment would fall below the baseline system cost. Critical to achieving these significant performance improvements were detailed nonlinear dynamic simulations coupling geothermal heat transfer, energy system operation, and realistic feedback control logic. These simulations identified essential design modifications and control strategy refinements that substantially reduced energy use, peak demand, and compressor shortcycling, thereby improving durability and reliability—issues that would have been significantly more expensive to resolve during operation. Additionally, the verification step highlighted sensitivities to key design parameters that could reduce initial investment by approximately $\$2$ million and reduce annual life-cycle costs more than $\$300,000$. We recommend adopting the PBD methodology for future feasibility studies and TEN pilot projects to gain valuable operational experience. Furthermore, we recommend that DoD invest in transferring and scaling the PBD methodology to other installations. This entails developing standardized computational frameworks and component libraries as well as training industry in conducting PBD. Such investments would enable rapid, robust, reliable, and cost-effective retrofits, supporting DoD’s ambitious energy system modernization goals.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Order conditions for nonlinearly partitioned Runge-Kutta methods

Recently, a new class of nonlinearly partitioned Runge–Kutta (NPRK) methods was proposed for nonlinearly partitioned systems of autonomous ordinary differential equations y' = F(y, y). The target class of problems are those in which different scales, stiffnesses, or physics are coupled in a nonlinear way, wherein the desired partition cannot be written in a classical additive or component-wise fashion. Here we use a rooted-tree analysis to derive full-order conditions for NPRKM methods, where M denotes the number of nonlinear partitions. Due to the nonlinear coupling and thereby the mixed product differentials, it turns out that the standard node-colored rooted tree analysis used in analyzing ODE integrators does not naturally apply. Instead we develop a new edge-colored rooted-tree framework to address the nonlinear coupling. The resulting order conditions are enumerated, are provided directly for up to fourth order with M = 2 and third order with M = 3, and are related to existing order conditions of additive and partitioned RK methods. We conclude with an example that shows how the nonlinear order conditions can be used to obtain an embedded estimate of the state-dependent nonlinear coupling strength in a dynamical system.

97 MATHEMATICS AND COMPUTING↗

Ultrafast nonlinear absorption of Haldane model quantum dots

We study theoretically the nonlinear absorbance of Haldane model quantum dots (QDs) placed in the field of an ultrashort and strong optical pulse. The absorbance strongly depends on the frequency of the pulse. When the frequency of the pulse is much less than the QD bandgap, the absorbance shows strong dependence on the pulse amplitude and, as a function of an internal phase of the Haldane model, the absorbance has maxima at intermediate values of the phase. When the frequency of the pulse becomes closer to, but still less than, the bandgap, the absorbance has a weak dependence on the pulse amplitude and, as a function of the internal phase, it has a maximum at the phase of 900 when the QD bandgap also has the smallest value. Furthermore, nonlinear electron dynamics in such QD systems changes from almost reversible one at small pulse frequencies to highly irreversible dynamics at large frequencies of the pulse.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Nonlinear Topological Photonics: Capturing Nonlinear Dynamics and Optical Thermodynamics

Combining multiple optical resonators or engineering dispersion of complex media has provided an effective method for demonstrating topological physics controlling photons in unprecedented ways such as unidirectional light propagation and spatially localized modes between an interface or on a corner. Further, adding nonlinear responses to those topological photonic systems has enabled achieving diverse phases of photons in both space and time, allowing for more functionalities in photonic devices that provide a new playground for studying dynamic features of nonlinear topological systems. However, most methods for describing nonlinear topological photonic systems rely on linear topological theories, making it challenging to accurately characterize the topology of nonlinear systems. Thus, substantial efforts have focused on rigorously describing nonlinear topological phases and developing effective tools to analyze nonlinear topological effects. Meanwhile, coupled multimode optical waveguides with nonlinear dynamic responses provide an excellent platform for the statistical description of photons, opening a new paradigm called “optical thermodynamics”. This review will introduce the basic concepts of nonlinear topological photonics and the recent development of theoretical approaches focusing on data-driven approaches for creating phase diagrams as well as the spectral localizer framework and the pseudospectrum method for understanding optical nonlinearities in topological systems. In addition, the new concept of optical thermodynamics will be introduced with some recent theoretical works.

Topological photonics↗

Uncertainty quantification of a physics-informed model based on sparse identification of a Thermal Energy Distribution System

Integrated energy systems (IES)s are crucial for enhancing the economy and efficiency of power generation sources (e.g., nuclear energy) necessary to unleash American energy dominance. These systems can be integrated with thermal energy storage (TES) and intermittent renewable energies to optimize overall energy use, peak-load regulation, and demand-side responses. However, the stabilization of energy generation, transport, and utilization introduces operational complexities that exceed the challenges of managing each sub-component individually. Currently, though IESs rely on human operators for efficiency and stability, reducing human error risk and enhancing performance through automation is highly desirable. Recent advances at Idaho National Laboratory have demonstrated successful control of the Thermal Energy Distributed System (TEDS). However, the automatic control system depends on a deterministic Sparse Identification of Nonlinear Dynamics with Control (SINDyC) model, which are trained based on simulation data from physics-based simulations. Because of uncertainties in physics-based simulation, SINDyC model results in large discrepancies against experimental data and cannot be reliably used in automatic control. In this paper, we present an innovative approach to address these discrepancies by quantifying uncertainties and developing a more robust model. We first generated trajectories by using first-principles physics codes to encapsulate the experiment. Next, we trained thousands of models by randomly sampling these trajectories. We then collapsed all those models into one probabilistic SINDyC by fitting a multivariate Gaussian distribution onto the resulting coefficient’s distribution. Despite its simplicity, our approach successfully produced 95% confidence intervals that captured the experimental trajectories. It even did so with a higher probability and better U-pooling score across six of the seven relevant quantities of interest (QoIs), as compared to other classical approaches. In conclusion, ongoing research is focusing on generating new experimental trajectories to validate this approach, and on employing Bayesian calibration to refine parametric uncertainties and guide future model development efforts.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Recurrent convolutional neural networks for modeling nonadiabatic dynamics of quantum-classical systems

Recurrent neural networks (RNNs) have recently been extensively applied to model the time evolution in fluid dynamics, weather predictions, and even chaotic systems due to their ability to capture temporal dependencies and sequential patterns in data. Here we present an RNN model based on convolutional neural networks for modeling the nonlinear nonadiabatic dynamics of hybrid quantum-classical systems. The dynamical evolution of the hybrid systems is governed by equations of motion for classical degrees of freedom and von Neumann equation for electrons. The Physics-Aware Recurrent Convolution (PARC) neural network structure incorporates a differentiator-integrator architecture that inductively models the spatiotemporal dynamics of generic physical systems. Here, we apply our RNN approach to learn the space-time evolution of a one-dimensional semiclassical Holstein model after an interaction quench. For shallow quenches (small changes in electron-lattice coupling), the deterministic dynamics can be accurately captured using a single-CNN-based recurrent network. In contrast, deep quenches induce chaotic evolution, making long-term trajectory prediction significantly more challenging. Nonetheless, we demonstrate that the PARC-CNN architecture can effectively learn the statistical climate of the Holstein model under deep-quench conditions.

Holstein model↗

Customizable wave tailoring nonlinear materials enabled by bilevel inverse design

Abstract Passive wave transformation via nonlinearity is ubiquitous in settings from acoustics to optics and electromagnetics. It is well known that different nonlinearities yield different effects on propagating signals, which raises the question of “what precise nonlinearity is the best for a given wave tailoring application?” In this work, considering a one-dimensional spring-mass chain connected by polynomial springs (a variant of the Fermi-Pasta-Ulam-Tsingou system), we introduce a bilevel inverse design method which couples the shape optimization of structures for tailored constitutive responses with reduced-order nonlinear dynamical inverse design. We apply it to two qualitatively distinct problems—minimization of peak transmitted kinetic energy from impact, and pulse shape transformation—demonstrating our method’s breadth of applicability. For the impact problem, we obtain two fundamental insights. First, small differences in nonlinearity can drastically change the dynamic response of the system, from severely under- to outperforming a comparative linear system. Second, the oft-used strategy of impact mitigation via “energy locking” bistability can be significantly outperformed by our optimal nonlinearity. We validate this case with impact experiments and find excellent agreement. This study establishes a framework for broader passive nonlinear mechanical wave tailoring material design, with applications to computing, signal processing, shock mitigation, and autonomous materials.

Science & Technology - Other Topics↗

Experimental Studies of Nonlinear Integrable Optics with Elliptic Potentials

The stable transport of charged particle beams is the core challenge in designing and operating accelerators. The current paradigm for transverse focusing is based on quadrupole and dipole elements, described as a linear integrable Hamiltonian by Courant and Snyder. The operational limits of high intensity accelerators are often determined by collective instabilities within the beam. Nonlinear elements may be added to suppress certain forms of these collective effects, but restrict the range of stable trajectories. These shortcomings motivate extending to a nonlinear integrable system which can offer the benefits of suppressing collective instabilities without limiting the stable trajectories. Danilov and Nagaitsev proposed a novel nonlinear integrable system with elliptical potentials, which could be implemented as magnetic elements in an accelerator. Such a system has been implemented for practical verification in the integrable optics test accelerator (IOTA), a small storage ring constructed for beam dynamics studies. Electron beam studies in IOTA treat the low-emittance beam as a macroparticle for detailed probing of the expected single particle dynamics. This dissertation details new measurements of the nonlinear integrable system relevant for practical implementation. Turn-by-turn measurements of the kicked beam responses are used for phase space reconstruction and analysis of the predicted nonlinear dynamics. The stability and aperture for various nonlinear configurations are measured with beam losses. Amplitude dependent detuning, a core figure of merit for suppressing instabilities, is measured and compared with high fidelity particle tracking simulations. The synchrotron radiation images of circulating beam allow direct measurements of the topology and lifetime in configurations where nonlinear focusing terms dominate.

Wieland, John [Michigan U.] (ORCID:000000028971852↗

Learning nonlinear operators in latent spaces for real-time predictions of complex dynamics in physical systems

Abstract Predicting complex dynamics in physical applications governed by partial differential equations in real-time is nearly impossible with traditional numerical simulations due to high computational cost. Neural operators offer a solution by approximating mappings between infinite-dimensional Banach spaces, yet their performance degrades with system size and complexity. We propose an approach for learning neural operators in latent spaces, facilitating real-time predictions for highly nonlinear and multiscale systems on high-dimensional domains. Our method utilizes the deep operator network architecture on a low-dimensional latent space to efficiently approximate underlying operators. Demonstrations on material fracture, fluid flow prediction, and climate modeling highlight superior prediction accuracy and computational efficiency compared to existing methods. Notably, our approach enables approximating large-scale atmospheric flows with millions of degrees, enhancing weather and climate forecasts. Here we show that the proposed approach enables real-time predictions that can facilitate decision-making for a wide range of applications in science and engineering.

97 MATHEMATICS AND COMPUTING↗

Equation-Free Coarse Control of Distributed Parameter Systems via Local Neural Operators

The control of high-dimensional distributed parameter systems (DPS) remains a challenge when explicit coarse-grained equations are unavailable. Classical equation-free (EF) approaches rely on fine-scale simulators treated as black-box timesteppers. However, repeated simulations for steady-state computation, linearization, and control design are often computationally prohibitive, or the microscopic timestepper may not even be available, leaving us with data as the only resource. We propose a data-driven alternative that uses local neural operators, trained on spatiotemporal microscopic/mesoscopic data, to obtain efficient short-time solution operators. These surrogates are employed within Krylov subspace methods to compute coarse steady and unsteady-states, while also providing Jacobian information in a matrix-free manner. Krylov-Arnoldi iterations then approximate the dominant eigenspectrum, yielding reduced models that capture the open-loop slow dynamics without explicit Jacobian assembly. Both discrete-time Linear Quadratic Regulator (dLQR) and pole-placement (PP) controllers are based on this reduced system and lifted back to the full nonlinear dynamics, thereby closing the feedback loop.

93B52, 93C20, 47N70, 65J15, 65M32, 68T07, 68T20, 6↗

Applying Particle Swarm Optimization and Extended Kalman Filtering to Model Kaplan Generation Dynamics for Hydropower Systems

Variable renewable generation is increasing the need for hydropower plants to provide fast and flexible grid support, which places new demands on plant-level dynamic models used for monitoring, control, and operational decision-making. This need is especially important for hydroelectric systems, where turbine and generator dynamics are strongly coupled, nonlinear, and time-varying, making accurate real-time representation difficult. To address this problem, this paper develops a digital twin (DT) framework for a synchronous generator–Kaplan turbine system using an explicit separation of slow turbine dynamics and fast generator dynamics. The turbine subsystem is represented by a six-coefficient model, whose parameters are identified offline using particle swarm optimization, while the generator subsystem is updated online through an extended Kalman filter for real-time state and parameter estimation. These models are integrated within a closed-loop simulation that includes a proportional–integral–derivative–double-derivative governor and excitation system, allowing the DT to track plant behavior under realistic operating conditions. Unlike prior studies that treat turbine and generator modeling separately or rely mainly on simulated inputs, the proposed framework is validated using real operational data from a hydropower plant. Results show that the DT reproduces terminal voltage, active power, and reactive power with a normalized root mean square error of approximately 5%. This hybrid offline–online formulation constitutes the main contribution of the work, providing an adaptive and practically deployable DT for hydropower systems with direct relevance to control improvement, performance monitoring, and grid-support applications under high renewable penetration.

13 HYDRO ENERGY↗

Accurate data-driven surrogates of dynamical systems for forward propagation of uncertainty

Stochastic collocation (SC) is a well-known non-intrusive method of constructing surrogate models for uncertainty quantification. In dynamical systems, SC is especially suited for full-field uncertainty propagation that characterizes the distributions of the high-dimensional solution fields of a model with stochastic input parameters. However, due to the highly nonlinear nature of the parameter-to-solution map in even the simplest dynamical systems, the constructed SC surrogates are often inaccurate. Here, this work presents an alternative approach, where we apply the SC approximation over the dynamics of the model, rather than the solution. By combining the data-driven sparse identification of nonlinear dynamics framework with SC, we construct dynamics surrogates and integrate them through time to construct the surrogate solutions. We demonstrate that the SC-over-dynamics framework leads to smaller errors, both in terms of the approximated system trajectories as well as the model state distributions, when compared against full-field SC applied to the solutions directly. We present numerical evidence of this improvement using three test problems: a chaotic ordinary differential equation, and two partial differential equations from solid mechanics.

42 ENGINEERING↗

Dynamics of McMillan mappings II. axially symmetric map

Here, in this article, we investigate the transverse dynamics of a single particle in a model integrable accelerator lattice, based on a McMillan axially-symmetric electron lens. Although the McMillan e-lens has been considered as a device potentially capable of mitigating collective space charge forces, some of its fundamental properties have not been described yet. The main goal of our work is to close this gap and understand the limitations and potentials of this device. It is worth mentioning that the McMillan axially symmetric map provides the first-order approximations of dynamics for a general linear lattice plus an arbitrary thin lens with motion separable in polar coordinates. Therefore, advancements in its understanding should give us a better picture of more generic and not necessarily integrable round beams. In the first part of the article, we classify all possible regimes with stable trajectories and find the canonical action-angle variables. This provides an evaluation of the dynamical aperture, Poincaré rotation numbers as functions of amplitudes, and thus determines the spread in nonlinear tunes. Also, we provide a parameterization of invariant curves, allowing for the immediate determination of the map image forward and backward in time. The second part investigates the particle dynamics as a function of system parameters. We show that there are three fundamentally different configurations of the accelerator optics causing different regimes of nonlinear oscillations. Each regime is considered in great detail, including the limiting cases of large and small amplitudes. In addition, we analyze the dynamics in Cartesian coordinates and provide a description of observable variables and corresponding spectra.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Multistage economic MPC for systems with a cyclic steady state: A gas network case study

Multistage model predictive control (MPC) provides a robust control strategy for dynamic systems with uncertainties and a setpoint tracking objective. Moreover, extending MPC to minimize an economic cost instead of tracking a pre-calculated optimal setpoint improves controller performance. This paper presents a novel multistage economic nonlinear model predictive control (E-NMPC) framework for dynamic systems operating under uncertainty, with specific application to natural gas transmission networks. A key innovation lies in the integration of cyclic steady-state (CSS) constraints within the multistage MPC formulation, enabling the controller to manage periodic operating conditions commonly observed in energy systems. A Lyapunov-based descent condition is enforced to ensure robust stability of the controller. The multistage economic MPC framework is validated on two gas pipeline case studies, where it successfully minimizes net energy consumption, respects operational constraints under uncertain demand profiles, and guides the network to optimal cyclic operation. The Lyapunov function remains bounded in both case studies, validating the robust stability of multistage E-NMPC.

03 NATURAL GAS↗

Reduced-action-integral approach for photon-photon interactions in vacuum

Electromagnetic waves propagating through vacuum can polarize virtual electron–positron pairs; this polarization, in turn, nonlinearly modifies their propagation. A semi-classical nonlinear wave equation describing the propagation is derived from the Euler–Heisenberg Lagrangian density, which captures vacuum polarization effects up to the one-loop level. In this article, we present a reduced-actionintegral approach that enables rapid modeling of nonlinear phenomena arising from the Euler– Heisenberg Lagrangian. Application of the variational principle to the reduced action provides equations of motion for familiar light-pulse parameters, such as spot size, phase, polarization, and phase-front curvature, without requiring full-field simulations. Three examples demonstrate the utility of the approach: phase modulation, birefringence, and frequency mixing.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Studying Open Quantum Systems Relevant to Chemistry on a Trapped-Ion Quantum Simulator (Final Technical Report)

This project advances the trapped-ion quantum simulator as a versatile platform for studying open quantum system phenomena. We aim to contribute to the emerging quantum simulation toolkits and enable simulation of nanoscale energy processes. Trapped-ion platforms offer unique capabilities: their vibrational motion can be precisely manipulated, measured, and coherently coupled to auxiliary qubits. The vibrational mode can function both as a highly sensitive quantum sensor and a programmable environment bath. Using this platform, we achieved three major outcomes. First, we demonstrated using the vibrational mode as an ultrasensitive probe for testing fundamental physics, including possible nonlinear quantum mechanics effects. Second, we established that these modes can act as controllable baths in which tunable noise and loss can enhance or modify energy-transfer dynamics, providing the experimental preparation toward studying mechanisms relevant to chemical reactions and light-harvesting systems. Third, by introducing controllable nonlinear gain and loss, we showed theoretically how simulations using trapped ions can model vibrationally-assisted energy transport in a non‐Hermitian quantum system comprising a chromophore dimer weakly coupled to a vibrational mode. Exploring the non‐Hermitian dynamics of the whole system including vibrations, we found that energy transfer accompanied by absorption of phonons from a vibrational mode can be significantly enhanced near an exceptional point. This theoretical work on simulation of energy transfer processes in driven non‐Hermitian quantum systems revealed an interesting novel path to study open quantum systems dynamics under conditions of gain and loss. We then further explored the benefits of controllable gain and loss with an experimental realization of quantum analogs of nonlinear oscillators, namely, the van der Pol oscillator. Here we observed mutual synchronization mediated by collective dissipation between two oscillators. In parallel, we explored related quantum networking protocols using the same trapped-ion platform, developing fast, high-fidelity schemes for distributing entanglement. Together, these achievements show that trapped-ion vibrational modes provide a highly programmable and high-fidelity platform for investigating complex dissipative quantum behavior, while enabling new approaches to remote quantum sensing, energy science, and nonlinear quantum dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Effects of Central Cut Width on the Dynamical Characteristics of Box Assembly Structure

An investigation into the Box Assembly with Removable Component (BARC) structure is conducted by utilizing computational simulations and experimental structural testing in order to determine the complex dynamical responses instigated by the central cut of the system. Because the dynamics of the BARC system is complex, this study focuses primarily on analyzing the behavior of the box assembly (BA) system. The investigation explores the dynamics of the BA system by varying the central cut widths, ranging from a cut as wide as 0.5” cut to a 0.25” cut system, as well as a 0.1” cut and a system with no cut at all. Experimental testing is performed on each system including a free vibration test using an impact hammer to excite and identify the dominant frequencies of each structure. This testing is followed by pseudo-random vibration tests and swept sinusoidal excitation tests to determine the nonlinear aspects of these systems, such as the possible existence of nonlinear softening, hardening, and/or damping. The results show that nonlinear softening and nonlinear damping are present in each system. The no-cut system demonstrated the highest peak frequencies throughout all the tests, being the most rigid structure. The 0.25” cut system was shown to have the highest peak frequencies among all the cut systems in both the finite elemenet analysis (FEA) and impact testing. This trend did not continue, though, in the random and harmonic testing, possibly due to the added stiffness of the test setup with the slip table and stinger. The results show the importance of accurately measuring the central cut width and how possible geometric uncertainties change the overall dynamical behaviors of complex systems, such as natural characteristics, nonlinear responses, coupling of modes, and oscillating amplitudes.

Padilla, Christopher (ORCID:000900033446732X)↗

Emergence of complex-regular eutectic patterns in Al–Ge: Observations from correlative nano-imaging

Eutectic solidification exemplifies nonequilibrium pattern formation, making it a well-studied moving boundary problem. Yet the mechanisms behind the formation of complex-regular microstructures – particularly in highly anisotropic systems with a significant volume fraction of a faceted phase – remain poorly understood. Our understanding of such systems is made complicated by the nonlinear interface kinetics and unique growth dynamics characteristic of faceted phases. To address these challenges, we investigate a model Al–Ge eutectic system, where the faceted Ge phase constitutes a substantial volume fraction (~0.35) and where the two solid phases arrange into so-called “fishbone” or “feather” complex-regular patterns. Using synchrotron-based x-ray nano-imaging and nanotomography with high spatial resolution (22 nm per pixel), we capture in real-time the evolution of the solid–liquid interfaces and the resulting three-dimensional microstructures in this faceted/non-faceted eutectic system. By integrating these observations with electron backscattered diffraction, we elucidate the crystallographic biases on the solidification process and the mechanisms driving the formation of such complex-regular microstructures. These findings inform a new growth model for irregular eutectics in (near-)symmetrical phase diagrams, offering insight on advanced microstructural design and processing strategies. More broadly, we demonstrate how interfacial curvature is generated in irregular eutectic alloys and how it depends on the volume fraction of the faceted phase.

36 MATERIALS SCIENCE↗