Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mathematical programming formulations”

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 55 records · Page 3

Three dimensional thermal pollution models. Volume 1: Review of mathematical formulations

A mathematical model package for thermal pollution analyses and prediction is presented. These models, intended as user's manuals, are three dimensional and time dependent using the primitive equation approach. Although they have sufficient generality for application at sites with diverse topographical features; they also present specific instructions regarding data preparation for program execution and sample problems. The mathematical formulation of these models is presented including assumptions, approximations, governing equations, boundary and initial conditions, numerical method of solution, and same results.

Lee, S. S.↗

The flow of plasma in the solar terrestrial environment

The overall goal of our NASA Theory Program was to study the coupling, time delays, and feedback mechanisms between the various regions of the solar-terrestrial system in a self-consistent, quantitative manner. To accomplish this goal, it will eventually be necessary to have time-dependent macroscopic models of the different regions of the solar-terrestrial system and we are continually working toward this goal. However, with the funding from this NASA program, we concentrated on the near-earth plasma environment, including the ionosphere, the plasmasphere, and the polar wind. In this area, we developed unique global models that allowed us to study the coupling between the different regions. These results are highlighted in the next section. Another important aspect of our NASA Theory Program concerned the effect that localized 'structure' had on the macroscopic flow in the ionosphere, plasmasphere, thermosphere, and polar wind. The localized structure can be created by structured magnetospheric inputs (i.e., structured plasma convection, particle precipitation or Birkland current patterns) or time variations in these input due to storms and substorms. Also, some of the plasma flows that we predicted with our macroscopic models could be unstable, and another one of our goals was to examine the stability of our predicted flows. Because time-dependent, three-dimensional numerical models of the solar-terrestrial environment generally require extensive computer resources, they are usually based on relatively simple mathematical formulations (i.e., simple MHD or hydrodynamic formulations). Therefore, another goal of our NASA Theory Program was to study the conditions under which various mathematical formulations can be applied to specific solar-terrestrial regions. This could involve a detailed comparison of kinetic, semi-kinetic, and hydrodynamic predictions for a given polar wind scenario or it could involve the comparison of a small-scale particle-in-cell (PIC) simulation of a plasma expansion event with a similar macroscopic expansion event. The different mathematical formulations have different strengths and weaknesses and a careful comparison of model predictions for similar geophysical situations provides insight into when the various models can be used with confidence.

Schunk, Robert W.↗

Open-source Tools for Solving Grid Optimization Problems: ARPA-e Benchmark Algorithm Overview [Slides]

This document contains the official formulation that will be used for evaluation in Challenge 2 of the Grid Optimization (GO) Competition. Minor changes may occur within the formulation. Entrants will be notified when a new version is released. Changes are not expected to be of a significance that would cause a change in approach for the Entrants. This formulation builds upon the Challenge 1 formulation published in ARPA-E DE-FOA-0001952. Entrants will be judged based on the current official Challenge 2 formulation posted on the GO Competition website (this document, which is subject to change), not the formulation posted in DE-FOA-0001952. Entrants are permitted and encouraged to use any alternative problem formulation and modeling convention within their own software (such as convex relaxation, decoupled power flow formulations, current-voltage formulations, etc.) in an attempt to produce an exact or approximate solution to this particular mathematical program. However, the judging of all submitted approaches must conform to the official formulation presented here.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Optimal Strategies for Hybrid Battery‐Storage Systems Design

As stationary hybrid energy‐storage systems (HESS) for power systems applications have recently drawn interest due to their enhanced performance and decreasing cost, developing systematic approaches for HESS design while considering controls is gaining traction. Herein, a method is presented to optimally design hybrid battery storage by proposing a mathematical modeling framework, formulated as a mixed integer linear programming model. The optimization is capable of handling multiple subsystems of batteries, considering their economic and technological performance. Decisions involve sizing of the batteries, optimal temporal and strategic dispatch to end uses, and energy sources for charging each battery. The applicability of the model is tested on four case studies for three battery chemistries representing distinct objectives: high‐power, high‐energy, and second life. Compared to traditionally designed battery storage with a homogeneous battery, optimally designed hybrid systems can save 12%–26% of system costs, depending on the nature of the dispatch profile. Findings point to design preference toward the second life battery supplemented with some high‐power or high‐energy battery capacity, or both. With the utilized electricity price structure, customers can experience approximately 10%–35% reduction in their bills.

Koleva, Mariya↗

Integrating Concentrating Solar Power Technologies into the Hybrid Optimization and Performance Platform (HOPP)

