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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Incorporation of market signals for the optimal design of post combustion carbon capture systems

Recent studies have shown that fossil generators equipped with post-combustion carbon capture (PCC) systems are needed to reduce the cost of deep decarbonization. Such generators need to be flexible and responsive to grid conditions, particularly in a high variable renewable energy (VRE) environment. In this work, we evaluate the net present value (NPV) of retrofitting an existing natural gas combined cycle (NGCC) unit with a flexible PCC system while incorporating market signals from a high VRE grid. We use our industrial partner’s NGCC configuration as representative of existing NGCC units and Svante’s rapid-temperature swing adsorption (TSA) for PCC. Because of its ability to rapidly startup/shutdown and ramp-up/ramp-down, the chosen capture technology is very attractive for load-following operations. For a given set of market signals, we formulate a two-stage stochastic multi-period optimization problem, under the price-taker assumption, to simultaneously optimize the design of the capture system and operation of the entire plant. Rigorous models for the NGCC unit, PCC system, and compression system are developed using commercial process simulators and validated with either plant or vendor data. For computational tractability, we develop surrogate/reduced-order models for use in the optimization problem. The surrogate model for the NGCC plant is constructed by linearizing the rigorous dynamic model at 75% load, while data-driven nonlinear surrogate models for the capture and compression systems are constructed using simulation data from the rigorous models. The optimization problem, formulated as a mixed integer bilinear program, is implemented in the IDAES® integrated platform and solved to global optimality using Gurobi 9.5. Using this formulation, we determine the profitability of retrofitting an existing NGCC unit with the chosen capture system for multiple regions in the U.S. under two scenarios with different carbon prices. Importantly, the results show that the optimal decision strongly depends on the region and on the carbon price, thereby demonstrating the importance of the inclusion of market signals in the design process.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

Multiscale simulations for multi-continuum Richards equations

In this paper, we study a multiscale method for simulating a dual-continuum unsaturated flow problem within complex heterogeneous fractured porous media. Mathematically, each of the dual continua is modeled by a multiscale Richards equation (for pressure head), and these equations are coupled to one another by transfer terms. On its own, Richards equation is already a nonlinear partial differential equation, and it is exceedingly difficult to solve numerically due to the extra nonlinear dependencies involving the soil water. To deal with multiple scales, our strategy is that starting from a microscopic scale, we upscale the coupled system of dual-continuum Richards equations via homogenization by the two-scale asymptotic expansion, to obtain a homogenized system, at an intermediate scale (level). Based on a hierarchical approach, the homogenization’s effective coefficients are computed through solving the arising cell problems. Furthermore, to tackle the nonlinearity, after time discretization, we use Picard iteration procedure for linearization of the homogenized Richards equations. At each Picard iteration, some degree of multiscale still remains from the intermediate level, so we utilize the generalized multiscale finite element method (GMsFEM) combining with a multi-continuum approach, to upscale the homogenized system to a macroscopic (coarse-grid) level. This scheme involves building uncoupled and coupled multiscale basis functions, which are used not only to construct coarse-grid solution approximation with high accuracy but also (with the coupled multiscale basis) to capture the interactions among continua. These prospects and convergence are demonstrated by several numerical results for the proposed method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

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↗
DOE OSTI2021

Simulation of adiabatic quantum computing for molecular ground states

Quantum computation promises to provide substantial speedups in many practical applications with a particularly exciting one being the simulation of quantum many-body systems. Adiabatic state preparation (ASP) is one way that quantum computers could recreate and simulate the ground state of a physical system. In this paper, we explore a novel approach for classically simulating the time dynamics of ASP with high accuracy and with only modest computational resources via an adaptive sampling configuration interaction scheme for truncating the Hilbert space to only the most important determinants. We verify that this truncation introduces negligible error and use this new approach to simulate ASP for sets of small molecular systems and Hubbard models. Furthermore, we examine two approaches to speeding up ASP when performed on quantum hardware: (i) using the complete active space configuration interaction (CASCI) wave function instead of the Hartree–Fock initial state and (ii) a nonlinear interpolation between the initial and target Hamiltonians. We find that starting with a CASCI wave function with a limited active space yields substantial speedups for many of the systems examined, while nonlinear interpolation does not. In additional, we observe interesting trends in the minimum gap location (based on the initial state) as well as how state preparation time can depend on certain molecular properties, such as the number of valence electrons. Importantly, we find that the required state preparation times do not show an immediate exponential wall that would preclude an efficient run of ASP on actual hardware.

