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,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

A Method for Evaluating Volt-VAR Optimization Field Demonstrations

In a regulated business environment a utility must be able to validate that deployed technologies provide quantifiable benefits to the end-use customers. For traditional technologies there are well established procedures for determining what benefits will be derived from the deployment. But for many emerging technologies procedures for determining benefits are less clear and completely absent in some cases. Volt-VAR Optimization is a technology that is being deployed across the nation, but there are still numerous discussions about potential benefits and how they are achieved. This paper will present a method for the evaluation, and quantification of benefits, for field deployments of Volt-VAR Optimization technologies. In addition to the basic methodology, the paper will present a summary of results, and observations, from two separate Volt-VAR Optimization field evaluations using the proposed method.

optimization↗
DOE OSTI2018

Integration of optimal operational dispatch and controller determined dynamics for microgrid survivability

We present that the reliability and resilience of the electrical power grids are essential to industry, economy and society. Microgrids that are able to island from the bulk electrical grid are one technology that may vastly improve electrical power service to customer loads. To achieve these improvements, an islanded microgrid should be able to operate through the loss of one of its generators without shedding electrical load. The loss of one generator will typically result in significant additional loads, including transient overloads, being placed on the remaining generators. There is also the possibility of additional generator tripping during these processes (i.e. cascading failures), which would likely result in the collapse of the microgrid. The novelty of our work consists in incorporating dynamic models of generator controllers into a microgrid optimal dispatch formulation with the ultimate goal to avoid operational failures and ensure the “survivability” of all-inverter microgrids to generator loss and transient overloads. The integration of generator and controller dynamics into the optimal dispatch formulation significantly increases the computational complexity. As we develop algorithms to restore speed of the optimization, our method can be readily implemented into a new operational strategy capable of an unprecedented level of reliability against generator contingencies. In addition, we quantitatively illustrate the effect of the survivability constraints on the microgrid operating costs and how the related trade-off between capital and operating costs should be taken into account at the design stage. Finally, the methods developed here also apply to the dispatch of off-grid microgrids.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2019

Coordination of Behind-the-Meter Energy Storage and Building Loads: Optimization with Deep Learning Models*

With the increasing penetration of renewable energy systems and energy storage systems in buildings, it is critical to optimize system operation. In this paper, we develop an optimization strategy to minimize the operation cost as well as maintain indoor thermal comfort, for a building integrated with battery and PV. A Recurrent Neural Network (RNN) is used to predict building thermal load and zone temperatures. A black-box optimization algorithm known as nonlinear optimization by mesh adaptive direct search (NOMADS) is employed in simulation to provide look-ahead optimal battery dispatch and zone temperature set-point schedules so that the operation cost is minimized. Field data collected from a medium sized office building at Pacific Northwest National Laboratory (PNNL) integrated with Photovoltaic (PV) system model and battery storage system model is used to demonstrate the proposed methodology. Compared with rule based methods, the optimization strategy obtained a lower cost of operation while satisfying comfort constraints.

Chen, Yimin↗
DOE OSTI2025

Rapid quantification of whole seed fatty acid amount, composition, and shape phenotypes from diverse oilseed species with large differences in seed size

Seed oils are widely used in the food, biofuel, and industrial feedstock industries, with their utility and value determined by total oil content and fatty acid composition. Current high throughput seed oil analysis methods either lack accuracy in total fatty acid profiling or require extensive labor for lipid extraction prior to derivatization to fatty acid methyl esters (FAME) and quantification by gas chromatography (GC). Alternatively, direct whole seed FAME production methods have been developed for the very small seeds in the model species Arabidopsis thaliana but these have generally not been adapted to larger seeds of most oilseed crops. High-throughput direct whole seed FAME production methods were optimized for seeds up to 5 mg each utilizing acid-catalyzed esterification. For the oilseed species Camelina sativa, Thlaspi avernse (pennycress), Cuphea viscosissima, and Brassica napus (var. Canola), the total seed fatty acid content and composition from direct seed esterification to FAME matched that of lipid extract derivatization demonstrating the accuracy of the methods. In combination with seed phenotyping using GridFree, this approach enabled the development of a rapid pipeline for simultaneous seed weight, count, size/shape phenotyping, and oil analysis. For the larger and tougher seeds produced by Limnanthes alba (Meadowfoam) and Cannabis sativa L. (hemp) the whole seed acid-based method proved insufficient, and prior laborious homogenization of seeds was required. Therefore, a rapid one-tube bead homogenization and base catalyzed-esterification method was developed. Base-derived fatty acid esterification cannot derivatize free fatty acids leading to slightly lower total seed fatty acid than acid-catalyzed methods, however the seed oil content and fatty acid composition that is valuable for screening large numbers of samples in research populations was accurately measured. New rapid whole seed fatty acid esterification and phenotyping protocols were developed to accurately assess oilseed lipid content. These methods are particularly valuable in oilseed research, breeding, and engineering applications where efficient analysis of large numbers of samples and accurate oil fatty acid profiling is essential. While having been developed for current and emerging oilseed crops, these methods also provide a foundation from which protocols might be established for new and emerging crop species.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

