Search NASA⌕ Search

SEARCH · Search NASA

Results for “optimization algorithm”

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

Area Coverage Planning with 3-axis Steerable, 2D Framing Sensors

Existing algorithms for Agile Earth Observing Satellites((Lemaitre et al. 2002)) were largely created for 1D line sensors that acquire images in linear swaths. However, imaging satellites increasingly use 2D framing sensors (cameras) that capture discrete rectangular images. We describe tiling step-stare approaches that are more suited to rectangular image footprints than are 1D swath-based algorithms. Optimal area planning for these 2D framing instruments is an NPcomplete problem and intractable for large areas, so we present four approximation algorithms. Strategies are compared against a prior 2D framing instrument algorithm (Knight 2014) in three computational experiments. The impact of observer agility on schedule makespan is examined. Makespans vary more as observer agility decreases toward a critical point, then vary less after the critical point, suggesting a possible problem phase transition.

Chien, Steve↗

The Computer Program LIAR for Beam Dynamics Calculations in Linear Accelerators

Linear accelerators are the central components of the proposed next generation of linear colliders. They need to provide acceleration of up to 750 GeV per beam while maintaining very small normalized emittances. Standard simulation programs, mainly developed for storage rings, do not meet the specific requirements for high energy linear accelerators. We present a new program LIAR ('LInear Accelerator Research code') that includes wakefield effects, a 6D coupled beam description, specific optimization algorithms and other advanced features. Its modular structure allows to use and to extend it easily for different purposes. The program is available for UNIX workstations and Windows PC's. It can be applied to a broad range of accelerators. We present examples of simulations for SLC and NLC.

43 PARTICLE ACCELERATORS↗

Optimal probe selection in diagnostic search

Probe selection (PS) in machine diagnosis is viewed as a collection of models that apply under specific conditions. This makes it possible for three polynomial-time optimal algorithms to be developed for simplified PS models that allow different probes to have different costs. The work is compared with the research of Simon and Kadane (1975), who developed a collection of models for optimal problem-solving search. The relationship between these models and the three newly developed algorithms for PS is explored. Two of the algorithms are unlike the ones discussed by Simon and Kadane. The third cannot be related to the problem-solving models.

Bhandari, Inderpal S.↗

An adaptive, multivariate partitioning algorithm for global optimization of nonconvex programs

We develop an adaptive, multivariate partitioning algorithm for solving nonconvex, Mixed-Integer Nonlinear Programs (MINLPs) with polynomial functions to global optimality. In particular, we present an iterative algorithm that exploits piecewise, convex relaxation approaches via disjunctive formulations to solve MINLPs that is different than conventional spatial branch-and-bound approaches. The algorithm partitions the domains of variables in an adaptive and non-uniform manner at every iteration to focus on productive areas of the search space. Furthermore, domain reduction techniques based on sequential, optimization-based bound-tightening and piecewise relaxation techniques, as a part of a presolve step, are integrated into the main algorithm. Finally, we demonstrate the effectiveness of the algorithm on well-known benchmark problems (including Pooling and Blending instances) from MINLPLib and compare our algorithm with state-of-the-art global optimization solvers. With our novel approach, we solve several large-scale instances, some of which are not solvable by state-of-the-art solvers. We also succeed in reducing the best known optimality gap for a hard, generalized pooling problem instance.

97 MATHEMATICS AND COMPUTING↗

Sequence of polyhedral relaxations for nonlinear univariate functions

Here, given a nonlinear, univariate, bounded, and differentiable function f(x), this article develops a sequence of Mixed Integer Linear Programming (MILP) and Linear Programming (LP) relaxations that converge to the graph of f(x) and its convex hull, respectively. Theoretical convergence of the sequence of relaxations to the graph of the function and its convex hull is established. For nonlinear non-convex optimization problems, the relaxations presented in this article can be used to construct tight MILP and LP relaxations. These MILP and the LP relaxations can also be used with MILP-based and spatial branch-and-bound based global optimization algorithms, respectively.

42 ENGINEERING↗

Chance-Constrained System of Systems Based Operation of Power Systems

In this paper, a chance-constrained system of systems (SoS) based decision-making approach is presented for stochastic scheduling of power systems encompassing active distribution grids. Based on the concept of SoS, the independent system operator (ISO) and distribution companies (DISCOs) are modeled as self-governing systems. These systems collaborate with each other to run the entire power system in a secure and economic manner. Each self-governing system accounts for its local reserve requirements and line flow constraints with respect to the uncertainties of load and renewable energy resources. A set of chance constraints are formulated to model the interactions between the ISO and DISCOs. The proposed model is solved by using analytical target cascading (ATC) method, a distributed optimization algorithm in which only a limited amount of information is exchanged between collaborative ISO and DISCOs. In this paper, a 6-bus and a modified IEEE 118-bus power systems are studied to show the effectiveness of the proposed algorithm.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Pilot guidance and display considerations for energy efficient flight profiles

