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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2011

Analytic Nonlinear Methods for Beam Optics

In recent years there have been major advances in the computation and use of high-order maps for the design, optimization and operation of beamlines. We will describe five practical examples for both linear and circular colliders. From least to most complex these examples will be: use of similarity transformations (FFTB design, SLC diagnosis); statistical maps (SSC smear and tune-shift); aberrations (SLC upgrade); resonance basis and nPB tracking (PEP-II design); and kick factorization (possibly LHC).

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Design and dispatch optimization of packaged ice storage systems within a connected community

The traditional implementation of cool thermal energy storage (CTES) must be reimagined within the context of a dynamic grid and smart buildings operating as connected communities. As most buildings do not operate central chillers or connect to district cooling loops, this necessitates a broader use of packaged CTES. Our objective is to begin answering the question of how such packaged CTES should be implemented within a connected community. We do so by presenting a simulation–optimization workflow employing building energy modeling software and a mixed-integer linear program to design and dispatch a packaged CTES technology to achieve minimum total annual cost. We demonstrate this methodology on a seven-building case study using current utility rates and find that total annual cooling energy costs can be reduced by 17.8% compared to baseline, after accounting for the cost of storage. We perform three parametric sensitivity studies to evaluate modeling assumptions and obtain the prioritization of storage procurement as a function of annualized life-cycle cost of storage. We find that a community optimization approach provides significantly different results than individual building optimizations and provides greater savings compared to baseline.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2015

Vision 20/20: Simultaneous CT-MRI — Next chapter of multimodality imaging

Multimodality imaging systems such as positron emission tomography-computed tomography (PET-CT) and MRI-PET are widely available, but a simultaneous CT-MRI instrument has not been developed. Synergies between independent modalities, e.g., CT, MRI, and PET/SPECT can be realized with image registration, but such postprocessing suffers from registration errors that can be avoided with synchronized data acquisition. The clinical potential of simultaneous CT-MRI is significant, especially in cardiovascular and oncologic applications where studies of the vulnerable plaque, response to cancer therapy, and kinetic and dynamic mechanisms of targeted agents are limited by current imaging technologies. The rationale, feasibility, and realization of simultaneous CT-MRI are described in this perspective paper. The enabling technologies include interior tomography, unique gantry designs, open magnet and RF sequences, and source and detector adaptation. Based on the experience with PET-CT, PET-MRI, and MRI-LINAC instrumentation where hardware innovation and performance optimization were instrumental to construct commercial systems, the authors provide top-level concepts for simultaneous CT-MRI to meet clinical requirements and new challenges. Simultaneous CT-MRI fills a major gap of modality coupling and represents a key step toward the so-called “omnitomography” defined as the integration of all relevant imaging modalities for systems biology and precision medicine.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Pyomo.GDP: an ecosystem for logic based modeling and optimization development

We present three core principles for engineering-oriented integrated modeling and optimization tool sets—intuitive modeling contexts, systematic computer-aided reformulations, and flexible solution strategies—and describe how new developments in Pyomo.GDP for Generalized Disjunctive Programming (GDP) advance this vision. We describe a new logical expression system implementation for Pyomo.GDP allowing for a more intuitive description of logical propositions. The logical expression system supports automated reformulation of these logical constraints to linear constraints. We also describe two new logic-based global optimization solver implementations built on Pyomo.GDP that exploit logical structure to avoid “zero-flow” numerical difficulties that arise in nonlinear network design problems when nodes or streams disappear. These new solvers also demonstrate the capability to link to external libraries for expanded functionality within an integrated implementation. We present these new solvers in the context of a flexible array of solution paths available to GDP models. Finally, we present results on a new library of GDP models demonstrating the value of multiple solution approaches.

42 ENGINEERING↗
DOE OSTI2019

Advanced Control Methods for Particle Accelerators (ACM4PA) 2019

