Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “optimization methods”

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,639 records · Page 91

Synthesis of Nickel-Doped Ceria Catalysts for Selective Acetylene Hydrogenation

Metallic nickel is known to be an active, but not a selective hydrogenation catalyst for conversion of alkynes to alkenes. On the other hand, nickel oxide is not active. Recently, we have demonstrated that nickel doped into ceria provides an inexpensive catalyst for selective hydrogenation of acetylene in the presence of ethylene. Here, we evaluate various synthesis methods to achieve optimal selective hydrogenation performance. We examined incipient wetness impregnation, coprecipitation, solution combustion, and sol-gel synthesis to study how the method of preparation affects catalytic structure and behavior. Sol-gel synthesis, coprecipitation, and solution combustion synthesis methods favor nickel incorporation into the ceria lattice, while incipient wetness impregnation creates segregated nickel species on the ceria surface. For hydrogenation of acetylene, these nickel surface species lead to poor ethylene selectivity due to ethane and oligomer formation. However, when nickel is incorporated into the ceria lattice, ethane formation is prevented even while achieving 100 % conversion of acetylene. Coke formation is also significantly reduced on these catalysts compared to conventional nanoparticle counterparts. Finally, we conclude that sol-gel synthesis provides the optimal method for creating a uniform dopant distribution within the high surface area ceria.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Encoding trade-offs and design toolkits in quantum algorithms for discrete optimization: coloring, routing, scheduling, and other problems

Challenging combinatorial optimization problems are ubiquitous in science and engineering. Several quantum methods for optimization have recently been developed, in different settings including both exact and approximate solvers. Addressing this field of research, this manuscript has three distinct purposes. First, we present an intuitive method for synthesizing and analyzing discrete (i.e., integer-based) optimization problems, wherein the problem and corresponding algorithmic primitives are expressed using a discrete quantum intermediate representation (DQIR) that is encoding-independent. This compact representation often allows for more efficient problem compilation, automated analyses of different encoding choices, easier interpretability, more complex runtime procedures, and richer programmability, as compared to previous approaches, which we demonstrate with a number of examples. Second, we perform numerical studies comparing several qubit encodings; the results exhibit a number of preliminary trends that help guide the choice of encoding for a particular set of hardware and a particular problem and algorithm. Our study includes problems related to graph coloring, the traveling salesperson problem, factory/machine scheduling, financial portfolio rebalancing, and integer linear programming. Third, we design low-depth graph-derived partial mixers (GDPMs) up to 16-level quantum variables, demonstrating that compact (binary) encodings are more amenable to QAOA than previously understood. We expect this toolkit of programming abstractions and low-level building blocks to aid in designing quantum algorithms for discrete combinatorial problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

SU-F-J-09: Radioactive Seed Localization for Breast Lumpectomy - Towards Optimization

Purpose: The aim of this work is to propose a method to optimize radioactive source localization (RSL) for non-palpable breast cancer surgery. RSL is commonly used as a guiding technique during surgery for excision of non-palpable tumors. A collimated hand-held detector is used to localize radioactive sources implanted in tumors. Incisions made by the surgeon are based on maximum observed detector counts, and tumors are subsequently resected based on an arbitrary estimate of the counts expected at the surgical margin boundary. This work focuses on building a framework to predict detector counts expected throughout the procedure to improve surgical margins. Methods: A gamma detection system called the Neoprobe GDS was used for this work. The probe consists of a cesium zinc telluride crystal and a collimator. For this work, an I-125 Best Medical model 2301 source was used. The source was placed in three different phantoms, a PMMA, a Breast (25%- glandular tissue/75%- adipose tissue) and a Breast (75-25) phantom with a backscatter thickness of 6 cm. Counts detected by the probe were recorded with varying amounts of phantom thicknesses placed on top of the source. A calibration curve was generated using MATLAB based on the counts recorded for the calibration dataset acquired with the PMMA phantom. Results: The observed detector counts data used as the validation set was accurately predicted to within ±3.2%, ±6.9%, ±8.4% for the PMMA, Breast (75-25), Breast (25–75) phantom respectively. The average difference between predicted and observed counts was −0.4%, 2.4%, 1.4% with a standard deviation of 1.2 %, 1.8%, 3.4% for the PMMA, Breast (75-25), Breast (25–75) phantom respectively. Conclusion: The results of this work provide a basis for characterization of a detector used for RSL. Counts were predicted to within ±9% for three different phantoms without the application of a density correction factor.

