Time-parallel optimal control of nonlinear dynamical systems
Explore the source record for details and available documents.
Ridzal, Denis [Sandia National Laboratories (SNL-N↗
SEARCH · Search NASA
Results for “Nonlinear 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.
Dynamically reconfigurable topological routing in nonlinear photonic systems
Explore the source record for details and available documents.
Wong, Stephan Manua [Sandia National Laboratories ↗
Detecting isolated resonance curves using fixed frequency voltage control tests
Isolated resonance curves, or isolas, are resonance branches of the harmonically forced system that exist separately from the main nonlinear forced response curve, leading to excessive vibrations. Traditional stepped or swept sine simulations and tests rely on continuation along the frequency parameter, typically resulting in a jump phenomenon along the primary resonance branch, prior to the disconnected isola. The main objective of this research is to propose an approach to identify isolated resonance curves by performing continuation along the input amplitude that initializes the response from a low-amplitude solution in the linear regime. Furthermore, this is achieved with the open-loop fixed frequency voltage control method that continues along the shaker voltage parameter and measures the so-called S-curves, which are theoretically a continuous solution branch that connect to the isola. The methodology is demonstrated on a fixture-wing-pylon assembly with a vibro-impact nonlinearity localized in a pylon subcomponent attachment. Multi-harmonic balance simulations are deployed to compute both the nonlinear forced response curves and S-curves to demonstrate the isola detection strategy on a reduced-order finite element model of the nonlinear system. Swept sine and fixed frequency voltage control tests are then conducted on the physical structure to demonstrate the isola detection experimentally, revealing the existence of the large amplitude vibrations that are undetected in the forces levels and frequencies measured with traditional frequency sweeping.
Characterization and Analytical Technique↗
A Tensor Network-Based Quantum Algorithm for the Nonlinear 1D Burgers' Equation
In this work, we implement a tensor network-based quantum algorithm to solve unsteady, nonlinear partial differential equations (PDEs). The challenge lies in how to effectively represent, encode, process, and evolve the nonlinear system of PDEs on quantum computers. We will discuss the new techniques using the compressible 1-dimensional (1D) Burgers' equation as an example, because it represents the fundamental nonlinear feature and yet removes certain complexity in physics, allowing us to focus on the design of quantum algorithms. Previous attempts to solve nonlinear PDEs in quantum computation have often involved storing multiple copies of solutions or employing linearizations. Neither is practical due to exponential scaling with evolution time or insufficient solution accuracy. Our framework is based on matrix product states (MPSs) and matrix product operators (MPOs). For example, the velocity field is represented by MPS, whereas the linear and nonlinear spatial differential terms of the velocity field are processed by MPOs. Our primary focus herein is to verify and validate the various tensor network components of the algorithm using solutions obtained by the classical algorithms on high performance computing (HPC) architectures. We use a classical time marching method to demonstrate the functionality of the tensor network operations to model the PDE and their robustness with the time evolution of the system. Our classical simulation results demonstrate the utility of tensor network-based operations in modeling nonlinear PDEs and highlight the necessity as well as potential advantages of using quantum simulations for these techniques.
Gopalakrishnan Meena, Murali [ORNL] (ORCID:0000000↗
Microscopic Theory of Nonlinear Hall Effect in Three-Dimensional Magnetic Systems
Abstract Nonlinear Hall effect (NLHE) has been detected in various of condensed matter systems. Unlike linear Hall effect, NLHE may exist in physical systems with broken inversion symmetry in crystals. On the other hand, real space spin texture may also break inversion symmetry and result in NLHE. We employ the Feynman diagrammatic technique to calculate non-linear Hall conductivity (NLHC) in three-dimensional magnetic systems. The results connect NLHE with the physical quantity of emergent electrodynamics which originates from the magnetic texture. The leading order contribution of NLHC, χabb , is proportional to the emergent toroidal moment T a e , which reflects how the spin textures wind in three dimensions.
Hou 侯, Wen-Tao 文涛↗
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↗
Optimal Control of SOEC-Based Hydrogen Production Systems for Demand Response Using Deep Reinforcement Learning in Smart Grids
Solid oxide electrolysis cell (SOEC) hydrogen production technology can range in size from small, appliance-size equipment to large-scale, central production facilities that can be tied directly to renewable or non-greenhouse-gas-emitting forms of electricity production, making it an ideal resource for demand response (DR). The SOEC hydrogen production system is a complex integrated system that encompasses fluid dynamics, electrical dynamics, and electrochemical and thermal dynamics, all of which involve non-linearity and non-convexity. Proper control of the SOEC hydrogen production system is crucial to enable its participation in the DR program. Here, to overcome the difficulty of designing an explicit control law for such nonlinear systems with nonconvex optimization features in DR applications, deep reinforcement learning (DRL) is explored to achieve the optimal control of the SOEC system for DR participation. Specifically, a twin delayed deterministic policy gradient (TD3) control framework is applied to achieve optimal response performance during DR events by considering power tracking error and hydrogen production efficiency with a suitable reward function. Two case studies with grid connections for tracking different DR commands were investigated. The first case study involved operating conditions reaching the boundaries, while the second involved operating conditions within the boundaries. The results showed that the proposed DRL-based control for SOEC can track the DR signal in a timely manner while maintaining high energy efficiency.
08 HYDROGEN↗
Efficient data-driven regression for reduced-order modeling of spatial pattern formation
We present an efficient data-driven regression approach for constructing reduced-order models (ROMs) of reaction-diffusion systems exhibiting pattern formation. The ROMs are learned non-intrusively from available training data of physically accurate numerical simulations. The method can be applied to general nonlinear systems through the use of polynomial model form, while not requiring knowledge of the underlying physical model, governing equations, or numerical solvers. The process of learning ROMs is posed as a low-cost least-squares problem in a reduced-order subspace identified via Proper Orthogonal Decomposition (POD). Numerical experiments on classical pattern-forming systems–including the Schnakenberg and Mimura–Tsujikawa models–demonstrate that higher-order surrogate models significantly improve prediction accuracy while maintaining low computational cost. The proposed method provides a flexible, non-intrusive model reduction framework, well suited for the analysis of complex spatio-temporal pattern formation phenomena.
Data-driven modeling↗
Towards a Quantum Algorithm for the Incompressible Nonlinear Navier-Stokes Equations
In this work, we present novel concepts for quantum algorithms to solve transient, nonlinear partial differential equations (PDEs). The challenge lies in how to effectively represent, encode, process, and evolve the nonlinear system of PDEs on quantum computers. We will discuss the new techniques using the incompressible Navier-Stokes equations as an example, because it represents the fundamental nonlinear feature and yet removes certain complexity in physics, allowing us to focus on the design of quantum algorithms. Previous attempts solving nonlinear PDEs in quantum computation have often involved storing multiple copies of solutions or employing linearizations. Neither is practical due to exponential scaling with evolution time or insufficient solution accuracy. We propose a new framework based on matrix product states (MPSs) and matrix product operators (MPOs), in addition to the Krylov subspace methods. For example, the solution variables of the Navier-Stokes equations are represented by MPSs, and the linear and nonlinear terms are processed by MPOs. The time evolution of the operators is attained by a fast-forwarding algorithm using Krylov subspace methods. Furthermore, we discuss various techniques for efficient encoding of MPSs, measurement reduction for MPOs, and use of tensor operations to treat multi-variate, multi-physics characteristics of Navier-Stokes.
Gopalakrishnan Meena, Murali [ORNL] (ORCID:0000000↗
Uncertainty propagation and sensitivity analysis for constrained optimization of nuclear waste vitrification
Abstract The vitrification of high‐level waste (HLW) by heating a mixture of glass‐forming chemicals (GFCs) with the waste can be improved using a constrained optimization problem. This study explores how different uncertainty propagation (UP) methods implemented with the optimization process can affect the glass formulation of nuclear waste glasses. UP is the effort of propagating uncertain inputs through a system to understand and quantify output distributions. Uncertainty intervals are crafted from output distributions to inform the optimization algorithm. UP is often implemented with Monte Carlo (MC) sampling for large nonlinear systems, which can be difficult to implement within a constrained optimization algorithm that requires derivative information. Other UP methods often used for optimization under uncertainty (OUU) can be designed to work within an established constrained optimization framework. Methods of UP are evaluated in this study including iterative sampling approaches, first‐order approximations, and surrogate modeling with machine learning (ML). A method of dimensional reduction based on global sensitivity analysis is introduced to support the UP methods for the large dimensionality of the problem. Analytical UP methods able to achieve similar optimums 10 times faster than the baseline MC approach, and produce 93.9% similar output distributions are reported.
12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
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↗
A Jacobian-free pseudo-arclength continuation method for phase transitions in inhomogeneous thermodynamic systems
Developing phase diagrams for inhomogeneous systems in thermodynamics is difficult, in part, due to the large phase space and the possibility of unstable and metastable solutions arising from first-order phase transitions. Pseudo-arclength continuation (PAC) is a method that allows one to trace out stable and unstable solutions of nonlinear systems. Typically, PAC utilizes the Jacobian in order to implement Newton (or quasi-Newton) steps. In this work, we present a Jacobian-free PAC method that is amenable to the usual workflows in inhomogeneous thermodynamics. We demonstrate our method in systems that have first-order phase transitions, including a novel example of polyelectrolyte complex coacervation in confinement, where multiple surface phase transitions occur and can overlap with one another.
Chemistry↗
Observation of Joule–Thomson photon-gas expansion
In recent years, a self-consistent optical thermodynamic framework has emerged that offers a systematic methodology to understand, harness and exploit the complex collective dynamics of multimode nonlinear systems. These developments now allow consideration of a series of longstanding problems in optics, including the prospect of funnelling the entire power flowing in a multimode system into its ground state, for which no methodology currently exists. Here, we demonstrate an all-optical Joule-Thomson expansion process mediated by photon-photon interactions whereby the temperature of the optical gas drops abruptly to zero. Our experiments in various configurations of coupled multicore nonlinear waveguide arrangements illustrate how light undergoing expansion-induced cooling can be channelled from arbitrary input states into the fundamental mode with near-unity efficiency. We show that the stability of the post-expansion state is ensured through an irreversible process of energy conversion. The all-optical thermodynamic phenomena explored in this study may enable innovative techniques where various uncorrelated but identical sources are merged into a unified spatially coherent state, offering a route for direct beam combining.
Kirsch, Marco S. (ORCID:0000000342798228)↗
Solid State Transformer Architecture and Control Compensation for Common Mode Currents
A high-altitude electromagnetic pulse (HEMP) or similar geomagnetic disturbance (GMD) has the potential to impact the operation of large-scale electric power grids. By introducing low-frequency common-mode (CM) currents, these events can degrade the performance of critical system components, such as large power transformers by introducing CM currents which can lead to magnetic saturation of the transformer core. In this work, a solid-state transformer (SST) is developed to replace susceptible equipment and improve grid resiliency by safely absorbing these CM disturbances. This device will be referred to as a common-mode solid-state transformer (CM-SST). An SST architecture based on a four-legged AC/DC converter is developed. This architecture enables active control of CM signals without disturbing the AC voltages or the real and reactive power delivery capabilities. A system-level model of this architecture is created, and time-domain simulations are performed to evaluate the SST’s performance in response to simulated CM disturbances. A control strategy for mitigating CM current is also investigated. Hamiltonian surface shaping and power flow control (HSSPFC) is used to design a nonlinear controller for the SST’s output inverter. The objectives of the controller are to suppress CM-induced AC current offsets and regulate AC currents to desired setpoints. Nonlinear system analysis is applied to design and validate the controller. Two cases are tested: (a) the proposed four-leg inverter and (b) a standard three-leg inverter. The results show that the proposed controller rapidly mitigates CM disturbances while maintaining high-quality AC current waveforms in the four-leg configuration. Finally, the hardware performance of an SST prototype is evaluated. In particular, the ability of the SST to safely redirect and absorb CM currents is demonstrated, showing how it can protect neighboring conventional transformers in the system. The study confirms that appropriate control laws allow the SST to protect both itself and adjacent transformers during a HEMP or GMD event.
24 POWER TRANSMISSION AND DISTRIBUTION↗
Convergence Criteria for Multiphysics Simulations
The behavior of engineered systems is often influenced by multiple physical phenomena, such as mechanical deformation, heat transfer, and chemical species transport and reactions. There are often strong interactions between these phenomena, and there is increasing interest in applying coupled-physics models to improve understanding of physical behavior under complex environmental conditions. Multiple simulation frameworks that facilitate coupled-physics simulations are in widespread use, and these employ a variety of techniques to account for interactions between those physics. Many frameworks solve the physics models independently and transfer results between them. Alternatively, a single monolithic system of equations for every physics model can be formed and solved. Each of these approaches has its benefits and drawbacks, and the optimal approach varies depending on the nature of the problem. The open-source MOOSE framework was developed targeting solution of large-scale multiphysics problems. Although it provides options for all these coupling approaches, its standard approach for multiphysics solutions is to form and solve a single monolithic system of equations containing the unknowns for all physics models. MOOSE provides a streamlined approach for users to define the solution variables, the terms in the partial differential equations pertaining to each variable, and interactions between solution variables. One aspect of the monolithic solution approach that can be problematic, however, is defining appropriate convergence criteria for the nonlinear system. A standard approach is to determine convergence is to simply take a norm of the residual vector corresponding to the full vector of unknowns. However, if the residual vector contains variables for multiple physics models, the magnitudes of those variables can differ significantly, and the variables can converge at significantly different rates from each other. It is important to ensure that the variables for each of the physics are converged, and also ensure that the convergence criteria are not excessively stringent in cases when there is little change in the solution. This talk presents representative multiphysics problems to highlight these issues, and shows strategies for convergence criteria in MOOSE that are robust for multiphysics models under a variety of conditions.
97 - MATHEMATICS AND COMPUTING↗
Stage-local partitioned two-step runge-kutta methods for large systems of ordinary differential equations
We introduce stage-local partitioned two-step Runge-Kutta methods are an extension of standard two-step Runge-Kutta methods, which are an alternative to the standard additive two-step Runge-Kutta methods currently existing in the literature. Furthermore, these new schemes are designed with an eye towards truly N-partitioned systems and leverage local stage approximations to make several computationally interesting approximations viable. Specifically, the focus on local stage approximations makes possible the construction of truly asynchronous schemes, in the parallel sense, possible. In addition, we show that an implicit-explicit approach to these schemes can lead to methods that require the inversion of only local nonlinear systems.
Applied Dynamical Systems↗
Complex Dependence of Calcite Crack Kinetics on Salinity: The Role of DLVO and Hydration Forces
Abstract Subcritical crack growth (SCG) plays an important role in many geological processes such as delayed earth rupture and rock weathering. The complex dependency of SCG on the in‐crack fluid chemistry, however, is still poorly understood. In this study, we utilize the newly developed surface force‐based fracture theory (SFFT) to elucidate the relative contributions of surface forces and solute transport to the crack growth kinetics of calcite in NaCl solutions. Expanding on Barenblatt's cohesive crack model, SFFT introduces an effective stress intensity at the crack tip that encompasses all the relevant intermolecular forces across the crack in addition to the external far‐field stresses. The nonlinear system of equations portraying the crack opening profile, the solute distribution in a propagating crack, and the crack growth velocity are numerically solved via an implicit scheme. After carefully calibrating the model for calcite‐water systems, the SFFT is used to predict the SCG response of calcite at different NaCl concentrations, based on various hypotheses. These predictions are then compared to existing SCG data from the literature. We demonstrate that the experimentally observed variation of SCG rate with NaCl concentration cannot be explained solely by DLVO forces (electrostatic and Van der Waals interactions). This can be remediated by introducing an exponentially decaying hydration force with a nonlinear, nonmonotonic dependence on NaCl concentration. Furthermore, we demonstrate that accounting for both diffusive and advective transport of ions is important in explaining the absence of a stage‐II SCG response for calcite in electrolyte solutions. Plain Language Summary Subcritical crack growth (SCG) refers to the slow propagation of cracks in materials under a stress below the threshold for catastrophic failure. SCG is a key process in many geological events, for example, delayed earth ruptures and rock weathering. New initiatives such as underground CO 2 and H 2 storage in carbonate reservoirs further call for better understanding of SCG in carbonate minerals subjected to varying fluid chemistry. This study examines the SCG of calcite, a key mineral found in carbonate rocks, intergranular cement in sandstones, and filling material in mineral veins and faults, determining their deformation and strength. A mathematical model is developed to describe how the crack opens and propagates, how solutes (like salts) distribute within the crack, and how the crack surfaces interact with each other. We used the model to predict calcite SCG in water at different salt concentrations and compared it with experimental data. Our results revealed that the hydration force is the dominating factor in determining the complex, non‐linear dependency of SCG on salinity. We also found that both the movement of ions by diffusion and by bulk water flow are crucial for explaining the SCG rates, especially when the cracks grow quickly. Key Points Surface Force‐Based Fracture Theory predicts the complex subcritical crack growth patterns of calcite crystals immersed in NaCl solutions Results highlight the dominant role of hydration forces in altering the fracture behavior of calcite compared to VdW and electric double‐layer forces Advective solute transport explains the absence of stages‐II and ‐III subcritical crack growth responses in solid‐liquid systems
DLVO↗
Probing Excited-State Dynamics of Transmon Ionization
The fidelity and quantum nondemolition character of the dispersive readout in circuit QED are limited by unwanted transitions to highly excited states at specific photon numbers in the readout resonator. This observation can be explained by multiphoton resonances between computational states and highly excited states in strongly driven nonlinear systems, analogous to multiphoton ionization in atoms and molecules. In this work, we utilize the multilevel nature of high-𝐸 𝐽 /𝐸 𝐶 transmons to probe the excited-state dynamics induced by strong drives during readout. With up to ten resolvable states, we quantify the critical photon number of ionization, the resulting state after ionization, and the fraction of the population transferred to highly excited states. Moreover, using pulse shaping to control the photon number in the readout resonator in the high-power regime, we tune the adiabaticity of the transition and verify that transmon ionization is a Landau-Zener-type transition. We further extend these methods to a typical transmon with 𝐸 𝐽 /𝐸 𝐶 ≈ 55 and probe the offset-charge dependence of ionization dynamics in a timed-resolved manner. Our experimental results agree well with the theoretical prediction from a semiclassical driven transmon model and may guide future exploration of strongly driven nonlinear oscillators.
cavity quantum electrodynamics↗