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 991 records · Page 55

Perform Life Cycle Energy and GHG Emission Analysis, Select Candidate Refrigerant(s), Developing Heat Pump by Using Low-GWP Refrigerants and 5 mm Diameter Tubes

Reducing environmental impacts from the energy conversion systems has been an important research topic due to recent severe global climate changes. In residential buildings, the space heating and cooling are main energy consumers, which are relying on vapor compression based heat pumps. To mitigate the environmental impacts of residential systems, the heat pumps using low-GWP refrigerants and 5mm-tube heat exchangers are proposed to replace the conventional 9 mm or 7 mm tube heat exchangers. This research uses a simulation-based design method to optimize the refrigerant circuitries of a 5-ton low-GWP refrigerant heat pumps whose indoor and outdoor heat exchangers adopt 5 mm tubes. Heat pump systems using four commercially available low-GWP refrigerants (R454A, R454C, ARM20A and ARM20B) are optimized and compared with the baseline R410A system. The results indicate that directly replacing 9 mm tubes into 5mm tubes in the baseline systems induces heat pump performance degradation; by optimizing the 5 mm tube heat pumps with added tube rows, the system refrigerant charge is reduced up to 33% compared with the baseline; the optimized low-GWP 5 mm tube heat pumps offer the same or higher performance under both cooling and heating mode; that under cooling mode, the EER improvement is up to 10.7% and under heating mode the COP improvement is up to 3.4%; and a smart flow control device to maintain counter-flow circuitry pattern on both operating modes can improve the heating performance of 5 mm tube optimized system by up to 3.0%. In an effort to better understand the environmental impacts of these systems, life cycle climate performance (LCCP) evaluation method is used to assess the baseline and optimized systems. LCCP considers the direct and indirect emissions of the system over the course of its lifetime from manufacturing to disposal. Our analyses indicate that the optimized systems using 5mm tubes and low-GWP refrigerants can reduce life-cycle CO 2 emission significantly. And improving system efficiency is the most effective way to reduce the emission of low-GWP refrigerant systems. This report shows that using 5 mm diameter tubes in heat exchangers is a feasible approach to develop heat pumps with low-GWP refrigerants. Using the proposed design approach can yield higher system efficiency with lower environmental impact.

42 ENGINEERING↗

Stochastic Least-Squares Petrov--Galerkin Method for Parameterized Linear Systems

We consider the numerical solution of parameterized linear systems where the system matrix, the solution, and the right-hand side are parameterized by a set of uncertain input parameters. We explore spectral methods in which the solutions are approximated in a chosen finite-dimensional subspace. It has been shown that the stochastic Galerkin projection technique fails to minimize any measure of the solution error [A. Mugler and H.-J. Starkloff, ESAIM Math. Model. Numer. Anal., 47 (2013), pp. 1237--1263]. As a remedy for this, we propose a novel stochatic least-squares Petrov--Galerkin (LSPG) method. The proposed method is optimal in the sense that it produces the solution that minimizes a weighted $\ell^2$-norm of the residual over all solutions in a given finite-dimensional subspace. Moreover, the method can be adapted to minimize the solution error in different weighted $\ell^2$-norms by simply applying a weighting function within the least-squares formulation. In addition, a goal-oriented seminorm induced by an output quantity of interest can be minimized by defining a weighting function as a linear functional of the solution. We establish optimality and error bounds for the proposed method, and extensive numerical experiments show that the weighted LSPG method outperforms other spectral methods in minimizing corresponding target weighted norms.

97 MATHEMATICS AND COMPUTING↗

Bayesian optimization for inverse calibration of expensive computer models: A case study for Johnson-Cook model in machining

Inverse model calibration for identifying the constitutive model parameters can be computationally demanding for expensive-to-evaluate simulation models. Here, this paper presents a modified Bayesian optimization (BO) method, denoted as BO-bound, that incorporates theoretical bounds on the quantity of interest. A case study for the inverse calibration of the Johnson Cook (J-C) flow stress model parameters is presented using machining (cutting) force data. The results show fast calibration of the five J-C parameters within 25 simulations. In general, the BO-bound method is applicable for inverse calibration of any expensive simulation models as well as optimization problems with known bounds.

Bayesian optimization↗

SU-G-BRC-02: A Novel Multi-Criteria Optimization Approach to Generate Deliverable Intensity-Modulated Radiation Therapy (IMRT) Treatment Plans