60 APPLIED LIFE SCIENCES↗

New technique for optimal smoothing of data

Recursive method for the optimal smoothing of data has numerical superiority and is more easily understood in terms of physical reasoning than earlier methods. Using a Kalman filter, the smoothing technique, applied to a nonlinear parameter identification problem, is useful in those situations where linearization about a reference solution is valid.

Fraser, D. C.↗

Perturbing engine performance measurements to determine optimal engine control settings

Methods and systems for optimizing a performance of a vehicle engine are provided. The method includes determining an initial value for a first engine control parameter based on one or more detected operating conditions of the vehicle engine, determining a value of an engine performance variable, and artificially perturbing the determined value of the engine performance variable. The initial value for the first engine control parameter is then adjusted based on the perturbed engine performance variable causing the engine performance variable to approach a target engine performance variable. Operation of the vehicle engine is controlled based on the adjusted initial value for the first engine control parameter. These acts are repeated until the engine performance variable approaches the target engine performance variable.

33 ADVANCED PROPULSION SYSTEMS↗

DiaMonD: An Integrated Multifaceted Approach to Mathematics at the Interfaces of Data, Models, and Decisions

For many problems in computational science and engineering, the so-called forward problem-solution of the underlying mathematical model to yield output quantities of interest, given input parameters--is difficult enough for frontier complex models, which are often characterized by coupled multiphysics and possibly stochastic behavior over a wide range of length and time scales. However, we face the need to move beyond the forward problem, to address what is often the ultimate goal of CS&E: decision-making. This requires us to tackle a spectrum of mathematical problems that subsume, and are thus even more difficult than, the forward problem. First, given experimental data, we wish to estimate unknown parameters characterizing a model of a natural or engineered system by solving an inverse problem. Second, we seek the optimal con figuration of the system (or experiments) by solving an optimal design problem. Third, the optimal operation of the system must be determined by solving an optimal control problem. And fourth, we must quantify uncertainties as they propagate through all of the preceding problems, from data to model inference to prediction and finally to optimal design and control. Optimal mathematical methods and algorithms for such end-to-end, data-to-decisions modeling and simulation of complex problems--requiring integrated solution of forward, inverse, optimization, and UQ problems for large, multiphysics, multiscale models--entails challenges of the highest order, while also presenting great opportunities for applied mathematics research. Unfortunately, research on forward, inverse, optimization, and UQ problems has all-too-often progressed in isolation. This has led to mathematical methods that perform well when considering the forward simulation in isolation, but are prohibitive or suboptimal or unstable when combined with other methods within the framework of inversion, optimization, or UQ. The reverse is also true: general-purpose methods developed within the inversion, optimization, and UQ fields often become prohibitive when applied to complex CS&E problems since they do not exploit their structure. DiaMonD is a multi-institutional DOE MMICC effort that aims to address the challenges of end-to-end, data-to-decisions modeling and simulation for complex CS&E problems in a unifi ed, systematic, and integrated fashion. Institutions involved are Colorado State, Florida State, Los Alamos, MIT, Oak Ridge, Stanford, and UT-Austin. The goals of DiaMonD are (1) to develop advanced mathematical methods and analysis for multimodel, multiphysics, multiscale model problems driven by frontier DOE applications, including those in subsurface energy and environmental flows, materials for energy storage and conversion, and ice sheet dynamics; (2) to create theory and algorithms for integrated inversion, optimization, and uncertainty quantification for these complex problems; and (3) to disseminate the philosophy of an integrated end-to-end, data-to-decisions approach to modeling and simulation of complex problems to the broader applied math and computational science communities through workshops and other outreach.

97 MATHEMATICS AND COMPUTING↗

Gamma-ray timing of redback PSR J2339–0533: Hints for gravitational quadrupole moment changes