Los Alamos is currently developing novel particle accelerator controls and diagnostics algorithms to enable higher quality beams with lower beam losses than is currently possible. The purpose of this workshop was to consider tuning and optimization challenges of a wide range of particle accelerators including linear proton accelerators such as the Los Alamos Neutron Science Center (LANSCE), rings such as the Advanced Photon Source (APS) synchrotron, free electron lasers (FEL) such as the Linac Coherent Light Source (LCLS) and LCLS-II, the European X-ray Free Electron Laser (EuXFEL), the Swiss FEL, and the planned MaRIE FEL, and plasma wake- eld accelerators such as FACET, FACET-II, and AWAKE at CERN. One major challenge is an the ability to quickly create very high quality, extremely intense, custom current and energy pro le beams while working with limited real time non-invasive diagnostics and utilizing time-varying uncertain initial beam distributions and accelerator components. Currently, a few individual accelerator labs have been developing and applying their own diagnostics tools and custom control and ML algorithms for au- tomated machine tuning and optimization. The goal of this workshop was to bring together a group of accelerator physicists and accelerator related control and ML experts in order to de ne which controls and diagnostics would be most useful for existing and future accelerators and to create a plan for developing a new family of algorithms that can be shared and maintained by the community.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Adapting Optimal Velocity Tracking Control to Account for WEC Constraints and Power-Take-Off Efficiencies: Preprint

Wave energy converters (WECs) come in many different forms, from point absorbers, oscillating water columns, to bulge wave devices. This paper focuses on the control of point absorber WECs, which typically have a narrow-banded frequency response and, therefore, control is well placed to improve the energy capture of such WECs. The acausal nature of the control problem means that theoretically optimal control is almost impossible to achieve inpractice; however optimal velocity tracking (OVT) offers a simple and robust approximation to optimal control that can achieve better power capture than passive linear damping methods, albeit with necessarily higher force demands. OVT is a form of impedance matching and the magnitude of the power-take off (PTO) force demand is often not linearly proportional to the WEC velocity, which can lead to PTO force and speed combinations far from the optimalPTO efficiency. The highly non-linear PTO force and speed to efficiency mapping can, without remedial measures, severely diminish the effectiveness of OVT techniques. In this paper, limits to the movement and force of the PTO are demonstrated, allowing OVT to be more easily applied in practice. The effect of PTO efficiency is explored, and a potential solution to the problem of adapting control to account for efficiency is presented. Both aspects of the work presented highlight the requirement for co-design of the WEC, PTO, and controller.

conversion efficiency↗
DOE OSTI2017

Superresolution Full-polarimetric Imaging for Radio Interferometry with Sparse Modeling

We propose a new technique for radio interferometry to obtain superresolution full-polarization images in all four Stokes parameters using sparse modeling. The proposed technique reconstructs the image in each Stokes parameter from the corresponding full-complex Stokes visibilities by utilizing two regularization functions: the ℓ {sub 1} norm and the total variation (TV) of the brightness distribution. As an application of this technique, we present simulated linear polarization observations of two physically motivated models of M87 with the Event Horizon Telescope. We confirm that ℓ {sub 1}+TV regularization can achieve an optimal resolution of ~25%–30% of the diffraction limit λ/D{sub max}, which is the nominal spatial resolution of a radio interferometer for both the total intensity (i.e., Stokes I) and linear polarizations (i.e., Stokes Q and U). This optimal resolution is better than that obtained from the widely used Cotton–Schwab CLEAN algorithm or from using ℓ {sub 1} or TV regularizations alone. Furthermore, we find that ℓ {sub 1}+TV regularization can achieve much better image fidelity in linear polarization than other techniques over a wide range of spatial scales, not only in the superresolution regime, but also on scales larger than the diffraction limit. Our results clearly demonstrate that sparse reconstruction is a useful choice for high-fidelity full-polarimetric interferometric imaging.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Reliable clock estimation using linear weighted fusion based on pairwise broadcast synchronization

Clock synchronization in wireless sensor networks (WSNs) has been studied extensively in recent years and many protocols are put forward based on the point of statistical signal processing, which is an effective way to optimize accuracy. However, the accuracy derived from the statistical data can be improved mainly by sufficient packets exchange, which will consume the limited power resources greatly. In this paper, a reliable clock estimation using linear weighted fusion based on pairwise broadcast synchronization is proposed to optimize sync accuracy without expending additional sync packets. As a contribution, a linear weighted fusion scheme for multiple clock deviations is constructed with the collaborative sensing of clock timestamp. And the fusion weight is defined by the covariance of sync errors for different clock deviations. Extensive simulation results show that the proposed approach can achieve better performance in terms of sync overhead and sync accuracy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Degeneracy engineering for classical and quantum annealing: A case study of sparse linear regression in collider physics