Kremenetski, Vladimir↗
DOE OSTI2024

Crosstalk Effects in Microwave SQUID Multiplexed TES Bolometer Readout

Transition-edge sensor (TES) bolometers are broadly used for background-limited astrophysical measurements from the far-infrared to mm-waves. Many planned future instruments require increasingly large detector arrays, but their scalability is limited by their cryogenic readout electronics. Microwave SQUID multiplexing offers a highly capable scaling solution through the use of inherently broadband circuitry, enabling readout of hundreds to thousands of channels per microwave line. As with any multiplexing technique, the channelization mechanism gives rise to electrical crosstalk which must be understood and controlled so as to not degrade the instrument sensitivity. Here, in this study, we explore implications relevant for TES bolometer array applications, focusing in particular on upcoming mm-wave observatories such as the Simons Observatory and AliCPT. We model the relative contributions of the various underlying crosstalk mechanisms, evaluate the difference between fixed-tone and tone-tracking readout systems, and discuss ways in which crosstalk nonlinearity will complicate on-sky measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

OC6 Phase Ib: Validation of the CFD predictions of difference-frequency wave excitation on a FOWT semisubmersible

During the previous OC5 11 project, state-of-the-art mid-fidelity engineering tools for floating wind systems were found to consistently underpredict the nonlinear, low-frequency responses of semisubmersible offshore wind platforms, leading to substantial errors in the structural loads because of low-frequency surge and pitch resonance. To examine this underprediction, a coordinated investigation with computational fluid dynamics (CFD) simulations and model-basin experiments was carried out. Both investigations involved a fixed and simplified OC5-DeepCwind semisubmersible in bichromatic waves. The wave excitations - especially the nonlinear, difference-frequency excitation - on the structure from the CFD simulations were compared to the experimental measurements for validation, with uncertainty analyses for both the experimental and the CFD results. Furthermore, the wave excitations on each column of the semisubmersible were measured separately in the experiment, allowing the validation of the CFD results to be done on a per-column basis. Overall, the CFD predictions of the difference-frequency excitations agree with the experimental measurements, suggesting the CFD solutions can be used as a reference for tuning and improving the engineering-level tools and can provide a means to better understand the underprediction at low frequencies.

16 TIDAL AND WAVE POWER↗
DOE OSTI2023

Performance characterization of freeform finished surfaces of potassium dihydrogen phosphate using fluid jet polishing with a nonaqueous slurry

Abstract Potassium dihydrogen phosphate (KDP) and its deuterated analog (DKDP) are unique nonlinear optical materials for high power laser systems. They are used widely for frequency conversion and polarization control by virtue of the ability to grow optical-quality crystals at apertures suitable for fusion-class laser systems. Existing methods for freeform figuring of KDP/DKDP optics do not produce surfaces with sufficient laser-induced–damage thresholds (LIDT’s) for operation in the ultraviolet portion of high-peak-power laser systems. In this work, we investigate fluid jet polishing (FJP) using a nonaqueous slurry as a sub-aperture finishing method for producing freeform KDP surfaces. This method was used to selectively polish surface areas to different depths on the same substrate with removals ranging from 0.16 μm to 5.13 μm. The finished surfaces demonstrated a slight increase in roughness as the removal depth increased along with a small number of fracture pits. Laser damage testing with 351 nm, 1 ns pulses demonstrated excellent surface damage thresholds, with the highest values in areas devoid of fracture pits. This work demonstrates, for the first time, a method that enables fabrication of a waveplate that provides tailored polarization randomization that can be scaled to meter-sized optics. Furthermore, this method is based on FJP technology that incorporates a nonaqueous slurry specially designed for use with KDP. This novel nonaqueous FJP process can be also used for figuring other types of materials that exhibit similar challenging inherent properties such as softness, brittleness, water-solubility, and temperature sensitivity.