As the world increases renewable energy deployment, there is an increasing interest in hybridizing various generation and storage technologies to maximize net benefit to the developer and/or off-taker. A particularly interesting combination of renewable technologies is concentrating solar power (CSP) with thermal energy storage (TES), photovoltaics (PV), and electrochemical battery energy storage (BESS). Due to the system complexity of CSP technology, it is difficult to evaluate the technological and financial performance of a CSP-PV hybrid system without detailed modeling of annual operations. To address this challenge, we have developed a modeling framework for evaluating the performance and financial viability of CSP systems hybridized with PV and battery technologies. This modeling effort incorporates CSP tower and trough systems into an existing modeling tool recently developed by NREL referred to as the Hybrid Optimization and Performance Platform (HOPP). This report outlines the modeling methodology as well as preliminary results from example case studies conducted using the model. The methodology describes: (i) the integration of CSP tower and troughs into HOPP using python interfaces to access System Advisor Model (SAM) underlining technology models, (ii) the mathematical formulation of the mixed integer linear program dispatch optimization model which optimizes operations of storage asset to either maximize system revenue or minimize operating cost while load following, (iii) the design analysis methods implemented within HOPP, and (iv) simulation clustering for the purposes of reducing computational expense. We exercise the model using a case study of a future scenario where we assume (i) CSP and PV technologies achieve the 2030 cost targets provided by the Solar Energy Technologies Office (SETO), (ii) battery costs reduce to the 2030 mid cost projection presented by NREL. Lastly, (iii) electricity prices for southern California in 2030 are provided by NREL's Cambium database, and (iv) a capacity payment of $150/kW-yr based on the system capacity factor during the to 100 net-load hours.

14 SOLAR ENERGY↗

Adaptive age replacement.

Loss structure models concerning age replacement policies, discussing mathematical aspects of maintenance with dynamic programming formulation

COST ESTIMATE↗

Subsonic panel method for the efficient analysis of multiple geometry perturbations

An accurate and efficient method was developed for the aerodynamic analysis of a series of arbitrary small geometry perturbations to a given baseline configuration. The method is appropriate for wing-fuselage configurations in incompressible potential flow. Mathematical formulations are presented for three computer programs that are employed. The first program is a conventional surface panel method for calculating the baseline singularity distribution. The second program calculates a partial derivative matrix. Each element of the matrix is the rate of change of singularity strength at one point with respect to a surface coordinate of a different point. For each baseline configuration, the calculated quantities from the first two programs establish an input file for the third. The third program calculates the surface pressure distribution and forces and moments for a series of geometry perturbations.

Bristow, D. R.↗

Capillary priming characteristics of a dual passage heat pipe in zero-g

Technical improvements of a long life heat rejection system, suitable for long duration high power missions, that can be constructed and deployed in orbit is discussed. A mathematical model is formulated and a computer program developed which describes the transient priming characteristics of a dual passage heat pipe. An experimental test package is described for flight in the KC-135 Zero-g Aircraft, to be used to verify the modeling predictions.

Peterson, G. P.↗

The flow of plasma in the solar terrestrial environment

The overall goal of our NASA Theory Program is to study the coupling, time delays, and feedback mechanisms between the various regions of the solar-terrestrial system in a self-consistent, quantitative, manner. To accomplish this goal, it will eventually be necessary to have time-dependent macroscopic models of the different regions of the solar-terrestrial system and we are continually working toward this goal. However, our immediate emphasis is on the near-earth plasma environment, including the ionosphere, the plasmasphere, and the polar wind. In this area, we have developed unique global models that allow us to study the coupling between the different regions. These results are highlighted. Another important aspect of our NASA Theory Program concerns the effect that localized structure has on the macroscopic flow in the ionosphere, plasmasphere, thermosphere and polar wind. The localized structure can be created by structured magnetospheric inputs (i.e., structured plasma convection, particle precipitation or Birkeland current patterns) or time variations in these inputs due to storms and substorms. Also, some of the plasma flows that we predict with our macroscopic models may be unstable. Another one of our goals is to examine the stability of our predicted flows. Because time-dependent three-dimensional numerical models of the solar-terrestrial environment generally require extensive computer resources, they are usually based on relatively simple mathematical formulations (i.e., simple MHD or hydrodynamic formulations). Therefore, another long-range goal of our NASA Theory Program is to study the conditions under which various mathematical formulations can be applied to specific solar-terrestrial regions. This may involve a detailed comparison of kinetic, semikinetic, and hydrodynamic predictions for a given polar wind scenario or it may involve the comparison of a small-scale particle-in-cell (PIC) simulation of a plasma expansion event with a similar macroscopic expansion event. The different mathematical formulations have different strengths and weaknesses and a careful comparison of model predictions for similar geophysical situations will provide insight into when the various models can be used with confidence.