Classical and quantum annealing are computing paradigms that have been proposed to solve a wide range of optimization problems. In this paper, we aim to enhance the performance of annealing algorithms by introducing the technique of degeneracy engineering, through which the relative degeneracy of the ground state is increased by modifying a subset of terms in the objective Hamiltonian. We illustrate this novel approach by applying it to the example of ℓ 0 -norm regularization for sparse linear regression, which is, in general, an NP-hard optimization problem. Specifically, we show how to cast ℓ 0 -norm regularization as a quadratic unconstrained binary optimization (QUBO) problem, suitable for implementation on annealing platforms. As a case study, we apply this QUBO formulation to energy flow polynomials in high-energy collider physics, finding that degeneracy engineering substantially improves the annealing performance. Furthermore, our results motivate the application of degeneracy engineering to a variety of regularized optimization problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Factorization Machine‐Based Active Learning for Functional Materials Design with Optimal Initial Data

The optimization of functional materials is important to enhance their properties, but their complex geometries pose great challenges to optimization. Data-driven algorithms efficiently navigate such complex design spaces by learning relationships between material structures and performance metrics to discover high-performance functional materials. Surrogate-based active learning, continually improving its surrogate model by iteratively including high-quality data points, has emerged as a cost-effective data-driven approach. Furthermore, it can be coupled with quantum computing to enhance optimization processes, especially when paired with a special form of surrogate model (i.e., quadratic unconstrained binary optimization), formulated by factorization machine (FM). However, current practices often overlook the variability in design space sizes when determining the initial data size for optimization. In this work, we investigate the optimal initial data sizes required for efficient convergence across various design space sizes. By employing averaged piecewise linear regression, we identify initiation points where convergence begins, highlighting the crucial role of employing adequate initial data in achieving efficient optimization. These results contribute to the efficient optimization of functional materials by ensuring faster convergence and reducing computational costs in FM-based active learning.

active learning↗
DOE OSTI2013

Optimal control and Galois theory

An important role is played in the solution of a class of optimal control problems by a certain special polynomial of degree 2(n−1) with integer coefficients. The linear independence of a family of k roots of this polynomial over the field Q implies the existence of a solution of the original problem with optimal control in the form of an irrational winding of a k-dimensional Clifford torus, which is passed in finite time. In the paper, we prove that for n≤15 one can take an arbitrary positive integer not exceeding [n/2] for k. The apparatus developed in the paper is applied to the systems of Chebyshev-Hermite polynomials and generalized Chebyshev-Laguerre polynomials. It is proved that for such polynomials of degree 2m every subsystem of [(m+1)/2] roots with pairwise distinct squares is linearly independent over the field Q. Bibliography: 11 titles.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Real-time dispatch optimization for concentrating solar power with thermal energy storage

Concentrating solar power (CSP) plants present a promising path towards utility-scale renewable energy. The power tower, or central receiver, configuration can achieve higher operating temperatures than other forms of CSP, and, like all forms of CSP, naturally pairs with comparatively inexpensive thermal energy storage, which allows CSP plants to dispatch electricity according to market price incentives and outside the hours of solar resource availability. Currently, CSP plants commonly include a steam Rankine power cycle and several heat exchange components to generate high-pressure steam using stored thermal energy. The efficiency of the steam Rankine cycle depends on the temperature of the plant's operating fluid, and so is a main concern of plant operators. However, the variable nature of the solar resource and the conservatism with which the receiver is operated prevent perfect control over the receiver outlet temperature. Therefore, during periods of solar variability, collection occurs at lower-than-design temperature. To support operator decisions in a real-time setting, we develop a revenue-maximizing non-convex mixed-integer, quadradically-constrained program which determines a dispatch schedule with sub-hourly time fidelity and considers temperature-dependent power cycle efficiency. The exact nonlinear formulation proves intractable for real-time decision support. Here we present exact and inexact techniques to improve problem tractability that include a hybrid nonlinear and linear formulation. Our approach admits solutions within approximately 3% of optimality, on average, within a five-minute time limit, demonstrating its usability for decision support in a real-time setting.

14 SOLAR ENERGY↗
DOE OSTI2015

Different binding modes of Cu and Pb vs. Cd, Ni, and Zn with the trihydroxamate siderophore desferrioxamine B at seawater ionic strength

