Search NASA⌕ Search

SEARCH · Search NASA

Results for “nonlinear programming”

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 217 records · Page 12

Modeling and Detection of Future Cyber-Enabled DSM Data Attacks

Demand-Side Management (DSM) is an essential tool to ensure power system reliability and stability. In future smart grids, certain portions of a customer’s load usage could be under the automatic control of a cyber-enabled DSM program, which selectively schedules loads as a function of electricity prices to improve power balance and grid stability. In this scenario, the security of DSM cyberinfrastructure will be critical as advanced metering infrastructure and communication systems are susceptible to cyber-attacks. Such attacks, in the form of false data injections, can manipulate customer load profiles and cause metering chaos and energy losses in the grid. The feedback mechanism between load management on the consumer side and dynamic price schemes employed by independent system operators can further exacerbate attacks. To study how this feedback mechanism may worsen attacks in future cyber-enabled DSM programs, we propose a novel mathematical framework for (i) modeling the nonlinear relationship between load management and real-time pricing, (ii) simulating residential load data and prices, (iii) creating cyber-attacks, and (iv) detecting said attacks. In this framework, we first develop time-series forecasts to model load demand and use them as inputs to an elasticity model for the price-demand relationship in the DSM loop. This work then investigates the behavior of such a feedback loop under intentional cyber-attacks. We simulate and examine load-price data under different DSM-participation levels with three types of random additive attacks: ramp, sudden, and point attacks. We conduct two investigations for the detection of DSM attacks. The first studies a supervised learning approach, with various classification models, and the second studies the performance of parametric and nonparametric change point detectors. Results conclude that higher amounts of DSM participation can exacerbate ramp and sudden attacks leading to better detection of such attacks, especially with supervised learning classifiers. We also find that nonparametric detection outperforms parametric for smaller user pools, and random point attacks are the hardest to detect with any method.

97 MATHEMATICS AND COMPUTING↗

First Measurements of Nonlinear Decoherence in the IOTA Ring

The Integrable Optics Test Accelerator (IOTA), at Fermi National Laboratory is aimed at testing nonlinear optics for the next generation of high intensity rings. Through use of a special magnetic element the ring is designed to induce a large tune spread with amplitude while maintaining integrable motion. This will allow for the suppression of instabilities in high-intensity beams without significant reduction in dynamic aperture. One important aspect of this is the nonlinear decoherence that occurs when a beam is injected off axis or receives a transverse kick while circulating in the ring. This decoherence has been studied in detail, with simulations, for protons in IOTA both with and without space-charge. However, it has yet to be demonstrated experimentally. During the first phase of the IOTA experimental program, the ring is operated with 100 MeV electrons, allowing for the study of nonlinear optics without the complications introduced by space charge. Here we present measurements taken during the IOTA commissioning, and an analysis of the results.

43 PARTICLE ACCELERATORS↗

Sequence of polyhedral relaxations for nonlinear univariate functions

Here, given a nonlinear, univariate, bounded, and differentiable function f(x), this article develops a sequence of Mixed Integer Linear Programming (MILP) and Linear Programming (LP) relaxations that converge to the graph of f(x) and its convex hull, respectively. Theoretical convergence of the sequence of relaxations to the graph of the function and its convex hull is established. For nonlinear non-convex optimization problems, the relaxations presented in this article can be used to construct tight MILP and LP relaxations. These MILP and the LP relaxations can also be used with MILP-based and spatial branch-and-bound based global optimization algorithms, respectively.

42 ENGINEERING↗

Bi-Level Volt-VAR Optimization to Coordinate Smart Inverters with Voltage Control Devices

Conservation voltage reduction (CVR) uses voltage control to reduce customer power demand and help improve the operational efficiency of the electric power distribution system. With the advent of advanced distribution management systems (ADMS), the underlying data and control-rich platform can be utilized to enable a centralized Volt-VAR optimization (VVO) scheme that includes smart inverter controls along with legacy voltage control devices such as voltage regulators and capacitor banks to ensure a maximum CVR benefit. The objective of this paper to present a VVO approach that controls system’s legacy voltage control devices and coordinates their operation with smart inverter control. An optimal power flow (OPF) formulation is proposed by developing linear and nonlinear power flow approximations for a three-phase unbalanced electric power distribution system. A bi-level VVO approach is proposed where, Level-1 optimizes the control of legacy devices and smart inverters using a linear approximate three-phase power flow. In Level-2, the control parameters for smart inverters are adjusted to obtain an optimal and feasible solution by solving the approximate nonlinear OPF model. Level-1 is modeled as a Mixed Integer Linear Program (MILP) while Level-2 as a Quadratically Constrained Quadratic Program (QCQP). The proposed approach is validated using IEEE test cases and the results demonstrate the applicability of the framework in achieving the CVR objective.

