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,747 records · Page 97
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Unleashed from constrained optimization: quantum computing for quantum chemistry employing generator coordinate inspired method

Hybrid quantum-classical approaches offer potential solutions to quantum chemistry problems, yet they often manifest as constrained optimization problems. Here, we explore the interconnection between constrained optimization and generalized eigenvalue problems through the Unitary Coupled Cluster (UCC) excitation generators. Inspired by the generator coordinate method, we employ these UCC excitation generators to construct non-orthogonal, overcomplete many-body bases, projecting the system Hamiltonian into an effective Hamiltonian, which bypasses issues such as barren plateaus that heuristic numerical minimizers often encountered in standard variational quantum eigensolver (VQE). Diverging from conventional quantum subspace expansion methods, we introduce an adaptive scheme that robustly constructs the many-body basis sets from a pool of the UCC excitation generators. This scheme supports the development of a hierarchical ADAPT quantum-classical strategy, enabling a balanced interplay between subspace expansion and ansatz optimization to address complex, strongly correlated quantum chemical systems cost-effectively, setting the stage for more advanced quantum simulations in chemistry.

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

Communication: Analytical optimal pulse shapes obtained with the aid of genetic algorithms: Controlling the photoisomerization yield of retinal

We recently proposed a Quantum Optimal Control (QOC) method constrained to build pulses from analytical pulse shapes [R. D. Guerrero et al., J. Chem. Phys. 143(12), 124108 (2015)]. This approach was applied to control the dissociation channel yields of the diatomic molecule KH, considering three potential energy curves and one degree of freedom. In this work, we utilized this methodology to study the strong field control of the cis-trans photoisomerization of 11-cis retinal. This more complex system was modeled with a Hamiltonian comprising two potential energy surfaces and two degrees of freedom. The resulting optimal pulse, made of 6 linearly chirped pulses, was capable of controlling the population of the trans isomer on the ground electronic surface for nearly 200 fs. The simplicity of the pulse generated with our QOC approach offers two clear advantages: a direct analysis of the sequence of events occurring during the driven dynamics, and its reproducibility in the laboratory with current laser technologies.

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

3D Printed Optical Quality Silica and Silica-Titania Glasses from Sol-Gel Feedstocks

We report a method for fabricating optical quality silica and silica-titania glasses by additive manufacturing, or 3D printing. Key to this success was the combination of sol-gel derived silica and silica-titania colloidal feedstocks, 3D direct ink writing (DIW) technology, and conventional glass thermal processing methods. Printable silica and silica-titania sol inks were prepared directly from molecular precursors by a simple one-pot method, which was optimized to yield viscous, shear-thinning colloidal suspensions with tuned rheology ideal for DIW. After printing, the parts were dried and sintered under optimized thermal conditions to ensure complete organic removal and uniform densification without crystallization. Characterizations of the 3Dprinted pure silica and silica-titania glasses show that they are equivalent to commercial optical fused silica (Corning® 7980) and silica-titania glasses (Corning ULE® 7972). More specifically they exhibit comparable chemical composition, SiO2 network structure, refractive index, dispersion, optical transmission, and coefficient of thermal expansion. 3D printed silica and silica-titania glasses also exhibited comparable polished surface roughness and meet refractive index homogeneity standards within range of commercial optical grade glasses. This method establishes 3D printing as a viable tool to create optical glasses with compositional and geometric configurations that are inaccessible by conventional optical fabrication methods.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Application of an automated machine learning-genetic algorithm (AutoML-GA) coupled with computational fluid dynamics simulations for rapid engine design optimization

In recent years, the use of machine learning-based surrogate models for computational fluid dynamics (CFD) simulations has emerged as a promising technique for reducing the computational cost associated with engine design optimization. However, such methods still suffer from drawbacks. One main disadvantage is that the default machine learning (ML) hyperparameters are often severely suboptimal for a given problem. This has often been addressed by manually trying out different hyperparameter settings, but this solution is ineffective in case of a high-dimensional hyperparameter space. Besides this problem, the amount of data needed for training is also not known a priori. In response to these issues that need to be addressed, the present work describes and validates an automated active learning approach, AutoML-GA, for surrogate-based optimization of internal combustion engines. In this approach, a Bayesian optimization technique is used to find the best machine learning hyperparameters based on an initial dataset obtained from a small number of CFD simulations. Subsequently, a genetic algorithm is employed to locate the design optimum on the ML surrogate surface. In the vicinity of the design optimum, the solution is refined by repeatedly running CFD simulations at the projected optima and adding the newly obtained data to the training dataset. It is demonstrated that AutoML-GA leads to a better optimum with a lower number of CFD simulations, compared to the use of default hyperparameters. The proposed framework offers the advantage of being a more hands-off approach that can be readily utilized by researchers and engineers in industry who do not have extensive machine learning expertise.