The solution speciation in seawater of divalent trace metals (Cd, Cu, Ni, Pb, Zn) is dominated by strong, ostensibly metal-specific organic ligands that may play important roles in microbial metal acquisition and/or detoxification processes. We compare the effective stabilities of these metal-organic complexes to the stabilities of their complexes with a model siderophore, desferrioxamine B (DFOB). While metal-DFOB complexation has been studied in various dilute but often moderately coordinating media, for the purpose of this investigation we measured the stability constants in a non-coordinating background electrolyte at seawater ionic strength (0.7 M NaClO4). Potentiometric titrations of single metals (M) were performed in the presence of ligand (L) at different M:L molar ratios, whereupon the stability constants of multiple complexes were simultaneously determined by non-linear regression of the titration curves with FITEQL, using the optimal binding mode for each metal. Cadmium, Ni, and Zn, like trivalent Fe, sequentially form a bi-, tetra-, and hexadentate complex with DFOB as pH increases, consistent with their coordination number of 6 and regular octahedral geometry. Copper has a Jahn-Teller-distorted square-bipyramidal geometry whereas the geometry of Pb is cryptic, involving a range of bond lengths. Supported by a thermodynamic argument, our data suggest that this impedes binding of the third hydroxamate group and that the hexadentate Cu-DFOB and Pb-DFOB complex identified in earlier reports may instead be a deprotonated tetradentate complex. Absence of the hexadentate complex promotes the formation of a dinuclear (bidentate-tetradentate) complex, M2HL2+, albeit not for Pb in 0.7 M NaCl, evidently due to extensive complexation with chloride. Stabilities of the hexadentate Ni-DFOB, Zn-DFOB, and the tetradentate Pb-DFOB complex are nearly equal, yet about 2 orders of magnitude higher and 4 orders of magnitude lower than those of the hexadentate Cd-DFOB and tetradentate Cu-DFOB complex, respectively. Linear free-energy relations defined by the rare earth elements are able to predict stabilities of the Cd, Zn, and one of the Pb complexes, but underestimate those of the Ni and Cu complexes. The comparison with metal-specific organic ligands detected in seawater yields fair agreement for three of the five metals, implying that they could be siderophore-like. The Cd- and Ni-specific ligands are much stronger and may contain quite different functional groups. Calculations with MINEQL incorporating our new stability constants indicate that very high DFOB concentrations would be required to match the extent of metal-organic complexation observed in seawater, however DFOB may well represent a much broader class of structurally related ligands.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

A variant nonmonotone smoothing algorithm with improved numerical results for large-scale LWCPs

This paper considers the linear weighted complementarity problem (denoted by LWCP) introduced by Potra (SIAM J Optim 22:1634–1654, 2012). Based on two weighted smoothing functions, we propose a new nonmonotone smoothing algorithm for solving the LWCP and establish its global and local quadratic convergence without the strict complementarity assumption. Compared to existing nonmonotone smoothing algorithms, the proposed algorithm solves the linear system only approximately which can save the computation work when one solves large-scale LWCPs. Moreover, the nonmonotone line search technique adopted in this paper includes the usual monotone line search and some existing nonmonotone line searches as special cases. Numerical results show that our algorithm is considerably efficient for solving large-scale LWCPs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Sensitivity Study for Forecasting Variables of WRF-Solar Using a Tangent Linear Approach

Integrating solar generation in recent years has highlighted the need for improved accuracy in predicting solar power. Confidence in solar power forecasting can be achieved by designing an ensemble that provides reliable probabilistic information for solar radiation with reduced uncertainty and error. Ideally, ensemble members are created through the optimized perturbation of the initial conditions in numerical weather prediction (NWP) models. Tangent linear models are capable of efficiently investigating the sensitivity of solar radiation to model input parameters because they do not require individual perturbation of each variable. This sensitivity study using tangent linear models provide us the capability to identify the right variables to perturb in an ensemble prediction system. In this study, we developed tangent linear models for WRF-Solar modules that directly impact the computation of solar radiation and the simulation of cloud formation and dissipation including the Fast All-sky Model for Solar Applications (FARMS), the Noah land surface model (LSM), the Thompson microphysics, the Mello-Yamada-Nakanishi-Niino (MYNN) boundary layer parameterization, and the Deng scheme for a shallow-convection parameterization. A sensitivity analysis was conducted under various scenarios based on satellite observations and model simulations from the National Solar Radiation Data Base (NSRDB) and WRF-Solar, respectively. Critical forecasting variables that are highly sensitive to the forecasting of global horizontal irradiance (GHI), direct normal irradiance (DNI), cloud mixing ratio, cloud tendency, cloud fraction, and sensible and latent heat fluxes were determined using the relevant WRF-Solar module. This study will be used as a guidance on future research leading to high-quality probabilistic solar forecasting. In this presentation, we discuss the validation of tangent linear approach for WRF-Solar modules and illustrate how the sensitivity results are valuable in the improvement of probabilistic solar prediction.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

