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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Rolling Horizon with K-Position Search Method for Strategic Deconfliction of Package Delivery UAS

In this research, the strategic deconfliction of unmanned aircraft systems for an urban package delivery environment with two depots and multiple drop-off locations is studied. This research aims to formulate a mathematical model to compute both the departure sequence and scheduled time of departure for each unmanned aircraft system at a depot, considering temporal constraints at en-route crossing waypoints and depots for strategic deconfliction. However, the problem formulation results in an NP-hard mixed-integer nonlinear programming problem for the global optimal solution, so instead, a "rolling horizon with𝑘-position search"heuristic method is developed. The simulation studies show that an increase in the value of𝑘(the parameter used to determine the size of the local neighborhood) reduces the average ground delay at the cost of an increase in the computation time for a given problem size. The study also shows an order of magnitude increase in the maximum number of flights scheduled with the integration of rolling horizon (time decomposition) compared to those without the integration of rolling horizon in the heuristic algorithm for a given computation time cut off.

UTM↗
DOE OSTI2015

Unit Commitment With Continuous-Time Generation and Ramping Trajectory Models

There is increasing evidence of shortage of ramping resources in the real-time operation of power systems. To explain and remedy this problem systematically, we take a novel look at the way the day-ahead unit commitment (UC) problem represents the information about load, generation and ramping constraints. We specifically investigate the approximation error made in mapping of the original problem, that would decide the continuous-time generation and ramping trajectories of the committed generating units, onto the discrete-time problem that is solved in practice. We first show that current practice amounts to approximating the trajectories with linear splines. We then offer a different representation through cubic splines that provides physically feasible schedules and increases the accuracy of the continuous-time generation and ramping trajectories by capturing sub-hourly variations and ramping of load in the day-ahead power system operation. The corresponding day-ahead UC model is formulated as an instance of mixed-integer linear programming (MILP), with the same number of binary variables as the traditional UC formulation. Numerical simulation over real load data from the California ISO demonstrate that the proposed UC model reduces the total day-ahead and real-time operation cost, and the number of events of ramping scarcity in the real-time operations.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1974

Power processing methodology

Discussion of the interim results of a program to investigate the feasibility of formulating a methodology for the modeling and analysis of aerospace electrical power processing systems. The object of the total program is to develop a flexible engineering tool which will allow the power processor designer to effectively and rapidly assess and analyze the tradeoffs available by providing, in one comprehensive program, a mathematical model, an analysis of expected performance, simulation, and a comparative evaluation with alternative designs. This requires an understanding of electrical power source characteristics and the effects of load control, protection, and total system interaction.

Fegley, K. A.↗
DOE OSTI2022

Residuals-based distributionally robust optimization with covariate information

We consider data-driven approaches that integrate a machine learning prediction model within distributionally robust optimization (DRO) given limited joint observations of uncertain parameters and covariates. Our framework is flexible in the sense that it can accommodate a variety of regression setups and DRO ambiguity sets. We investigate asymptotic and finite sample properties of solutions obtained using Wasserstein, sample robust optimization, and phi-divergence-based ambiguity sets within our DRO formulations, and explore cross-validation approaches for sizing these ambiguity sets. Through numerical experiments, we validate our theoretical results, study the effectiveness of our approaches for sizing ambiguity sets, and illustrate the benefits of our DRO formulations in the limited data regime even when the prediction model is misspecified.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Strategic Offering for Wind Power Producers Considering Energy and Flexible Ramping Products

The increasing deployments of renewable generation methods, such as wind, affects the flexibility of electric power system operating due to their inherent variability and uncertainty. To mitigate this, power systems need flexible resources. This paper investigates the potential for wind power to provide flexible ramping products in the real-time market, an additional value stream to the energy it provides. The proposed model for wind power’s strategic offering is formulated as a bi-level optimization problem with wind profit maximization at the upper level and the independent system operator’s economic dispatch—considering both the energy balance and the flexible ramping requirement to counter uncertainty—at the lower level. This bi-level model is converted to a mathematical program with equilibrium constraints (MPEC) by recasting the lower level problem with its Karush–Kuhn–Tucker optimality conditions. Then, through strong duality theory and the big-M method, the MPEC model is converted to a mixed-integer linear programming model. The opportunity cost and the price for wind power-providing ramping products are analyzed. Numerical examples based on a 5-bus network are presented to verify the proposed model and concept.