Purpose: To present a novel multi-criteria optimization (MCO) solution approach that generates treatment plans with deliverable apertures using column generation. Methods: We demonstrate our method with 10 locally advanced head-and-neck cancer cases retrospectively. In our MCO formulation, we defined an objective function for each structure in the treatment volume. This resulted in 9 objective functions, including 3 distinct objectives for primary target volume, high-risk and low-risk target volumes, 5 objectives for each of the organs-at-risk (OARs) (two parotid glands, spinal cord, brain stem and oral cavity), and one for the non-target non-OAR normal tissue. Conditional value-at-risk (CVaR) constraints were utilized to ensure at least certain fraction of the target volumes receiving the prescription doses. To directly generate deliverable plans, column generation algorithm was embedded within our MCO approach for aperture shape generation. Final dose distributions for all plans were generated using a Monte Carlo kernel-superposition dose calculation. We compared the MCO plans with the clinical plans, which were created by clinicians. Results: At least 95% target coverage was achieved by both MCO plans and clinical plans. However, the average conformity indices of clinical plans and the MCO plans were 1.95 and 1.35, respectively (31% reduction, p<0.01). Compared to the conventional clinical plan, the proposed MCO method achieved average reductions in left parotid mean dose of 5% (p=0.06), right parotid mean dose of 18% (p<0.01), oral cavity mean dose of 21% (p=0.03), spinal cord maximum dose of 20% (p<0.01), brain stem maximum dose of 61% (p<0.01), and normal tissue maximum dose of 5% (p<0.01), respectively. Conclusion: We demonstrated that the proposed MCO method was able to obtain deliverable IMRT treatment plans while achieving significant improvements in dosimetric plan quality.

60 APPLIED LIFE SCIENCES↗

Effects of Ink Formulation on the Structure and Performance of PGM-Free Catalyst Layer in PEMFCs

Platinum group metal (PGM) catalysts are the major electrocatalysts for oxygen reduction reaction (ORR) in the polymer electrolyte membrane fuel cells (PEMFCs). The cost becomes unaffordable if the PEMFC is in massive application. The PGM-Free catalyst shows very promising activity in rotation disk electrode (RDE) testing. The half-wave potential could reach 0.91 V versus standard hydrogen electrode (SHE). However, in a membrane electrode assembly (MEA), the performance of PGM-Free catalysts is not good enough to replace the PGM catalysts. Since the PGM-free catalysts are so different from the PGM catalysts in terms of catalytic activity, stability, surface conditions, particle size, etc., the fabrication of PGM-Free catalyst MEA cannot simply copy the method of making PGM MEA. Here we proposed a novel method of fabricating PGM-Free catalyst MEA, so that the intrinsic catalyst activity from RDE can be translated into MEA performance. The method is based on the catalyst coated membrane (CCM) method using optimized ionomer to carbon (I/C) ratio and solvent mixture of catalyst ink. Using this method, the PGM-free catalyst MEA achieved the current density 44.9 mA cm -2 at 0.9 V iR-free in H 2 /O 2 and 150 mA cm -2 at 0.8 V in H 2 /air, which surpassed the performance targets of US Department of Energy (DOE)for PGM-Free catalyst MEA. The property (solvent composition, dispersion of catalyst and ionomer in an ink), structure (pore structure) and the MEA performance have been characterized using, mercury intrusion porosimetry (MIP), MEA testing. A property-structure-performance relationship has been established.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fast yaw optimization for wind plant wake steering using Boolean yaw angles

Abstract. In wind plants, turbines can be yawed into the wind to steer their wakes away from downstream turbines and achieve an overall increase in plant power. Mathematical optimization is typically used to determine the best yaw angles at which to operate the turbines in a plant. In this paper, we present a new heuristic to rapidly determine the yaw angles in a wind plant. In this method, we define the turbine yaw angles as Boolean – either yawed at a predefined angle or nonyawed – as opposed to the typical methods of defining yaw angles as continuous or with fine discretizations. We then optimize which turbines should be yawed with an algorithm that sweeps through the turbines from the most upstream to the most downstream. We demonstrate that our new Boolean optimization method can find turbine yaw angles that perform well compared to a traditionally used gradient-based optimizer for which the yaw angles are defined as continuous. There is less than 0.6 % difference in the optimized power between the two optimization methods for randomly placed turbine layouts and less than a 0.6 % difference in the optimal annual energy production between the two optimization methods for a real wind farm. Additionally, we show that our new method is much more computationally efficient than the traditional method. For plants with nonzero optimal yaw angles, our new method is generally able to solve for the turbine yaw angles 50–150 times faster, and in some extreme cases up to 500 times faster, than the traditional method.

