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,783 records · Page 99
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Quantum Orbital Minimization Method for Excited States Calculation on a Quantum Computer

Herein we propose a quantum-classical hybrid variational algorithm, the quantum orbital minimization method (qOMM), for obtaining the ground state and low-lying excited states of a Hermitian operator. Given parametrized ansatz circuits representing eigenstates, qOMM implements quantum circuits to represent the objective function in the orbital minimization method and adopts a classical optimizer to minimize the objective function with respect to the parameters in ansatz circuits. The objective function has an orthogonality constraint implicitly embedded, which allows qOMM to apply a different ansatz circuit to each input reference state. We carry out numerical simulations that seek to find excited states of H 2 , LiH, and a toy model consisting of four hydrogen atoms arranged in a square lattice in the STO-3G basis with UCCSD ansatz circuits. Comparing the numerical results with existing excited states methods, qOMM is less prone to getting stuck in local minima and can achieve convergence with more shallow ansatz circuits.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Manufacturing of Dysprosium-Iron Alloys by Electrolysis in Fluoride-Based Electrolytes: Oxide Solubility Determinations

Electrolytic production of light rare earth elements and alloys takes place in a fluoride-based electrolyte using rare earth oxides as raw material. The optimization of this method, mainly in terms of the energy efficiency and environmental impact control, is rather challenging. Anode effects, evolution of fluorine-containing compounds, and side cathode reactions could largely be minimized by a good control of the amount of rare earth oxide species dissolved in the fluoride-based electrolyte and their dissolution rate. The oxide content of the fluoride melts REF{sub 3}-LiF (RE = Nd, Dy) at different compositions and temperatures were experimentally determined by carbothermal analysis of melt samples. The highest solubility values of oxide species, added as Dy{sub 2}O{sub 3} and Dy{sub 2}(CO{sub 3}){sub 3}, were obtained to be of ca. 3 wt pct (expressed as Dy{sub 2}O{sub 3}) in the case of the equimolar DyF{sub 3}-LiF melt at 1323 K (1050 °C). The oxide saturation values increased with the amount of REF{sub 3} present in the molten bath and the working temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Energetics and mechanisms for the acetonyl radical + O 2 reaction: An important system for atmospheric and combustion chemistry

The acetonyl radical (•CH 2 COCH 3 ) is relevant to atmospheric and combustion chemistry due to its prevalence in many important reaction mechanisms. One such reaction mechanism is the decomposition of Criegee intermediates in the atmosphere that can produce acetonyl radical and OH. In order to understand the fate of the acetonyl radical in these environments and to create more accurate kinetics models, we have examined the reaction system of the acetonyl radical with O 2 using highly reliable theoretical methods. Structures were optimized using coupled cluster theory with singles, doubles, and perturbative triples [CCSD(T)] with an atomic natural orbital (ANO0) basis set. Energetics were computed to chemical accuracy using the focal point approach involving perturbative treatment of quadruple excitations [CCSDT(Q)] and basis sets as large as cc-pV5Z. The addition of O 2 to the acetonyl radical produces the acetonylperoxy radical, and multireference computations on this reaction suggest it to be barrierless. Additionally, no submerged pathways were found for the unimolecular isomerization of the acetonylperoxy radical. Besides dissociation to reactants, the lowest energy pathway available for the acetonylperoxy radical is a 1-5 H shift from the methyl group to the peroxy group through a transition state that is 3.3 kcal mol -1 higher in energy than acetonyl radical + O 2 . The ultimate products from this pathway are the enol tautomer of the acetonyl radical along with O 2 . Multiple pathways that lead to OH formation are considered; however, all of these pathways are predicted to be energetically inaccessible, except at high temperatures.

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

Optimizing Service Restoration in Distribution Systems With Uncertain Repair Time and Demand

This paper proposes a novel method to co-optimize the distribution system operation and repair crew routing for outage restoration after extreme weather events. A two-stage stochastic mixed integer linear program is developed. The first stage is to dispatch the repair crews to the damaged components. The second stage is distribution system restoration using distributed generators, and reconfiguration. We consider demand uncertainty in terms of a truncated normal forecast error distribution, and model the uncertainty of the repair time using a lognormal distribution. A new decomposition approach, combined with the progressive hedging algorithm, is developed for solving large-scale outage management problems in an effective and timely manner. Furthermore, the proposed method is validated on modified IEEE 34- and 8500-bus distribution test systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2017

