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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

The geostatistical approach for structural and stratigraphic framework analysis of offshore NW Bonaparte Basin, Australia

Geostatistics or statistical approach is based on the studies of temporal and spatial trend, which depend upon spatial relationships to model known information of variable(s) at unsampled locations. The statistical technique known as kriging was used for petrophycial and facies analysis, which help to assume spatial relationship to model the geological continuity between the known data and the unknown to produce a single best guess of the unknown. Kriging is also known as optimal interpolation technique, which facilitate to generate best linear unbiased estimation of each horizon. The idea is to construct a numerical model of the lithofacies and rock properties that honor available data and further integrate with interpreting seismic sections, techtonostratigraphy chart with sea level curve (short term) and regional tectonics of the study area to find the structural and stratigraphic growth history of the NW Bonaparte Basin. By using kriging technique the models were built which help to estimate different parameters like horizons, facies, and porosities in the study area. The variograms were used to determine for identification of spatial relationship between data which help to find the depositional history of the North West (NW) Bonaparte Basin.

58 GEOSCIENCES↗
DOE OSTI2017

Magnetic field production via the Weibel instability in interpenetrating plasma flows

Here, many astrophysical systems are effectively “collisionless,” that is, the mean free path for collisions between particles is much longer than the size of the system. The absence of particle collisions does not preclude shock formation, however, as shocks can be the result of plasma instabilities that generate and amplify electromagnetic fields. The magnetic fields required for shock formation may either be initially present, for example, in supernova remnants or young galaxies, or they may be self-generated in systems such as gamma-ray bursts (GRBs). In the case of GRB outflows, the Weibel instability is a candidate mechanism for the generation of sufficiently strong magnetic fields to produce shocks. In experiments on the OMEGA Laser, we have demonstrated a quasi-collisionless system that is optimized for the study of the non-linear phase of Weibel instability growth. Using a proton probe to directly image electromagnetic fields, we measure Weibel-generated magnetic fields that grow in opposing, initially unmagnetized plasma flows. The collisionality of the system is determined from coherent Thomson scattering measurements, and the data are compared to similar measurements of a fully collisionless system. The strong, persistent Weibel growth observed here serves as a diagnostic for exploring large-scale magnetic field amplification and the microphysics present in the collisional–collisionless transition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

HyKKT: a hybrid direct-iterative method for solving KKT linear systems

Here, we propose a solution strategy for the large indefinite linear systems arising in interior methods for nonlinear optimization. The method is suitable for implementation on hardware accelerators such as graphical processing units (GPUs). The current gold standard for sparse indefinite systems is the LBLT factorization where L is a lower triangular matrix and B is 1×1 or 2×2 block diagonal. However, this requires pivoting, which substantially increases communication cost and degrades performance on GPUs. Our approach solves a large indefinite system by solving multiple smaller positive definite systems, using an iterative solver on the Schur complement and an inner direct solve (via Cholesky factorization) within each iteration. Cholesky is stable without pivoting, thereby reducing communication and allowing reuse of the symbolic factorization. We demonstrate the practicality of our approach on large optimal power flow problems and show that it can efficiently utilize GPUs and outperform LBL T factorization of the full system.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

The high-power helicon program at DIII-D: gearing up for first experiments

Helicon current drive, also called fast wave current drive in the lower hybrid range of frequencies, has long been regarded as a promising current drive tool for reactor grade plasmas. A newly installed MW-level system at DIII-D will be the first test of this technology in reactor-relevant plasmas, in the sense that full single-pass absorption is expected. A 30-module traveling wave antenna has been installed and optimized in-vessel in early 2020. The linear electromagnetic characteristics of the unloaded module array have been extensively tested both on the bench and in the vessel at instrumentation power levels. Excellent performance has been achieved, ~2% reflected power and ~1.5% dissipated power per module in air, in a 10 MHz band around 476 MHz. Stripline feeds on both ends of the antenna allow either co or counter current drive. The installation of a 1.2 MW klystron and associated high-power electronics was completed in Fall 2020. Commissioning of the antenna is ongoing. An important goal of this experiment is to validate the helicon current drive physics basis using an extensive set of new and upgraded diagnostics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Status of ADRIANO R&D in T1015 Collaboration

The physics program for future High Energy and High Intensity experiments requires an energy resolution of the calorimetric component of detectors at limits of traditional techniques and an excellent particle identification. The novel ADRIANO technology (A Dualreadout Integrally Active Non-segmented Option), currently under development at Fermilab, is showing excellent performance on those respects. Results from detailed Monte Carlo studies on the performance with respect to energy resolution, linear response and transverse containment and a preliminary optimization of the layout are presented. A baseline configuration is chosen with an estimated energy resolution of σ(E)/E ≈ 30%/√E , to support an extensive R&D program recently started by T1015 Collaboration at Fermilab. Furthermore, preliminary results from several test beams at the Fermilab Test Beam Facility (FTBF) of a ~ 1λI prototype are presented. Future prospects with ultra-heavy glass are, also, summarized.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
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↗