We present the results of precision gamma-ray timing measurements of the binary millisecond pulsar PSR J2339–0533, an irradiating system of the “redback” type, using data from the Fermi Large Area Telescope. We describe an optimized analysis method to determine a long-term phase-coherent timing solution spanning more than six years, including a measured eccentricity of the binary orbit and constraints on the proper motion of the system. A major result of this timing analysis is the discovery of an extreme variation of the nominal 4.6 hr orbital period P{sub orb} over time, showing alternating epochs of decrease and increase. We inferred a cyclic modulation of P{sub orb} with an approximate cycle duration of 4.2 yr and a modulation amplitude of ΔP{sub orb}/P{sub orb}=2.3×10{sup −7}. Considering different possible physical causes, the observed orbital-period modulation most likely results from a variable gravitational quadrupole moment of the companion star due to cyclic magnetic activity in its convective zone.

79 ASTRONOMY AND ASTROPHYSICS↗

DEVELOPMENT AND APPLICATION OF RISK ANALYSIS TOOLKIT FOR PLANT RESOURCE OPTIMIZATION

This paper presents the development of methods and tools that are being designed to optimize plant operations (e.g., maintenance/replacement schedules and optimal maintenance postures for plant components) in a manner that is more cost effective than current approaches and makes better use of available component health and cost data. These methods include both data- and model-based optimization methods. Model-based optimization methods directly include reliability and cost models to determine an optimal plant operational strategy. We consider gradient-based and evolutionary (based on genetic algorithms) optimization methods. The second class of methods target more specific use cases (e.g., project schedule optimization) and are not based on reliability models directly, but they require specific component reliability and cost data. This class of methods is based on variants of the knapsack problem with an aim to determine an optimal project schedule that maximizes the overall NPV. This paper also presents multi-objective methods designed to identify an optimal maintenance posture based on a Pareto frontier analysis. Rather than dictating the “right” tradeoff (i.e., identify the absolute best posture), we show how it is possible to perform a trade space exploration approach (i.e., identify value and costs of several postures and let the analysis account for desired value and cost metrics). This is performed by identifying maintenance postures that maximize value (e.g., system availability) and minimize operational costs, i.e., the Pareto frontier in a value-cost trade space. For all these methods we present detailed applicative examples that show their validity from a decision-making perspective.

97 - MATHEMATICS AND COMPUTING↗

The application of quadratic optimal cooperative control synthesis to a CH-47 helicopter

A control-system design method, Quadratic Optimal Cooperative Control Synthesis (CCS), is applied to the design of a Stability and Control Augmentation Systems (SCAS). The CCS design method is different from other design methods in that it does not require detailed a priori design criteria, but instead relies on an explicit optimal pilot-model to create desired performance. The design model, which was developed previously for fixed-wing aircraft, is simplified and modified for application to a Boeing Vertol CH-47 helicopter. Two SCAS designs are developed using the CCS design methodology. The resulting CCS designs are then compared with designs obtained using classical/frequency-domain methods and Linear Quadratic Regulator (LQR) theory in a piloted fixed-base simulation. Results indicate that the CCS method, with slight modifications, can be used to produce controller designs which compare favorably with the frequency-domain approach.

Townsend, Barbara K.↗

Optimization of human factors methods within the design cycle - 178

Amec Foster Wheeler developed and implemented the Human Factors Engineering program to analyze the design of a large-scale electrical component readying for installation in a nuclear facility. The HFE program applied the program elements from NUREG-0711 throughout the detailed design phase of the project. This paper provides: (1) a review of the Human Factors Engineering (HFE) program methodology that was implemented, and (2) a summary of the integrated approach to the design, verification, and validation of the Human-Machine Interface (HMI) display and mimic panel. In general, it was found that a collaborative approach between the design team, equipment vendor, HFE, and the end-users of the component created a more streamlined approach to the design and allowed for user feedback to be implemented throughout the progression of the design. However, the trade-off to this approach is that an integrated validation is likely required once all design aspects are complete, to ensure the modification is acceptable in the context of the overall system. (authors)

42 ENGINEERING↗

DiaMonD: An Integrated Multifaceted Approach to Mathematics at the Interfaces of Data, Models, and Decisions

