Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear optimization”

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 1,675 records · Page 93
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Detecting deviations from Gaussianity in high-redshift CMB lensing maps

ABSTRACT While the probability density function of the cosmic microwave background (CMB) convergence field approximately follows a Gaussian distribution, primordial non-Gaussianities, and small contributions from structures at low redshifts make the overall distribution slightly non-Gaussian. Some of the late-time component can be modelled using the distribution of galaxies and subtracted off from the original CMB lensing map to produce a map of matter distribution at high redshifts. Using this high-redshift mass map, we are able to directly study the early phases of structure formation. In this work, we forecast the detectability of signatures of non-Gaussianity due to non-linear structure formation at z > 1.2. Assuming the optimal case of no systematics, we find that it is challenging to detect such signatures in current surveys, but future experiments such as the deep field of CMB-S4 will be able to make detections of ∼7σ.

Astronomy & Astrophysics↗
DOE OSTI2021

Sequential Perturbation-based FDIA Detection using DERs for Unbalanced Distribution System

The power distribution system with its reliance on automated monitoring and control infrastructure makes it a complicated cyber-physical system and vulnerable to various cyber-attacks. This paper proposes a moving target defense inspired false data injection attack (FDIA) detection mechanism for an unbalanced distribution system. In the proposed grid diagnosis service framework, the distribution system operator judiciously manipulates power outputs of a subset of inverter-based distribution energy resources (DERs) to create secret low magnitude voltage perturbations which are inconsequential to the normal operation of the grid. A mixed-integer-linear-programming algorithm is developed to select the optimal set of DERs that can create a voltage perturbation signal of the required magnitude at each sensor location at a minimum cost. Then, a sequential detector is applied that detects for the FDIA as measurements are received from individual sensors. The performance of the proposed perturbation-based FDIA detection framework is demonstrated via simulation of the IEEE 123 bus test system.

Attack Detection↗
DOE OSTI2019

Community Microgrid Planning Considering Building Thermal Dynamics

In this paper, the building thermal dynamic characteristics are introduced in the community microgrid (MG) planning model. The proposed planning model is formulated as a mixed integer linear programming (MILP) which seeks to determine the optimal deployment strategy for various distributed energy resources (DER). The objective is to minimize the annualized cost of community MG, including investment cost in DER, operation cost for dispatchable fuel-generators (DFG), energy storage system (ESS) degradation cost, energy purchasing and peak demand charge at PCC, customer discomfort cost due to the room temperature deviation and load curtailment cost. Given the slow thermal dynamic characteristics of buildings, the heating, ventilation and air-conditioning (HVAC) system in the proposed model is treated as a demand side management (DSM) component, whose dispatch commands are provided by the central MG controller. Numerical results based on a community MG comprising of 20 residential buildings demonstrate the effectiveness of the proposed planning model and the benefits of introducing the building thermal dynamic model.

Zhang, Xiaohu↗
DOE OSTI2025

Improving Resiliency in Planning MW-Scale Medium and Heavy Duty EV Charging Stations Considering TSCOTS Optimization

Electrification of heavy-duty (HD) vehicles marks an important milestone and technical challenge in the electric vehicle (EV) industry and the public grid. However, implementing EV charging at this scale will necessitate that traditional truck stops be updated with EV charging infrastructure that could represent 10's of MW in electricity consumption. Furthermore, as the transportation sector is represented as critical infrastructure, supporting resiliency considerations in EV charging infrastructure will be critical. This paper proposes an optimization-based approach for optimally designing a MW-scale microgrid charging network. This approach transforms conventional designed truck stops into a reliable HDEV charging stations capable of overnight slow charging and 30-minute to 1 hour fast charging. Using a mixed-integer linear program formulation blending capacity planning and reliability constraints, an optimal network configuration can be solved for a proposed EV charging station that includes photovoltaic and battery energy storage capabilities.