Acceleration of the Particle Swarm Optimization for Peierls–Nabarro modeling of dislocations in conventional and high-entropy alloys

Dislocations are among the most important defects in determining the mechanical properties of both conventional alloys and high-entropy alloys. The Peierls-Nabarro model supplies an efficient pathway to their geometries and mobility. The difficulty in solving the integro-differential Peierls-Nabarro equation is how to effectively avoid the local minima in the energy landscape of a dislocation core. Among the other methods to optimize the dislocation core structures, we choose the algorithm of Particle Swarm Optimization, an algorithm that simulates the social behaviors of organisms. By employing more particles (bigger swarm) and more iterative steps (allowing them to explore for longer time), the local minima can be effectively avoided. But this would require more computational cost. The advantage of this algorithm is that it is readily parallelized in modern high computing architecture. We demonstrate the performance of our parallelized algorithm scales linearly with the number of employed cores.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Optimizing the Accelerated Recursive Doubling Algorithm for Block Tridiagonal Systems of Equations

The need to solve block tridiagonal systems with hundreds or thousands of right-hand sides for the same block tridiagonal matrix is common in a variety of disciplines. To meet this need, the Accelerated Recursive Doubling Algorithm was developed. After a right-hand side independent phase, the algorithm allows for the quick, online calculation of solutions for different right-hand sides. In this work, we present methods to optimize the Accelerated Recursive Doubling Algorithm in memory usage and computation time in a hybrid parallelization model. The right-hand side independent phase of the naïve implementation takes ≥ 11/3 the amount of memory required to store the tridiagonal matrix, while our implementation reduces the fraction to ≈ 5/3 . The right-hand side dependent phase of the naïve implementation takes ≥ 6 times the amount of memory required to store the right-hand side, while our implementation reduces the fraction to ≈ 3. The computation time for the independent phase is reduced to ≈ 2/3 times that of the naïve implementation, while the computation time for the dependent phase is reduced to ≈ 5/9 . With increasing numbers of shared-memory threads q on every distributed processing element, we have O(q) theoretical speedup.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Adaptive sampling quasi-Newton methods for zeroth-order stochastic optimization

Here, we consider unconstrained stochastic optimization problems with no available gradient information. Such problems arise in settings from derivative-free simulation optimization to reinforcement learning. We propose an adaptive sampling quasi-Newton method where we estimate the gradients using finite differences of stochastic function evaluations within a common random number framework. We develop modified versions of a norm test and an inner product quasi-Newton test to control the sample sizes used in the stochastic approximations and provide global convergence results to the neighborhood of a locally optimal solution. We present numerical experiments on simulation optimization problems to illustrate the performance of the proposed algorithm. When compared with classical zeroth-order stochastic gradient methods, we observe that our strategies of adapting the sample sizes significantly improve performance in terms of the number of stochastic function evaluations required.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Anisotropic thermal transport property of defect-free GaN

