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

Bridging Equipment Reliability Data and Robust Decisions in a Plant Operation Context

In order to reduce operation and maintenance (O&M) costs, nuclear power plants (NPPs) are moving from corrective and periodic maintenance to predictive maintenance strategies. Such transition requires changes on the data that needs to be retrieved and on the type of decision processes to be employed. Advanced monitoring and data analysis technologies are essential to support predictive strategies. They can in fact provide precise information about health of a component, track its degradation trends, and provide information of its expected failure time. With such information, maintenance operations for a component can be performed right before its expected failure time. This dynamic context of O&M operations requires new methods to analyze data, propagate component health information from the component to the system level, and optimize plant resources. In this respect, the risk informed asset management (RIAM) project has been tasked to develop and test this new class of methods into a risk analytics toolset. This toolset consists of data analytics tools coupled with reliability methods designed to manage plant assets and performances in a predictive maintenance context. This report shows the latest improvements on such development and the initial testing of our methods on the three main research areas that the RIAM project is focusing on. These areas are the following: equipment reliability data analytics, system reliability modeling, and plant resources optimization methods. We show how the methods developed in these areas can support predictive maintenance strategies by: 1) analyzing equipment reliability data (either in numeric and textual form), 2) assessing component and system health through an innovative margin-based reliability approach, and 3) identifying the most critical components and set optimal maintenance schedule based on plant economic and operational constraints.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Optimization of one-step in situ transesterification method for accurate quantification of EPA in Nannochloropsis gaditana

Microalgae are a valuable source of lipid feedstocks for biodiesel and valuable omega-3 fatty acids. Nannochloropsis gaditana has emerged as a promising producer of eicosapentaenoic acid (EPA) due to its fast growth rate and high EPA content. In the present study, the fatty acid profile of Nannochloropsis gaditana was found to be naturally high in EPA and devoid of docosahexaenoic acid (DHA), thereby providing an opportunity to maximize the efficacy of EPA production. Using an optimized one-step in situ transesterification method (methanol:biomass = 90 mL/g; HCl 5% by vol.; 70 °C; 1.5 h), the maximum fatty acid methyl ester (FAME) yield of Nannochloropsis gaditana cultivated under rich condition was quantified as 10.04% ± 0.08% by weight with EPA-yields as high as 4.02% ± 0.17% based on dry biomass. The total FAME and EPA yields were 1.58- and 1.23-fold higher separately than that obtained using conventional two-step method (solvent system: methanol and chloroform). Furthermore, this one-step in situ method provides a fast and simple method to measure fatty acid methyl ester (FAME) yields and could serve as a promising method to generate eicosapentaenoic acid methyl ester from microalgae.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2013

Large-Scale Optimization for Bayesian Inference in Complex Systems

The SAGUARO (Scalable Algorithms for Groundwater Uncertainty Analysis and Robust Optimization) Project focused on the development of scalable numerical algorithms for large-scale Bayesian inversion in complex systems that capitalize on advances in large-scale simulation-based optimization and inversion methods. The project was a collaborative effort among MIT, the University of Texas at Austin, Georgia Institute of Technology, and Sandia National Laboratories. The research was directed in three complementary areas: efficient approximations of the Hessian operator, reductions in complexity of forward simulations via stochastic spectral approximations and model reduction, and employing large-scale optimization concepts to accelerate sampling. The MIT₋Sandia component of the SAGUARO Project addressed the intractability of conventional sampling methods for large-scale statistical inverse problems by devising reduced-order models that are faithful to the full-order model over a wide range of parameter values; sampling then employs the reduced model rather than the full model, resulting in very large computational savings. Results indicate little effect on the computed posterior distribution. On the other hand, in the Texas₋Georgia Tech component of the project, we retain the full-order model, but exploit inverse problem structure (adjoint-based gradients and partial Hessian information of the parameter-to-observation map) to implicitly extract lower dimensional information on the posterior distribution; this greatly speeds up sampling methods, so that fewer sampling points are needed. We can think of these two approaches as "reduce then sample" and "sample then reduce." In fact, these two approaches are complementary, and can be used in conjunction with each other. Moreover, they both exploit deterministic inverse problem structure, in the form of adjoint-based gradient and Hessian information of the underlying parameter-to-observation map, to achieve their speedups.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