17 WIND ENERGY↗

Optimization of the Jastrow factor using the random-phase approximation and a similarity-transformed Hamiltonian: Application to band-structure calculation for some semiconductors and insulators

Based on the random-phase approximation and the transcorrelated (TC) method, we optimize the Jastrow factor together with one-electron orbitals in the Slater determinant in the correlated wave function with a new scheme for periodic systems. The TC method is one of the promising wave function theories for first-principles electronic structure calculation, where the many-body wave function is approximated as a product of a Slater determinant and a Jastrow factor, and the Hamiltonian is similarity-transformed by the Jastrow factor. Using this similarity-transformed Hamiltonian, we can optimize the one-electron orbitals without evaluating 3N-dimensional integrations for the N-electron system. In contrast, optimization of the Jastrow factor within the framework of the TC method is computationally much more expensive and has not been performed for solid-state calculations before. In this study, we also benefit from the similarity-transformation in optimizing the Jastrow factor. Our optimization scheme is tested in applications to some solids from narrow-gap semiconductors to wide-gap insulators, and it is verified that the band gap of a wide-gap insulator and the lattice constants of some solids are improved by this optimization with reasonable computational cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Role of structural anisotropy in geometric birefringence of high-birefringence index-guiding PCFs: analysis by moment of inertia

We have studied the role of structural anisotropy on the geometric birefringence of four typical high-birefringence (HB) index-guiding photonic crystal fibers (PCFs) with the same air-silica structure and circular, ellipse, and rectangle air holes, respectively. The normalized difference between the moment of inertia ΔI in x and y direction of the PCF structure is utilized to quantificationally describe the structural anisotropy. We demonstrate that the phase birefringence B increases monotonously when the normalized ΔI rises in different cases of PCFs structures. In order to obtain high birefringence B, it can be an effective method to optimize the parameters of the structure to increase the normalized ΔI. This work can be not only a reference of founding the essence of geometric birefringence and proposing the optimal design of HB-IG-PCFs, but also an enlightenment for bringing the different method (compared to the finite element) into the investigation of PCFs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

An optimal guarding scheme for thermal conductivity measurement using a guarded cut-bar technique, part 1 experimental study

In the guarded cut-bar technique, a guard surrounding the measured sample and reference (meter) bars is temperature controlled to carefully regulate heat losses from the sample and reference bars. Guarding is typically carried out by matching the temperature profiles between the guard and the test stack of sample and meter bars. Problems arise in matching the profiles, especially when the thermal conductivitiesof the meter bars and of the sample differ, as is usually the case. In a previous numerical study, the applied guarding condition (guard temperature profile) was found to be an important factor in measurement accuracy. Different from the linear-matched or isothermal schemes recommended in literature, the optimal guarding condition is dependent on the system geometry and thermal conductivity ratio of sample to meter bar. To validate the numerical results, an experimental study was performed to investigate the resulting error under different guarding conditions using stainless steel 304 as both the sample and meter bars. The optimal guarding condition was further verified on a certified reference material, pyroceram 9606, and 99.95% pure iron whose thermal conductivities are much smaller and much larger, respectively, than that of the stainless steel meter bars. Additionally, measurements are performed using three different inert gases to show the effect of the insulation effective thermal conductivity on measurement error, revealing low conductivity, argon gas, gives the lowest error sensitivity when deviating from the optimal condition. The result of this study provides a general guideline for the specific measurement method and for methods requiring optimal guarding or insulation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Co-design optimization of combined heat and power-based microgrids

With the emergent need for clean and reliable energy resources, hybrid energy systems, such as the microgrid, are widely adopted in the United States. A microgrid can consist of various distributed energy resources, for instance, combined heat and power (CHP) systems. Here, the CHP module is a distributed cogeneration technology that produces electricity and recaptures heat generated as a by-product. It is an energy-efficient technology converting heat that would otherwise be wasted to valuable thermal energy. For an optimal system configuration, this study develops a novel co-design optimization framework for CHP-based cogeneration microgrids. The framework provides the stakeholder with a method to optimize investments and attain resilient operations. The proposed co-design framework has a mixed integer programming (MIP) model that outputs decisions for both plant designs and operating controls. The microgrid considered in this study contains six components: the CHP, boiler, heat recovery unit, thermal storage system, power storage system, and photovoltaic plant. After solving the MIP model, the optimal design parameters of each component can be found to minimize the total installation cost of all components in the microgrid. Furthermore, the online costs from energy production, operation, maintenance, machine startup, and disruption-induced unsatisfied loads are minimized by solving the optimal control decisions for operations. Case studies based on designing a CHP-based microgrid with empirical data are conducted. Moreover, we consider both nominal and disruptive operational scenarios to validate the performance of the proposed co-design framework in terms of a cost-effective, resilient system.