Non-equilibrium molecular dynamics (MD) simulation is performed to calculate the thermal conductivity of defect-free GaN along three high-symmetry directions. It is found that the thermal conductivity along [001] direction is about 25% higher than that along [100] or [120] direction. The calculated phonon dispersion relation and iso-energy surface from lattice dynamics show that the difference of the sound speeds among the three high-symmetry directions is quite small for the same mode. However, the variation of phonon irradiation with direction is qualitatively consistent with that of the calculated thermal conductivity. Our results indicate that the anisotropic thermal conductivity may partly result from the phonons in the low-symmetry region of the first Brillouin zone due to phonon focus effects, even though the elastic properties along the three high-symmetry directions are nearly isotropic. Thus, the phonon irradiation is able to better describe the property of thermal conductivity as compared to the commonly used phonon dispersion relation. The present investigations uncover the physical origin of the anisotropic thermal conductivity in defect-free GaN, which would provide an important guide for optimizing the thermal management of GaN-based device.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Profitability, Scalability and Optimized Siting for Waste-to-Energy Supply Chains: A Case Study on Hydrothermal Liquefaction

Waste-to-Energy (WtE) conversion technologies have the potential to simultaneously reduce waste volumes, treatment-related emissions, and carbon intensity while producing valuable energy services. However, questions remain regarding under which conditions emerging WtE technologies are feasible to deploy. In this project, we develop an optimized siting method for WtE technologies to assess cost-effective WtE processor locations, throughput scales, and profitability. We apply the model using Hydrothermal Liquefaction (HTL) as an example WtE technology. We develop calibrated capital and operating expense cost curves based on literature data to study the techno-economic characteristics of WtE supply chains. Prior to solving the siting optimization, we analyze the model parameters to describe organic waste management system behavior under four different HTL deployment configurations: Co-Processing, where HTL plants send biocrude to existing conventional refineries via assumed pipelines; Co-Location, where biorefineries are co-located with HTL plants; and two alternative Standalone Biorefining cases, where distributed HTL plants transport biocrude intermediate via either assumed pipelines or trucking to centralized biorefineries. We also define and apply several feedstock gate fee calculations, representing different revenue distribution options, to investigate how bioproduct value could impact WtE supply chain economics in terms of waste producer and processor profit or cost reduction.

09 BIOMASS FUELS↗
DOE OSTI2014

First principles molecular dynamics without self-consistent field optimization

We present a first principles molecular dynamics approach that is based on time-reversible extended Lagrangian Born-Oppenheimer molecular dynamics [A. M. N. Niklasson, Phys. Rev. Lett. 100, 123004 (2008)] in the limit of vanishing self-consistent field optimization. The optimization-free dynamics keeps the computational cost to a minimum and typically provides molecular trajectories that closely follow the exact Born-Oppenheimer potential energy surface. Only one single diagonalization and Hamiltonian (or Fockian) construction are required in each integration time step. The proposed dynamics is derived for a general free-energy potential surface valid at finite electronic temperatures within hybrid density functional theory. Even in the event of irregular functional behavior that may cause a dynamical instability, the optimization-free limit represents a natural starting guess for force calculations that may require a more elaborate iterative electronic ground state optimization. Our optimization-free dynamics thus represents a flexible theoretical framework for a broad and general class of ab initio molecular dynamics simulations.

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

Ground-state structures of Hafnium clusters

Hafnium (Hf) is a very large tetra-valence d-block element which is able to form relatively long covalent bond. Researchers are interested to search for substitution to silicon in the semi-conductor industry. We attempt to obtain the ground-state structures of small Hf clusters at both empirical and density-functional theory (DFT) levels. For calculations at the empirical level, charge-optimized many-body functional potential (COMB) is used. The lowest-energy structures are obtained via a novel global-minimum search algorithm known as parallel tempering Monte-Carlo Basin-Hopping and Genetic Algorithm (PTMBHGA). The virtue of using COMB potential for Hf cluster calculation lies in the fact that by including the charge optimization at the valence shells, we can encourage the formation of proper bond hybridization, and thus getting the correct bond order. The obtained structures are further optimized using DFT to ensure a close proximity to the ground-state.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

PuReMD-GPU: A reactive molecular dynamics simulation package for GPUs