SU-E-I-77: A Noise Reduction Technique for Energy-Resolved Photon-Counting Detectors

Purpose: Finding the optimal energy threshold setting for an energy-resolved photon-counting detector has an important impact on the maximization of contrast-to-noise-ratio (CNR). We introduce a noise reduction method to enhance CNR by reducing the noise in each energy bin without altering the average gray levels in the projection and image domains. Methods: We simulated a four bin energy-resolved photon-counting detector based on Si with a 10 mm depth of interaction. TASMIP algorithm was used to simulate a spectrum of 65 kVp with 2.7 mm Al filter. A 13 mm PMMA phantom with hydroxyapatite and iodine at different concentrations (100, 200 and 300 mg/ml for HA, and 2, 4, and 8 mg/ml for Iodine) was used. Projection-based and Image-based energy weighting methods were used to generate weighted images. A reference low noise image was used for noise reduction purposes. A Gaussian-like weighting function which computes the similarity between pixels of interest was calculated from the reference image and implemented on a pixel by pixel basis for the noisy images. Results: CNR improvement compared to different methods (Charge-Integrated, Photon-Counting and Energy-Weighting) and after noise reduction was highly task-dependent. The CNR improvement with respect to the Charge-Integrated CNR for hydroxyapatite and iodine were 1.8 and 1.5, respectively. In each of the energy bins, the noise was reduced by approximately factor of two without altering their respective average gray levels. Conclusion: The proposed noise reduction technique for energy-resolved photon-counting detectors can significantly reduce image noise. This technique can be used as a compliment to the current energy-weighting methods in CNR optimization.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

An inexact semismooth Newton method with application to adaptive randomized sketching for dynamic optimization

In many applications, one can only access the inexact gradients and inexact hessian times vector products. Thus it is essential to consider algorithms that can handle such inexact quantities with a guaranteed convergence to solution. An inexact adaptive and provably convergent semismooth Newton method is considered to solve constrained optimization problems. In particular, dynamic optimization problems, which are known to be highly expensive, are the focus. A memory efficient semismooth Newton algorithm is introduced for these problems. The source of efficiency and inexactness is the randomized matrix sketching. Further, applications to optimization problems constrained by partial differential equations are also considered.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

Quantum-mechanical predictions of electron-induced ionization cross sections of DNA components

Ionization of biomolecules remains still today rarely investigated on both the experimental and the theoretical sides. In this context, the present work appears as one of the first quantum mechanical approaches providing a multi-differential description of the electron-induced ionization process of the main DNA components for impact energies ranging from the target ionization threshold up to about 10 keV. The cross section calculations are here performed within the 1st Born approximation framework in which the ejected electron is described by a Coulomb wave whereas the incident and the scattered electrons are both described by a plane wave. The biological targets of interest, namely, the DNA nucleobases and the sugar-phosphate backbone, are here described by means of the GAUSSIAN 09 system using the restricted Hartree-Fock method with geometry optimization. The theoretical predictions also obtained have shown a reasonable agreement with the experimental total ionization cross sections while huge discrepancies have been pointed out with existing theoretical models, mainly developed within a semi-classical framework.

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

Drop-on-demand sample delivery for studying biocatalysts in action at X-ray free-electron lasers