Viral phylogenomics using an alignment-free method: A three-step approach to determine optimal length of k-mer

The development of rapid, economical genome sequencing has shed new light on the classification of viruses. As of October 2016, the National Center for Biotechnology Information (NCBI) database contained >2 million viral genome sequences and a reference set of ~4000 viral genome sequences that cover a wide range of known viral families. Whole-genome sequences can be used to improve viral classification and provide insight into the viral tree of life . However, due to the lack of evolutionary conservation amongst diverse viruses, it is not feasible to build a viral tree of life using traditional phylogenetic methods based on conserved proteins. In this study, we used an alignment-free method that uses k-mers as genomic features for a large-scale comparison of complete viral genomes available in RefSeq. To determine the optimal feature length, k (an essential step in constructing a meaningful dendrogram), we designed a comprehensive strategy that combines three approaches: (1) cumulative relative entropy, (2) average number of common features among genomes, and (3) the Shannon diversity index. This strategy was used to determine k for all 3,905 complete viral genomes in RefSeq. Lastly, the resulting dendrogram shows consistency with the viral taxonomy of the ICTV and the Baltimore classification of viruses.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

On the ground state calculation of a many-body system using a self-consistent basis and quasi-Monte Carlo: An application to water hexamer

Given a quantum many-body system, the Self-Consistent Phonons (SCP) method provides an optimal harmonic approximation by minimizing the free energy. In particular, the SCP estimate for the vibrational ground state (zero temperature) appears to be surprisingly accurate. We explore the possibility of going beyond the SCP approximation by considering the system Hamiltonian evaluated in the harmonic eigenbasis of the SCP Hamiltonian. It appears that the SCP ground state is already uncoupled to all singly- and doubly-excited basis functions. So, in order to improve the SCP result at least triply-excited states must be included, which then reduces the error in the ground state estimate substantially. For a multidimensional system two numerical challenges arise, namely, evaluation of the potential energy matrix elements in the harmonic basis, and handling and diagonalizing the resulting Hamiltonian matrix, whose size grows rapidly with the dimensionality of the system. Using the example of water hexamer we demonstrate that such calculation is feasible, i.e., constructing and diagonalizing the Hamiltonian matrix in a triply-excited SCP basis, without any additional assumptions or approximations. Our results indicate particularly that the ground state energy differences between different isomers (e.g., cage and prism) of water hexamer are already quite accurate within the SCP approximation.

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

Factorizations of hyperpower family of iterative methods via least squares approach

A new method for finding proper factorizations of the hyperpower family of iterative methods for computing generalized inverses is proposed. The method is based on numerical optimization of an appropriate least squares problem. An iterative method in a non-normalized form for finding the generalized outer inverses with prescribed range and null space of any given complex matrix is proposed and considered. Computational complexity and convergence of the method are studied. Several numerical experiments are also presented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Impacts of annealing treatment on the microstructure of U-Mo monolithic fuel plates

The objective of this paper is to analyze the effect of different heat treatments on the fuel microstructure and texture in laboratory scale monolithic U-Mo fuel plates. Such analyses are relevant as they could possibly inform future fabrication methods for the optimization of microstructure and of material properties to ultimately improve fuel performance during irradiation. For this reason, detail characterization including Energy x-ray dispersive spectroscopy (EDS) and electron backscattered diffraction (EBSD) techniques were applied on U-Mo plates fabricated with different heat treatments to understand the influence of fabrication on fuel microstructure. Here, U-Mo fresh fuel specimens with or without annealing treatments (homogenization and stress release annealing) after rolling were analyzed in this work. In this study for the first time EBSD was systematically used to analyze the heat treatments influence on texture present in fresh fuel monolithic U-Mo fuel plates. Such changes were connected to other known microstructural and chemical changes. In this work minimization of molybdenum concentration variation and gamma phase decomposition of the U-Mo fuel core after homogenization was observed. Also, an increase in the interaction layer of the U-Mo fuel core with the Zr interlayer diffusion barrier was observed with the UZr 2 growing up to 1 µm in thickness. Such microstructural changes were aligned to the changes in texture and grain structure. Indeed, common texture observed in body-centered cubic (BCC) metal after cold rolling (e.g., α, ξ, γ fibers) were minimized by the annealing/homogenization process. Moreover, grain structure was influenced by the fabrication route with elongated grains formed in the rolled sample and equiaxed grains found after annealing. Such features are relevant since the formation of microstructural features formed during fuel fabrication can impact fuel performance in reactor.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2015