Two computer programs are applied to energy efficient flight operation in order to minimize aircraft operating costs. One algorithm (OPTIM) computes vertical flight profiles which optimize direct operating costs, including fuel and time costs, for an aircraft flying over a fixed range and with a fixed time-of-arrival. The second program (TRAGEN) simulates an aircraft steered to fly along a specified vertical flight trajectory, in order to examine fuel and cost penalties involved in flying nonoptimal trajectories. Constraints such as air traffic control procedures, and atmospheric and weight conditions are considered and supported by graphs and diagrams. The use of the algorithms as preflight planning tools is discussed, emphasizing OPTIM's future application for on-board energy management. Finally, research questions concerning pilot guidance and display considerations for advanced energy/flight management systems are addressed.

Samms, K. H.↗

The role of under-determined approximations in engineering and science application

There is currently a great deal of interest in using response surfaces in the optimization of aircraft performance. The objective function and/or constraint equations involved in these optimization problems may come from numerous disciplines such as structures, aerodynamics, environmental engineering, etc. In each of these disciplines, the mathematical complexity of the governing equations usually dictates that numerical results be obtained from large computer programs such as a finite element method program. Thus, when performing optimization studies, response surfaces are a convenient way of transferring information from the various disciplines to the optimization algorithm as opposed to bringing all the sundry computer programs together in a massive computer code. Response surfaces offer another advantage in the optimization of aircraft structures. A characteristic of these types of optimization problems is that evaluation of the objective function and response equations (referred to as a functional evaluation) can be very expensive in a computational sense. Because of the computational expense in obtaining functional evaluations, the present study was undertaken to investigate under-determinined approximations. An under-determined approximation is one in which there are fewer training pairs (pieces of information about a function) than there are undetermined parameters (coefficients or weights) associated with the approximation. Both polynomial approximations and neural net approximations were examined. Three main example problems were investigated: (1) a function of one design variable was considered; (2) a function of two design variables was considered; and (3) a 35 bar truss with 4 design variables was considered.

Carpenter, William C.↗

Structural design utilizing updated, approximate sensitivity derivatives

A method to improve the computational efficiency of structural optimization algorithms is investigated. In this method, the calculations of 'exact' sensitivity derivatives of constraint functions are performed only at selected iterations during the optimization process. The sensitivity derivatives utilized within other iterations are approximate derivatives which are calculated using an inexpensive derivative update formula. Optimization results are presented for an analytic optimization problem (i.e., one having simple polynomial expressions for the objective and constraint functions) and for two structural optimization problems. The structural optimization results indicate that up to a factor of three improvement in computation time is possible when using the updated sensitivity derivatives.

Scotti, Stephen J.↗

Removing cosmic-ray hits from multiorbit HST Wide Field Camera images

We present an optimized algorithm that removes cosmic rays ('CRs') from multiorbit Hubble Space Telescope (HST) Wide Field/Planetary Camera ('WF/PC') images. It computes the image noise in every iteration from the WF/PC CCD equation. This includes all known sources of random and systematic calibration errors. We test this algorithm on WF/PC stacks of 2-12 orbits as a function of the number of available orbits and the formal Poissonian sigma-clipping level. We find that the algorithm needs greater than or equal 4 WF/PC exposures to locate the minimal sky signal (which is noticeably affected by CRs), with an optimal clipping level at 2-2.5 x sigma(sub Poisson). We analyze the CR flux detected on multiorbit 'CR stacks,' which are constructed by subtracting the best CR filtered images from the unfiltered 8-12 orbit average. We use an automated object finder to determine the surface density of CRS as a function of the apparent magnitude (or ADU flux) they would have generated in the images had they not been removed. The power law slope of the CR 'counts' (gamma approximately = 0.6 for N(m) m(exp gamma)) is steeper than that of the faint galaxy counts down to V approximately = 28 mag. The CR counts show a drop off between 28 less than or approximately V less than or approximately 30 mag (the latter is our formal 2 sigma point source sensitivity without spherical aberration). This prevents the CR sky integral from diverging, and is likely due to a real cutoff in the CR energy distribution below approximately 11 ADU per orbit. The integral CR surface density is less than or approximately 10(exp 8)/sq. deg, and their sky signal is V approximately = 25.5-27.0 mag/sq. arcsec, or 3%-13% of our NEP sky background (V = 23.3 mag/sq. arcsec), and well above the EBL integral of the deepest galaxy counts (B(sub J) approximately = 28.0 mag/sq. arcsec). We conclude that faint CRs will always contribute to the sky signal in the deepest WF/PC images. Since WFPC2 has approximately 2.7x lower read noise and a thicker CCD, this will result in more CR detections than in WF/PC, potentially affecting approximately 10%-20% of the pixels in multiorbit WFPC2 data cubes.