42 ENGINEERING↗

WE-AB-207A-09: Optimization of the Design of a Moving Blocker for Cone-Beam CT Scatter Correction: Experimental Evaluation

Purpose: A moving blocker based strategy has shown promising results for scatter correction in cone-beam computed tomography (CBCT). Different geometry designs and moving speeds of the blocker affect its performance in image reconstruction accuracy. The goal of this work is to optimize the geometric design and moving speed of the moving blocker system through experimental evaluations. Methods: An Elekta Synergy XVI system and an anthropomorphic pelvis phantom CIRS 801-P were used for our experiment. A blocker consisting of lead strips was inserted between the x-ray source and the phantom moving back and forth along rotation axis to measure the scatter signal. Accoriding to our Monte Carlo simulation results, three blockers were used, which have the same lead strip width 3.2mm and different gap between neighboring lead strips, 3.2, 6.4 and 9.6mm. For each blocker, three moving speeds were evaluated, 10, 20 and 30 pixels per projection (on the detector plane). Scatter signal in the unblocked region was estimated by cubic B-spline based interpolation from the blocked region. CBCT image was reconstructed by a total variation (TV) based algebraic iterative reconstruction (ART) algorithm from the partially blocked projection data. Reconstruction accuracy in each condition is quantified as CT number error of region of interest (ROI) by comparing to a CBCT reconstructed image from analytically simulated unblocked and scatter free projection data. Results: Highest reconstruction accuracy is achieved when the blocker width is 3.2 mm, the gap between neighboring lead strips is 9.6 mm and the moving speed is 20 pixels per projection. RMSE of the CT number of ROIs can be reduced from 436 to 27. Conclusions: Image reconstruction accuracy is greatly affected by the geometry design of the blocker. The moving speed does not have a very strong effect on reconstruction result if it is over 20 pixels per projection.

60 APPLIED LIFE SCIENCES↗

Cleaner Energy for Cars

Representing the Center for Catalytic Hydrocarbon Functionalization (CCHF), this document is one of the entries in the Ten Hundred and One Word Challenge. As part of the challenge, the 46 Energy Frontier Research Centers were invited to represent their science in images, cartoons, photos, words and original paintings, but any descriptions or words could only use the 1000 most commonly used words in the English language, with the addition of one word important to each of the EFRCs and the mission of DOE: energy. The mission of CCHF is to develop, validate, and optimize new methods to rearrange the bonds of hydrocarbons, implement enzymatic strategies into synthetic systems, and design optimal environments for catalysts that can be used to reversibly functionalize hydrocarbons, especially for more efficient use of natural gas including low temperature conversion to liquid fuels.

08 HYDROGEN↗

Characterization and Clearance of m/s Sigyn - 16412

M/s Sigyn was a 2044 ton ship that transported contaminated material between 1982 and 2013. Most of the transports were with spent nuclear fuel between Swedish nuclear power plants and the Swedish fuel repository, Clab. In addition to spent fuel, m/s Sigyn also transported big contaminated components, from nuclear power plants in Europe, to Studsvik for analysis or decontamination and melting. She also transported components from Swedish nuclear power plants to final disposal of short lived waste, SFR. Even though all transported items have been cleaned, checked and wrapped with plastics, some activity ended up as contamination onboard. In December 2013 it was decided that m/s Sigyn was subject for decommissioning due to her age. This Paper describes the method in which the activity onboard m/s Sigyn was characterized and clearance levels were calculated. This method for radiological characterization and clearance has recently been developed in Sweden, the method is efficient in that it does not need a big pre characterization survey. Only the nuclide distribution and history of the object, is needed to start the final survey. In the future there will be further development of the method, making it possible to reduce the number of measurements even more. This will further optimize the method and reduce the cost for final surveys and clearance. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Machine Learning for Fairness-Aware Load Shedding: A Real-Time Solution via Identifying Binding Constraints: Preprint

