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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Data-Centric Approach to Capture Non-Polynomial Nonlinear Dynamics

We propose an analytical construction of observable functions in the extended dynamic mode decomposition (EDMD) algorithm. EDMD is a numerical method for approximating the spectral properties of the Koopman operator. The choice of observable functions is fundamental for applying EDMD to nonlinear problems arising in systems and control. Existing methods either start from a set of dictionary functions and look for the subset that best fits the underlying nonlinear dynamics or rely on machine learning algorithms to “learn” observable functions. Conversely, in this paper, we start from the dynamical system model and lift it through the Lie derivatives, rendering it into a polynomial form. This proposed transformation into a polynomial form is exact and provides an adequate set of observable functions. The strength of the proposed approach is its applicability to a broader class of nonlinear dynamical systems, particularly those with nonpolynomial functions and compositions thereof. Moreover, it retains the physical interpretability of the underlying dynamical system and can be readily integrated into existing numerical libraries. We demonstrate the proposed approach with an application to electric power systems. The modeled system consists of a single generator connected to an infinite bus, where nonlinear terms include sine and cosine functions. The results demonstrate the effectiveness of the proposed procedure in off-attractor nonlinear dynamics for estimation and prediction; the observable functions obtained from the proposed construction outperform methods that use dictionary functions comprising monomials or radial basis functions.

extended dynamic mode decomposition↗
DOE OSTI2025

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)↗
DOE OSTI2022

Operability of a Power System with Synchronous Condensers and Grid-Following Inverters

Growing shares of inverter-based resources gener-ally correlate with a reduction in inertia as synchronous genera-tors are displaced. Along the path to high shares of inverter-based resources, in particular with only conventional grid-following inverter controls, a proffered solution is the use of synchronous condensers as a source of inertia to maintain the frequency-power balance relationship. An outstanding question is whether only grid-following inverters and synchronous condensers yield a viable power system; i.e., all frequency response is derived from inverters. A validated electromagnetic transient model of the Maui system, with many nonlinear elements such as load-shedding, line tripping, and inverter ride through criteria disabled, is used to investigate the stability of such a system. Two types of perturbations were applied, a 15% generation loss and a fault event, and it was found that the system remains stable.

grid-following inverters↗
DOE OSTI2025

Flat nonlinear optics with intersubband polaritonic metasurfaces

Nonlinear intersubband polaritonic metasurfaces produce some of the strongest second- and third-order nonlinear optical responses reported for condensed matter systems at infrared frequencies. These metasurfaces are fabricated as two-dimensional arrays of nanoresonators from multi-quantum-well semiconductor heterostructures, designed to produce strong nonlinear responses associated with intersubband transitions. By optimally coupling the optical modes of the nanoresonators to vertically polarized intersubband transitions in semiconductor heterostructures, one can boost the nonlinear response associated with intersubband transitions, make intersubband transitions interact with free-space radiation at normal incidence, and hence produce optically thin flat nonlinear optical elements compatible with free-space optical setups. As a result of the strong nonlinear response in these systems, significant nonlinear conversion efficiencies (>0.1 %) can be attained in deeply subwavelength optical films using modest pumping intensities of only 10–100 kW/cm 2 . Subwavelength metasurface thickness relaxes phase-matching constraints limiting the operation of bulk nonlinear crystals. Furthermore, the amplitude and phase of the nonlinear optical response in intersubband polaritonic metasurfaces can be tailored for a specific pump wavelength and a nonlinear process of interest through the co-optimization of quantum engineering of electron states in semiconductor heterostructures and photonic engineering of the metasurface nanoresonators design. Additionally, an applied voltage can dynamically control the amplitude and phase of the nonlinear optical response at a nanoresonator level. Here, we review the current state of the art in this rapidly expanding field, focusing on nonlinear processes supporting second-harmonic generation, saturable absorption, and optical power limiting.

flat optics↗
DOE OSTI2021

Network-Level Optimization for Unbalanced Power Distribution System: Approximation and Relaxation