A solution framework for linear PDE-constrained mixed-integer problems

Abstract We present a general numerical solution method for control problems with state variables defined by a linear PDE over a finite set of binary or continuous control variables. We show empirically that a naive approach that applies a numerical discretization scheme to the PDEs to derive constraints for a mixed-integer linear program (MILP) leads to systems that are too large to be solved with state-of-the-art solvers for MILPs, especially if we desire an accurate approximation of the state variables. Our framework comprises two techniques to mitigate the rise of computation times with increasing discretization level: First, the linear system is solved for a basis of the control space in a preprocessing step. Second, certain constraints are just imposed on demand via the IBM ILOG CPLEX feature of a lazy constraint callback. These techniques are compared with an approach where the relations obtained by the discretization of the continuous constraints are directly included in the MILP. We demonstrate our approach on two examples: modeling of the spread of wildfire and the mitigation of water contamination. In both examples the computational results demonstrate that the solution time is significantly reduced by our methods. In particular, the dependence of the computation time on the size of the spatial discretization of the PDE is significantly reduced.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Optimization of the beam extraction systems for the Linac4 H{sup −} ion source

The development of the Linac 4 and its integration into CERN’s acceleration complex is part of the foreseen luminosity upgrade of the Large Hadron Collider (LHC). The goal is to inject a 160 MeV H{sup −} beam into the CERN PS Booster (PSB) in order to increase the beam brightness by a factor of 2 compared to the 50 MeV proton linac, Linac 2, that is currently in operation. The requirements for the ion source are a 45 keV H{sup −} beam of 80 mA intensity, 2 Hz repetition rate and 0.5 ms pulse length within a normalized rms-emittance of 0.25 mm· mrad. The previously installed beam extraction system has been designed for an H{sup −} ion beam intensity of 20 mA produced by an RF-volume source with an electron to H{sup −} ratio of up to 50. For the required intensity upgrades of the Linac4 ion source, a new beam extraction system is being produced and tested; it is optimized for a cesiated surface RF-source with a nominal beam current of 40 mA and an electron to H{sup −} ratio of 4. The simulations, based on the IBSIMU code, are presented. At the Brookhaven National Laboratory (BNL), a peak beam current of more than 100 mA was demonstrated with a magnetron H{sup −} source at an energy of 35 keV and a repetition rate of 2 Hz. A new extraction system is required to operate at an energy of 45 keV; simulation of a two stage extraction system dedicated to the magnetron is presented.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Controlling and Optimizing Photoinduced Charge Transfer across Ultrathin Silica Separation Membrane with Embedded Molecular Wires for Artificial Photosynthesis

Ultrathin amorphous silica membranes with embedded organic molecular wires (oligo( p -phenylenevinylene), three aryl units) provide chemical separation of incompatible catalytic environments of CO 2 reduction and H 2 O oxidation while maintaining electronic and protonic coupling between them. For an efficient nanoscale artificial photosystem, important performance criteria are high rate and directionality of charge flow. Here, the visible-light-induced charge flow from an anchored Ru bipyridyl light absorber across the silica nanomembrane to Co 3 O 4 water oxidation catalyst is quantitatively evaluated by photocurrent measurements. Charge transfer rates increase linearly with wire density, with 5 nm -2 identified as an optimal target. Accurate measurement of wire and light absorber densities is accomplished by the polarized FT-IRRAS method. Guided by density functional theory (DFT) calculations, four wire derivatives featuring electron-donating (methoxy) and -withdrawing groups (sulfonate, perfluorophenyl) with highest occupied molecular orbital (HOMO) potentials ranging from 1.48 to 0.64 V vs NHE were synthesized and photocurrents evaluated. Charge transfer rates increase sharply with increasing driving force for hole transfer from the excited light absorber to the embedded wire, followed by a decrease as the HOMO potential of the wire moves beyond the Co 3 O 4 valence band level toward more negative values, pointing to an optimal wire HOMO potential around 1.3 V vs NHE. Comparison with photocurrents of samples without nanomembrane indicates that silica layers with optimized wires are able to approach undiminished electron flux at typical solar intensities. Combined with the established high proton conductivity and small-molecule blocking property, the charge transfer measurements demonstrate that oxidation and reduction catalysis can be efficiently integrated on the nanoscale under separation by an ultrathin silica membrane.

photocurrent measurements↗