36 MATERIALS SCIENCE↗
DOE OSTI2025

DAmodel: hierarchical Bayesian modelling of DA white dwarfs for spectrophotometric calibration

We use hierarchical Bayesian modelling to calibrate a network of 32 all-sky faint DA white dwarf (DA WD) spectrophotometric standards (⁠16.5 < V , 19.5⁠) alongside three CALSPEC standards, from 912 Å to 32 μm. The framework is the first of its kind to jointly infer photometric zero points and WD parameters (surface gravity log g⁠, effective temperature T eff ⁠, extinction A V ⁠, dust relation parameter R V ) by simultaneously modelling both photometric and spectroscopic data. We model panchromatic Hubble Space Telescope Wide Field Camera 3 (HST/WFC3) UVIS and IR photometry, HST/STIS UV spectroscopy, and ground-based optical spectroscopy to sub-per cent precision. Photometric residuals for the sample are the lowest yet yielding < 0.004 mag RMS on average from the UV to the NIR, achieved by jointly inferring time-dependent changes in system sensitivity and WFC3/IR count-rate nonlinearity. Our GPU-accelerated implementation enables efficient sampling via Hamiltonian Monte Carlo, critical for exploring the high-dimensional posterior space. The hierarchical nature of the model enables population analysis of intrinsic WD and dust parameters. Inferred spectral energy distributions from this model will be essential for calibrating the James Webb Space Telescope as well as next-generation surveys, including Vera Rubin Observatory’s Legacy Survey of Space and Time and the Nancy Grace Roman Space Telescope.

methods: statistical↗
DOE OSTI2025

Simulating quantum-classical interfaces via the Lindblad master equation

In hybrid quantum systems, the interface between quantum and classical domains is essential for the generation, control, and measurement of quantum states. Quantum-classical interfaces (QCIs) are ubiquitous in devices such as optical modulators, quantum sensors, and signal processors, where classical signals influence quantum dynamics. In this paper, we employ the Lindblad master equation to simulate the evolution of a quantum system interacting with a classical control system. Our model captures both linear and nonlinear interactions by incorporating first- and second-order susceptibilities, and it quantifies the influence of externally applied control parameters on decoherence and state evolution. As an illustrative example, we analyze an optical modulator and demonstrate how variations in material response and drive conditions affect photon statistics, coherence, and phase-space distributions. In conclusion, the findings offer a path to an all-encompassing model for understanding and optimizing QCIs, with wide-ranging implications for the performance, design, and robustness of next-generation quantum devices.

Quantum engineering↗
DOE OSTI2021

Exact results for nonlinear Drude weights in the spin-1/2 XXZ chain

Nonlinear Drude weight (NLDW) is a generalization of the linear Drude weight, which characterizes the nonlinear transport in quantum many-body systems. We investigate these weights for the spin-1/2 XXZ chain in the critical regime. The effects of the Dzyaloshinskii-Moriya interaction and an external magnetic field are also studied. Solving the Bethe equations numerically, we obtain these weights for very large system sizes and identify parameter regimes where the weights diverge in the thermodynamic limit. These divergences appear in all the orders studied in this Letter and can be regarded as a generic feature of the NLDWs. We study the origin of these divergences and reveal that they result from nonanalytic finite-size corrections to the ground-state energy. As a result, we compute closed-form expressions for several weights in the thermodynamic limit and find excellent agreement with the numerical results.

1-dimensional systems↗
DOE OSTI2025

Deep Reinforcement Learning Based Control of Wind Turbines for Fast Frequency Response