17 WIND ENERGY↗
DOE OSTI2016

A Piecewise Deterministic Markov Toy Model for Traffic/Maintenance and Associated Hamilton–Jacobi Integrodifferential Systems on Networks

We study optimal control problems in infinite horizon whxen the dynamics belong to a specific class of piecewise deterministic Markov processes constrained to star-shaped networks (corresponding to a toy traffic model). We adapt the results in Soner (SIAM J Control Optim 24(6):1110–1122, 1986) to prove the regularity of the value function and the dynamic programming principle. Extending the networks and Krylov’s “shaking the coefficients” method, we prove that the value function can be seen as the solution to a linearized optimization problem set on a convenient set of probability measures. The approach relies entirely on viscosity arguments. As a by-product, the dual formulation guarantees that the value function is the pointwise supremum over regular subsolutions of the associated Hamilton–Jacobi integrodifferential system. This ensures that the value function satisfies Perron’s preconization for the (unique) candidate to viscosity solution.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1988

Benchmark solutions for the galactic ion transport equations with spatial and energy coupling

In order to anticipate future space shielding requirements, NASA has initiated an effort to formulate computational methods to simulate radiation effects in space. As part of the program, numerical transport algorithms have been developed for the deterministic Boltzman equation describing galactic cosmic ray (GCR) interactions with matter. It thus becomes necessary to assess the accuracy of proposed deterministic algorithms. For this reason, analytical benchmark solutions to mathematically tractable galactic cosmic ray equations have recently been obtained. Even though these problems involve simplifying assumptions of the associated physics, they still contain the essential features of the basic transport processes. The solutions obtained are features of the basic transport processes. The solutions obtained are compared to results from numerical algorithms in order to ensure proper coding and to provide a measure of the accuracy of the numerical methods used in the algorithm. For the first time, mathematical methods have been applied to the galactic ion transport (GIT) equations in the straight ahead approximation with constant nuclear properties. The approach utilizes a Laplace transforms inversion yielding a closed form benchmark solution which is also computationally efficient.

Ganapol, Barry D.↗
DOE OSTI2019

Robust and High-Waste Loaded Glass Formulations for Hanford Low-Activity Waste - 19050

The Hanford Low-Activity Waste (LAW) vitrification facility will receive LAW batches with abruptly changing compositions that span relatively broad ranges. In particular, the concentration ratios of sodium-to-sulfate, -chloride, -chromate, -phosphate, -potassium, and -aluminate create challenges in glass formulation. Glass compositions that simultaneously satisfy several processability, product quality, and composition constraints with sufficient confidence must be formulated for each LAW batch in near-realtime to enable desired throughput. A method for developing glass compositions while maximizing the waste loading and the process robustness was developed. This method compiled the glass composition-property data collected over nearly two decades, primarily by Pacific Northwest National Laboratory, The Catholic University of America, and Savannah River National Laboratory, in support of both the Hanford Tank Waste Treatment and Immobilization Plant project and the U.S. Department of Energy Office of River Protection Enhanced Waste Glass project. These data were used to develop models that link the melter feed composition to the resulting glass properties. Mathematical constructs were developed to optimize glass compositions for each waste batch in real-time that will simultaneously satisfy processability, product quality, and composition constraints with sufficient statistical confidence. Property prediction and various process uncertainties were reflected in the formulation. The formulation approach directly lends itself to adjustment of glass formulation to achieve program goals. This presentation will summarize the formulation approach and supporting data. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Update on NNL TSL Validation Efforts [Slides]

This presentation covers how a new sub-thermal neutron transmission capability at RPI LINAC is now in use to help validate TSL evaluations. Additionally, NNL has contributed multiple material evaluations to ENDF/B-VIII.1. Lastly, this presentation covers new methods to improve TSL physics and accuracy are being formulated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2021

A Study into Validating a Coupled Method of Characteristics and Direct Simulation Monte Carlo Method Against Empirical Data

