Engineering Papers⌕ Search

SEARCH · Engineering Papers

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.

At least 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Magnetic Energy Conversion in Magnetohydrodynamics: Curvature Relaxation and Perpendicular Expansion of Magnetic Fields

Abstract The mechanisms and pathways of magnetic energy conversion are an important subject for many laboratory, space, and astrophysical systems. Here, we present a perspective on magnetic energy conversion in magnetohydrodynamics through magnetic field curvature relaxation (CR) and perpendicular expansion (PE) due to magnetic pressure gradients, and quantify their relative importance in two representative cases, namely 3D magnetic reconnection and 3D kink-driven instability in an astrophysical jet. We find that the CR and PE processes have different temporal and spatial evolutions in these systems. The relative importance of the two processes tends to reverse as the system enters the nonlinear stage from the instability growth stage. Overall, the two processes make comparable contributions to magnetic energy conversion, with the PE process somewhat stronger than the CR process. We further explore how these energy conversion terms can be related to particle energization in these systems.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

An antiferromagnetic diode effect in even-layered MnBi 2 Te 4

In a p–n junction, the separation of positive and negative charges leads to diode transport, in which charge flows in only one direction. Non-centrosymmetric polar conductors are intrinsic diodes that could be of use in the development of nonlinear applications. Such systems have recently been extended to non-centrosymmetric superconductors, and the superconducting diode effect has been observed. Here, in this paper, we report an antiferromagnetic diode effect in a centrosymmetric crystal without directional charge separation. We observed large second-harmonic transport in a nonlinear electronic device enabled by the compensated antiferromagnetic state of even-layered MnBi 2 Te 4 . We show that this antiferromagnetic diode effect can be used to create in-plane field-effect transistors and microwave-energy-harvesting devices. We also show that electrical sum-frequency generation can be used as a tool to detect nonlinear responses in quantum materials.

electronic properties and materials↗
DOE OSTI2026

Lyapunov-Based Iterative Learning of Regions of Attraction for Autonomous Systems

This paper proposes a novel algorithm for estimating the region of attraction of equilibrium points for nonlinear discrete-time autonomous systems. The method iteratively expands an initial estimate of the region of attraction by constructing unions of sublevel sets of learned functions parametrized as neural networks. Unlike conventional techniques that rely on a single global Lyapunov function, the proposed approach provides a collection of local Lyapunov-like functions, enabling richer representations and potentially larger region of attraction estimates. These functions are trained using sampled state-space data, and their Lipschitz continuity ensures that desirable properties extend beyond the training samples. The devised strategy is tested via numerical simulations, demonstrating the effectiveness of the proposed approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

On the Impact of High-Order Harmonic Generation in Electrical Distribution Systems

The modern power grid has seen a rise in the integration of non-linear loads, presenting a significant concern for operators. These loads introduce unwanted harmonics, leading to potential issues such as overheating and improper functioning of circuit breakers. In pursuing a more sustainable grid, the adoption of electric vehicles (EVs) and photovoltaic (PV) systems in residential networks has increased. Understanding and examining the effects of high-order harmonic frequencies beyond $1.5$ kHz is crucial to understanding their impact on the operation and planning of electrical distribution systems under varying nonlinear loading conditions. This study investigates a diverse set of critical power electronic loads within a household modeled using PSCAD/EMTdc, analyzing their unique harmonic spectra. This information is utilized to run the time-series harmonic analysis program in OpenDSS on a modified IEEE 34 bus test system model. The impact of high-order harmonics is quantified using metrics that evaluate total harmonic distortion (THD), transformer harmonic-driven eddy current loss component, and propagation of harmonics from the source to the substation transformer.

Peerzada, Aaqib A. [BATTELLE (PACIFIC NW LAB)]↗
DOE OSTI2023

Experimental Performance of a Nonlinear Control Strategy to Regulate Temperature of a High-Temperature Solar Reactor