Timely and effective load shedding in power systems is critical for maintaining supply-demand balance and preventing cascading blackouts. To eliminate load shedding bias against specific regions in the system, optimization-based methods are uniquely positioned to help balance between economic and fairness considerations. However, the resulting optimization problem involves complex constraints, which can be time-consuming to solve and thus cannot meet the real-time requirements of load shedding. To tackle this challenge, in this paper we present an efficient machine learning algorithm to enable millisecond-level computation for the optimization-based load shedding problem. Numerical studies on both a 3-bus toy example and a realistic RTS-GMLC system have demonstrated the validity and efficiency of the proposed algorithm for delivering fairness-aware and real-time load shedding decisions.

97 MATHEMATICS AND COMPUTING↗

Automatic Detection of Clear-Sky Periods Using Ground and Satellite Based Solar Resource Data

Solar resource availability and variability are important aspects of monitoring performance of photovoltaic installations. For example, recent degradation studies have highlighted the importance of considering cloud cover when calculation degradation rates. With this in mind, we present a method for optimizing clear sky detection algorithms given only modeled clear sky irradiance and ground-measured irradiance values. This method is tested on global horizontal irradiance (GHI) data from ground collectors at six sites across the US and was trained against clear sky classifications determined from satellite data. Thirty models were optimized on each individual site at GHI data frequencies of 1, 5, 10, 15, and 30 minutes. The models had an average F 0.5 score of 0.945 +/- .021 on a holdout test set. In comparison, the un-optimized clear sky detection algorithm produced F 0.5 score that averaged to 0.707 +/- 0.187.

14 SOLAR ENERGY↗

Fast leaf-fitting with generalized underdose/overdose constraints for real-time MLC tracking

Purpose: Real-time multileaf collimator (MLC) tracking is a promising approach to the management of intrafractional tumor motion during thoracic and abdominal radiotherapy. MLC tracking is typically performed in two steps: transforming a planned MLC aperture in response to patient motion and refitting the leaves to the newly generated aperture. One of the challenges of this approach is the inability to faithfully reproduce the desired motion-adapted aperture. This work presents an optimization-based framework with which to solve this leaf-fitting problem in real-time. Methods: This optimization framework is designed to facilitate the determination of leaf positions in real-time while accounting for the trade-off between coverage of the PTV and avoidance of organs at risk (OARs). Derived within this framework, an algorithm is presented that can account for general linear transformations of the planned MLC aperture, particularly 3D translations and in-plane rotations. This algorithm, together with algorithms presented in Sawant et al. [“Management of three-dimensional intrafraction motion through real-time DMLC tracking,” Med. Phys. 35, 2050–2061 (2008)] and Ruan and Keall [Presented at the 2011 IEEE Power Engineering and Automation Conference (PEAM) (2011) (unpublished)], was applied to apertures derived from eight lung intensity modulated radiotherapy plans subjected to six-degree-of-freedom motion traces acquired from lung cancer patients using the kilovoltage intrafraction monitoring system developed at the University of Sydney. A quality-of-fit metric was defined, and each algorithm was evaluated in terms of quality-of-fit and computation time. Results: This algorithm is shown to perform leaf-fittings of apertures, each with 80 leaf pairs, in 0.226 ms on average as compared to 0.082 and 64.2 ms for the algorithms of Sawant et al., Ruan, and Keall, respectively. The algorithm shows approximately 12% improvement in quality-of-fit over the Sawant et al. approach, while performing comparably to Ruan and Keall. Conclusions: This work improves upon the quality of the Sawant et al. approach, but does so without sacrificing run-time performance. In addition, using this framework allows for complex leaf-fitting strategies that can be used to account for PTV/OAR trade-off during real-time MLC tracking.

61 RADIATION PROTECTION AND DOSIMETRY↗

sMF-BO-2CoGP: A sequential multi-fidelity constrained Bayesian optimization framework for design applications

Bayesian optimization (BO) is an effective surrogate-based method that has been widely used to optimize simulation-based applications. While the traditional Bayesian optimization approach only applies to single-fidelity models, many realistic applications provide multiple levels of fidelity with various levels of computational complexity and predictive capability. In this work, we propose a multi-fidelity Bayesian optimization method for design applications with both known and unknown constraints. The proposed framework, called sMF-BO-2CoGP, is built on a multi-level CoKriging method to predict the objective function. An external binary classifier, which we approximate using a separate CoKriging model, is used to distinguish between feasible and infeasible regions. Finally, the sMF-BO-2CoGP method is demonstrated using a series of analytical examples and a flip-chip application for design optimization to minimize the deformation due to warping under thermal loading conditions.

42 ENGINEERING↗