The following will outline the methodology and results of validating a coupled Method of Characteristics (MOC) and Direct Simulation Monte Carlo (DSMC) method. This research focused specifically on modeling plume impingement, induced by Reaction Control System (RCS) thrusters that flew on the National Aeronautics and Space Administration’s (NASA’s) space shuttle Discovery. For each simulation, the continuum portion of the RCS thruster was simulated using MOC for solving hyperbolic Partial Differential Equations (PDEs) and computed with the NASA code, Reacting and Multi-phase Program (RAMP). The solution was then implemented as a starting condition into the NASA DSMC code, Direct Simulation and Monte Carlo Analysis Code (DAC). Typically, DSMC models rely on code-to-code validation for fidelity. The significance of this research is in its ability to validate its models against empirical data. Prior to computing solutions for these simulations, the mesh size and structure were optimized and many variants of DSMC input parameters were iterated on in order to acquire a reliable, mesh-independent, fully optimized numerical solution. This research will discuss the mathematical formulation of MOC for nozzle flow and DSMC for rarefied gases. Additionally, it will provide an explanation of how to implement these mathematical concepts into the two solvers: RAMP and DAC. Ultimately, this research will demonstrate that the overall process illustrated produces results in good agreement with empirical data. As a consequence, the methodology presented is granted an increased level of confidence and will greatly contribute to the aerospace industry and its effort in understanding and predicting rarefied flow fields.

Direct Simulation Monte Carlo Analysis Code↗
DOE OSTI2015

A Characterization of the Reflected Quasipotential

Recent interest in the reflected quasipotential comes from the queueing theory literature, specifically the analysis of so-called (b,A,D) reflected Brownian motion where it is the large deviation rate function for the stationary distribution. Our purpose here is to characterize the reflected quasipotential in terms of a first-order Hamilton–Jacobi equation. Using conventional dynamic programming ideas, along with a complementarity problem formulation of the effect of the Skorokhod map on absolutely continuous paths, we will derive necessary conditions in the form of viscosity-sense boundary conditions. It turns out that even with these boundary conditions solutions are not unique. Thus a unique characterization needs to refer to some additional property of V(⋅). We establish such a characterization in two dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Modeling and optimization of steady flow of natural gas and hydrogen mixtures in pipeline networks

Here, we extend the canonical problems of simulation and optimization of steady-state gas flows in pipeline networks with compressors to the transport of mixtures of highly heterogeneous gases injected throughout a network. Our study is motivated by proposed projects to blend hydrogen generated using clean energy into existing natural gas pipeline systems as part of efforts to reduce the reliance of energy systems on fossil fuels. Flow in a pipe is related to endpoint pressures by a basic Weymouth equation model, with an ideal gas equation of state, where the wave speed depends on the hydrogen concentration. At vertices, in addition to mass balance, we also consider mixing of incoming flows of varying hydrogen concentrations. The problems of interest are the heterogeneous gas flow simulation (HGFS), which determines system pressures and flows given fixed boundary conditions and compressor settings, as well as the heterogeneous gas flow optimization (HGFO), which extremizes an objective by determining optimal boundary conditions and compressor settings. We examine conditions for uniqueness of solutions to the HGFS, as well as compare and contrast mixed-integer and continuous nonlinear programming formulations for the HGFO. We develop computational methods to solve both problems, and examine their performance using four test networks of increasing complexity.

08 HYDROGEN↗
DOE OSTI2022

Quantization: History and problems

In this work, I explore the concept of quantization as a mapping from classical phase space functions to quantum operators. I discuss the early history of this notion of quantization with emphasis on the works of Schrödinger and Dirac, and how quantization fit into their overall understanding of quantum theory in the 1920's. Dirac, in particular, proposed a quantization map which should satisfy certain properties, including the property that quantum commutators should be related to classical Poisson brackets in a particular way. However, in 1946, Groenewold proved that Dirac's mapping was inconsistent, making the problem of defining a rigorous quantization map more elusive than originally expected. This result, known as the Groenewold-Van Hove theorem, is not often discussed in physics texts, but here I will give an account of the theorem and what it means for potential ``corrections” to Dirac's scheme. Other proposals for quantization have arisen over the years, the first major one being that of Weyl in 1927, which was later developed by many, including Groenewold, and which has since become known as Weyl Quantization in the mathematical literature. Another, known as Geometric Quantization, formulates quantization in differential-geometric terms by appealing to the character of classical phase spaces as symplectic manifolds; this approach began with the work of Souriau, Kostant, and Kirillov in the 1960's. I will describe these proposals for quantization and comment on their relation to Dirac's original program. Along the way, the problem of operator ordering and of quantizing in curvilinear coordinates will be described, since these are natural questions that immediately present themselves when thinking about quantization.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Pyomo v5.0