Windhorst, Rogier A.↗

Parallelization of Program to Optimize Simulated Trajectories (POST3D)

This paper describes the parallelization of the Program to Optimize Simulated Trajectories (POST3D). POST3D uses a gradient-based optimization algorithm that reaches an optimum design point by moving from one design point to the next. The gradient calculations required to complete the optimization process, dominate the computational time and have been parallelized using a Single Program Multiple Data (SPMD) on a distributed memory NUMA (non-uniform memory access) architecture. The Origin2000 was used for the tests presented.

Hammond, Dana P.↗

TH-C-BRD-12: Robust Intensity Modulated Proton Therapy Plan Can Eliminate Junction Shifts for Craniospinal Irradiation

Purpose: The passive scattering proton therapy (PSPT) technique is the commonly used radiotherapy technique for craniospinal irradiation (CSI). However, PSPT involves many numbers of junction shifts applied over the course of treatment to reduce the cold and hot regions caused by field mismatching. In this work, we introduced a robust planning approach to develop an optimal and clinical efficient techniques for CSI using intensity modulated proton therapy (IMPT) so that junction shifts can essentially be eliminated. Methods: The intra-fractional uncertainty, in which two overlapping fields shift in the opposite directions along the craniospinal axis, are incorporated into the robust optimization algorithm. Treatment plans with junction sizes 3,5,10,15,20,25 cm were designed and compared with the plan designed using the non-robust optimization. Robustness of the plans were evaluated based on dose profiles along the craniospinal axis for the plans applying 3 mm intra-fractional shift. The dose intra-fraction variations (DIV) at the junction are used to evaluate the robustness of the plans. Results: The DIVs are 7.9%, 6.3%, 5.0%, 3.8%, 2.8% and 2.2%, for the robustly optimized plans with junction sizes 3,5,10,15,20,25 cm. The DIV are 10% for the non-robustly optimized plans with junction size 25 cm. The dose profiles along the craniospinal axis exhibit gradual and tapered dose distribution. Using DIVs less than 5% as maximum acceptable intrafractional variation, the overlapping region can be reduced to 10 cm, leading to potential reduced number of the fields. The DIVs are less than 5% for 5 mm intra-fractional shifts with junction size 25 cm, leading to potential no-junction-shift for CSI using IMPT. Conclusion: This work is the first report of the robust optimization on CSI based on IMPT. We demonstrate that robust optimization can lead to much efficient carniospinal irradiation by eliminating the junction shifts.

60 APPLIED LIFE SCIENCES↗

Performance seeking control excitation mode

Flight testing of the performance seeking control (PSC) excitation mode was successfully completed at NASA Dryden on the F-15 highly integrated digital electronic control (HIDEC) aircraft. Although the excitation mode was not one of the original objectives of the PSC program, it was rapidly prototyped and implemented into the architecture of the PSC algorithm, allowing valuable and timely research data to be gathered. The primary flight test objective was to investigate the feasibility of a future measurement-based performance optimization algorithm. This future algorithm, called AdAPT, which stands for adaptive aircraft performance technology, generates and applies excitation inputs to selected control effectors. Fourier transformations are used to convert measured response and control effector data into frequency domain models which are mapped into state space models using multiterm frequency matching. Formal optimization principles are applied to produce an integrated, performance optimal effector suite. The key technical challenge of the measurement-based approach is the identification of the gradient of the performance index to the selected control effector. This concern was addressed by the excitation mode flight test. The AdAPT feasibility study utilized the PSC excitation mode to apply separate sinusoidal excitation trims to the controls - one aircraft, inlet first ramp (cowl), and one engine, throat area. Aircraft control and response data were recorded using on-board instrumentation and analyzed post-flight. Sensor noise characteristics, axial acceleration performance gradients, and repeatability were determined. Results were compared to pilot comments to assess the ride quality. Flight test results indicate that performance gradients were identified at all flight conditions, sensor noise levels were acceptable at the frequencies of interest, and excitations were generally not sensed by the pilot.

Schkolnik, Gerard↗

Microhole drilling technology utilizing a golden section search algorithm