Owoyele, Opeoluwa↗
DOE OSTI2016

Real-time Tumor Oxygenation Changes After Single High-dose Radiation Therapy in Orthotopic and Subcutaneous Lung Cancer in Mice: Clinical Implication for Stereotactic Ablative Radiation Therapy Schedule Optimization

Purpose: To investigate the serial changes of tumor hypoxia in response to single high-dose irradiation by various clinical and preclinical methods to propose an optimal fractionation schedule for stereotactic ablative radiation therapy. Methods and Materials: Syngeneic Lewis lung carcinomas were grown either orthotopically or subcutaneously in C57BL/6 mice and irradiated with a single dose of 15 Gy to mimic stereotactic ablative radiation therapy used in the clinic. Serial [{sup 18}F]-misonidazole (F-MISO) positron emission tomography (PET) imaging, pimonidazole fluorescence-activated cell sorting analyses, hypoxia-responsive element-driven bioluminescence, and Hoechst 33342 perfusion were performed before irradiation (day −1), at 6 hours (day 0), and 2 (day 2) and 6 (day 6) days after irradiation for both subcutaneous and orthotopic lung tumors. For F-MISO, the tumor/brain ratio was analyzed. Results: Hypoxic signals were too low to quantitate for orthotopic tumors using F-MISO PET or hypoxia-responsive element-driven bioluminescence imaging. In subcutaneous tumors, the maximum tumor/brain ratio was 2.87 ± 0.483 at day −1, 1.67 ± 0.116 at day 0, 2.92 ± 0.334 at day 2, and 2.13 ± 0.385 at day 6, indicating that tumor hypoxia was decreased immediately after irradiation and had returned to the pretreatment levels at day 2, followed by a slight decrease by day 6 after radiation. Pimonidazole analysis also revealed similar patterns. Using Hoechst 33342 vascular perfusion dye, CD31, and cleaved caspase 3 co-immunostaining, we found a rapid and transient vascular collapse, which might have resulted in poor intratumor perfusion of F-MISO PET tracer or pimonidazole delivered at day 0, leading to decreased hypoxic signals at day 0 by PET or pimonidazole analyses. Conclusions: We found tumor hypoxia levels decreased immediately after delivery of a single dose of 15 Gy and had returned to the pretreatment levels 2 days after irradiation and had decreased slightly by day 6. Our results indicate that single high-dose irradiation can produce a rapid, but reversible, vascular collapse in tumors.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2019

O(N) Quantum Molecular Dynamics (QMD) Algorithms, Code Requirements, Release Notes, and Library Deployment

The Exascale Computing Project (ECP) Co-design Center for Particle Applications (CoPA) is addressing the challenges for particle applications to run on upcoming exascale computing architectures. This includes short- and long-range molecular dynamics (MD), particle-in-cell (PIC), and quantum molecular dynamics (QMD). This document addresses $\mathcal{O}(N)$ QMD. The goal of this CoPA team is to develop $\mathcal{O}(N)$ linear scaling methods for QMD optimized for new and next generation architectures. This will provide the capability to simulate larger systems and run for longer timeframes than is possible today. Focus is initially on the density matrix construction. This requires efficient matrix formats and linear algebra operations for distributed architectures. A significant reduction in the wallclock time for an MD timestep is necessary to achieve practical large scale QMD. Equally challenging is to develop scalable distributed implementations. This requires the exploration of efficient data formats (dense, sparse, etc.), data decomposition and distribution across compute nodes using block- and graph-based methods. Traditional programming models such as MPI/OpenMP, hybrid MPI/OpenMP/Accelerator, as well as asynchronous task-based models will be explored to achieve maximum parallelism. Parallelism at the level of a matrix operation, solver, and/or application needs to be explored to gain the most benefit on new and emerging architectures. This report serves as a “living document” of the additions and enhancements throughout the ECP CoPA project for $\mathcal{O}(N)$ QMD. The following sections address the applications we are working with, the quantum chemistry libraries (BML and PROGRESS), and our co-design exploration through the use of a proxy application.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Machine learning-based surrogate modeling for data-driven optimization: a comparison of subset selection for regression techniques