We present an efficient and highly accurate GP-GPU implementation of our community code, PuReMD, for reactive molecular dynamics simulations using the ReaxFF force field. PuReMD and its incorporation into LAMMPS (Reax/C) is used by a large number of research groups worldwide for simulating diverse systems ranging from biomembranes to explosives (RDX) at atomistic level of detail. The sub-femtosecond time-steps associated with ReaxFF strongly motivate significant improvements to per-timestep simulation time through effective use of GPUs. This paper presents, in detail, the design and implementation of PuReMD-GPU, which enables ReaxFF simulations on GPUs, as well as various performance optimization techniques we developed to obtain high performance on state-of-the-art hardware. Comprehensive experiments on model systems (bulk water and amorphous silica) are presented to quantify the performance improvements achieved by PuReMD-GPU and to verify its accuracy. In particular, our experiments show up to 16× improvement in runtime compared to our highly optimized CPU-only single-core ReaxFF implementation. PuReMD-GPU is a unique production code, and is currently available on request from the authors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Similarity Metric for Data Optimization and Efficient Training of Reactive Machine Learning Force Fields for Hydrocarbon Radiolysis

Radiolysis is a common approach to sterilize polymers, chemically modify them for upcycling, and accelerate their decomposition for recycling purposes. Reactive molecular dynamics (MD) simulations provide a powerful tool to generate atomic-level trajectories of the reactive processes and quantify radiolytic chemical degradation pathways. For this, machine learning (ML) surrogate models for reactive force fields with quantum mechanical accuracy are now widely used, which require ML training data sets that can provide information on atomic environments for target chemical systems. However, radiolysis chemistry can be highly complex and diverse, which poses significant challenges for generating training data to parametrize ML models. In this regard, we developed a method for optimizing the training data set using a cosine similarity metric to help guide training set selection for radiolysis of polyethylene, a model hydrocarbon polymer, as well as to enhance the transferability of our reactive ML force field (MLFF) to a variety of molecular and polymeric systems. Our approach performs atom-by-atom comparisons between local atomic environments to pinpoint important data points associated with rare and localized events, such as radiolysis damage within structures. We apply this approach to train the Chebyshev Interaction Model for Efficient Simulation (ChIMES) MLFF model, which expresses the atomic interaction potentials in terms of linear combinations of many-body Chebyshev polynomials. We first show that our method can reduce our training set size by ∼70% while improving overall accuracy compared to more standard MD model fitting approaches. We then validate our optimum model against diverse hydrocarbon simulation data, including simple alkanes and systems with unsaturated carbon bonds, over a wide range of thermodynamic conditions. Finally, we use our ChIMES model to perform MD simulations of radiolytic damage with large-scale systems that help avoid system size effects. Overall, our approach yields an MD force field that retains most of the accuracy of the underlying quantum method while yielding many orders of improvement in computational efficiency. In conclusion, our efforts will have impact on future hydrocarbon polymer radiolysis studies, where the chemical details of the polymer–radiation interactions can have a strong effect on the resulting products observed in experiments.

Hydrocarbons↗
DOE OSTI2014

A staggered-grid finite-difference scheme optimized in the time–space domain for modeling scalar-wave propagation in geophysical problems

For modeling scalar-wave propagation in geophysical problems using finite-difference schemes, optimizing the coefficients of the finite-difference operators can reduce numerical dispersion. Most optimized finite-difference schemes for modeling seismic-wave propagation suppress only spatial but not temporal dispersion errors. We develop a novel optimized finite-difference scheme for numerical scalar-wave modeling to control dispersion errors not only in space but also in time. Our optimized scheme is based on a new stencil that contains a few more grid points than the standard stencil. We design an objective function for minimizing relative errors of phase velocities of waves propagating in all directions within a given range of wavenumbers. Dispersion analysis and numerical examples demonstrate that our optimized finite-difference scheme is computationally up to 2.5 times faster than the optimized schemes using the standard stencil to achieve the similar modeling accuracy for a given 2D or 3D problem. Compared with the high-order finite-difference scheme using the same new stencil, our optimized scheme reduces 50 percent of the computational cost to achieve the similar modeling accuracy. This new optimized finite-difference scheme is particularly useful for large-scale 3D scalar-wave modeling and inversion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Optimization of a superconducting adiabatic radio frequency neutron resonant spin flipper

The neutron radio frequency (RF) spin flipper is one of the key components to implement neutron resonance spin echo (NRSE) for high resolution measurements of dynamics in materials. Various types of RF flippers have been built and implemented including adiabatic and non-adiabatic RF flippers. For a pulsed neutron source, adiabatic RF flipper is more favorable due to its wide operation bandwidth. However, the optimization criteria of an adiabatic RF flipper is still lacking. Also they are difficult to be operated in bootstrap mode, which can be use to minimize the phase aberration and increase the achievable resolution. In the report by Dadisman et al. (2019), an adiabatic RF flipper has been built for the first time with superconducting technologies, which allows bootstrap mode. Here in this article, the method of optimizing the device and the detailed design will be discussed for each component of the device.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