Abstract Despite the significant potential of solar thermochemical process technology for storing solar energy as solid-state solar fuel, several challenges have made its industrial application difficult. It is important to note that solar energy has a transient nature that causes instability and reduces process efficiency. Therefore, it is crucial to implement a robust control system to regulate the process temperature and tackle the shortage of incoming solar energy during cloudy weather. In our previous works, different model-based control strategies were developed namely a proportional integral derivative controller (PID) with gain scheduling and adaptive model predictive control (MPC). These methods were tested numerically to regulate the temperature inside a high-temperature tubular solar reactor. In this work, the proposed control strategies were experimentally tested under various operation conditions. The controllers were challenged to track different setpoints (500 °C, 1000 °C, and 1450 °C) with different amounts of gas/particle flowrates. Additionally, the flow controller was tested to regulate the reactor temperature under a cloudy weather scenario. The ultimate goal was to produce 5 kg of reduced solar fuel magnesium manganese oxide (MgMn2O4) successfully, and the controllers were able to track the required process temperature and reject disturbances despite the system's strong nonlinearity. The experimental results showed a maximum error in the temperature setpoint of less than 0.5% (6 °C), and the MPC controller demonstrated superior performance in reducing the control effort and rejecting disturbances.

Energy & Fuels↗
DOE OSTI2023

Neutrino fast flavor instability in three dimensions for a neutron star merger

The flavor evolution of neutrinos in core collapse supernovae and neutron star mergers is a critically important unsolved problem in astrophysics. Following the electron flavor evolution of the neutrino system is essential for calculating the thermodynamics of compact objects as well as the chemical elements they produce. Accurately accounting for flavor transformation in these environments is challenging for a number of reasons, including the large number of neutrinos involved, the small spatial scale of the oscillation, and the nonlinearity of the system. We take a step in addressing these issues by presenting a method which describes the neutrino fields in terms of angular moments. We apply our moment method to neutron star merger conditions and show it simulates fast flavor neutrino transformation in a region where this phenomenon is expected to occur. By comparing with particle-in-cell calculations we show that the moment method is able to capture the three phases of growth, saturation, and decoherence, and correctly predicts the lengthscale of the fastest growing fluctuations in the neutrino field.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Revisiting two thiophosphate compounds constituting d 0 transition metal HfP 2 S 6 and d 10 transition metal α-Ag 4 P 2 S 6 as multifunctional materials for combining second harmonic generation response and photocurrent response

Two acentric thiophosphate compounds, HfP 2 S 6 and α-Ag 4 P 2 S 6 , are revisited and studied as infrared nonlinear optical materials. HfP 2 S 6 and α-Ag 4 P 2 S 6 were structurally characterized without any property measurements. Here, in this work, HfP 2 S 6 and α-Ag 4 P 2 S 6 were synthesized via high temperature salt flux reactions. Low-temperature polymorph acentric α-Ag 4 P 2 S 6 was purified and grown as mm-sized crystals with the aid of AgBr flux. The AgBr flux was revealed to play an important role in stabilizing the acentric α-Ag 4 P 2 S 6 . The acentric α-Ag 4 P 2 S 6 transforms to centrosymmetric β-Ag 4 P 2 S 6 at 850(5) K, which is revealed by differential scanning calorimetry (DSC) analysis and powder X-ray diffraction experiments. HfP 2 S 6 and α-Ag 4 P 2 S 6 are discovered by UV-vis spectrum measurements as indirect bandgap semiconductors with bandgaps of 2.2(1) eV and 2.5(1) eV, respectively, which is supported by DFT calculations and TB-LMTO-ASA calculations. The bonding pictures of α-Ag 4 P 2 S 6 were studied by crystal orbital Hamilton population calculations (COHP) coupled with electron localization function (ELF) analysis. DFT calculations predict that HfP 2 S 6 and α-Ag 4 P 2 S 6 would exhibit different optical performances regardless of being constructed from identical [P 2 S 6 ] motifs. HfP 2 S 6 exhibits a low second harmonic generation (SHG) response, ~0.21 × AGS (for the sample of particle size of 25 μm). α-Ag 4 P 2 S 6 possesses moderate SHG response, ~0.61 × AGS (for the sample of particle size of 225 μm) coupled with a high laser damage threshold (LDT) of ~3.2 × AGS. Characteristics of high ambient stability, moderate bandgap and SHG response, type-I phase-matching capability, and high LDT together with the easy growth of large crystals make α-Ag 4 P 2 S 6 attractive for future infrared nonlinear optical applications. Photocurrent measurements found that α-Ag 4 P 2 S 6 and β-Ag 4 P 2 S 6 have high photocurrent response, 165 nA cm -2 and 135 nA cm -2 , respectively. α-Ag 4 P 2 S 6 is a new multifunctional material of the ternary Ag–P–S system, which combines nonlinear optical (NLO) properties and photocurrent response.

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