Optimization of simulation-based or data-driven systems is a challenging task, which has attracted significant attention in the recent literature. A very efficient approach for optimizing systems without analytical expressions is through fitting surrogate models. Due to their increased flexibility, nonlinear interpolating functions, such as radial basis functions and Kriging, have been predominantly used as surrogates for data-driven optimization; however, these methods lead to complex nonconvex formulations. Alternatively, commonly used regression-based surrogates lead to simpler formulations, but they are less flexible and inaccurate if the form is not known a priori. In this work, we investigate the efficiency of subset selection regression techniques for developing surrogate functions that balance both accuracy and complexity. Subset selection creates sparse regression models by selecting only a subset of original features, which are linearly combined to generate a diverse set of surrogate models. Five different subset selection techniques are compared with commonly used nonlinear interpolating surrogate functions with respect to optimization solution accuracy, computation time, sampling requirements, and model sparsity. Furthermore, our results indicate that subset selection-based regression functions exhibit promising performance when the dimensionality is low, while interpolation performs better for higher dimensional problems.

42 ENGINEERING↗
DOE OSTI2018

Selection of Spectral Variables and Improvement of the Accuracy of Calibration of Temperature by Projection onto Latent Structures Using the Fluorescence Spectra of Yb{sup 3+}:CaF{sub 2}

The efficiency of interval methods for selecting variables and the divisive method of interval optimization by a genetic algorithm is considered with the purpose to increase the accuracy of temperature calibration by projection onto latent structures. The fluorescence spectra of Yb{sup 3+}:CaF{sub 2} recorded in the 880–1120 nm band with a resolution of approximately 0.2 nm in the temperature range from 66 to 150°C in 2°C increments are used as an example. The best result for the root-mean-square error of temperature prediction in the test sample (0.45°C) was obtained by interval projection onto latent structures by a combination of moving windows. The use of methods for selecting spectral variables improved more than two times the accuracy of temperature calibration.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Large-scale quasi-Newton trust-region methods with low-dimensional linear equality constraints

We propose two limited-memory BFGS (L-BFGS) trust-region methods for large-scale optimization with linear equality constraints. Here, the methods are intended for problems where the number of equality constraints is small. By exploiting the structure of the quasi-Newton compact representation, both proposed methods solve the trust-region subproblems nearly exactly, even for large problems. We derive theoretical global convergence results of the proposed algorithms, and compare their numerical effectiveness and performance on a variety of large-scale problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

Linear Energy Transfer-Guided Optimization in Intensity Modulated Proton Therapy: Feasibility Study and Clinical Potential

Purpose: To investigate the feasibility and potential clinical benefit of linear energy transfer (LET) guided plan optimization in intensity modulated proton therapy (IMPT). Methods and Materials: A multicriteria optimization (MCO) module was used to generate a series of Pareto-optimal IMPT base plans (BPs), corresponding to defined objectives, for 5 patients with head-and-neck cancer and 2 with pancreatic cancer. A Monte Carlo platform was used to calculate dose and LET distributions for each BP. A custom-designed MCO navigation module allowed the user to interpolate between BPs to produce deliverable Pareto-optimal solutions. Differences among the BPs were evaluated for each patient, based on dose–volume and LET–volume histograms and 3-dimensional distributions. An LET-based relative biological effectiveness (RBE) model was used to evaluate the potential clinical benefit when navigating the space of Pareto-optimal BPs. Results: The mean LET values for the target varied up to 30% among the BPs for the head-and-neck patients and up to 14% for the pancreatic cancer patients. Variations were more prominent in organs at risk (OARs), where mean LET values differed by a factor of up to 2 among the BPs for the same patient. An inverse relation between dose and LET distributions for the OARs was typically observed. Accounting for LET-dependent variable RBE values, a potential improvement on RBE-weighted dose of up to 40%, averaged over several structures under study, was noticed during MCO navigation. Conclusions: We present a novel strategy for optimizing proton therapy to maximize dose-averaged LET in tumor targets while simultaneously minimizing dose-averaged LET in normal tissue structures. MCO BPs show substantial LET variations, leading to potentially significant differences in RBE-weighted doses. Pareto-surface navigation, using both dose and LET distributions for guidance, provides the means for evaluating a large variety of deliverable plans and aids in identifying the clinically optimal solution.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2016