The nonlinear programming (NLP) problem to solve distribution-level optimal power flow (D-OPF) poses convergence issues and does not scale well for unbalanced distribution systems. The existing scalable D-OPF algorithms either use approximations that are not valid for an unbalanced power distribution system, or apply relaxation techniques to the nonlinear power flow equations that do not guarantee a feasible power flow solution. In this paper, we propose scalable D-OPF algorithms that simultaneously achieve optimal and feasible solutions by solving multiple iterations of approximate, or relaxed, D-OPF subproblems of low complexity. The first algorithm is based on a successive linear approximation of the nonlinear power flow equations around the current operating point, where the D-OPF solution is obtained by solving multiple iterations of a linear programming (LP) problem. The second algorithm is based on the relaxation of the nonlinear power flow equations as conic constraints together with directional constraints, which achieves optimal and feasible solutions over multiple iterations of a second-order cone programming (SOCP) problem. Finally, it is demonstrated that the proposed algorithms are able to reach an optimal and feasible solution while significantly reducing the computation time as compared to an equivalent NLPD-OPF model for the same distribution system.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Massively parallel and universal approximation of nonlinear functions using diffractive processors

Nonlinear computation is essential for a wide range of information processing tasks, yet implementing nonlinear functions using optical systems remains a challenge due to the weak and power-intensive nature of optical nonlinearities. Overcoming this limitation without relying on nonlinear optical materials could unlock unprecedented opportunities for ultrafast and parallel optical computing systems. Here, we demonstrate that large-scale nonlinear computation can be performed using linear optics through optimized diffractive processors composed of passive phase-only surfaces. In this framework, the input variables of nonlinear functions are encoded into the phase of an optical wavefront—e.g., via a spatial light modulator (SLM)—and transformed by an optimized diffractive structure with spatially varying point-spread functions to yield output intensities that approximate a large set of unique nonlinear functions–all in parallel. We provide proof establishing that this architecture serves as a universal function approximator for an arbitrary set of bandlimited nonlinear functions, also covering wavelength-multiplexed nonlinear functions as well as multi-variate and complex-valued functions that are all-optically cascadable. Our analysis also indicates the successful approximation of typical nonlinear activation functions commonly used in neural networks, including the sigmoid, tanh, ReLU (rectified linear unit), and softplus. We numerically demonstrate the parallel computation of one million distinct nonlinear functions, accurately executed at wavelength-scale spatial density at the output of a diffractive optical processor. Furthermore, we experimentally validated this framework using in situ optical learning and approximated 35 unique nonlinear functions in a single shot using a compact setup consisting of an SLM and an image sensor. These results establish diffractive optical processors as a scalable platform for massively parallel universal nonlinear function approximation, paving the way for new capabilities in analog optical computing based on linear materials.