X-ray crystallography at X-ray free-electron laser (XFEL) sources is a powerful method for studying macromolecules at biologically relevant temperatures. Moreover, when combined with complementary techniques like X-ray emission spectroscopy (XES), both global structures and chemical properties of metalloenzymes can be obtained concurrently, providing new insights into the interplay between the protein structure/dynamics and chemistry at an active site. However, implementing such a multimodal approach can be compromised by conflicting requirements to optimize each individual method. In particular, the method used for sample delivery greatly impacts the data quality. We present here a new, robust way of delivering controlled sample amounts on demand using acoustic droplet ejection coupled with a conveyor belt drive that is optimized for crystallography and spectroscopy measurements of photochemical and chemical reactions over a wide range of time scales. Studies with photosystem II, the phytochrome photoreceptor, and ribonucleotide reductase R2 illustrate the power and versatility of this method.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Optimization of Random Phase Approximation Calculations for Improved Energies of Molecules, Solids, and Surfaces

We present an optimized random phase approximation method (optRPA26) that significantly improves upon conventional RPA through an optimized choice of reference orbitals and energy components, rather than a modification of the RPA correlation functional itself. The method employs an empirically constructed hybrid functional to generate DFT orbitals to evaluate the RPA correlation energy, which is then scaled by a constant. Comprehensive benchmarks across molecules, bulk solids, and surface systems demonstrate that optRPA26 consistently achieves high accuracy, with mean absolute errors of 0.05 eV for W4-11-RE reaction energies, 0.07 eV for cohesive energies, 0.09 eV for metal oxide formation energies, 0.11–0.12 eV for adsorption of small molecules on metals, and 0.06 eV for adsorption on oxides. In addition, optRPA26 correctly captures phase stability in metal oxides and magnetic metals. The optRPA26 approach can be run using standard RPA implementations, highlighting its potential as a general-purpose reference method that can accurately capture covalent, ionic, metallic, and van der Waals bonding in molecules, solids, and interfaces.

Adsorption↗
DOE OSTI2023

Differentiable Neural Architecture, Mixed Precision and Accelerator Co-Search

Quantization, effective Neural Network architecture, and efficient accelerator hardware are three important design paradigms to maximize accuracy and efficiency. Mixed Precision Quantization is a process of assigning different precision to different Neural Network layers for optimized inference. Neural Architecture Search (NAS) is a process of automatically designing the neural network for a task and can also be extended to search for the precision of each weight and activation matrix. In this paper, we develop the following three methods: (i) Fast Differentiable Hardware-aware Mixed Precision Quantization Search method to find optimal precision, (ii) Joint Differentiable hardware-aware Architecture and Mixed Precision Quantization Co-search, (iii) Joint Accelerator, Architecture, and Precision triple co-search to find best possibilities in all the three worlds. We demonstrate the effectiveness of our proposed methods targeting Bitfusion accelerator by searching mixed precision models on MobilenetV2. We achieve better accuracy-latency trade-off models than the manually designed and previously proposed search methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Optimization of micropillar sequences for fluid flow sculpting

Inertial fluid flow deformation around pillars in a microchannel is a new method for controlling fluid flow. Sequences of pillars have been shown to produce a rich phase space with a wide variety of flow transformations. Previous work has successfully demonstrated manual design of pillar sequences to achieve desired transformations of the flow cross section, with experimental validation. However, such a method is not ideal for seeking out complex sculpted shapes as the search space quickly becomes too large for efficient manual discovery. We explore fast, automated optimization methods to solve this problem. We formulate the inertial flow physics in microchannels with different micropillar configurations as a set of state transition matrix operations. These state transition matrices are constructed from experimentally validated streamtraces for a fixed channel length per pillar. This facilitates modeling the effect of a sequence of micropillars as nested matrix-matrix products, which have very efficient numerical implementations. With this new forward model, arbitrary micropillar sequences can be rapidly simulated with various inlet configurations, allowing optimization routines quick access to a large search space. We integrate this framework with the genetic algorithm and showcase its applicability by designing micropillar sequences for various useful transformations. We computationally discover micropillar sequences for complex transformations that are substantially shorter than manually designed sequences. We also determine sequences for novel transformations that were difficult to manually design. Finally, we experimentally validate these computational designs by fabricating devices and comparing predictions with the results from confocal microscopy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Petascale solvers for anisotropic PDEs in atmospheric modelling on GPU clusters