Integration of TomoPy and the ASTRA toolbox for advanced processing and reconstruction of tomographic synchrotron data

The processing of tomographic synchrotron data requires advanced and efficient software to be able to produce accurate results in reasonable time. In this paper, the integration of two software toolboxes, TomoPy and the ASTRA toolbox, which, together, provide a powerful framework for processing tomographic data, is presented. The integration combines the advantages of both toolboxes, such as the user-friendliness and CPU-efficient methods of TomoPy and the flexibility and optimized GPU-based reconstruction methods of the ASTRA toolbox. It is shown that both toolboxes can be easily installed and used together, requiring only minor changes to existing TomoPy scripts. Furthermore, it is shown that the efficient GPU-based reconstruction methods of the ASTRA toolbox can significantly decrease the time needed to reconstruct large datasets, and that advanced reconstruction methods can improve reconstruction quality compared with TomoPy's standard reconstruction method.

Pelt, Daniël M.↗
DOE OSTI2016

WE-AB-BRA-07: Operating Room Quality Assurance (ORQA) for Spine Surgery Using Known-Component 3D-2D Image Registration

Purpose: Intraoperative x-ray radiography/fluoroscopy is commonly used to qualitatively assess delivery of surgical devices (e.g., spine pedicle screws) but can fail to reliably detect suboptimal placement (e.g., breach of adjacent critical structures). We present a method wherein prior knowledge of the patient and surgical components is leveraged to match preoperative CT and intraoperative radiographs for quantitative assessment of 3D pose. The method presents a new means of operating room quantitative quality assurance (ORQA) that could improve quality and safety, and reduce the frequency of revision surgeries. Methods: The algorithm (known-component registration, KC-Reg) uses patient-specific preoperative CT and parametrically defined surgical component models within a robust 3D-2D registration method to iteratively optimize gradient similarity using the covariance matrix adaptation evolution strategy. Advances from previous work address key challenges to clinical translation: i) absolving the need for offline geometric calibration of the C-arm; and ii) solving multiple component bodies simultaneously, thereby allowing QA in a single step (e.g., spinal construct with 4–20 screws), rather than sequential QA of each component. Performance was tested in a spine phantom with 10 pedicle screws, and first results from clinical studies are reported. Results: Phantom experiments demonstrated median target registration error (TRE) of (1.0±0.3) mm at the screw tip and (0.7°±0.4°) in angulation. The simultaneous multi-body registration approach improved TRE from the previous (sequential) method by 42%, reduced outliers, and fits into the natural workflow. Initial application of KC-Reg in clinical data shows TRE of (2.5±4.5) mm and (4.7°±0.5°). Conclusion: The KC-Reg algorithm offers a potentially valuable method for quantitative QA of the surgical product, using radiographic systems that are already within the surgical arsenal. For spine surgery, the method offers a near-real-time independent check on the quality of surgical product, facilitating immediate revision if necessary and potentially avoiding postoperative morbidity and/or revision surgery. Gerhard Kleinszig and Sebastian Vogt are employees of Siemens Healthcare.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Utilization of Radiation Therapy in Norway After the Implementation of The National Cancer Plan—A National, Population-Based Study

Purpose: To estimate actual utilization rates of radiation therapy (RT) in Norway, describe time trends (1997-2010), and compare these estimates with corresponding optimal RT rates. Methods and Materials: Data from the population-based Cancer Registry of Norway was used to identify all patients diagnosed with cancer and/or treated by RT for cancer in 1997-2010. Radiation therapy utilization rates (RURs) were calculated as (1) the proportion of incident cancer cases who received RT at least once within 1 year of diagnosis (RUR{sub 1Y}); and (2) the proportion who received RT within 5 years of diagnosis (RUR{sub 5Y}). The number of RT treatment courses per incident cancer case (TCI) was also calculated for all cancer sites combined. The actual RURs were compared with corresponding Australian and Canadian epidemiologic- and evidence-based model estimates and criterion-based benchmark estimates of optimal RURs. The TCIs were compared with TCI estimates from the 1997 Norwegian/National Cancer Plan (NCP). Joinpoint regression was used to identify changes in trends and to estimate annual percentage change (APC) in actual RUR{sub 1Y} and actual TCI. Results: The actual RUR{sub 5Y} (all sites) increased significantly to 29% in 2005 but still differed markedly from the Australian epidemiologic- and evidence-based model estimate of 48%. With the exception of RUR{sub 5Y} for breast cancer and RUR{sub 1Y} for lung cancers, all actual RURs were markedly lower than optimal RUR estimates. The actual TCI increased significantly during the study period, reaching 42.5% in 2010, but was still lower than the 54% recommended in the NCP. The trend for RUR{sub 1Y} (all sites) and TCI changed significantly, with the annual percentage change being largest during the first part of the study period. Conclusions: Utilization rates of RT in Norway increased after the NCP was implemented and RT capacity was increased, but they still seem to be lower than optimal levels.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2023