Rahman, Md Sadman Sakib [University of California,↗
DOE OSTI2025

Towards a Robust Adaptive Digital Twin for Fusion Applications

The development of a digital twin system for fusion applications is essential for enhancing the prediction, analysis, and optimization of complex plasma processes. Machine learning (ML), particularly deep learning has demonstrated strong capabilities in modeling such highly nonlinear and intricate systems. However, two critical challenges limit the deployment of deep learning-based digital twins: Uncertainty Quantification (UQ) and data drift. UQ is vital for ensuring trustworthy predictions, especially in decision-support scenarios. Additionally, data-driven models are often sensitive to changes in the underlying data distribution, such as shot-to-shot variations in fusion experiments, which can lead to performance degradation over time. To address these challenges, we are developing an uncertainty-aware, adaptive digital twin framework. Our approach incorporates deep learning models enhanced with Gaussian Process approximations for predictive uncertainty estimation, coupled with an online learning mechanism that enables continuous model adaptation to new experimental data. This adaptive capability allows the data driven models to respond effectively to evolving plasma behaviors and equipment conditions. Specifically, to mitigate the effects of shot-to-shot drift, our system updates itself incrementally as new data becomes available, improving both robustness and fidelity. Our vision is to evolve this data driven model into a self-sustaining digital twin system that leverages UQ based feedback to continuously refine itself and potentially support real-time decision making. This presentation will cover a brief background on uncertainty quantification for ML, our ongoing effort on development of UQ capabilities for ML, our data science pipeline from data collection to model development and analysis and online learning framework for modeling coil deflection at DIII-D. I will also briefly touch upon opportunities and challenges in development of digital twin framework.

Sammuli, Brian [General Atomics]↗
DOE OSTI2023

Three-dimensional broadband nonlinear phased array imaging

System include an ultrasonic transducer configured to couple to a nondestructive testing (NDT) sample and configured to produce and direct an ultrasonic probe wave at a selected frequency into a subsurface region of the NDT sample, a 3D laser scanning vibrometer configured to direct a detection beam in a scan area on a surface of the NDT sample and to receive a return beam from the scan area, and to detect, based on the return beam, a 3D motion of the surface across a wideband frequency range, and a processor, and a memory configured with instructions that, when executed by the processor, cause the processor to produce sub-surface image data of the NDT sample at multiple harmonics of the selected frequency in the wideband frequency range based on the detected 3D surface motion, wherein the sub-surface image data describes a nonlinear defect response produced in the NDT sample by interaction of the ultrasonic probe wave with the subsurface region.

Remillieux, Marcel↗
DOE OSTI2021

Spatially co-registered wide-field nonlinear optical imaging of living and complex biosystems in a total internal reflection geometry

Nonlinear optical microscopy that leverages an objective based total internal reflection (TIR) excitation scheme is an attractive means for rapid, wide-field imaging with enhanced surface sensitivity. Through select combinations of distinct modalities, one can, in principle, access complementary chemical and structural information for various chemical species near interfaces. Here, we report a successful implementation of such a wide-field nonlinear optical microscope system, which combines coherent anti-Stokes Raman scattering (CARS), two-photon fluorescence (TPF), second harmonic generation (SHG), and sum frequency generation (SFG) modalities on the same platform. The intense optical fields needed to drive these high order nonlinear optical processes are achieved through the use of femtosecond pulsed light in combination with the intrinsic field confinement induced by TIR over a large field of view. The performance of our multimodal microscope was first assessed through the experimental determination of its chemical fidelity, intensity and polarization dependences, and spatial resolution using a set of well-defined model systems. Subsequently, its unique capabilities were validated through imaging complex biological systems, including Hydrangea quercifolia pollen grains and Pantoea sp. YR343 bacterial cells. Specifically, the spatial distribution of different molecular groups in the former was visualized via vibrational contrast mechanisms of CARS, whereas co-registered TPF imaging allowed the identification of spatially localized intrinsic fluorophores. We further demonstrate the feasibility of our microscope for wide-field CARS imaging on live cells through independent characterization of cell viability using spatially co-registered TPF imaging. Overall, this approach to TIR enabled wide-field imaging is expected to provide new insights into bacterial strains and their interactions with other species in the rhizosphere in a time-resolved and chemically selective manner.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Insights into ENSO Diversity from an Intermediate Coupled Model. Part II: Role of Nonlinear Dynamics and Stochastic Forcing

Abstract In this study, we investigate how a single leading linear El Niño–Southern Oscillation (ENSO) mode, as studied in Part I, leads to the irregular coexistence of central Pacific (CP) and eastern Pacific (EP) ENSO, a phenomenon known as ENSO spatiotemporal diversity. This diversity is fundamentally generated by deterministic nonlinear pathways to chaos via the period-doubling route and, more prevailingly, the subharmonic resonance route with the presence of a seasonally varying basic state. When residing in the weakly nonlinear regime, the coupled system sustains a weak periodic oscillation with a mixed CP/EP pattern as captured by the linear ENSO mode. With a stronger nonlinearity effect, the ENSO behavior experiences a period-doubling bifurcation. The single ENSO orbit splits into coexisting CP-like and EP-like ENSO orbits. A sequence of period-doubling bifurcation results in an aperiodic oscillation featuring irregular CP and EP ENSO occurrences. The overlapping of subharmonic resonances between ENSO and the seasonal cycle allows this ENSO irregularity and diversity to be more readily excited. In the strongly nonlinear regime, the coupled system is dominated by regular EP ENSO. The deterministic ENSO spatiotemporal diversity is thus confined to a relatively narrow range corresponding to a moderately unstable ENSO mode. Stochastic forcing broadens this range and allows ENSO diversity to occur when the ENSO mode is weakly subcritical. A close relationship among a weakened mean zonal temperature gradient, stronger ENSO activity, and more (fewer) occurrences of EP (CP) ENSO is noted, indicating that ENSO–mean state interaction may yield ENSO regime modulations on the multidecadal time scale.

Geng, Licheng↗
DOE OSTI2021

On the mechanism of ionization oscillations in Hall thrusters

Low-frequency ionization oscillations involving plasma and neutral density (breathing modes) are the most violent perturbations in Hall thrusters for electric propulsion. Because of its simplicity, the zero-dimensional (0D) predator–prey model of two nonlinearly coupled ordinary differential equations for plasma and neutral density has often been used for the characterization of such oscillations and scaling estimates. In this work, we investigate the properties of its continuum analog, the one-dimensional (1D) system of two nonlinearly coupled equations in partial derivatives (PDEs) for plasma and neutral density. This is a more general model, of which the standard 0D predator–prey model is a special limit case. We show that the 1D model is stable and does not show any oscillations for the boundary conditions relevant to Hall thrusters and the uniform ion velocity. We then propose a reduced 1D model based on two coupled PDEs for plasma and neutral densities that is unstable and exhibit oscillations if the ion velocity profile with the near-the-anode back-flow (toward the anode) region is used. Comparisons of the reduced model with the predictions of the full model that takes into account the self-consistent plasma response show that the main properties of the breathing mode are well captured. In particular, it is shown that the frequency of the breathing mode oscillations is weakly dependent on the final ion velocity but shows a strong correlation with the width of the ion back-flow region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A quantum Monte Carlo study of systems with effective core potentials and node nonlinearities

In this report we study beryllium dihydride (BeH2) and acetylene (C2H2) molecules using real-space diffusion Monte Carlo (DMC) method. The molecules serve as perhaps the simplest prototypes that illustrate the difficulties with biases in the fixed-node DMC calculations that might appear with the use of effective core potentials (ECPs) or other nonlocal operators. This is especially relevant for the recently introduced correlation consistent ECPs (ccECPs) for 2s2p elements. Corresponding ccECPs exhibit deeper potential functions due to higher fidelity to all-electron counterparts, which could lead to larger local energy fluctuations. We point out that the difficulties stem from issues that are straightforward to address by upgrades of basis sets, use of T-moves for nonlocal terms, inclusion of a few configurations into the trial function and similar. The resulting accuracy corresponds to the ccECP target (chemical accuracy) and it is in consistent agreement with independent correlated calculations. Further possibilities for upgrading the reliability of the DMC algorithm and considerations for better adapted and more robust Jastrow factors are discussed as well.

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

Extended real-time voltage instability identification method based on synchronized phasor measurements

This paper presents an extended adaptive approach designed to accurately estimate the Thévenin equivalent parameters using phasor measurements at a given bus for measurements lying in any of four quadrants of the PQ-plane. The improvement is achieved by using a new condition to properly update the estimated parameters after an initial guess. Based on an adaptive philosophy, the proposed approach can correctly account for the intrinsic nonlinearities of a power system, can provide a real-time estimation of Thévenin parameters, and does not require network topology knowledge. The method is validated using the Kundur 2-area system, showing estimation improvements compared to the current adaptive approach and the classical recursive least-squares method. The proposed approach is able to estimate both sides of the system with respect to the measurement bus. In addition, a data-driven voltage stability index is developed. To illustrate the performance of the proposed approach in a larger power system, a voltage stability assessment is carried out on the IEEE 39-bus system, considering the action of overexcitation limiters of generators and nonlinear loads. The proposed approach is suitable for applications that require an accurate Thévenin equivalent estimation in real-time, such as for voltage stability assessment. The new approach provides a reliable tool for the system operators to make proper and timely decisions.

42 ENGINEERING↗
DOE OSTI2022

Nonlinear Optimization and the Modeling of Energy Systems [Slides]

Energy delivery systems are critical for the function of modern society. (Up to) continental-scale engineered systems move energy from source points to consumers. These systems are increasingly complex and interconnected.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

A 1D–0D–3D coupled model for simulating blood flow and transport processes in breast tissue

In this work, we present mixed dimensional models for simulating blood flow and transport processes in breast tissue and the vascular tree supplying it. These processes are considered, to start from the aortic inlet to the capillaries and tissue of the breast. Large variations in biophysical properties and flow conditions exist in this system necessitating the use of different flow models for different geometries and flow regimes. In total, we consider four different model types. First, a system of 1D nonlinear hyperbolic partial differential equations (PDEs) is considered to simulate blood flow in larger arteries with highly elastic vessel walls. Second, we assign 1D linearized hyperbolic PDEs to model the smaller arteries with stiffer vessel walls. The third model type consists of ODE systems (0D models). It is used to model the arterioles and peripheral circulation. Finally, homogenized 3D porous media models are considered to simulate flow and transport in capillaries and tissue within the breast volume. Sink terms are used to account for the influence of the venous and lymphatic systems. Combining the four model types, we obtain two different 1D–0D–3D coupled models for simulating blood flow and transport processes: The first model results in a fully coupled 1D–0D–3D model covering the complete path from the aorta to the breast combining a generic arterial network with a patient specific breast network and geometry. The second model is a reduced one based on the separation of the generic and patient specific parts. The information from a calibrated fully coupled model is used as inflow condition for the patient specific sub-model allowing a significant computational cost reduction. Several numerical experiments are conducted to calibrate the generic model parameters and to demonstrate realistic flow simulations compared to existing data on blood flow in the human breast and vascular system. Moreover, we use two different breast vasculature and tissue data sets to illustrate the robustness of our reduced sub-model approach.

multiscale blood flow models↗
DOE OSTI2026

Latent space dynamics identification for interface tracking with application to shock-induced pore collapse

Capturing sharp, evolving interfaces remains a central challenge in reduced-order modeling, especially when data is limited and the system exhibits localized nonlinearities or discontinuities. Here, we propose LaSDI-IT (Latent Space Dynamics Identification for Interface Tracking), a data-driven framework that combines low-dimensional latent dynamics learning with explicit interface-aware encoding to enable accurate and efficient modeling of physical systems involving moving material boundaries. At the core of LaSDI-IT is a revised autoencoder architecture that jointly reconstructs the physical field and an indicator function representing material regions or phases, allowing the model to track complex interface evolution without requiring detailed physical models or mesh adaptation. The latent dynamics are learned through linear regression in the encoded space and generalized across parameter regimes using Gaussian process interpolation with greedy sampling. We demonstrate LaSDI-IT on the problem of shock-induced pore collapse in high explosives, a process characterized by sharp temperature gradients and dynamically deforming pore geometries. The method achieves relative prediction errors below 9% across the parameter space, accurately recovers key quantities of interest such as pore area and hot spot formation, and matches the performance of dense training with only half the data. This latent dynamics prediction was 10 6 times faster than the conventional high-fidelity simulation, proving its utility for multi-query applications. These results highlight LaSDI-IT as a general, data-efficient framework for modeling discontinuity-rich systems in computational physics, with potential applications in multiphase flows, fracture mechanics, and phase change problems.

Gaussian process↗
DOE OSTI2020

Electronic Population Transfer via Impulsive Stimulated X-Ray Raman Scattering with Attosecond Soft-X-Ray Pulses

Free electron lasers provide a source of X-ray pulses short enough and intense enough to drive nonlinearities in molecular systems. Impulsive interactions driven by these X-ray pulses provide a way to create and probe valence electron motions with high temporal and spatial resolution. Observing these electronic motions is crucial to understand the role of electronic coherence in chemical processes. A simple nonlinear technique for probing electronic motion, impulsive stimulated X-ray Raman scattering (ISXRS), involves a single impulsive interaction to produce a coherent superposition of electronic states. In this work, we demonstrate electronic population transfer via ISXRS using broad bandwidth (5.5 eV full width at half maximum) attosecond X-ray pulses produced by the Linac Coherent Light Source (LCLS). The impulsive excitation is resonantly enhanced by the oxygen 1s → 2π* resonance of nitric oxide (NO), and excited state neutral molecules are probed with a time-delayed UV laser pulse.

74 ATOMIC AND MOLECULAR PHYSICS↗