Schunk, Robert W.↗

The flow of plasma in the solar terrestrial environment

The overall goal of our NASA Theory Program is to study the coupling, time delays, and feedback mechanisms between the various regions of the solar-terrestrial system in a self-consistent, quantitative manner. To accomplish this goal, it will eventually be necessary to have time-dependent macroscopic models of the different regions of the solar-terrestrial system and we are continually working toward this goal. However, our immediate emphasis is on the near-earth plasma environment, including the ionosphere, the plasmasphere, and the polar wind. In this area, we have developed unique global models that allow us to study the coupling between the different regions. Another important aspect of our NASA Theory Program concerns the effect that localized structure has on the macroscopic flow in the ionosphere, plasmasphere, thermosphere, and polar wind. The localized structure can be created by structured magnetospheric inputs (i.e., structured plasma convection, particle precipitation or Birkeland current patterns) or time variations in these inputs due to storms and substorms. Also, some of the plasma flows that we predict with our macroscopic models may be unstable, and another one of our goals is to examine the stability of our predicted flows. Because time-dependent, three-dimensional numerical models of the solar-terrestrial environment generally require extensive computer resources, they are usually based on relatively simple mathematical formulations (i.e., simple MHD or hydrodynamic formulation). Therefore, another long-range goal of our NASA Theory Program is to study the conditions under which various mathematical formulations can be applied to specific solar-terrestrial regions. This may involve a detailed comparison of kinetic, semikinetic, and hydrodynamic predictions for a given polar wind scenario or it may involve the comparison of a small-scale particle-in-cell (PIC) simulation of a plasma expansion event with a similar macroscopic expansion event. The different mathematical formulations have different strengths and weaknesses and a careful comparison of model predictions for similar geophysical situations will provide insight into when the various models can be used with confidence.

Schunk, R. W.↗

Simple intuitive models of programming

It is hypothesized that mathematical models of the programming process can be formulated to gauge the sensitivities of that process to various given parameters, and that such models can be calibrated on an empirical basis and used as guides toward maximizing productivity, documentation quality, and programming reliability. Three oversimplified models are presented as illustrations.

Tausworthe, R. C.↗

Mixed Integer Nonlinear Programming Approaches to Enhance Resiliency and Response Strategies in Critical Infrastructure