For many problems in computational science and engineering, the so-called forward problem—solution of the underlying mathematical model to yield output quantities of interest, given input parameters—is difficult enough for frontier complex models, which are often characterized by coupled multiphysics and possibly stochastic behavior over a wide range of length and time scales. However, we face the need to move beyond the forward problem, to address what is often the ultimate goal of CS&E: decision making. This requires us to tackle a spectrum of mathematical problems that subsume, and are thus even more difficult than, the forward problem. First, given experimental data, we wish to estimate unknown parameters characterizing a model of a natural or engineered system by solving an inverse problem. Second, we seek the optimal configuration of the system (or experiments) by solving an optimal design problem. Third, the optimal operation of the system must be determined by solving an optimal control problem. And fourth, we must quantify uncertainties as they propagate through all of the preceding problems, from data to model inference to prediction and finally to optimal design and control. Optimal mathematical methods and algorithms for such end-to-end, data-to-decisions modeling and simulation of complex problems—requiring integrated solution of forward, inverse, optimization, and UQ problems for large, multiphysics, multiscale models—entails challenges of the highest order, while also presenting great opportunities for applied mathematics research. Unfortunately, research on forward, inverse, optimization, and UQ problems has all-too-often progressed in isolation. This has led to mathematical methods that perform well when considering the forward simulation in isolation, but are prohibitive or suboptimal or unstable when combined with other methods within the framework of inversion, optimization, or UQ. The reverse is also true: general-purpose methods developed within the inversion, optimization, and UQ fields often become prohibitive when applied to complex CS&E problems since they do not exploit their structure.

58 GEOSCIENCES↗

Temperature-Dependent Evolution of the Oxidation States of Cobalt and Platinum in Co 1– x Pt x Clusters under H 2 and CO + H 2 Atmospheres

Co1-xPtx clusters of 2.9-nm size with a range of atomically precise Pt/Co atomic ratios (x = 0, 0.25, 0.5, 0.75, 1) were synthesized using the mass-selected low-energy cluster beam deposition (LECBD) technique and soft-landed onto an amorphous alumina thin film prepared by atomic layer deposition (ALD). Utilizing ex situ X-ray photoemission spectroscopy (XPS), the oxidation state of the as-made clusters supported on Al2O3 was determined after both a 1-h-long exposure to air and aging for several weeks while exposed to air. Next, the aged duster samples were characterized by grazing-incidence X-ray absorption spectroscopy (GIXAS) and then pretreated with diluted hydrogen and further exposed to the mixture of diluted CO and H-2 up to 225 degrees C at atmospheric pressure, and the temperature-dependent evolutions of the particle size/shape and the oxidation states of the individual metal components within the dusters were monitored using in situ grazing-incidence small-angle X-ray scattering and X-ray absorption spectroscopy (GISAXS/GIXAS). The changes in the oxidation states of Co and Pt exhibited a nonlinear dependence on the Pt/Co atomic ratio of the dusters. For example, a low Pt/Co ratio (x <= 0.5) facilitates the formation of Co(OH)(2), whereas a high Pt/Co ratio (x = 0.75) stabilizes the Co3O4 composition instead through the formation of a Co-Pt core-shell structure where the platinum shell inhibits the reduction of cobalt in the core of the Co1-xPtx alloy dusters. The obtained results indicate methods for optimizing the composition and structure of binary alloy clusters for catalysis.

Yang, Bing↗

Shaping Diffraction-Grating Grooves to Optimize Efficiency