Jha, Rahul Ranjan↗

Cultural Perspectives of the 2016 Nonlinear Mechanics and Dynamics Research Institute

This report analyzes the results of a study on culture and its capability to influence research. The study occurred during the 2016 Nonlinear Mechanics and Dynamics Summer Research Institute, a six-week research program sponsored by Sandia National Laboratories and the University of New Mexico consisting of 27 graduate students participating in ten different projects. Two separate surveys were administered at the beginning and end of the Institute, in addition to interviews and observation, in order to study the effects of various cultural factors on engineering processes and maintaining professional interactions. The results of this study indicate that cultural differences are not a significant barrier to engineering progress and most cultural issues are minor. A variety of cultures instead provide new perspectives, advancing universal understanding.

42 ENGINEERING↗

Resilience Project Implementation Memo for the Town of Stowe, Vermont

This Resilience Project Implementation Process Memo was developed with and for the Town of Stowe Electric Department (SED) under the DOE C2C Expert Match Program. It is intended to outline the steps of a process, albeit nonlinear and iterative, to plan, design, and implement community resilience hubs.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Final Report for Statistical Nonlinear Optics of High Energy Density Plasmas: The Physics of Multiple Crossing Laser Beams

The various implementations of the STUD pulse program (spike trains of uneven duration and delay) for LPI (laser-plasma instability) control were studied in depth, and novel regimes were found. How to generate STUD pulses with large time-bandwidth products, how to measure their optical scattering signatures, and how to experimentally demonstrate their usefulness were explored. Theoretical and numerical studies were conducted on Stimulated Brillouin Scattering (SBS) and Crossed Beam Energy Transfer (CBET) including statistical models. We established how LPI can be tamed and gain democratized in space and time. Implementing STUD pulses on NIF was also studied. Future high rep rate lasers and fast diagnostics will aid in the adoption of the whole STUD pulse program for LPI control in High Energy Density Plasmas (HEDP).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Spatiotemporal control of laser intensity using differentiable programming

Optical techniques for spatiotemporal control can produce laser pulses with custom amplitude, phase, or polarization structure. In nonlinear optics and plasma physics, the use of structured pulses typically follows a forward design approach, in which the efficacy of a known structure is analyzed for a particular application. Inverse approaches, in contrast, enable the discovery of new structures with the potential for superior performance. Here, an implementation of the unidirectional pulse propagation equation that supports automatic differentiation is combined with gradient-based optimization to design structured pulses with features that are advantageous for a range of nonlinear optical and plasma-based applications: (1) a longitudinally uniform intensity over an extended region, (2) a superluminal intensity peak that travels many Rayleigh ranges with constant duration, spot size, and amplitude, and (3) a laser pulse that ionizes a gas to form a uniform column of plasma. In the final case, optimizing the full spatiotemporal structure improves the performance by a factor of 15 compared to optimizing only spatial or only temporal structure, highlighting the advantage of spatiotemporal control.

automatic differentiation↗

Applications of the Dulmage–Mendelsohn decomposition for debugging nonlinear optimization problems

Nonlinear modeling and optimization is a valuable tool for aiding decisions by engineering practitioners, but programming an optimization problem based on a complex electrical, mechanical, or chemical process is a time-consuming and error-prone activity. Therefore, there is a need for model analysis and debugging tools that can detect and diagnose modeling errors. One such tool is the Dulmage–Mendelsohn decomposition, which identifies structurally under- and over-determined subsets in systems of equations and variables by partitioning the bipartite graph of the system. This work provides the necessary background to understand the Dulmage–Mendelsohn decomposition and its application to the analysis of nonlinear optimization problems, demonstrates its use in diagnosing a variety of modeling errors, and introduces software implementations for analyzing nonlinear optimization problems in the Pyomo and JuMP algebraic modeling languages.