A fundamental issue in microhole drilling is that delivering high weight-on-bit (WOB), high torque rotational horsepower to a conventional drill bit does not scale down to the hole sizes necessary to realize the envisioned cost savings An optimization algorithm called a golden section search (GSS) was used to systematically identify the preferred WOB for a given set of conditions. This research focused on implementing and evaluating two low WOB drilling technologies for microhole drilling: - Laser-assisted mechanical drill, which was tested in the laboratory - Lightly modified commercial off-the-shelf (COTS) percussive hammer, which was tested in a limited field test. Data were collected for microhole GSS using WOB optimization via simulation as well as at the Blue Canyon Dome Site in Socorro, NM. Information on the attached files and folders are as follows: - the .tdms files are LabView data files, which can be opened within Excel using a .tdms add-in or using a Matlab .tdms converter - the .tdms_index files are part of the .tdms file structure - sampling rate, column headers, and length data within the .tdms files follow SOP when utilizing Excel and/or Matlab as described above

15 GEOTHERMAL ENERGY↗

Using a derivative-free optimization method for multiple solutions of inverse transport problems

Identifying unknown components of an object that emits radiation is an important problem for national and global security. Radiation signatures measured from an object of interest can be used to infer object parameter values that are not known. This problem is called an inverse transport problem. An inverse transport problem may have multiple solutions and the most widely used approach for its solution is an iterative optimization method. This paper proposes a stochastic derivative-free global optimization algorithm to find multiple solutions of inverse transport problems. The algorithm is an extension of a multilevel single linkage (MLSL) method where a mesh adaptive direct search (MADS) algorithm is incorporated into the local phase. Furthermore, numerical test cases using uncollided fluxes of discrete gamma-ray lines are presented to show the performance of this new algorithm.

42 ENGINEERING↗

Trajectory Modulated Arc Therapy: A Fully Dynamic Delivery With Synchronized Couch and Gantry Motion Significantly Improves Dosimetric Indices Correlated With Poor Cosmesis in Accelerated Partial Breast Irradiation

Purpose: To develop planning and delivery capabilities for linear accelerator–based nonisocentric trajectory modulated arc therapy (TMAT) and to evaluate the benefit of TMAT for accelerated partial breast irradiation (APBI) with the patient in prone position. Methods and Materials: An optimization algorithm for volumetrically modulated arc therapy (VMAT) was generalized to allow for user-defined nonisocentric TMAT trajectories combining couch rotations and translations. After optimization, XML scripts were automatically generated to program and subsequently deliver the TMAT plans. For 10 breast patients in the prone position, TMAT and 6-field noncoplanar intensity modulated radiation therapy (IMRT) plans were generated under equivalent objectives and constraints. These plans were compared with regard to whole breast tissue volume receiving more than 100%, 80%, 50%, and 20% of the prescription dose. Results: For TMAT APBI, nonisocentric collision-free horizontal arcs with large angular span (251.5 ± 7.9°) were optimized and delivered with delivery time of ∼4.5 minutes. Percentage changes of whole breast tissue volume receiving more than 100%, 80%, 50%, and 20% of the prescription dose for TMAT relative to IMRT were −10.81% ± 6.91%, −27.81% ± 7.39%, −14.82% ± 9.67%, and 39.40% ± 10.53% (P≤.01). Conclusions: This is a first demonstration of end-to-end planning and delivery implementation of a fully dynamic APBI TMAT. Compared with IMRT, TMAT resulted in marked reduction of the breast tissue volume irradiated at high doses.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Low-depth Clifford circuits approximately solve MaxCut

We introduce a quantum-inspired approximation algorithm for MaxCut based on low-depth Clifford circuits. We start by showing that the solution unitaries found by the adaptive quantum approximation optimization algorithm (ADAPT-QAOA) for the MaxCut problem on weighted fully connected graphs are (almost) Clifford circuits. Motivated by this observation, we devise an approximation algorithm for MaxCut, ADAPT-Clifford, that searches through the Clifford manifold by combining a minimal set of generating elements of the Clifford group. Our algorithm finds an approximate solution of MaxCut on an N -vertex graph by building a depth O ( N ) Clifford circuit. The algorithm has runtime complexity O ( N 2 ) and O ( N 3 ) for sparse and dense graphs, respectively, and space complexity O ( N 2 ) , with improved solution quality achieved at the expense of more demanding runtimes. We implement ADAPT-Clifford and characterize its performance on graphs with positive and signed weights. The case of signed weights is illustrated with the paradigmatic Sherrington-Kirkpatrick model, for which our algorithm finds solutions with ground-state mean energy density corresponding to ∼ 94 % of the Parisi value in the thermodynamic limit. The case of positive weights is investigated by comparing the cut found by ADAPT-Clifford with the cut found with the Goemans-Williamson (GW) algorithm. For both sparse and dense instances we provide copious evidence that, up to hundreds of nodes, ADAPT-Clifford finds cuts of lower energy than GW. Published by the American Physical Society 2024

Muñoz-Arias, Manuel H. (ORCID:000000025711029X)↗