Prospects for the methods of radionuclide production

In the present report, methods of radionuclide production for the nuclear-medicine purposes are described. In a budget approach, the application of low-energy accelerators is especially advantageous. Intense flux of bremsstrahlung at electron accelerators or high-current cyclotron beams of alpha particles must supply a great yield for many isotopes. The choice of a target material and of the projectile energy provides enough variation for concrete species formation. The innovating procedures are here proposed for optimizing of methods, for instance, application of the noble-gas target for production and transport of activities. The known and new variants of the “generator” scheme are discussed. Many isotopes are listed as promising in the context of the therapeutic and theragnostic applications. Among them are isotopes/isomers emitting soft radiation for the selective and careful body treatment, also the positron emitters for PET, and the halogen and alkali-metal species convenient for chemical separation.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2022

An optimized approach and inflation media for obtaining complimentary mass spectrometry-based omics data from human lung tissue

Human disease states are biomolecularly multifaceted and can span across phenotypic states, therefore it is important to understand diseases on all levels, across cell types, and within and across microanatomical tissue compartments. To obtain an accurate and representative view of the molecular landscape within human lungs, this fragile tissue must be inflated and embedded to maintain spatial fidelity of the location of molecules and minimize molecular degradation for molecular imaging experiments. Here, we evaluated agarose inflation and carboxymethyl cellulose embedding media and determined effective tissue preparation protocols for performing bulk and spatial mass spectrometry-based omics measurements. Mass spectrometry imaging methods were optimized to boost the number of annotatable molecules in agarose inflated lung samples. This optimized protocol permitted the observation of unique lipid distributions within several airway regions in the lung tissue block. Laser capture microdissection of these airway regions followed by high-resolution proteomic analysis allowed us to begin linking the lipidome with the proteome in a spatially resolved manner, where we observed proteins with high abundance specifically localized to the airway regions. We also compared our mass spectrometry results to lung tissue samples preserved using two other inflation/embedding media, but we identified several pitfalls with the sample preparation steps using this preservation method. Overall, we demonstrated the versatility of the inflation method, and we can start to reveal how the metabolome, lipidome, and proteome are connected spatially in human lungs and across disease states through a variety of different experiments.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

SU-E-T-551: Monitor Unit Optimization in Stereotactic Body Radiation Therapy for Stage I Lung Cancer

Purpose: The study aims to reduce the monitor units (MUs) in the stereotactic body radiation therapy (SBRT) treatment for lung cancer by adjusting the optimizing parameters. Methods: Fourteen patients suffered from stage I Non-Small Cell Lung Cancer (NSCLC) were enrolled. Three groups of parameters were adjusted to investigate their effects on MU numbers and organs at risk (OARs) sparing: (1) the upper objective of planning target volume (UOPTV); (2) strength setting in the MU constraining objective; (3) max MU setting in the MU constraining objective. Results: We found that the parameters in the optimizer influenced the MU numbers in a priority, strength and max MU dependent manner. MU numbers showed a decreasing trend with the UOPTV increasing. MU numbers with low, medium and high priority for the UOPTV were 428±54, 312±48 and 258±31 MU/Gy, respectively. High priority for UOPTV also spared the heart, cord and lung while maintaining comparable PTV coverage than the low and medium priority group. It was observed that MU numbers tended to decrease with the strength increasing and max MU setting decreasing. With maximum strength, the MU numbers reached its minimum while maintaining comparable or improved dose to the normal tissues. It was also found that the MU numbers continued to decline at 85% and 75% max MU setting but no longer to decrease at 50% and 25%. Combined with high priority for UOPTV and MU constraining objectives, the MU numbers can be decreased as low as 223±26 MU/Gy. Conclusion:: The priority of UOPTV, MU constraining objective in the optimizer impact on the MU numbers in SBRT treatment for lung cancer. Giving high priority to the UOPTV, setting the strength to maximum value and the max MU to 50% in the MU objective achieves the lowest MU numbers while maintaining comparable or improved OAR sparing.

62 RADIOLOGY AND NUCLEAR MEDICINE↗