Mixed Delay/Nondelay Embeddings Based Neuromorphic Computing with Patterned Nanomagnet Arrays

Patterned nanomagnet arrays (PNAs) have been shown to exhibit a strong geometrically frustrated dipole interaction. Some PNAs have also shown emergent domain wall dynamics. Previous works have demonstrated methods to physically probe these magnetization dynamics of PNAs to realize neuromorphic reservoir systems that exhibit chaotic dynamical behavior and high-dimensional nonlinearity. These PNA reservoir systems from prior works leverage echo state properties and linear/nonlinear short-term memory of component reservoir nodes to map and preserve the dynamical information of the input time-series data into nondelay spatial embeddings. Such mappings enable these PNA reservoir systems to imitate and predict/forecast the input time series data. However, these prior PNA reservoir systems are based solely on the nondelay spatial embeddings obtained at component reservoir nodes. As a result, they require a massive number of component reservoir nodes, or a very large spatial embedding (i.e., high-dimensional spatial embedding) per reservoir node, or both, to achieve acceptable imitation and prediction accuracy. These requirements reduce the practical feasibility of such PNA reservoir systems. To address this shortcoming, we present a mixed delay/nondelay embeddings-based PNA reservoir system. Our system uses a single PNA reservoir node with the ability to obtain a mixture of delay/nondelay embeddings of the dynamical information of the time-series data applied at the input of a single PNA reservoir node. Our analysis shows that when these mixed delay/nondelay embeddings are used to train a perceptron at the output layer, our reservoir system outperforms existing PNA-based reservoir systems for the imitation of NARMA 2, NARMA 5, NARMA 7, and NARMA 10 time series data, and for the short-term and long-term prediction of the Mackey Glass time series data.

Ti, Changpeng↗
DOE OSTI2022

A modified model parametrization algorithm for solving a special type of heat and mass transfer systems

A new method for solving nonlinear heat and mass transfer design tasks was considered. Systems using the Number of Transfer Units (NTU) method are a special type of mathematical model of heat and mass exchangers. It was observed, that the NTU models in a form of differential-algebraic equations (DAEs) cannot be directly solved with higher values of NTU. The requirements for consistent initial conditions, as well as numerical limitations of DAEs solvers, result, that the solution to the considered design problems that cannot be obtained by a classical direct shooting procedure. To overcome the presented difficulties, the αDAE model optimization algorithm was adjusted for solving NTU-based models. The new approach consists of 3 main steps: 1) task discretization by a multiple-shooting approach, 2) design an appropriate function $f_{NTU}$(α) to effectively influence the variability of the state variables described by dynamical relations, 3) the iterative numerical optimization algorithm for the new parametrized system. Moreover, computations can be performed by a chosen numerical optimization approach, which can be communicated with an available outer procedure for solving differential-algebraic equations. The presented algorithm was implemented and applied to solve the design task with the NTU model of a counter-flow exchanger. Here, the new approach was used to modify the system dynamics to influence the difficulty of the considered problem. Finally, the presented method enabled failure-free numerical computations for the higher values of the NTU parameter.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A nonlinear journey from structural phase transitions to quantum annealing