Noncoplanar Beam Angle Class Solutions to Replace Time-Consuming Patient-Specific Beam Angle Optimization in Robotic Prostate Stereotactic Body Radiation Therapy

Purpose: To investigate development of a recipe for the creation of a beam angle class solution (CS) for noncoplanar prostate stereotactic body radiation therapy to replace time-consuming individualized beam angle selection (iBAS) without significant loss in plan quality, using the in-house “Erasmus-iCycle” optimizer for fully automated beam profile optimization and iBAS. Methods and Materials: For 30 patients, Erasmus-iCycle was first used to generate 15-, 20-, and 25-beam iBAS plans for a CyberKnife equipped with a multileaf collimator. With these plans, 6 recipes for creation of beam angle CSs were investigated. Plans of 10 patients were used to create CSs based on the recipes, and the other 20 to independently test them. For these tests, Erasmus-iCycle was also used to generate intensity modulated radiation therapy plans for the fixed CS beam setups. Results: Of the tested recipes for CS creation, only 1 resulted in 15-, 20-, and 25-beam noncoplanar CSs without plan deterioration compared with iBAS. For the patient group, mean differences in rectum D{sub 1cc}, V{sub 60GyEq}, V{sub 40GyEq}, and D{sub mean} between 25-beam CS plans and 25-beam plans generated with iBAS were 0.2 ± 0.4 Gy, 0.1% ± 0.2%, 0.2% ± 0.3%, and 0.1 ± 0.2 Gy, respectively. Differences between 15- and 20-beam CS and iBAS plans were also negligible. Plan quality for CS plans relative to iBAS plans was also preserved when narrower planning target volume margins were arranged and when planning target volume dose inhomogeneity was decreased. Using a CS instead of iBAS reduced the computation time by a factor of 14 to 25, mainly depending on beam number, without loss in plan quality. Conclusions: A recipe for creation of robust beam angle CSs for robotic prostate stereotactic body radiation therapy has been developed. Compared with iBAS, computation times decreased by a factor 14 to 25. The use of a CS may avoid long planning times without losses in plan quality.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2023

A Model Calibration Method for Grid-Forming Inverters Using Iterative Bayesian Optimization

As inverter-based resources (IBRs) are rapidly deployed, especially at the distribution and microgrid levels, the need to include their accurate representations in power systems models increases. With well-calibrated models, IBR-interactions can be examined, preventing any stability and operational issues, especially in islanded systems. This paper proposes a method to calibrate generic inverter models using Bayesian optimization, but with a parameter-grouping approach to improve speed and accuracy. This approach is illustrated for a grid-forming inverter using synthetic data and field measurements from a simple IBRbased microgrid. Tests show that the calibration algorithm has modest computation cost and is robust to noisy measurements..

Biswas, Shuchismita↗
DOE OSTI2015

Recent progress of quantum annealing

We review the recent progress of quantum annealing. Quantum annealing was proposed as a method to solve generic optimization problems. Recently a Canadian company has drawn a great deal of attention, as it has commercialized a quantum computer based on quantum annealing. Although the performance of quantum annealing is not sufficiently understood, it is likely that quantum annealing will be a practical method both on a conventional computer and on a quantum computer.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Automated characterization of perceptual quality of clinical chest radiographs: Validation and calibration to observer preference