42 ENGINEERING↗

A unified funnel restoration SQP algorithm

We consider nonlinearly constrained optimization problems and discuss a generic double-loop framework consisting of basic algorithmic ingredients that unifies a broad range of nonlinear optimization solvers. This framework has been implemented in the open-source solver Uno, a Swiss Army knife-like C++ optimization framework that unifies many nonlinearly constrained nonconvex optimization solvers. We illustrate the framework with a sequential quadratic programming (SQP) algorithm that maintains an acceptable upper bound on the constraint violation, called a funnel, that is monotonically decreased to control the feasibility of the iterates. Infeasible quadratic subproblems are handled by a feasibility restoration strategy. Globalization is controlled by a line search or a trust-region method. We prove global convergence of the trust-region funnel SQP method, building on known results from filter methods. We implement the algorithm in Uno, and we provide extensive test results for the trust-region line-search funnel SQP on small CUTEst instances.

Kiessling, David [Katholieke Univ. Leuven, Heverle↗

A symbolic framework to obtain mid-fidelity models of flexible multibody systems with application to horizontal-axis wind turbines

Abstract. The article presents a symbolic framework (also called computer algebra program) that is used to obtain, in symbolic mathematical form, the linear and nonlinear equations of motion of a mid-fidelity multibody system including rigid and flexible bodies. Our approach is based on Kane's method and a nonlinear shape function representation for flexible bodies. The shape function approach does not represent the state of the art for flexible multibody dynamics but is an effective trade-off to obtain mid-fidelity models with few degrees of freedom, taking advantage of the separation of space and time. The method yields compact symbolic equations of motion with implicit account of the constraints. The general and automatic framework facilitates the creation and manipulation of models with various levels of complexity by adding or removing degrees of freedom. The symbolic treatment allows for analytical gradients and linearized equations of motion. The linear and nonlinear equations can be exported to Python code or dedicated software. There are multiple applications, such as time domain simulation, stability analyses, frequency domain analyses, advanced controller design, state observers, and digital twins. In this article, we describe the method we used to systematically generate the equations of motion of multibody systems and present the implementation of the framework using the Python package SymPy. We apply the framework to generate illustrative land-based and offshore wind turbine models. We compare our results with OpenFAST simulations and discuss the advantages and limitations of the method. The Python implementation is provided as an open-source project.

Branlard, Emmanuel (ORCID:0000000277506128)↗

Improved Measurements of Nonlinear Integrable Optics Invariants at IOTA

Nonlinear integrable optics (NIO) are a promising novel approach at improving the stability of high intensity beams. Implementations of NIO based on specialized magnetic elements are being tested at the Integrable Optics Test Accelerator (IOTA) at Fermilab. One method of verifying proper implementation of these solutions is by measuring the analytic invariants predicted by theory. The initial measurements of nonlinear invariants were performed during IOTA run in 2019/20, however the covid-19 pandemic prevented the full-scale experimental program from being completed. Several important improvements were implemented in IOTA for the 2022/23 run, including the operation at higher beam energy of 150 MeV, improved optics control, and chromaticity correction. This report presents the results of improved measurements of nonlinear invariants.

43 PARTICLE ACCELERATORS↗

Improved measurements of nonlinear integrable optics at IOTA

Nonlinear integrable optics (NIO) are a promising novel approach at improving the stability of high intensity beams. Implementations of NIO based on specialized magnetic elements are being tested at the Integrable Optics Test Accelerator (IOTA) at Fermilab. One method of verifying proper implementation of these solutions is by measuring the analytic invariants predicted by theory. The initial measurements of nonlinear invariants were performed during IOTA run in 2019/20, however the covid-19 pandemic prevented the full-scale experimental program from being completed. Several important improvements were implemented in IOTA for the 2022/23 run, including the operation at higher beam energy of 150 MeV, improved optics control, and chromaticity correction. This report presents on the improved calibrations of the NIO for nonlinear invariant measurements.

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↗

2024 Incubator Final Report: Active Hybrid Mooring

This Incubator project has defined and developed a representation of an active hybrid mooring (AHM) system in Simulink to assess whether the speed and accuracy of AHM is sufficient to warrant further investigation in a larger project. We first built a representative mooring system for the VolturnUS-S 15MW floating offshore wind turbine (FOWT) in 1,000-meter water depth to provide baseline loads and displacements needed in the AHM system, with the results showing off-the-shelf actuators would be practical. We then defined model configurations for the numerical and physical substructures, using standalone MoorDyn v2 with a novel input strategy and chained hydraulic actuators represented in Simulink Simscape, respectively. The remainder of the project has focused on verifying the existing models for speed and accuracy against OpenFAST simulations, increasing their complexity to enable six degree of freedom (6DOF) operation, and connecting the numerical and physical substructures into a single holistic model. The verification phase of the project has shown promising results, though further refinement is needed to resolve accuracy discrepancies and for robust performance across a wider range of metocean conditions. The numerical model we developed is computationally stable and operates faster than real time using a time step of 0.01 seconds or shorter, indicating the model is likely fast enough for our AHM design to operate successfully. Additionally, the Simscape model of a hydraulic actuator operates with less than 0.3% relative error when actuating mooring loads from an OpenFAST simulation of the VolturnUS-S 1,000-meter mooring system. The aforementined refinements will be performed in the next year as a Continuing Incubator, with key remaining tasks include programming the predictor-corrector loop to minimize movement error, testing against more complex nonlinear wave cases, and expanding to enable loads from connected shared mooring lines.

17 WIND ENERGY↗

Iteration-based Linearized Distribution-level Locational Marginal Price for Three-phase Unbalanced Distribution Systems

Distributed energy resources (DERs) are rocking the utilities’ business landscape. It calls for competitive market environments that incentivize DERs to form maximum operating efficiency. Among proposed pricing schemes, distribution-level locational marginal price (DLMP) is effective in signaling the marginal generation cost differences driven by energy losses and network constraints. It can be derived from a distribution-level optimal power flow (OPF) framework, as it essentially presents the sensitivity of optimized generation cost towards incremental loads. However, due to the high resistance-to-inductance ratio and unbalanced characteristics of distribution networks, computational affordable DLMPs are highly challenged. This article provides a linear-approximated DLMP that can be solved efficiently and generalized to account for reactive power flow, three-phase unbalanced loads and meshed network structure. The successive linear programming technique is introduced to enhance the model accuracy. Case studies on an IEEE 123-Bus system validate its accuracy against a nonlinear benchmark and capability in offering proper incentives.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Solution of monotone complementarity and general convex programming problems using a modified potential reduction interior point method

We present a homogeneous algorithm equipped with a modified potential function for the monotone complementarity problem. We show that this potential function is reduced by at least a constant amount if a scaled Lipschitz condition (SLC) is satisfied. A practical algorithm based on this potential function is implemented in a software package named iOptimize. The implementation in iOptimize maintains global linear and polynomial time convergence properties, while achieving practical performance. It either successfully solves the problem, or concludes that the SLC is not satisfied. When compared with the mature software package MOSEK (barrier solver version 6.0.0.106), iOptimize solves convex quadratic programming problems, convex quadratically constrained quadratic programming problems, and general convex programming problems in fewer iterations. Moreover, several problems for which MOSEK fails are solved to optimality. In addition, we also find that iOptimize detects infeasibility more reliably than the general nonlinear solvers Ipopt (version 3.9.2) and Knitro (version 8.0).

97 MATHEMATICS AND COMPUTING↗

Model Identification of Inverter Nonlinear Control Dynamics

To comprehensively predict in silico potential disruptive dynamic performance of converter-based generation and develop solutions, the application of electromagnetic transient programs is desirable. However, these accurate simulations highly depend on inverter models of actual field implementations, which are often unavailable because they are proprietary. To address this problem, we propose a time-domain model identification method that can be used to determine the control models of inverters in continuous time from multiple experimental studies. The focus is on identifying the dynamic representation of controls because they dominate the dynamic response of renewable assets. To illustrate the applicability of the technique, we consider a subclass of inverter controllers for autonomous grid applications that rely on nonlinear oscillators.

24 POWER TRANSMISSION AND DISTRIBUTION↗