Ponce, Moises [University of Tennessee, Knoxville ↗
DOE OSTI2018

Restoration using Distributed Energy Resources for Resilient Power Distribution Networks

In this paper, we propose a restorative method to minimize the unserved demand in a power distribution network where some islands are formed due to multiple line outages. The restoration is achieved by tie-line switching and by solving a mixed-integer non-linear problem. Due to combinatorial issues of the optimization problem, a method is proposed which restores the system by setting up microgrids boundaries within the islands and connecting them through tie-lines for exchange of active power in between them. The optimal tie-line to interconnect the microgrids is decided through information exchange between the microgrid operators. Case studies on IEEE 33-bus distribution network is presented to illustrate the results and validate the performance of the proposed method.

Sarkar, Ranadhir↗
DOE OSTI2013

Architectures and Control of Submodule Integrated DC-DC Converters for Photovoltaic Applications

This paper describes photovoltaic (PV) module architectures with parallel-connected submodule-integrated dc-dc converters (subMICs) that improve efficiency of energy capture in the presence of partial shading or other mismatch conditions. The subMICs are bidirectional isolated dc-dc converters capable of injecting or subtracting currents to balance the module substring voltages. When no mismatches are present, the subMICs are simply shut down, resulting in zero insertion losses. It is shown that the objective of minimum subMIC power processing can be solved as a linear programming problem. A simple close-to-optimal distributed control approach is presented that allows autonomous subMIC control without the need for a central controller or any communication among the subMICs. Furthermore, the proposed control approach is well suited for an isolated-port architecture, which yields additional practical advantages including reduced subMIC power and voltage ratings. The architectures and the control approach are validated by simulations and experimental results using three bidirectional flyback subMICs attached to a standard 180-W, 72-cell PV module, yielding greater than 98% module-level power processing efficiency for a mismatch less than 25%.

Olalla, C↗
DOE OSTI2017

A Geometric Approach to Aggregate Flexibility Modeling of Thermostatically Controlled Loads

Coordinated aggregation of a large population of thermostatically controlled loads (TCLs) presents a great potential to provide various ancillary services to the grid. One of the key challenges in coordination and control of TCLs is developing a simple and portable model to accurately extract their aggregate flexibility. In this paper, we propose a novel geometric approach to model the aggregate flexibility of TCLs. We show that the set of admissible power profiles of an individual TCL is a polytope, and their aggregate flexibility is the Minkowski sum of the individual polytopes. In order to represent their aggregate flexibility in an intuitive way and achieve a tractable approximation, we develop optimization-based algorithms to approximate the polytopes by the homothets of a given convex set. As a special application, this set is chosen as a virtual battery model and the corresponding optimal approximation problems are solved efficiently by equivalent linear programs. Numerical results show that our algorithms yield significant improvement in characterizing the aggregate flexibility over existing modeling methods. Finally, we also conduct case studies to demonstrate the efficacy of our approaches by coordinating TCLs to track a frequency regulation signal from the Pennsylvania-New Jersey-Maryland (PJM) Interconnection.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Sparse measurement medical CT reconstruction using multi-fused block matching denoising priors

A major challenge for medical X-ray CT imaging is reducing the number of X-ray projections to lower radiation dosage and reduce scan times without compromising image quality. However these under-determined inverse imaging problems rely on the formulation of an expressive prior model to constrain the solution space while remaining computationally tractable. Traditional analytical reconstruction methods like Filtered Back Projection (FBP) often fail with sparse measurements, producing artifacts due to their reliance on the Shannon-Nyquist Sampling Theorem. Consensus Equilibrium, which is a generalization of Plug and Play, is a recent advancement in Model-Based Iterative Reconstruction (MBIR), has facilitated the use of multiple denoisers are prior models in an optimization free framework to capture complex, non-linear prior information. However, 3D prior modelling in a Plug and Play approach for volumetric image reconstruction requires long processing time due to high computing requirement. Instead of directly using a 3D prior, this work proposes a BM3D Multi Slice Fusion (BM3D-MSF) prior that uses multiple 2D image denoisers fused to act as a fully 3D prior model in Plug and Play reconstruction approach. Our approach does not require training and are thus able to circumvent ethical issues related with patient training data and are readily deployable in varying noise and measurement sparsity levels. In addition, reconstruction with the BM3D-MSF prior achieves similar reconstruction image quality as fully 3D image priors, but with significantly reduced computational complexity. We test our method on clinical CT data and demonstrate that our approach improves reconstructed image quality.

Hossain, Maliha [ORNL]↗
DOE OSTI2023

Tools for analysis of optimization models

The software is code for analyzing, debugging, and solving optimization models. The code implements several algorithms and provides convenient APIs to apply these algorithms to optimization models. The algorithms are primarily based on graph theory. They compute well-known partitions of graphs, and use these partitions to provide the user information about their optimization model, including diagnosing certain types of structural modeling errors. Some algorithms may be used as a subroutine to solve an optimization problem, and may call an optimization solver to facilitate this. Well-known linear algebra routines are called to provide further diagnostics. The software contains no data and no models other than toy models used for testing.

Parker, Robert↗
DOE OSTI2016

Hydrodynamics of suspensions of passive and active rigid particles: a rigid multiblob approach

We develop a rigid multiblob method for numerically solving the mobility problem for suspensions of passive and active rigid particles of complex shape in Stokes flow in unconfined, partially confined, and fully confined geometries. As in a number of existing methods, we discretize rigid bodies using a collection of minimally resolved spherical blobs constrained to move as a rigid body, to arrive at a potentially large linear system of equations for the unknown Lagrange multipliers and rigid-body motions. Here we develop a block-diagonal preconditioner for this linear system and show that a standard Krylov solver converges in a modest number of iterations that is essentially independent of the number of particles. Key to the efficiency of the method is a technique for fast computation of the product of the blob-blob mobility matrix and a vector. For unbounded suspensions, we rely on existing analytical expressions for the Rotne-Prager-Yamakawa tensor combined with a fast multipole method (FMM) to obtain linear scaling in the number of particles. For suspensions sedimented against a single no-slip boundary, we use a direct summation on a graphical processing unit (GPU), which gives quadratic asymptotic scaling with the number of particles. For fully confined domains, such as periodic suspensions or suspensions confined in slit and square channels, we extend a recently developed rigid-body immersed boundary method by B. Kallemov, A. P. S. Bhalla, B. E. Griffith, and A. Donev (Commun. Appl. Math. Comput. Sci. 11 (2016), no. 1, 79-141) to suspensions of freely moving passive or active rigid particles at zero Reynolds number. We demonstrate that the iterative solver for the coupled fluid and rigid-body equations converges in a bounded number of iterations regardless of the system size. In our approach, each iteration only requires a few cycles of a geometric multigrid solver for the Poisson equation, and an application of the block-diagonal preconditioner, leading to linear scaling with the number of particles. We optimize a number of parameters in the iterative solvers and apply our method to a variety of benchmark problems to carefully assess the accuracy of the rigid multiblob approach as a function of the resolution. We also model the dynamics of colloidal particles studied in recent experiments, such as passive boomerangs in a slit channel, as well as a pair of non-Brownian active nanorods sedimented against a wall.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2013

Embedded Sensors and Controls to Improve Component Performance and Reliability - System Dynamics Modeling and Control System Design

This report documents the current status of the modeling, control design, and embedded control research for the magnetic bearing canned rotor pump being used as a demonstration platform for deeply integrating instrumentation and controls (I&C) into nuclear power plant components. This pump is a highly inter-connected thermo/electro/mechanical system that requires an active control system to operate. Magnetic bearings are inherently unstable system and without active, moment by moment control, the rotor would contact fixed surfaces in the pump causing physical damage. This report details the modeling of the pump rotordynamics, fluid forces, electromagnetic properties of the protective cans, active magnetic bearings, power electronics, and interactions between different dynamical models. The system stability of the unforced and controlled rotor are investigated analytically. Additionally, controllers are designed using proportional derivative (PD) control, proportional integral derivative (PID) control, voltage control, and linear quadratic regulator (LQR) control. Finally, a design optimization problem that joins the electrical, mechanical, magnetic, and control system design into one problem to balance the opposing needs of various design criteria using the embedded system approach is presented.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

Nalu's Linear System Assembly using Tpetra

The Nalu Exascale Wind application assembles linear systems using data structures provided by the Tpetra package in Trilinos. This note describes the initialization and assembly process. The purpose of this note is to help Nalu developers and maintainers to understand the code surrounding linear system assembly, in order to facilitate debugging, optimizations, and maintenance.

17 WIND ENERGY↗
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 OSTI2016

CRYSTAL CHEMISTRY OF THREE-COMPONENT WHITE DWARFS AND NEUTRON STAR CRUSTS: PHASE STABILITY, PHASE STRATIFICATION, AND PHYSICAL PROPERTIES

A systematic search for multicomponent crystal structures is carried out for five different ternary systems of nuclei in a polarizable background of electrons, representative of accreted neutron star crusts and some white dwarfs. Candidate structures are “bred” by a genetic algorithm and optimized at constant pressure under the assumption of linear response (Thomas–Fermi) charge screening. Subsequent phase equilibria calculations reveal eight distinct crystal structures in the T = 0 bulk phase diagrams, five of which are complicated multinary structures not previously predicted in the context of compact object astrophysics. Frequent instances of geometrically similar but compositionally distinct phases give insight into structural preferences of systems with pairwise Yukawa interactions, including and extending to the regime of low-density colloidal suspensions made in a laboratory. As an application of these main results, we self-consistently couple the phase stability problem to the equations for a self-gravitating, hydrostatically stable white dwarf, with fixed overall composition. To our knowledge, this is the first attempt to incorporate complex multinary phases into the equilibrium phase-layering diagram and mass–radius-composition dependence, both of which are reported for He–C–O and C–O–Ne white dwarfs. Finite thickness interfacial phases (“interphases”) show up at the boundaries between single-component body-centered cubic (bcc) crystalline regions, some of which have lower lattice symmetry than cubic. A second application—quasi-static settling of heavy nuclei in white dwarfs—builds on our equilibrium phase-layering method. Tests of this nonequilibrium method reveal extra phases that play the role of transient host phases for the settling species.

79 ASTRONOMY AND ASTROPHYSICS↗