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,711 records · Page 95
All sourcesNASA NTRSDOE OSTI
NASA NTRS2022

Multi-Robot Assembly Scheduling for the Lunar Crater Radio Telescope on the Far-Side of the Moon

The Lunar Crater Radio Telescope (LCRT) is a pro- posed ultra-long-wavelength radio telescope to be constructed on the far side of the moon. The proposed telescope will be constructed by deploying a 1km wire mesh in a 3-5km crater using a team of wall-climbing DuAxel robots. In this work, we consider the problem of generating minimum-time assembly sequences for LCRT, using realistic models of travel speed and lighting. Specifically, we pose the assembly sequencing problem as a mixed-integer linear program (MILP), which we solve to global optimality using commercial solvers. We present methods for modeling time-varying travel and assembly times, based on variable lighting conditions (including crater shadowing), and show how such time-varying parameters can be incorporated into the MILP. Finally, we present numerical studies of our method, showing how makespan varies with the number of assembly robots.

Schwager, Mac↗
DOE OSTI2024

Evaluating Nationwide Supply Chain for Circularity of PET and Olefin Plastics

PET (Polyethylene Terephthalate, #1) and olefin plastics including HDPE (High Density Polyethylene, #2), LDPE/LLDPE (Low Density/Linear Low-Density Polyethylene, #4) and PP (Polypropylene, #5) together comprise nearly 80% of the U.S. plastic market. Due to their high market share and extensive application in packaging these polymers have better potential for circularity than other polymer types. To understand the potential of future scenarios with higher recycling rates, supply chain scenarios for PET and olefin plastic packaging need to be analyzed with increased availability and collection of plastics. Currently there exists a knowledge gap to understand the circular supply chain on a national level. Performing a nationwide analysis introduces certain challenges such as variability in the plastic mix of the recycling stream in different municipalities, lack of collection of some types of plastics, regional differences in cost and estimating a national supply curve based on state level access rate and participation rate. Currently, a limited number of Plastic Reclaimers recycle the nation’s collected plastic involving transportation over long distances surpassing state boundaries. Modeling a nationwide scenario instead of regional/state based scenario will facilitate transportation beyond state boundaries for the development of a circular economy. A Mixed-integer Linear Programming model was developed, to identify the optimal location and capacity of the Material Recovery Facilities (MRFs) nationwide subject to maximizing the profit margin of the industrial entities within the model. We considered three different scenarios including the base case scenario (S1) collecting 2.27 million metric tons/year as well as two additional scenarios with available plastic supply of 1.8 times (S2) and 2.2 times (S3) of the base case scenario. The model identified 166, 275 and 319 counties as potential MRF locations for scenarios S1, S2 and S3 respectively. Compared to the existing number of counties having MRFs in the US, the model results indicated a reduction of 38% of MRFs for S1. For S2 and S3, the number of counties with MRFs increased 3% and 20% respectively compared to counties currently with MRFs. The average profit remains constant between $180-$181/ton regardless of the increased plastic collection.

54 - ENVIRONMENTAL SCIENCES/GLOBAL CLIMATE CHANGE ↗
DOE OSTI2019

Online Primal-Dual Methods With Measurement Feedback for Time-Varying Convex Optimization

This paper addresses the design and analysis of feedback-based online algorithms to control systems or networked systems based on performance objectives and engineering constraints that may evolve over time. The emerging timevarying convex optimization formalism is leveraged to model optimal operational trajectories of the systems, as well as explicit local and network-level operational constraints. Departing from existing batch and feed-forward optimization approaches, the design of the algorithms capitalizes on an online implementation of primal-dual projected-gradient methods; the gradient steps are, however, suitably modified to accommodate feedback from the system in the form of measurements - hence, the term 'online optimization with feedback.' By virtue of this approach, the resultant algorithms can cope with model mismatches in the algebraic representation of the system states and outputs, they avoid pervasive measurements of exogenous inputs, and they naturally lend themselves to a distributed implementation. Under suitable assumptions, analytical convergence claims are established in terms of dynamic regret. Furthermore, when the synthesis of the feedback-based online algorithms is based on a regularized Lagrangian function, Q-linear convergence to solutions of the time-varying optimization problem is shown.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1988

Optimization and guidance of penetration landing trajectories in a windshear

The optimization and guidance of penetration landing trajectories in a windshear are considered. It is assumed that the aircraft is controlled by the angle of attack and the power setting. For the optimal trajectory, the performance index being minimized measures the deviation of the flight trajectory from the nominal trajectory. In turn, the nominal trajectory includes two parts: the approach part (nominal glide slope constant) and the flare part (nominal glide slope varying linearly with the horizontal distance). Numerical results show that the optimal trajectory deviates somewhat from the nominal trajectory in the shear region. A guidance scheme is developed to approximate the optimal trajectory. The angle of attack is determined by the windshear intensity, the absolute path inclination, and the glide slope angle, while the power setting is determined by the windshear intensity and the velocity. Numerical results indicate that the guidance trajectory is close to the optimal trajectory.

Miele, A.↗
DOE OSTI2026

A Linear Programming Approach to Backtracking for Single-Axis Trackers on Rolling Terrain

In this article, we present a computationally efficient method for determining optimal backtracking rotations for single-axis solar trackers on nonuniform terrain. The method allows for ganged tracking, mechanical rotation constraints, uneven row spacing, and arbitrary maximum allowable shaded fractions (to enable “fractional backtracking”). As with previous 2-D approaches, the method is suitable for terrain that varies in the transverse direction with respect to the rotation axis of the trackers. The novelty of the method lies in formulating the problem of shade avoidance as a linear problem, which is achieved by using the row interception width as the optimization variable instead of rotation angles. Formulating backtracking as a linear problem enables the use of extremely efficient linear programming algorithms, making the method highly scalable, requiring less than 1 min to compute optimal rotation schedules for hundreds of trackers. It also produces more effective backtracking rotations, reducing the frequency of shading by 4× and improving system energy output by 1%–2%.

Optimization↗
NASA NTRS1994

Parameter estimation of the vibrational model for the SCOLE experimental facility

The objective of this study is to experimentally determine an empirical model of the vibrational dynamics of the Spacecraft COntrol Laboratory Experiment (SCOLE) facility. The first two flexible modes of this test article are identified using a linear least-square identification procedure and the data utilized for this procedure are obtained by exciting the structure from a quiescent state with torque wheels. The time history data of rate gyro sensors and accelerometers due to excitation and after excitation in terms of free-decay are used in the parameter estimation of the vibrational model. The free-decay portion of the data is analyzed using the Discrete Fourier transform to determine the optimal model order to use in modelling the response. Linear least-square analysis is then used to select the parameters that best fit the output of an Autoregressive (AR) model to the data. The control effectiveness of the torque wheels is then determined using the excitation portion of the test data, again using linear least squares.

Crotts, B. D.↗
DOE OSTI2020

Optimal Equilibrium Selection of Price-Maker Agents in Performance-Based Regulation Market

This paper analyzes the oligopolistic equilibrium of multiple price-maker agents in performance-based regulation (PBR) markets. In these markets, there are price-maker agents representing some frequency regulation (FR) providers and a number of independent price-taker FR providers. An equilibrium problem with equilibrium constraints (EPEC) model is employed in this paper to study the equilibria of a PBR market in the presence of price-maker agents and pricetaker FR providers. Due to the incorporation of the FR providers' dynamics, the proposed model is reformulated as a mixed-integer linear programming (MILP) problem over innovative mathematical techniques. An optimal equilibrium point is also selected for the market, where none of the agents is unique deviator and the dynamic performance of power system is improved simultaneously. The effectiveness of proposed optimal equilibrium point is evaluated in the numerical results section by comparing the outputs with conventional optimal dispatches of the FR providers.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2014

SU-E-T-197: Helical Cranial-Spinal Treatments with a Linear Accelerator

Purpose: Craniospinal irradiation (CSI) of systemic disease requires a high level of beam intensity modulation to reduce dose to bone marrow and other critical structures. Current helical delivery machines can take 30 minutes or more of beam-on time to complete these treatments. This pilot study aims to test the feasibility of performing helical treatments with a conventional linear accelerator using longitudinal couch travel during multiple gantry revolutions. Methods: The VMAT optimization package of the Eclipse 10.0 treatment planning system was used to optimize pseudo-helical CSI plans of 5 clinical patient scans. Each gantry revolution was divided into three 120° arcs with each isocenter shifted longitudinally. Treatments requiring more than the maximum 10 arcs used multiple plans with each plan after the first being optimized including the dose of the others (Figure 1). The beam pitch was varied between 0.2 and 0.9 (couch speed 5- 20cm/revolution and field width of 22cm) and dose-volume histograms of critical organs were compared to tomotherapy plans. Results: Viable pseudo-helical plans were achieved using Eclipse. Decreasing the pitch from 0.9 to 0.2 lowered the maximum lens dose by 40%, the mean bone marrow dose by 2.1% and the maximum esophagus dose by 17.5%. (Figure 2). Linac-based helical plans showed dose results comparable to tomotherapy delivery for both target coverage and critical organ sparing, with the D50 of bone marrow and esophagus respectively 12% and 31% lower in the helical linear accelerator plan (Figure 3). Total mean beam-on time for the linear accelerator plan was 8.3 minutes, 54% faster than the tomotherapy average for the same plans. Conclusions: This pilot study has demonstrated the feasibility of planning pseudo-helical treatments for CSI targets using a conventional linac and dynamic couch movement, and supports the ongoing development of true helical optimization and delivery.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Maximizing dynamic range and performance of anatase TiO 2 ECRAM through structure and programming

Here, in this study, we investigate the structure-dependent modulation characteristics of all-solid-state three-terminal electrochemical random-access memory (ECRAM) based on an anatase Li x TiO 2 channel. By directly comparing “asymmetric” and “symmetric” ECRAM device architectures, we reveal significant insight into the impact of a non-zero gate-drain open-circuit voltage and its influence on voltage vs. current-controlled gating. We also explore the impact of potentiation/depression write parameters on the symmetry, linearity, and dynamic range of the device response. Together, initial results from optimizing structure and programming approaches yielded unprecedented G max /G min ratios of >1,000 for ECRAM and hundreds of tunable memory states with excellent linearity and symmetry. Simulations based on these ECRAM devices further illustrate the promise of this analog memory technology, achieving near 2% classification error in the MNIST digit recognition benchmark for a range of training parameters compared to a theoretical best of 1.66% and outperforming other device models extracted from the literature.

AIHWKit↗
NASA NTRS1983

Optimal processing of satellite-derived magnetic anomaly data

It is shown how the concept of the power spectrum can be extended to two-dimensional spatial power spectra and how it can be used in determining optimal data processing methods for satellite-derived magnetic anomaly data and planning missions to obtain such data. The analysis techniques are applied to the data set and data-processing procedure described by Mayhew et al. (1980), a study that treats magnetic anomaly data for Australia and the surrounding ocean obtained by the polar orbit POGO series satellites. It is shown that the data-processing method used by Mayhew et al. is approximately equivalent to an invariant two-dimensional linear filter and that it is reasonably close to optimal with respect to accuracy, although some possible improvements are suggested. However, as is usual when filtering data, some real 'signal' is unavoidably removed along with the 'noise' resulting in errors that can be quite large. A method for reducing these errors by using additional data from a medium inclination orbit satellite (for example, 60 deg inclination) is proposed.

Mcleod, M. G.↗
DOE OSTI2023

Multi-output multilevel best linear unbiased estimators via semidefinite programming

Multifidelity forward uncertainty quantification (UQ) problems often involve multiple quantities of interest and heterogeneous models (e.g., different grids, equations, dimensions, physics, surrogate and reduced-order models). While computational efficiency is key in this context, multi-output strategies in multilevel/multifidelity methods are either sub-optimal or non-existent. In this paper we extend multilevel best linear unbiased estimators (MLBLUE) to multi-output forward UQ problems and we present new semidefinite programming formulations for their optimal setup. Not only do these formulations yield the optimal number of samples required, but also the optimal selection of low-fidelity models to use. While existing MLBLUE approaches are single-output only and require a non-trivial nonlinear optimization procedure, the new multi-output formulations can be solved reliably and efficiently. Here, we demonstrate the efficacy of the new methods and formulations in practical UQ problems with model heterogeneity.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1978

Invariants of optimal minimal-order observer-based compensators

The relations characterizing the optimal minimal-order observer-based compensators for a linear time-invariant multivariable system with a random initial state (or, equivalently, known initial state and white plant-driving noise) have been reported by Miller (1973). In this note we establish in general that the plant transfer function uniquely determines the optimal compensator transfer function, and that this characterizes precisely the degrees of freedom in the compensator design; computational implications of this result are indicated.

Blanvillain, P. J.↗
DOE OSTI2021

Transient Efficiency Flexibility and Reliability Optimization of Coal-Fired Power Plants: Model-Predictive Control Library Development (Report)

This document pertains to the reporting requirements of DOE contract FE-0031767. The document covers the development of a model predictive control (MPC) library for implementing MPC for a general dynamic system. The library is implemented in a standardized manner in Matlab/Simulink, where core functions on model prediction, linearization and formulation and solution of a quadratic programming (QP) optimization problem is done in the core library - independent of the specific application. The user can provide the application-specific dynamic model in continuous and discrete time, to rapidly implement and test the MPC performance in a desktop simulation. The MPC optimization objective and constraints are also easily configured via an Excel file to allow iterative refinement as needed. Finally, the MPC library enables a rapid deployment to a target environment through auto C-code generation and containerization. The MPC library works seamlessly with the model based estimation (MBE) library to obtain the overall output feedback control solution. In this program, the reduced order model (ROM) of a coal-fired power plant (CFPP) is used to implement and test the MPC solution.

01 COAL, LIGNITE, AND PEAT↗
NASA NTRS1979

Linear feature selection with applications

Several ways in which feature selection techniques were used in LACIE are discussed. In all cases, the methods require some a priori information and assumptions; in most, the classification procedure (Bayes optimal) was chosen in advance. The transformations used for dimensionality reduction are linear, that is, the variables in feature space are always linear combinations of the original measurements. Several numerically tractable criteria developed for LACIE, which provide information about the probability of misclassification, are discussed. Recent results on linear feature selection techniques are included. Their use in LACIE is discussed. Related open questions are mentioned.

Decell, H. P., Jr.↗
DOE OSTI2014

An optimal guarding scheme for thermal conductivity measurement using a guarded cut-bar technique, part 1 experimental study

In the guarded cut-bar technique, a guard surrounding the measured sample and reference (meter) bars is temperature controlled to carefully regulate heat losses from the sample and reference bars. Guarding is typically carried out by matching the temperature profiles between the guard and the test stack of sample and meter bars. Problems arise in matching the profiles, especially when the thermal conductivitiesof the meter bars and of the sample differ, as is usually the case. In a previous numerical study, the applied guarding condition (guard temperature profile) was found to be an important factor in measurement accuracy. Different from the linear-matched or isothermal schemes recommended in literature, the optimal guarding condition is dependent on the system geometry and thermal conductivity ratio of sample to meter bar. To validate the numerical results, an experimental study was performed to investigate the resulting error under different guarding conditions using stainless steel 304 as both the sample and meter bars. The optimal guarding condition was further verified on a certified reference material, pyroceram 9606, and 99.95% pure iron whose thermal conductivities are much smaller and much larger, respectively, than that of the stainless steel meter bars. Additionally, measurements are performed using three different inert gases to show the effect of the insulation effective thermal conductivity on measurement error, revealing low conductivity, argon gas, gives the lowest error sensitivity when deviating from the optimal condition. The result of this study provides a general guideline for the specific measurement method and for methods requiring optimal guarding or insulation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2017

Hydrodynamic Design Optimization and Wave Tank Testing of a Self-Reacting Two-Body Wave Energy Converter

As worldwide energy consumption continues to increase, so does the demand for renewable energy sources. The total available wave energy resource for the United States alone is 2,640 TWh/yr; nearly two thirds of the 4,000 TWh of electricity used in the United States each year. It is estimated that nearly half of that available energy is recoverable through wave energy conversion techniques. In this thesis, a two-body ‘point absorber’ type wave energy converter with a mechanical power-takeoff is investigated. The two-body wave energy converter extracts energy through the relative motion of a floating buoy and a neutrally buoyant submerged body. Using a linear frequency-domain model, analytical solutions of the optimal power and the corresponding power-takeoff components are derived for the two-body wave energy converter. Using these solutions, a case study is conducted to investigate the influence of the submerged body size on the absorbed power of the device in regular and irregular waves. Here it is found that an optimal mass ratio between the submerged body and floating buoy exists where the device will achieve resonance. Furthermore, a case study to investigate the influence of the submerged body shape on the absorbed power is conducted using a time-domain numerical model. Here it is found that the submerged body should be designed to reduce the effects of drag, but to maintain relatively large hydrodynamic added mass and excitation force. To validate the analytical and numerical models, a 1/30th scale model of a two-body wave energy converter is tested in a wave tank. The results of the wave tank tests show that the two-body wave energy converter can absorb nearly twice the energy of a single-body ‘point absorber’ type wave energy converter.

16 TIDAL AND WAVE POWER↗
DOE OSTI2016

Development and clinical evaluation of an ionization chamber array with 3.5 mm pixel pitch for quality assurance in advanced radiotherapy techniques

Purpose: To characterize a new air vented ionization chamber technology, suitable to build detector arrays with small pixel pitch and independence of sensitivity on dose per pulse. Methods: The prototype under test is a linear array of air vented ionization chambers, consisting of 80 pixels with 3.5 mm pixel pitch distance and a sensitive volume of about 4 mm{sup 3}. The detector has been characterized with {sup 60}Co radiation and MV x rays from different linear accelerators (with flattened and unflattened beam qualities). Sensitivity dependence on dose per pulse has been evaluated under MV x rays by changing both the source to detector distance and the beam quality. Bias voltage has been varied in order to evaluate the charge collection efficiency in the most critical conditions. Relative dose profiles have been measured for both flattened and unflattened distributions with different field sizes. The reference detectors were a commercial array of ionization chambers and an amorphous silicon flat panel in direct conversion configuration. Profiles of dose distribution have been measured also with intensity modulated radiation therapy (IMRT), stereotactic radiosurgery (SRS), and volumetric modulated arc therapy (VMAT) patient plans. Comparison has been done with a commercial diode array and with Gafchromic EBT3 films. Results: Repeatability and stability under continuous gamma irradiation are within 0.3%, in spite of low active volume and sensitivity (∼200 pC/Gy). Deviation from linearity is in the range [0.3%, −0.9%] for a dose of at least 20 cGy, while a worsening of linearity is observed below 10 cGy. Charge collection efficiency with 2.67 mGy/pulse is higher than 99%, leading to a ±0.9% sensitivity change in the range 0.09–2.67 mGy/pulse (covering all flattened and unflattened beam qualities). Tissue to phantom ratios show an agreement within 0.6% with the reference detector up to 34 cm depth. For field sizes in the range 2 × 2 to 15 × 15 cm{sup 2}, the output factors are in agreement with a thimble chamber within 2%, while with 25 × 25 cm{sup 2} field size, an underestimation of 4.0% was found. Agreement of field and penumbra width measurements with the flat panel is of the order of 1 mm down to 1 × 1 cm{sup 2} field size. Flatness and symmetry values measured with the 1D array and the reference detectors are comparable, and differences are always smaller than 1%. Angular dependence of the detector, when compared to measurements taken with a cylindrical chamber in the same phantom, is as large as 16%. This includes inhomogeneity and asymmetry of the design, which during plan verification are accounted for by the treatment planning system (TPS). The detector is capable to reproduce the dose distributions of IMRT and VMAT plans with a maximum deviation from TPS of 3.0% in the target region. In the case of VMAT and SRS plans, an average (maximum) deviation of the order of 1% (4%) from films has been measured. Conclusions: The investigated technology appears to be useful both for Linac QA and patient plan verification, especially in treatments with steep dose gradients and nonuniform dose rates such as VMAT and SRS. Major limitations of the present prototype are the linearity at low dose, which can be solved by optimizing the readout electronics, and the underestimation of output factors with large field sizes. The latter problem is presently not completely understood and will require further investigations.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Optimal Design of Food Packaging Considering Waste Management Technologies to Achieve Circular Economy

Plastic packaging plays a fundamental role in the food industry, avoiding food waste and facilitating food access. The increasing plastic production and the lack of appropriate plastic waste management technologies represent a threat to the environmental and human welfare. Therefore, there is an urgent need to identify sustainable packaging solutions. Circular economy (CE) promotes reducing waste and increasing recycling practices to achieve sustainability. In this work, we propose a CE framework based on multi-objective optimization, considering both economic and environmental impacts, to identify optimal packaging designs and waste management technologies. Using mixed-integer linear programming (MILP), techno-economic analysis (TEA), and life cycle assessment (LCA), this work aims to build the first steps in packaging design, informing about the best packaging alternatives and the optimal technology or technologies to process packaging waste. For the economic analysis, we consider the minimum increase in price (MIP) when adding recycling to the cost of each packaging solution, while for the environmental analysis, the greenhouse gas emissions impact was considered. A case study on ground coffee packaging is used to illustrate the proposed framework. The results demonstrate that the multilayer bag option is the most convenient when considering both the chosen economic and environmental impacts.

Life Cycle Analysis↗