A method of shaping diffraction-grating grooves to optimize the spectral efficiency, spectral range, and image quality of a spectral imaging instrument is under development. The method is based on the use of an advanced design algorithm to determine the possibly complex shape of grooves needed to obtain a desired efficiency-versus-wavelength response (see figure). Then electron- beam fabrication techniques are used to realize the required groove shape. The method could be used, for example, to make the spectral efficiency of the grating in a given wavelength range proportional to the inverse of the spectral efficiency of a photodetector array so that the overall spectral efficiency of the combination of the grating and the photodetector array would be flat. The method has thus far been applied to one-dimensional gratings only, but in principle, it is also applicable to two-dimensional gratings. The algorithm involves calculations in the spatial-frequency domain. The spatial-frequency spectrum of a grating is represented as a diffraction-order spectral-peak-width function multiplied by an efficiency function for a single grating groove. This representation affords computational efficiency and accuracy by making it possible to consider only the response from one grating groove (one period of the grating), instead of from the whole grating area, in determining the response from the entire grating. This combination of efficiency and accuracy is crucial for future extensions of the algorithm to two-dimensional designs and to designs in which polarization must also be taken into account. The algorithm begins with the definition of target values of relative efficiency that represent the desired spectral response of the grating in certain spectral frequencies calculated from the diffraction order and wavelength. The grating period is divided into a number of cells - typically, 100. The phase contribution from each cell is determined from the phase of the incident electromagnetic wave and the height of the grating surface in the cell. The total contribution from all cells to each target value is then calculated. Then a method known to specialists as the optimum-rotation-angle method is used to adjust the height of each cell so that the total response from all cells is optimized. The computation is iterative and continues until the desired response is obtained. In the event that the desired response is unphysical, the algorithm nevertheless strives to generate a grating-grove profile for which the response approximates the desired one as closely as possible.

Backlund, John↗

Optimization with artificial neural network systems - A mapping principle and a comparison to gradient based methods

General formulae for mapping optimization problems into systems of ordinary differential equations associated with artificial neural networks are presented. A comparison is made to optimization using gradient-search methods. The performance measure is the settling time from an initial state to a target state. A simple analytical example illustrates a situation where dynamical systems representing artificial neural network methods would settle faster than those representing gradient-search. Settling time was investigated for a more complicated optimization problem using computer simulations. The problem was a simplified version of a problem in medical imaging: determining loci of cerebral activity from electromagnetic measurements at the scalp. The simulations showed that gradient based systems typically settled 50 to 100 times faster than systems based on current neural network optimization methods.

Leong, Harrison Monfook↗

Delivery validation of an automated modulated electron radiotherapy plan

Purpose: Modulated electron radiation therapy (MERT) represents an active area of interest that offers the potential to improve healthy tissue sparing in treatment of certain cancer cases. Challenges remain however in accurate beamlet dose calculation, plan optimization, collimation method, and delivery accuracy. In this work, the authors investigate the accuracy and efficiency of an end-to-end MERT plan and automated delivery method. Methods: Treatment planning was initiated on a previously treated whole breast irradiation case including an electron boost. All dose calculations were performed using Monte Carlo methods and beam weights were determined using a research-based treatment planning system capable of inverse optimization. The plan was delivered to radiochromic film placed in a water equivalent phantom for verification, using an automated motorized tertiary collimator. Results: The automated delivery, which covered four electron energies, 196 subfields, and 6183 total MU was completed in 25.8 min, including 6.2 min of beam-on time. The remainder of the delivery time was spent on collimator leaf motion and the automated interfacing with the accelerator in service mode. Comparison of the planned and delivered film dose gave 3%/3mm gamma pass rates of 62.1%, 99.8%, 97.8%, 98.3%, and 98.7% for the 9, 12, 16, and 20 MeV, and combined energy deliveries, respectively. Delivery was also performed with a MapCHECK device and resulted in 3%/3 mm gamma pass rates of 88.8%, 86.1%, 89.4%, and 94.8% for the 9, 12, 16, and 20 MeV energies, respectively. Conclusions: Results of the authors’ study showed that an accurate delivery utilizing an add-on tertiary electron collimator is possible using Monte Carlo calculated plans and inverse optimization, which brings MERT closer to becoming a viable option for physicians in treating superficial malignancies.

60 APPLIED LIFE SCIENCES↗

Four-body trajectory optimization

The two methods which are suitable for use in a 4-body trajectory optimization program are both multiconic methods. They include an approach due to Wilson (1970) and to Byrnes and Hooper (1970) and a procedure developed by Stumpff and Weiss (1968). The various steps in a trajectory optimization program are discussed, giving attention to variable step integration, the correction of errors by quadrature formulas, questions of two-impulse transfer, three-impulse transfer, and two examples which illustrate the implementation of the computational approaches.

Pu, C. L.↗