Motivated by an exact mapping between equilibrium properties of a one-dimensional chain of quantum Ising spins in a transverse field (the transverse field Ising (TFI) model) and a two-dimensional classical array of particles in double-well potentials (the “$\phi$ 4 model”) with weak inter-chain coupling, we explore connections between the driven variants of the two systems. Here, we argue that coupling between the fundamental topological solitary waves in the form of kinks between neighboring chains in the classical $\phi$ 4 system is the analog of the competing effect of the transverse field on spin flips in the quantum TFI model. As an example application, we mimic simplified measurement protocols in a closed quantum model system by studying the classical $\phi$ 4 model subjected to periodic perturbations. This reveals memory/loss of memory and coherence/decoherence regimes, whose quantum analogs are essential in annealing phenomena. In particular, we examine regimes where the topological excitations control the thermal equilibration following perturbations. This paves the way for further explorations of the analogy between lower-dimensional linear quantum and higher-dimensional classical nonlinear systems.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Lyapunov-Based Iterative Learning of the Region of Attraction for Autonomous Systems

This presentation introduces a novel algorithm for estimating the region of attraction of equilibrium points for nonlinear discrete-time autonomous systems. The method iteratively expands an initial estimate of the region of attraction by constructing unions of sublevel sets of learned functions parametrized as neural networks. Unlike conventional techniques that rely on a single global Lyapunov function, the proposed approach provides a collection of local Lyapunov-like functions, enabling richer representations and potentially larger region of attraction estimates. These functions are trained using sampled state-space data, and their Lipschitz continuity ensures that desirable properties extend beyond the training samples. The devised strategy is tested via numerical simulations, demonstrating the effectiveness of the proposed approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Insights on the Influence of the Central-Cut Width of the Box Assembly with Removable Component on Its Dynamical Responses

This investigation focuses on the dynamical effects caused by varying the central-cut width within the Box Assembly with Removable Component (BARC) system. The central-cut widths included in this study are a 0.5″ cut, a 0.25″ cut, a thin 0.1″ cut, and a structure that did not have a cut at all. Finite element analysis was conducted to determine the mode shapes and natural frequencies of each of the BARC structures. Structural dynamics experiments were run to examine the effects of the central-cut width on the dynamical responses and nonlinear characteristics of the BARC system. Free vibration testing with an impact hammer was carried out to excite the system and extract the dominant frequencies and directions of the significant responses. A pseudorandom vibration test that allows for the qualitative determination of any nonlinear behavior within the system was performed. This type of behavior can include nonlinear softening, nonlinear hardening, and the most common, nonlinear damping due to the presence of several bolted-joint connections and the possible activation of geometric and inertia nonlinearities. To quantitatively investigate the impacts of the central-cut width on the dynamics of the system, swept sinusoidal testing was conducted. It is determined that almost all systems with central cuts demonstrate the presence of nonlinear softening, but at times, nonlinear hardening trends are seen, particularly in the 0.1″ cut and no-cut systems when testing harmonically. Each of the central-cut systems displays nonlinear damping, with the amount of damping generally increasing as the central cut decreases in size. The effect of the central cut of the BARC system on the mode-switching ability of the system is negligible; however, mode switching takes place when comparing the central-cut configurations to the no-cut one. These results show the significance of accurately measuring the central-cut width and how geometric uncertainty may change the dynamical responses and nonlinear properties of the system.

Padilla, Christopher (ORCID:000900033446732X)↗
DOE OSTI2024

Practical and Optimal Sequential Bayesian Experimental Design for Complex Systems Incorporating Human Experimenter Preferences (Final Scientific/Technical Report)

Experiments are indispensable for developing models of complex systems. Carefully designed experiments can provide substantial savings for these expensive data-acquisition opportunities. However, designs based on heuristics are often suboptimal for systems with multiphysics, nonlinear dynamics, and uncertain and noisy environments. Optimal experimental design, while leveraging predictive models, seeks to systematically quantify and maximize the value of experiments. In this project, we focused on the design of multiple experiments, where current approaches are largely suboptimal: batch-design does not adapt to new data acquired during the experiment campaign (no feedback), and greedy/myopic design ignores future dynamics and consequences (no lookahead). We developed the mathematical framework and computational methods for sequential optimal experimental design (sOED) for complex systems. We enabled tractable model-based sOED in a rigorous manner through novel algorithms based on reinforcement learning, and investigated the effects of human experimenters on the design process. Our methods are fully Bayesian, able to quantify and update uncertainty in a principled manner. The traits aimed by our approach—mathematical rigor and optimality, human effects and uncertainty quantification, computational practicality—are crucial for elevating the standards of artificial intelligence (AI) to support decision-making in scientific domains, and contribute toward trust and realistic adoption of AI in experimental design practice.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