Purpose: The authors previously proposed an image-based technique [Y. Lin et al. Med. Phys. 39, 7019–7031 (2012)] to assess the perceptual quality of clinical chest radiographs. In this study, an observer study was designed and conducted to validate the output of the program against rankings by expert radiologists and to establish the ranges of the output values that reflect the acceptable image appearance so the program output can be used for image quality optimization and tracking. Methods: Using an IRB-approved protocol, 2500 clinical chest radiographs (PA/AP) were collected from our clinical operation. The images were processed through our perceptual quality assessment program to measure their appearance in terms of ten metrics of perceptual image quality: lung gray level, lung detail, lung noise, rib–lung contrast, rib sharpness, mediastinum detail, mediastinum noise, mediastinum alignment, subdiaphragm–lung contrast, and subdiaphragm area. From the results, for each targeted appearance attribute/metric, 18 images were selected such that the images presented a relatively constant appearance with respect to all metrics except the targeted one. The images were then incorporated into a graphical user interface, which displayed them into three panels of six in a random order. Using a DICOM calibrated diagnostic display workstation and under low ambient lighting conditions, each of five participating attending chest radiologists was tasked to spatially order the images based only on the targeted appearance attribute regardless of the other qualities. Once ordered, the observer also indicated the range of image appearances that he/she considered clinically acceptable. The observer data were analyzed in terms of the correlations between the observer and algorithmic rankings and interobserver variability. An observer-averaged acceptable image appearance was also statistically derived for each quality attribute based on the collected individual acceptable ranges. Results: The observer study indicated that, for each image quality attribute, the averaged observer ranking strongly correlated with the algorithmic ranking (linear correlation coefficient R > 0.92), with highest correlation (R = 1) for lung gray level and the lowest (R = 0.92) for mediastinum noise. There was a strong concordance between the observers in terms of their rankings (i.e., Kendall’s tau agreement > 0.84). The observers also generally indicated similar tolerance and preference levels in terms of acceptable ranges, as 85% of the values were close to the overall tolerance or preference levels and the differences were smaller than 0.15. Conclusions: The observer study indicates that the previously proposed technique provides a robust reflection of the perceptual image quality in clinical images. The results established the range of algorithmic outputs for each metric that can be used to quantitatively assess and qualify the appearance quality of clinical chest radiographs.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Physics-guided logistic classification for tool life modeling and process parameter optimization in machining

This paper describes a physics-guided logistic classification method for tool life modeling and process parameter optimization in machining. Tool life is modeled using a classification method since the exact tool life cannot be measured in a typical production environment where tool wear can only be directly measured when the tool is replaced. Here, in this study, laboratory tool wear experiments are used to simulate tool wear data normally collected during part production. Two states are defined: tool not worn (class 0) and tool worn (class 1). The non-linear reduction in tool life with cutting speed is modeled by applying a logarithmic transformation to the inputs for the logistic classification model. A method for interpretability of the logistic model coefficients is provided by comparison with the empirical Taylor tool life model. The method is validated using tool wear experiments for milling. Results show that the physics-guided logistic classification method can predict tool life using limited datasets. A method for pre-process optimization of machining parameters using a probabilistic machining cost model is presented. The proposed method offers a robust and practical approach to tool life modeling and process parameter optimization in a production environment.

Machine learning↗
DOE OSTI2017

Optimized pulses for Raman excitation through the continuum: Verification using the multiconfigurational time-dependent Hartree-Fock method

Here, we have verified a mechanism for Raman excitation of atoms through continuum levels previously obtained by quantum optimal control using the multiconfigurational time-dependent Hartree-Fock (MCTDHF) method. For the optimal control, which requires running multiple propagations to determine the optimal pulse sequence, we used the computationally inexpensive time-dependent configuration interaction singles (TDCIS) method. TDCIS captures all of the necessary correlation of the desired processes but assumes that ionization pathways reached via double excitations are not present. MCTDHF includes these pathways and all multiparticle correlations in a set of time-dependent orbitals. The mechanism that was determined to be optimal in the Raman excitation of the Ne 1s 2 2s 2 2p 5 3p 1 valence state via the metastable 1s 2 2s 1 2p 6 3p 1 resonance state involves a sequential resonance-valence excitation. First, a long pump pulse excites the core-hole state, and then a shorter Stokes pulse transfers the population to the valence state. This process represents the first step in a multidimensional x-ray spectroscopy scheme that will provide a local probe of valence electronic correlations. Although at the optimal pulse intensities at the TDCIS level of theory the MCTDHF method predicts multiple ionization or excitation ionization of the atom, at slightly lower intensities (reduced by a factor of about 4) the TDCIS mechanism is shown to hold qualitatively. Quantitatively, the MCTDHF populations are reduced from the TDCIS calculations by a factor of 4.

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

Bounding the moment deficit rate on crustal faults using geodetic data: Methods