Pyomo supports the formulation and analysis of mathematical models for complex optimization applications. This capability is commonly associated with algebraic modeling languages (AMLs), which support the description and analysis of mathematical models with a high-level language. Although most AMLs are implemented in custom modeling languages, Pyomo's modeling objects are embedded within Python, a full- featured high-level programming language that contains a rich set of supporting libraries.

Woodruff, David↗
DOE OSTI2026

Classical-Quantum Algorithm for Solving Stochastic Programs

Stochastic programming provides a rigorous mathematical framework for making decisions under uncertainty in a risk-aware manner. Two-stage stochastic programming is, perhaps, the simplest form of this framework. Here the first-stage variables represent decisions that must be made "here and now" in the face of uncertainty, while the second-stage variables are decisions made after uncertain events. However, the broad adoption of stochastic programming has been hindered by computational challenges caused by the two-stage stochastic programming formulation which requires solving an ensemble of optimization problems. Using quantum amplitude estimation (QAE), quantum computers have shown the theoretic ability to compute expectations with Monte-Carlo methods with quadratically fewer samples than classical methods. In this work, we present a quantum algorithm for computing the expectation term using QAE for given first-stage decisions. Further, we detail methods of computing gradient information from the quantum calculation enabling the application of classical gradient-based optimization techniques. The result is a classical-quantum hybrid method of solving two-stage stochastic programs. These techniques are demonstrated with computational experiments based an engineering optimization problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Pyomo

Pyomo supports the formulation and analysis of mathematical models for complex optimization applications. This capability is commonly associated with algebraic modeling languages (AMLs), which support the description and analysis of mathematical models with a high-level language. Although most AMLs are implemented in custom modeling languages, Pyomo's modeling objects are embedded within Python, a full-featured high-level programming language that contains a rich set of supporting libraries.

Source record↗
DOE OSTI2021

A Tutorial for Using an Open-Source Solver for the Regional Energy Deployment System (ReEDS) Model

ReEDS is a publicly available model developed at NREL that can be used to analyze the potential evolution of the U.S. electric power system into the future. ReEDS is formulated as a linear program, written in GAMS, and solved using a linear programming solver (e.g., CPLEX). This tutorial provides context and understanding for the implications of using an open-source solver for ReEDS.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Data-driven optimization of mixed-integer bi-level multi-follower integrated planning and scheduling problems under demand uncertainty

The coordination of interconnected elements across the different layers of the supply chain is essential for all industrial processes and the key to optimal decision-making. Yet, the modeling and optimization of such interdependent systems are still burdensome. Here we address the simultaneous modeling and optimization of medium-term planning and short-term scheduling problems under demand uncertainty using mixed-integer bi-level multi-follower programming and data-driven optimization. Bi-level multi-follower programs model the natural hierarchy between different layers of supply chain management holistically, while scenario analysis and data-driven optimization allow us to retrieve the guaranteed feasible solutions of the integrated formulation under various demand considerations. We address the data-driven optimization of this challenging class of problems using the DOMINO framework, which was initially developed to solve single-leader single-follower bi-level optimization problems to guaranteed feasibility. This framework is extended to solve single-leader multi-follower stochastic formulations and its performance is characterized by well-known single and multi-product process scheduling case studies. Through our data-driven algorithmic approach, we present guaranteed feasible solutions to linear and nonlinear mixed-integer bi-level formulations of simultaneous planning and scheduling problems and further characterize the effects of the scheduling level complexity on the solution performance, which spans over several hundred continuous and binary variables, and thousands of constraints.

42 ENGINEERING↗