Motivation. Critical infrastructures are large, complex engineered systems that must be operated robustly under abnormal conditions resulting from natural hazards or intentional acts. For example, electric power systems must be robust to line faults, water utilities must rapidly mitigate contamination incidents, and computing networks must adapt to adversarial intrusions to protect critical information. Problem. It is difficult for decision-makers within resiliency analysis in critical infrastructure to optimize designs and develop effective response strategies that can account for uncertainties. Facing incomplete information and the sheer scope that a natural hazard or attack vector may incorporate, response can be ineffective without reliable, scalable decision support tools. These problems are intrinsically nonlinear and involve discrete decisions, and unfortunately, existing off-the-shelf mathematical programming methods cannot support optimization-based decision-making of these nonlinear at scale. Method/Approach and Results. This project emphasized development of fundamental optimization strategies that supported real- time mitigation and response for critical infrastructures. In particular, the project developed multi- tree approaches based on piecewise outer-approximations for solution of mixed-integer nonlinear programming (MINLP) problems. These techniques alternate between an MILP or MISOCP relaxation to obtain a lower bound and candidate discrete solutions and an NLP subproblem to obtain upper bounds. Using tailored relaxations based on problem structure, these methods were used to solve several key applications in resilience and response of critical infrastructure. This work resulted in two open-source, copyrighted software packages: CORAMIN (https://github.com/Coramin/Coramin) -- an object-oriented mathematical programming framework that supports tailored multi-tree algorithms for solution of large- scale mixed-integer nonlinear programming; and EGRET (Electrical Grid Research and Engineering Toolkit) (https://github.com/grid- parity-exchange/Egret) -- a declarative mathematical programming framework built upon CORAMIN and Pyomo for formulation and solution of resilience and operations problems in power grid systems. Furthermore, these tools resulted in several important published results, including the following: The first known global optimization approach that could solve the unit-commitment problem with nonlinear power flow constraints on medium-sized test problems; Improved parallel optimization-based bounds tightening and strengthening of relaxations of AC power flow constraints; and, Optimization-based approaches for improved grid resilience and use of demand response to improve grid resilience with reduction in capital requirements. Result Implications. This project developed first-of-a-kind algorithms for decision-making in critical infrastructure resilience operations and planning, as well as a next-generation toolkit for MINLP researchers. These approaches leveraged high-performance computing architectures to solve some of the largest, most challenging nonlinear discrete optimization problems to global optimality, and these successes were captured in open-source software to enable optimization-based decision-making, and efficient solution of MINLP formulations for electric power transmission grids.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Primal and dual formulations of sequential gradient-restoration algorithms for trajectory optimization problems

One of the most effective first-order algorithms for solving trajectory optimization problems is the sequential gradient-restoration algorithm (SGRA). Originally developed in the primal formulation, this algorithm is extended to incorporate a dual formulation. Both the primal formulation and the dual formulation involve a sequence of two-phase cycles, each cycle including a gradient phase and a restoration phase. In turn, each iteration of the gradient phase and the restoration phase requires the solution of an auxiliary minimization problem (AMP). In the primal formulation, the AMP is solved with respect to the variations of the state, the control, and the parameter. In the dual formulation, the AMP is solved with respect to the Lagrange multipliers. A characteristic of the dual formulation is that the AMPs associated with the gradient phase and the restoration phase of SGRA can be reduced to mathematical programming problems involving a finite number of parameters as unknowns. A comparison of the primal formulation and the dual formulation is presented. The comparison is done in terms of several trajectory optimization problems having current aerospace interest.

Miele, A.↗

Continued research on selected parameters to minimize community annoyance from airport noise

A mathematical model of the annoyance created at an airport by aircraft operations is developed. The model incorporates population distribution considerations around an airport and the annoyance caused by aircraft noise. The objective function of this model corresponds to seeking to minimize total population annoyance created by all aircraft operations in a 24-hour period. Several factors are included in this model as constraint bounded. Demand for flight services is incorporated by including lower bounds on the number of operations by type of aircraft, runway used and time period. Also upper bounds on the number of operations for each runway are included. The mathematical model as formulated is recognized as corresponding to a nonlinear integer mathematical programming problem.

Frair, L.↗

Analysis of dynamic systems with DAP4H computer program

Dynamic Analysis Program, FORTRAN 4 Level H /DAP4H/, developed from 27 subprograms, features liberal use of the subroutines, subprograms, and skeletonization to minimize programming effort in formulating models of new systems and components. It formulates mathematical models of complex mechanical, pneumatic, and hydraulic dynamic systems.

Absalom, J. G.↗

A design procedure for a tension-wire stiffened truss-column

A deployable, tension wire stiffened, truss column configuration was considered for space structure applications. An analytical procedure, developed for design of the truss column and exercised in numerical studies, was based on equivalent beam stiffness coefficients in the classical analysis for an initially imperfect beam column. Failure constraints were formulated to be used in a combined weight/strength and nonlinear mathematical programming automated design procedure to determine the minimum mass column for a particular combination of design load and length. Numerical studies gave the mass characteristics of the truss column for broad ranges of load and length. Comparisons of the truss column with a baseline tubular column used a special structural efficiency parameter for this class of columns.

Greene, W. H.↗

A centralized-based method to determine the local voltage control strategies of distributed generator operation in active distribution networks

The increasing penetration of distributed generators (DGs) exacerbates the risk of voltage violations in active distribution networks (ADNs). The var capacity provided by DG inverters is a potential solution for voltage regulation. The conventional centralized control strategies limited by computation and communication burdens are difficult to meet the requirement of rapid voltage control. However, the local control strategies based on real-time measurements have a fast response to the frequent fluctuations of DG outputs. The performance of these local controllers depends on the tuning of the control parameters. This paper proposes a centralized-based method to determine the local voltage control strategies for DGs. A centralized parameter tuning model of control curves is built, in which and curves of DG inverters are mathematically formulated based on piecewise linearization. The original mixed-integer nonlinear programming (MINLP) model is converted into an effectively solved mixed-integer second-order cone programming (MISOCP) model using convex relaxation. The potential benefits of DG inverters are explored to regulate both reactive and curtailed active power, based on local voltage measurements. Case studies on a modified PG&E 69-node distribution system are carried out to verify the effectiveness of the proposed method. Results show that power losses of ADNs are significantly reduced and voltage profiles are also improved.

24 POWER TRANSMISSION AND DISTRIBUTION↗