Memory bound applications such as solvers for large sparse systems of equations remain a challenge for GPUs. Fast solvers should be based on numerically efficient algorithms and implemented such that global memory access is minimised. To solve systems with trillions ( O ( 10 12 ) ) unknowns the code has to make efficient use of several million individual processor cores on large GPU clusters. We describe the multi-GPU implementation of two algorithmically optimal iterative solvers for anisotropic PDEs which are encountered in (semi-) implicit time stepping procedures in atmospheric modelling. In this application the condition number is large but independent of the grid resolution and both methods are asymptotically optimal, albeit with different absolute performance. In particular, an important constant in the discretisation is the CFL number; only the multigrid solver is robust to changes in this constant. We parallelise the solvers and adapt them to the specific features of GPU architectures, paying particular attention to efficient global memory access. We achieve a performance of up to 0.78 PFLOPs when solving an equation with 0.55 · 10 12 unknowns on 16384 GPUs; this corresponds to about 3% of the theoretical peak performance of the machine and we use more than 40% of the peak memory bandwidth with a Conjugate Gradient (CG) solver. Although the other solver, a geometric multigrid algorithm, has a slightly worse performance in terms of FLOPs per second, overall it is faster as it needs less iterations to converge; the multigrid algorithm can solve a linear PDE with half a trillion unknowns in about one second.

Computer Science↗
DOE OSTI2024

Machine learning enhanced characterization and optimization of photonic cured MAPbI 3 for efficient perovskite solar cells

Photonic curing (PC) can facilitate high-speed perovskite solar cell (PSC) manufacturing because it uses high-intensity light pulses to crystallize perovskite films in milliseconds. However, optimizing PC conditions is challenging due to its many variables, and using power conversion efficiency (PCE) as the optimization metric is both time-consuming and labor-intensive. This work presents a machine learning (ML) approach to optimize PC conditions for fabricating methylammonium lead iodide (MAPbI 3 ) films by quantitatively comparing their ultraviolet-visible (UV-vis) absorbance spectra to thermal annealed (TA) films using four similarity metrics. We perform Bayesian optimization coupled with Gaussian process regression (BO-GP) to minimize the similarity metrics. Refining PC conditions using active learning based on BO-GP models, we achieve a PC MAPbI3 film with an absorbance spectrum closely matching a TA reference film, which is further verified by its crystalline and morphological properties. Thus, we demonstrate that the UV-vis absorption spectrum can accurately proxy film quality. Additionally, we use an AI-based segmentation model for a more efficient grain size analysis. However, when we use the optimized PC condition to fabricate PSCs, we find that interaction between MAPbI 3 and the hole transport layer (HTL) during PC critically degrades the PSC performance. By adding a buffer layer between the HTL and MAPbI 3 , the optimized PC PSCs produce a champion PCE of 11.8%, comparable to the TA reference of 11.7%. Using UV-vis similarity metrics instead of device PCE as the objective in our BO-GP method accelerates the optimization of PC processing conditions for MAPbI 3 films.

14 SOLAR ENERGY↗
DOE OSTI2024

Method Validation Summary for L16.2 AD-ISO-0015 for ISO-17025 / NFAC Applications

Methods for the analysis of semi-volatile organic compounds are employed at Savannah River National Laboratory (SRNL) for routine and non-routine samples from tank waste, process control, waste acceptance, and a wide array of other process and research samples. A method has been developed by SRNL, based on existing methods for semi-volatile analysis, for the analysis of nitroaromatic high explosives by gas chromatography / mass spectrometry in soils and sediments. Relative to past work performed at SRNL on nitroaromatic explosives, this developed, optimized, and validated method can achieve lower limits of detection and quantitation, greater precision, lower bias, higher linearity, and accuracy across a greater linear range.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

INTEGRATION OF DATA ANALYTICS WITH SYSTEM HEALTH PROGRAMS