Here, the geodetically derived interseismic moment deficit rate (MDR) provides a first-order constraint on earthquake potential and can play an important role in seismic hazard assessment, but quantifying uncertainty in MDR is a challenging problem that has not been fully addressed. We establish criteria for reliable MDR estimators, evaluate existing methods for determining the probability density of MDR, and propose and evaluate new methods. Geodetic measurements moderately far from the fault provide tighter constraints on MDR than those nearby. Previously used methods can fail catastrophically under predictable circumstances. The bootstrap method works well with strong data constraints on MDR, but can be strongly biased when network geometry is poor. We propose two new methods: the Constrained Optimization Bounding Estimator (COBE) assumes uniform priors on slip rate (from geologic information) and MDR, and can be shown through synthetic tests to be a useful, albeit conservative estimator; the Constrained Optimization Bounding Linear Estimator (COBLE) is the corresponding linear estimator with Gaussian priors rather than point-wise bounds on slip rates. COBE matches COBLE with strong data constraints on MDR. We compare results from COBE and COBLE to previously published results for the interseismic MDR at Parkfield, on the San Andreas Fault, and find similar results; thus, the apparent discrepancy between MDR and the total moment release (seismic and afterslip) in the 2004 Parkfield earthquake remains.

58 GEOSCIENCES↗
DOE OSTI2017

Trade-Offs in Thin Film Solar Cells with Layered Chalcostibite Photovoltaic Absorbers

Discovery of novel semiconducting materials is needed for solar energy conversion and other optoelectronic applications. However, emerging low-dimensional solar absorbers often have unconventional crystal structures and unusual combinations of optical absorption and electrical transport properties, which considerably slows down the research and development progress. Here, the effect of stronger absorption and weaker carrier collection of 2D-like absorber materials are studied using a high-throughput combinatorial experimental approach, complemented by advanced characterization and computations. It is found that the photoexcited charge carrier collection in CuSbSe 2 solar cells is enhanced by drift in an electric field, addressing a different absorption/collection balance. The resulting drift solar cells efficiency is <5% due to inherent J SC / V OC trade-off, suggesting that improved carrier diffusion and better contacts are needed to further increase the CuSbSe 2 performance. Furthermore, this study also illustrates the advantages of high-throughput experimental methods for fast optimization of the optoelectronic devices based on emerging low-dimensional semiconductor materials.

14 SOLAR ENERGY↗
DOE OSTI2017

Grain Boundary Engineering for Achieving High Thermoelectric Performance in n-Type Skutterudites

Grain or phase boundaries play a critical role in the carrier and phonon transport in bulk thermoelectric materials. Previous investigations about controlling boundaries primarily focused on the reducing grain size or forming nanoinclusions. Herein, liquid phase compaction method is first used to fabricate the Yb–filled CoSb 3 with excess Sb content, which shows the typical feature of low–angle grain boundaries with dense dislocation arrays. Seebeck coefficients show a dramatic increase via energy filtering effect through dislocation arrays with little deterioration on the carrier mobility, which significantly enhances the power factor over a broad temperature range with a high room–temperature value around 47 μW cm –2 K –1 . Simultaneously, the lattice thermal conductivity could be further suppressed via scattering phonons via dense dislocation scattering. As a result, the highest average figure of merit ZT of ≈1.08 from 300 to 850 K could be realized, comparable to the best reported result of single or triple–filled Skutterudites. Furthermore, this work clearly points out that low–angle grain boundaries fabricated by liquid phase compaction method could concurrently optimize the electrical and thermal transport properties leading to an obvious enhancement of both power factor and ZT .

36 MATERIALS SCIENCE↗
DOE OSTI2020

Measuring the dynamic response of the thylakoid architecture in plant leaves by electron microscopy

The performance of the photosynthesis machinery in plants, including light harvesting, electron transport, and protein repair is controlled by structural changes of the thylakoid membrane system inside the chloroplasts. In particular, the structure of the stacked grana area of thylakoid membranes is highly dynamic, changing in response to different environmental cues such as light intensity. For example, the aqueous thylakoid lumen enclosed by thylakoid membranes in grana has been documented to swell in the presence of light. However, light-induced alteration of the stromal gap in the stacked grana (partition gap) and of the unstacked stroma lamellae have not been well characterized. Light-induced changes of the entire thylakoid membrane system including the lumen in both stacked and unstacked domains as well as the partition gap are presented here, and the functional implications are discussed. This structural analysis was made possible by development of a robust semi-automated image analysis method combined with optimized plant tissue fixation techniques for transmission electron microscopy generating quantitative structural results for the analysis of thylakoid ultrastructure.

59 BASIC BIOLOGICAL SCIENCES↗