First experimental report of nonlinear collimator in the SuperKEKB

SuperKEKB collimator system protects the superconducting final focusing magnets (QC1 and QC2) from quenching caused by beam halo particles. For effective shielding, the minimum aperture of vertical collimators is set smaller than the physical aperture of QC1. However, such collimators become dominant sources of beam coupling impedance, which can limit the maximum beam current due to impedance-driven instabilities. To address this issue, a nonlinear collimator (NLC) was introduced in February 2024. NLC is designed to reduce both beam backgrounds and collimator impedance, thereby suppressing beam instabilities caused by short-range wakefields. In this paper, we describe the results on impedance effects, beam background reduction, and beam lifetime when the nonlinear collimation system is turned on and off. These results validate the NLC concept and offer a promising approach for future high-luminosity colliders.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

A forced Korteweg–de Vries model for nonlinear mixing of oscillations in a dusty plasma

Nonlinear mixing of oscillations in a dusty plasma due to the harmonic time varying modulation of a nonlinear compressional oscillation is analyzed using a simple mathematical model consisting of a forced Korteweg-de Vries equation. An exact analytic solution of this equation is found to exhibit nonlinear mixing in the system. Here, the model solution can be usefully employed to predict the existence of nonlinear mixing of oscillations in a two-dimensional dusty plasma system of a particular experimental configuration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Nonlinear analysis and vibro-impact characteristics of a shaft-bearing assembly

Here this study investigates the nonlinear frequency response of a shaft-bearing assembly with vibro-impacts occurring at the bearing clearances. The formation of nonlinear behavior as system parameters change is examined, along with the effects of asymmetries in the nominal, inherently symmetric system. The primary effect of increasing the forcing magnitude or decreasing the contact gap sizes is the formation of grazing-induced chaotic solution branches occurring over a wide frequency range near each system resonance. The system's nominal setup has very hard contact stiffness and shows no evidence of isolas or superharmonic resonances over the frequency ranges of interest. Moderate contact stiffnesses cause symmetry breaking and introduce superharmonic resonance branches of primary resonances. Even if some primary resonances are not present due to the system's inherent symmetry, their superharmonic resonances still manifest. Branches of quasiperiodic isolas (isolated resonance branches) are also discovered, along with a cloud of isolas near a high-frequency resonance. Parameter asymmetries are found to produce a few significant changes in behavior: asymmetric linear stiffness, contact stiffness, and gap size could affect the behavior of primary resonant frequencies and isolas.

42 ENGINEERING↗
DOE OSTI2026

Suppression of ion temperature gradient modes by Alfvén activity above a drive threshold in DIII-D

A recent study demonstrates that the suppression of ion temperature gradient (ITG) modes can occur during the nonlinear evolution of toroidicity-induced Alfvén eigenmodes (TAEs) in the DIII-D tokamak [Du et al., Phys. Rev. Lett. 135, 265101 (2025)]. In that work, ITG was suppressed by the formation of a narrow, TAE-induced shear flow layer, whose shearing rate exceeded the ITG decorrelation rate. Here, the shear flow arises from an imbalance between Reynolds and Maxwell stress forces, as the TAE departs from the conventional shear Alfvén wave polarization. This follow-up paper systematically identifies the plasma conditions required for robust ITG suppression through a series of comparative experiments. The results show that TAEs routinely suppress ITG turbulence in plasmas with higher local safety factor (q), elevated fast-ion beta, and larger populations of fast ions on passing orbits, i.e., conditions consistent with the substantial TAE drive. Database analysis further reveals the existence of a threshold in TAE drive for ITG suppression. That is, once the drive exceeds a critical value within a favorable q window, the system undergoes the nonlinear bifurcation process, characterized by a sharp increase in TAE saturated amplitudes at a nearly fixed fast ion drive, suppression of ITG turbulence, and the formation of an internal thermal transport barrier at the localized radii.

Du, X. D. [General Atomics, San Diego, CA (United ↗