Industry equipment reliability and asset management programs are essential elements that help ensure the safe and economical operation of nuclear power plants. The effectiveness of these programs is addressed in several industry developed and regulatory programs. However, these programs have proven to be labor intensive and expensive. There is an opportunity to significantly enhance the collection, analysis, and use of this information to provide more cost-effective plant operation. Additionally, there is an acute industry need to leverage advanced technology to reduce costs and improve operational effectiveness. The goal of this paper is to provide effective and efficient analytical methods and tools to support risk-informed decisions for the equipment reliability and asset management programs at nuclear power plants. This is accomplished by creating a direct bridge between component health/lifecycle data and decision making (e.g., maintenance scheduling and project prioritization). Here we are supporting typical system engineer decisions regarding maintenance activity scheduling and component ageing management. This is performed in a risk-informed context where herein the term “risk” is broadly constructed to include both plant reliability and economics. This framework combines data analytics tools to analyze equipment reliability data with risk-informed methods designed to support system engineer decisions (e.g., maintenance and replacement schedules, optimal maintenance posture) in a customizable workflow. A challenge is that the structure of this workflow strongly depends on the decision that needs to be made, the type of data available, and the constraints that need to be considered. Current methods are designed to provide specific answers to specific problems; however, these methods might prove to be inadequate even when problem settings slightly change (e.g., different types of requirements, additional dependencies between system reliability and economics). We tackled this challenge by designing framework in a flexible and modular fashion such that the user can assemble and customize his/her own workflow that integrates SSC economic lifecycle models (e.g., maintenance and replacement costs), system reliability models, and optimization methods.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Integrating novel stellarator single-stage optimization algorithms to design the Columbia stellarator experiment

Abstract The Columbia Stellarator eXperiment (CSX), currently being designed at Columbia University, aims to test theoretical predictions related to QA plasma behavior, and to pioneer the construction of an optimized stellarator using three-dimensional, non-insulated high-temperature superconducting (NI-HTS) coils. The magnetic configuration is generated by a combination of two circular planar poloidal field (PF) coils and two 3D-shaped interlinked (IL) coils, with the possibility to add windowpane coils to enhance shaping and experimental flexibility. The PF coils and vacuum vessel are repurposed from the former Columbia Non-Neutral Torus experiment, while the IL coils will be custom-wound in-house using NI-HTS tapes. To obtain a plasma shape that meets the physics objectives with a limited number of coils, novel single-stage optimization techniques are employed, optimizing both the plasma and coils concurrently, in particular targeting a tight aspect ratio QA plasma and minimized strain on the HTS tape. Despite the increased complexity due to the expanded degrees of freedom, these methods successfully identify optimized plasma geometries that can be realized by coils meeting engineering specifications. This paper discusses the derivation of the constraints and objectives specific to CSX, and describe how two recently developed single-stage optimization methodologies are applied to the design of CSX. A set of selected configurations for CSX is then described in detail.

Baillod, A. (ORCID:0000000303529180)↗
DOE OSTI2014

SU-E-T-234: Modulated Photon Radiotherapy (XMRT):The Impact of Incorporating Energy Modulation Into Intensity Modulated Radiotherapy (IMRT) Optimization

Purpose: To develop a new radiotherapy plan optimization technique that, for a given organ geometry, will find the optimal photon beam energies and fluences to produce a desirable dose distribution. This new modulated (both in energy and fluence) photon radiotherapy (XMRT) was compared with intensity modulated radiotherapy (IMRT) for a simple organ geometry. Methods: The XMRT optimization was formulated using a linear programming approach where the objective function is the mean dose to the healthy organs and dose-point constraints were assigned to each organ of interest. The organ geometry consisted of a target, two organs at risk (OARs), and normal tissue. A seven-equispaced-coplanar beam arrangement was used. For conventional IMRT, only 6 MV beams were available, while XMRT was optimized using 6 and 18 MV beams. A prescribed dose (PD) of 72 GY was assigned to the target, with upper and lower bounds of 110% and 95% of the PD, respectively. Both OARs were assigned a maximum dose of 64 Gy, while the normal tissue was assigned a maximum dose of 66 Gy. A numerical solver, Gurobi, generated solutions for the XMRT and IMRT problems. The dose-volume histograms from IMRT and XMRT solutions were compared. Results: The maximum, minimum, mean, and homogeneity of the dose to the target were comparable between IMRT and XMRT. Though IMRT had improved dose conformity relative to XMRT, XMRT reduced the mean dose to both OARs by more than 1 Gy. For normal tissue, an increase of 5 Gy in mean dose and 27 percent in integral dose was seen for IMRT relative to XMRT. Conclusion: This work demonstrates the benefits of simultaneously modulating photon beam energy and fluence using our XMRT approach in a given phantom geometry. While target coverage was comparable, dose to healthy structures was reduced using XMRT.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2016