In order to fulfill vital auxiliary grid services, such as load regulation, spin and non-spin reserve provision, and frequency support during emergencies, there is often a requirement for certain wind farms to operate in de-loaded modes. Leveraging the swift response capabilities of wind farms, this study demonstrates that reserving power in de-loaded modes can significantly enhance power grid stability and reliability during system contingencies. Controlling wind farms optimally for frequency support is intricate due to the nonlinearity of models and controllers and the complexity of wind farm interactions with power systems. Here, to address this challenge, this paper introduces a novel approach that integrates wind turbines into reinforcement learning-based solutions for frequency response. This innovative methodology utilizes the state-of-the-art reinforcement learning algorithm known as the surrogate-gradient-based evolutionary strategy. The proposed learning-based algorithm provides continuous control of wind farm output to rapidly stabilize system frequency and prevent unnecessary trips of under-frequency load shedding relays. To facilitate efficient training, parallel computing techniques are employed. The proposed methodology is evaluated on a modified IEEE-39 bus system, and simulation results reveal its efficacy in reliably supporting power system frequency and preventing the need for unnecessary load shedding.

Gao, Wei [Argonne National Laboratory (ANL), Argon↗
DOE OSTI2021

Convex Optimization of Integrated Power-Gas Energy Flow Model With Applications to Probabilistic Energy Flow

Energy flow calculation is a fundamental problem of the integrated power and gas system (IPGS) operation and planning. However, the nonlinear gas flow model introduces major challenges to the energy flow calculation. In this paper, we propose a tractably convex optimization model to solve the energy flow problem in IPGSs. It is demonstrated that the proposed optimization model has the same optimal solution as the original nonlinear steady energy flow model. Also, piecewise linearization is adopted to tightly linearize the nonlinear objective function of the model, which transforms the formulated convex optimization into a linear program one. Thus, the computation complexity of the proposed energy flow model is significantly reduced as compared with the existing methods. In addition, the proposed model can be extended to probabilistic energy flow estimation. Extensive case studies are conducted to validate the effectiveness of the proposed energy flow model using three IPGSs.

42 ENGINEERING↗
DOE OSTI2024

Heterosynaptic plasticity in memristive and memcapacitive lipid bilayers: A snapshot review

Synaptic plasticity refers to activity-dependent synaptic strengthening or weakening between neurons. It is usually associated with homosynaptic plasticity, which refers to a synaptic junction controlled by interactions between specific neurons. Heterosynaptic plasticity, on the other hand, lacks this specificity. It involves much larger populations of synapses and neurons and can be associated with changes in synaptic strength due to nonlocal alterations in the ambient electrochemical environment. Here, this paper presents specific examples demonstrating how variations in the ambient electrochemical environment of lipid membranes can impact the nonlinear dynamical behaviors of memristive and memcapacitive systems in droplet interface bilayers (DIBs). Examples include the use of pH as a modulatory factor that alters the voltage-dependent memristive behavior of alamethicin ion channels in DIB lipid bilayers, and the discovery of long-term potentiation (LTP) in a lipid bilayer-only system after application of electrical stimulation protocols.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Advanced Features and Recent Developments in the WEC-Sim Open-Source Design Tool: Preprint

WEC-Sim is an open-source tool developed in Matlab/Simulink to simulate wave energy converters. WEC-Sim has been developed jointly by the National Renewable Energy Laboratory (NREL) and Sandia National Laboratories (SNL) since its initial release in 2014. WEC-Sim simulations are performed in the time domain based on the radiation and diffraction method using hydrodynamics coefficients derived from boundary element method (BEM)-based frequency-domain potential flow solvers (e.g., WAMIT, NEMOH, Capytaine, or ANSYS-AQWA). Within this level of modeling fidelity, WEC-Sim can handle floating body hydrodynamics, mechanical and electrical power generation methods, advanced control implementation, mooring systems, and other unique applications such as desalination. WEC-Sim functionalities are developed using Simulink blocks and Matlab functions. WEC-Sim is capabilities include simulation of different wave conditions, constraints and PTOs, body dynamics, mooring systems, among other advanced features such as nonlinear hydrodynamics, nonlinear PTOs, body-to-body interactions, etc.

hydrodynamics modeling↗
DOE OSTI2023

Analytical Large-Signal Modeling of Inverter-Based Microgrids With Koopman Operator Theory for Autonomous Control

The microgrid (MG) plays a crucial role in the energy transition, but its nonlinearity presents a significant challenge for large-signal power systems studies in the electromagnetic transient (EMT) time scale. In this paper, we develop a large-signal linear MG model that considers the detailed dynamics of the primary and zero-control levels based on the Koopman operator (KO) theory. Firstly, a set of observable functions is carefully designed to capture the nonlinear dynamics of the MG. The corresponding linear KO is then analytically derived based on these observables, resulting in the linear representation of the original nonlinear MG with observables as the new coordinate. The influence of external input on the system dynamics is also considered during the derivation, enabling control of the MG. We solve the voltage control problem using the traditional linear quadratic integrator (LQI) method to demonstrate that textbook linear control techniques can accurately control the original nonlinear MG via the developed KO-linearized MG model. Our proposed KO linearization method is generic and can be easily extended for different control objectives and MG structures using our analytical derivation procedure. Here, we validate the effectiveness of our methodology through various case studies.

25 ENERGY STORAGE↗
DOE OSTI2022

A Realistic Theory of Quantum Measurement

Abstract We propose that the ontic understanding of quantum mechanics can be extended to a fully realistic theory that describes the evolution of the wavefunction at all times, including during a measurement. In such an approach the wave equation should reduce to the standard wave equation when there is no measurement, and describe state reduction when the system is measured. The general wave equation must be nonlinear and nonlocal, and we require it to be time-symmetric; consequently, this approach is not a new interpretation but a new theory. The wave equation is an integrodifferential equation (IDE). The time symmetry requirement leads to a retrocausal approach, in which the wave equation is solved subject to initial and final conditions to determine history at intermediate times. We propose that different outcomes from (apparently) identically prepared experiments may result from uncontrolled parameters; both the nonlocality and the retrocausality of the theory imply that Bell’s Theorem cannot rule out such “hidden variables.” Beginning with Hamilton’s principle, we demonstrate the construction of such a theory by replacing the action with a functional designed to give rise to a nonlinear, nonlocal IDE as the wave equation. This IDE reduces to the standard wave equation (a differential equation) in the absence of a measurement, but exhibits state reduction to a single eigenvalue when the system interacts with another system with the properties of a measurement apparatus. We demonstrate several desirable features of this theory; for other properties we indicate their plausibility and possible avenues to a proof.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Nonlinear Photovoltaic Degradation Rates: Modeling and Comparison Against Conventional Methods

Although common practice for estimating photovoltaic (PV) degradation rate (RD) assumes a linear behavior, field data have shown that degradation rates are frequently nonlinear. This article presents a new methodology to detect and calculate nonlinear RD based on PV performance time-series from nine different systems over an eight-year period. Prior to performing the analysis and in order to adjust model parameters to reflect actual PV operation, synthetic datasets were utilized for calibration purposes. A change-point analysis is then applied to detect changes in the slopes of PV trends, which are extracted from constructed performance ratio (PR) time-series. Once the number and location of change points is found, the ordinary least squares method is applied to the different segments to compute the corresponding rates. The obtained results verified that the extracted trends from the PR time-series may not always be linear and therefore, “nonconventional” models need to be applied. All thin-film technologies demonstrated nonlinear behavior whereas nonlinearity detected in the crystalline silicon systems is thought to be due to a maintenance event. A comparative analysis between the new methodology and other conventional methods demonstrated levelized cost of energy differences of up to 6.14%, highlighting the importance of considering nonlinear degradation behavior

14 SOLAR ENERGY↗