Transient Stability Preventive Control via Tuning the Parameters of Virtual Synchronous Generators

This paper presents an optimal preventive control (OPC) method to improve the power system transient stability via tuning the transient parameter of virtual synchronous generators. The novelty of this work is that we formulate the preventive control as an optimization problem so that inverter parameters can be adjusted at the pre-contingency stage. A reinforcement learning (RL)-driven method is proposed to solve the OPC problem with the fault energy-based reward function. An ANDES-based RL environment is also developed. Versatile functions included in the proposed environment have been presented in this paper. The proposed OPC formulation, the RLdriven method, and the fault energy-based reward function are verified on several standard test systems.

Inverter-based resources, virtual synchronous gene↗
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↗
DOE OSTI2023

Operational control with set-points tuning - application to mobile robots

This paper proposes a novel method for the optimal tuning of set points for multiple-layered control system structure widely seen in robotics and other complex industrial processes composed of a number of subsystems. The terminal sliding mode control (SMC) is used as the low-level control strategy to ensure the stability of subsystems. When uncertainties exist, it can be shown that the deteriorated system performance will be improved by the outer loop with set points tuning. For this purpose, the learning of the new set point is designed to compensate for the effects caused by uncertainties during the system operation. At the same time, the system is proven to stay with the original set point when the compensation is introduced. A practical application to a holonomic mobile robot system is given to illustrate the presented method. Desired results have been obtained.

42 ENGINEERING↗
DOE OSTI2024

Accurate numerical simulations of open quantum systems using spectral tensor trains

Decoherence between qubits is a major bottleneck in quantum computations. Decoherence results from intrinsic quantum and thermal fluctuations as well as noise in the external fields that perform the measurement and preparation processes. With prescribed colored noise spectra for intrinsic and extrinsic noise, we present a numerical method, Quantum Accelerated Stochastic Propagator Evaluation (Q-ASPEN), to solve the time-dependent noise-averaged reduced density matrix in the presence of intrinsic and extrinsic noise. Q-ASPEN is arbitrarily accurate and can be applied to provide estimates for the resources needed to error-correct quantum computations. We employ spectral tensor trains, which combine the advantages of tensor networks and pseudospectral methods, as a variational ansatz to the quantum relaxation problem and optimize the ansatz using methods typically used to train neural networks. Here, the spectral tensor trains in Q-ASPEN make accurate calculations with tens of quantum levels feasible. We present benchmarks for Q-ASPEN on the spin-boson model in the presence of intrinsic noise and on a quantum chain of up to 32 sites in the presence of extrinsic noise. In our benchmark, the memory cost of Q-ASPEN scales as a low-order polynomial in the size of the system once the number of system states surpasses the number of basis functions used in the spectral expansion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Optimal Control Strategy With Efficiency and Reliability Improvement for Offshore DC Microgrids

Offshore microgrids, due to their remote location and lack of external energy support, face significant challenges in wide-range load operation and maintenance. Consequently, efficiency and reliability are critical concerns for converters in offshore dc microgrids. This article presents an optimal control strategy aimed at enhancing both efficiency and reliability. A normalized nonlinear relationship between power loss and thermal stress of a paralleled converter is first established. Based on this, a dual-objective optimization function with an active weight function as well as a system overall performance index is established. The active weight function dynamically adjusts the control priority based on converter efficiency and switching device thermal stress. Then, the optimal power-sharing strategy is derived by the Lagrange multiplier method with the proposed optimal function. Additionally, to accommodate a wide load range, an optimal selection strategy for operating converter combinations is proposed, requiring only low-bandwidth communication. Experiment verification is given to validate the effectiveness of the proposed control strategy. The experiment results demonstrate that the proposed control strategy can improve the overall performance of offshore microgrids by optimizing efficiency and reliability.

24 POWER TRANSMISSION AND DISTRIBUTION↗