TH-EF-BRB-04: 4π Dynamic Conformal Arc Therapy Dynamic Conformal Arc Therapy (DCAT) for SBRT

Purpose: To develop an efficient and effective trajectory optimization methodology for 4π dynamic conformal arc treatment (4π DCAT) with synchronized gantry and couch motion; and to investigate potential clinical benefits for stereotactic body radiation therapy (SBRT) to breast, lung, liver and spine tumors. Methods: The entire optimization framework for 4π DCAT inverse planning consists of two parts: 1) integer programming algorithm and 2) particle swarm optimization (PSO) algorithm. The integer programming is designed to find an optimal solution for arc delivery trajectory with both couch and gantry rotation, while PSO minimize a non-convex objective function based on the selected trajectory and dose-volume constraints. In this study, control point interaction is explicitly taken into account. Beam trajectory was modeled as a series of control points connected together to form a deliverable path. With linear treatment planning objectives, a mixed-integer program (MIP) was formulated. Under mild assumptions, the MIP is tractable. Assigning monitor units to control points along the path can be integrated into the model and done by PSO. The developed 4π DCAT inverse planning strategy is evaluated on SBRT cases and compared to clinically treated plans. Results: The resultant dose distribution of this technique was evaluated between 3D conformal treatment plan generated by Pinnacle treatment planning system and 4π DCAT on a lung SBRT patient case. Both plans share the same scale of MU, 3038 and 2822 correspondingly to 3D conformal plan and 4π DCAT. The mean doses for most of OARs were greatly reduced at 32% (cord), 70% (esophagus), 2.8% (lung) and 42.4% (stomach). Conclusion: Initial results in this study show the proposed 4π DCAT treatment technique can achieve better OAR sparing and lower MUs, which indicates that the developed technique is promising for high dose SBRT to reduce the risk of secondary cancer.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

A higher‐order unsplit 2D direct Eulerian finite volume method for two‐material compressible flows based on the MOOD paradigms

SUMMARY A higher‐order unsplit multi‐dimensional discretization of the diffuse interface model for two‐material compressible flows proposed by R. Saurel, F. Petitpas and R. A. Berry in 2009 is developed. The proposed higher‐order method is based on the concepts of the Multidimensional Optimal Order Detection (MOOD) method introduced in three recent papers for single‐material flows. The first‐order unsplit multi‐dimensional Finite Volume discretization presented by SPB serves as foundation for the development of the higher‐order unlimited schemes. Specific detection criteria along with a novel decrementing algorithm for the MOOD method are designed in order to deal with the complexity of multi‐material flows. Numerically, we compare errors and computational times on several 1D problems (stringent shock tube and cavitation problems) computed on 2D meshes with the second‐ and fourth‐order MOOD methods using a classical MUSCL method as reference. Several simulations of a 2D shocked R22 bubble in the air are also presented on Cartesian and unstructured meshes with the second‐ and fourth‐order MOOD methods, and qualitative comparisons confirm the conclusions obtained with 1D problems. These numerical results demonstrate the robustness of the MOOD approach and the interest of using more than second‐order methods even for locally singular solutions of complex physics models. Copyright © 2014 John Wiley & Sons, Ltd.

Diot, S.↗