Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “code optimization”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

iGeoT v1.0

iGeoT v1.0 automates into a stand-alone computer program the multicomponent chemical geothermometry code GeoT v2.1 and the numerical optimization engine of iTOUGH2. iGeoT allows for optimizations of GeoT runs using multiple water chemical analyses. The underlying geothermometry method is that previously developed by Reed and Spycher (1984) using the computed saturation indices of multiple minerals. GeoT automatically reconstitutes deep fluid compositions and estimates reservoir temperature from statistical evaluation of computed mineral saturation indices. The output include estimated temperatures following various statistical methods and their range of uncertainty. The computer program also outputs temperatures predicted using classical geothermometers. Input parameters include water composition and, optionally, gas composition, fraction of gas discharge, dilution/evaporation factor (in case of boiling or mixing with dilute waters), and end-member water composition (in case of mixing with other non-dilute waters). The dilution/evaporation factor, fraction of gas discharge, and concentration of various aqueous and gas species can be automatically estimated by numerical optimization. This program was designed for geothermal applications.

Spycher, Nicolas↗
DOE OSTI2020

An advanced membrane-based HMX model and its application in utilizing sensible and latent heat removed in space cooling

This paper presents a comprehensive simulation model for advanced membrane-based heat and mass exchanger (HMX) components which account for three separate fluid streams with complex flow patterns. The model was implemented with various liquid desiccants and their properties (including LiCl, CaCl and ionic liquid (IL)). Preliminary model validation is presented using data from a laboratory prototype absorber. The paper also presents implementation of the model into a sorption system modeling platform and analyzes potential benefit analysis of the HMX component in utilizing the sensible and latent heat removed in space cooling to heat water. The developed model and codes are expected to enable detailed optimization studies and identify novel configurations for the three-fluid HMX maximum performance and cost effectiveness.

Gao, Zhiming↗
DOE OSTI2022

Core Design of the Holos-Quad Microreactor

The Holos-Quad micro-reactor concept, developed by HolosGen LLC, is equipped with a 22 MWt (Mega-Watt thermal) core and an integral power conversion system converting the core thermal energy into approximately 10 MWe (Mega-Watt electric). This design can be configured to support a wide range of applications. It is a very innovative high-temperature gas-cooled reactor concept using TRI-structural ISOtropic particle fuel (TRISO) distributed in graphite hexagonal blocks, cooled with helium in a direct Brayton cycle independently executed by four Subcritical Power Modules (SPMs) fitted into a hardened 40-foot container whose dimensions are in compliance with ISO shipping containers requirements. In FY2019 HolosGen LLC was awarded by the Department of Energy Advanced Research Project Agency-Energy (DOE ARPA-E) under the MEITNER funding program. As part of the MEITNER award, the Argonne National Laboratory (ANL) contributed expertise through two specialized teams: The “Design Team” and the “Resource Team”. The Design Team was dedicated to validate feasibility of the Holos-Quad core and to optimize its core design through neutronics analyses. The Resource Team was dedicated to feasibility verification via high-fidelity codes of Holos-Quad thermal-hydraulic, heat transfer, shielding, and structural aspects. This report summarizes the activities conducted by ANL Design Team. A rigorous design approach based on multi-criteria optimization and code-to-code comparison involving stochastic and high-fidelity deterministic solutions was developed and employed at several evolutionary stages of the Holos-Quad design. Several generations of the Holos-Quad core were designed within this project before converging to the current full-scale Gen 2+ design that is detailed in this report. Figure EA-1 illustrates a cross-sectional view of Gen 2+ Holos-Quad core configuration, and Figure EA-2 provides a simplified perspective view of 1-of-4 SPMs. The Holos-Quad uses four thermal-hydraulically independent SPMs locked into stationary positions during power operation, surrounded by BeO reflector and structural component fully comprised within the dimensional constraints represented by traditional ISO containers. One of the benefits of this approach is to enable transportation of each SPM promptly after irradiation in shielded containers. The core is designed to operate for approximately 8 full-power years while the reactivity controls and power conversion system enable load-following operations. The reactivity controls are represented by independent, diversified, and redundant reactivity control systems based on control drums and redundant sets of shutdown rods. The high-fidelity simulation tools were used to assess detailed power and flux distributions of the three-dimensional full-core or quarter-core of the Gen 2+ configuration. Single-physics and multi-physics simulations of the neutronics code PROTEUS and the thermal-hydraulic code System Analysis Module (SAM) were performed to analyze the Holos design configurations with detailed high-fidelity solutions. The design work performed confirmed feasibility of the Holos-Quad concept, provided realistic design description for detailed design of the operational system, and identified several core design improvements to be further considered for future reactor development activities.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

An OpenMP GPU-offload implementation of a non-equilibrium solidification cellular automata model for additive manufacturing

Here, in this paper, performance strategies on GPU-based HPC platforms of a cellular automata (CA) simulation code for non-equilibrium solidification, including nucleation, grain growth, solute partitioning and transport for the metal additive manufacturing (AM) process are investigated using OpenMP 4.5. To accurately report the speed-up for multicore CPUs and GPUs, a rigorous performance analysis employed optimizations appropriate for both CPU-only code (baseline) and GPU offload codes for an isothermal test problem. The performance results on Summit at the Oak Ridge Leadership Computing Facility indicate that using a precomputed list of interface cells significantly decreased the wall-clock time on GPUs. The speedup due to GPU acceleration was evaluated for a full Summit node and measured to be 1.8X when comparing a 6 MPI tasks run with 6 GPUs versus 36 MPI tasks on the CPU only. That speed-up was found to be 7.9X when comparing 6 MPI tasks with 6 GPUs versus the 6 MPI tasks running on the CPU only. Performance measurements showed that system total time is almost constant for runs with more than 96 MPI tasks (or GPUs), indicating that the GPU-accelerated code showed an excellent weak scaling performance. Finally, a rapid directional solidification problem was considered to demonstrate the CA code capability on Summit. It was found that a mesh size of at least 0.05 μm is recommended for the AM-like simulations in order to obtain accurate elongated grain microstructure and elongated subgrain features, which are in qualitative good agreement with experimental data. The results presented in this study indicate that the performance strategies on GPU-based HPC platforms for the CA code are appropriate for novel HPC exascale platforms.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Land-based wind turbines with flexible rail-transportable blades – Part 2: 3D finite element design optimization of the rotor blades

Abstract. Increasing growth in land-based wind turbine blades to enable higher machine capacities and capacity factors is creating challenges in design, manufacturing, logistics, and operation. Enabling further blade growth will require technology innovation. An emerging solution to overcome logistics constraints is to segment the blades spanwise and chordwise, which is effective, but the additional field-assembled joints result in added mass and loads, as well as increased reliability concerns in operation. An alternative to this methodology is to design slender flexible blades that can be shipped on rail lines by flexing during transport. However, the increased flexibility is challenging to accommodate with a typical glass-fiber, upwind design. In a two-part paper series, several design options are evaluated to enable slender flexible blades: downwind machines, optimized carbon fiber, and active aerodynamic controls. Part 1 presents the system-level optimization of the rotor variants as compared to conventional and segmented baselines, with a low-fidelity representation of the blades. The present work, Part 2, supplements the system-level optimization in Part 1 with high-fidelity blade structural optimization to ensure that the designs are at feasible optima with respect to material strength and fatigue limits, as well as global stability and structural dynamics constraints. To accommodate the requirements of the design process, a new version of the Numerical Manufacturing And Design (NuMAD) code has been developed and released. The code now supports laminate-level blade optimization and an interface to the International Energy Agency Wind Task 37 blade ontology. Transporting long, flexible blades via controlled flapwise bending is found to be a viable approach for blades of up to 100 m. The results confirm that blade mass can be substantially reduced by going either to a downwind design or to a highly coned and tilted upwind design. A discussion of active and inactive constraints consisting of material rupture, fatigue damage, buckling, deflection, and resonant frequencies is presented. An analysis of driving load cases revealed that the downwind designs are dominated by loads from sudden, abrupt events like gusts rather than fatigue. Finally, an analysis of carbon fiber spar caps for downwind machines finds that, compared to typical carbon fibers, the use of a new heavy-tow carbon fiber in the spar caps is found to yield between 9 % and 13 % cost savings.

17 WIND ENERGY↗
DOE OSTI2014

Rapid Optimization Library

ROL provides interfaces to and implementations of algorithms for gradient-based unconstrained and constrained optimization. ROL can be used to optimize the response of any client simulation code that evaluates scalar-valued response functions. If the client code can provide gradient information for the response function, ROL will take advantage of it, resulting in faster runtimes. ROL's interfaces are matrix-free, in other words ROL only uses evaluations of scalar-valued and vector-valued functions. ROL can be used to solve optimal design problems and inverse problems based on a variety of simulation software.

Kouri, DrewKouri↗
DOE OSTI2014

Rapid Optimization Library

ROL provides interfaces to and implementations of algorithms for gradient-based unconstrained and constrained optimization. ROL can be used to optimize the response of any client simulation code that evaluates scalar-valued response functions. If the client code can provide gradient information for the response function, ROL will take advantage of it, resulting in faster runtimes. ROL's interfaces are matrix-free, in other words ROL only uses evaluations of scalar-valued and vector-valued functions. ROL can be used to solve optimal design problems and inverse problems based on a variety of simulation software.

Denis Rldzal, Drew Kouri↗
DOE OSTI2018

Finite element analysis of compliant cladding and base metal systems

This report describes a roadmap towards a design method for high temperature nuclear cladded components that does not require long term testing of clad materials. Molten Salt Reactors (MSRs) will operate using corrosive fuel or coolant salts. The existing Class A metallic materials qualified for ASME Boiler and Pressure Vessel Code Section III, Division 5 design are not optimal for corrosion resistance when exposed to MSR salts. Corrosion could then severely limit the service life of MSRs. One option to alleviate this could be to Code qualify new corrosion resistant materials. However, this is a lengthy and expensive process and MSR developers have a near term need for a design solution that enables reactor designs with long service lives. An alternative would be to use cladded components – overlay some existing Class A material with a thin layer of some corrosion resistant material. This report describes the outlines of a design strategy that would enable this approach. The strategy relies on approximate design analysis methods for certain types of clad/base systems that avoids requiring long-term mechanical properties for the clad materials. This would allow the use of cladded components in the near term, without requiring long testing programs. However, it will limit the types of clad materials that can be used. This report describes the design framework, demonstrates the effectiveness of the approximate analysis methods, and confirms the approach using full scale finite element analysis of prototypical cladded high temperature components. It also describes a roadmap towards a complete design method and surveys available clad materials and bonding methods for suitably, considering only mechanical and not environmental effects.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

Conceptual moderator studies for the Spallation Neutron Source short-pulse second target station

Candidate moderator configurations for a short-pulse second target station (STS) at the Oak Ridge National Laboratory Spallation Neutron Source (SNS) have been identified using a global optimizer framework built around the MCNPX particle transport code. Neutron brightness metrics were selected as the figure-of-merit. We assumed that STS would use one out of six proton pulses produced by an SNS accelerator upgraded to operate at 1.3 GeV proton energy, 2.8 MW power and 60 Hz repetition rate. The simulations indicate that the peak brightness can be increased by a factor of 5 and 2.5 on a per proton pulse basis compared to the SNS first target station for both coupled and decoupled para-hydrogen moderators, respectively. Additional increases by factors of 3 and 2 were demonstrated for coupled and decoupled moderators, respectively, by reducing the area of neutron emission from 100 × 100 mm{sup 2} to 20 × 20 mm{sup 2}. This increase in brightness has the potential to translate to an increase of beam intensity at the instruments’ sample positions even though the total neutron emission of the smaller moderator is less than that of the larger. This is especially true for instruments with small samples (beam dimensions). The increased fluxes in the STS moderators come at accelerated poison and de-coupler burnout and higher radiation-induced material damage rates per unit power, which overall translate into lower moderator lifetimes. A first effort was undertaken to group decoupled moderators into a cluster collectively positioning them at the peak neutron production zone in the target and having a three-port neutron emission scheme that complements that of a cylindrical coupled moderator.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2013

Analysis and optimization of gyrokinetic toroidal simulations on homogenous and heterogenous platforms

The Gyrokinetic Toroidal Code (GTC) uses the particle-in-cell method to efficiently simulate plasma microturbulence. This paper presents novel analysis and optimization techniques to enhance the performance of GTC on large-scale machines. We introduce cell access analysis to better manage locality vs. synchronization tradeoffs on CPU and GPU-based architectures. Finally, our optimized hybrid parallel implementation of GTC uses MPI, OpenMP, and NVIDIA CUDA, achieves up to a 2× speedup over the reference Fortran version on multiple parallel systems, and scales efficiently to tens of thousands of cores.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Wind Farm Layout Optimization with Loads Considerations

The objective of this paper is to improve the annual energy production of a wind farm by optimizing the layout of a wind farm, while considering fatigue loads on turbines. In this paper, the loads are estimated using the edgewise bending moment computed using CCBlade, a steady-state blade element momentum code. The edgewise bending moment is then used to calculate fatigue damage using Miner's rule. The fatigue damage is used to constrain the layout optimization problem. We show that our method can predict blade root damage with similar trends to damage calculated with other methods, such as a complex, computationally expensive large-eddy simulation and unsteady aeroelastic code. We also optimize wind farm layouts and show that, for a simple problem with two wind directions and ten turbines, the fatigue damage can be constrained without sacrificing wind farm power production.

17 WIND ENERGY↗
DOE OSTI2014

Optimized aperiodic multilayer structures for use as narrow-angular absorbers

In this paper, we investigate aperiodic multilayer structures for use as narrow-angular absorbers. The layer thicknesses and materials are optimized using a genetic global optimization algorithm coupled to a transfer matrix code to maximize the angular selectivity in the absorptance at a single or multiple wavelengths. We first consider structures composed of alternating layers of tungsten and silicon or silica, and find that it is not possible to achieve angular selectivity in the absorptance with such structures. We next consider structures composed of alternating layers of silicon and silica, and show that when optimized they exhibit high angular selectivity in absorptance. In addition, as the angular selectivity in absorptance increases, the wavelength range of high angular selectivity also decreases. Optimizing the material composition of the multilayer structures, in addition to optimizing the layer thicknesses, leads to marginal improvement in angular selectivity. Finally, we show that by optimizing the absorptance of the multilayer structures at multiple wavelengths, we can obtain structures exhibiting almost perfect absorptance at normal incidence and narrow angular width in absorptance at these wavelengths. Similar to the structures optimized at a single wavelength, the wavelength range of high angularly selective absorptance is narrow.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Numerical Analysis of 2-D and 3-D MHD Flows Relevant to Fusion Applications

Here, the analysis of many fusion applications such as liquid-metal blankets requires application of computational fluid dynamics (CFD) methods for electrically conductive liquids in geometrically complex regions and in the presence of a strong magnetic field. A current state of the art general purpose CFD code allows modeling of the flow in complex geometric regions, with simultaneous conjugated heat transfer analysis in liquid and surrounding solid parts. Together with a magnetohydrodynamics (MHD) capability, the general purpose CFD code will be a valuable tool for the design and optimization of fusion devices. This paper describes an introduction of MHD capability into the general purpose CFD code CFX, part of the ANSYS Workbench. The code was adapted for MHD problems using a magnetic induction approach. CFX allows introduction of user-defined variables using transport or Poisson equations. For MHD adaptation of the code three additional transport equations were introduced for the components of the magnetic field, in addition to the Poisson equation for electric potential. The Lorentz force is included in the momentum transport equation as a source term. Fusion applications usually involve very strong magnetic fields, with values of the Hartmann number of up to tens of thousands. In this situation a system of MHD equations become very rigid with very large source terms and very strong variable gradients. To increase system robustness, special measures were introduced during the iterative convergence process, such as linearization using source coefficient for momentum equations. The MHD implementation in general purpose CFD code was tested against benchmarks, specifically selected for liquid-metal blanket applications. Results of numerical simulations using present implementation closely match analytical solutions for a Hartmann number of up to 1500 for a 2-D laminar flow in the duct of square cross section, with conducting and nonconducting walls. Results for a 3-D test case are also included.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2011

Commnity Petascale Project for Accelerator Science And Simulation: Advancing Computational Science for Future Accelerators And Accelerator Technologies

The design and performance optimization of particle accelerators are essential for the success of the DOE scientific program in the next decade. Particle accelerators are very complex systems whose accurate description involves a large number of degrees of freedom and requires the inclusion of many physics processes. Building on the success of the SciDAC-1 Accelerator Science and Technology project, the SciDAC-2 Community Petascale Project for Accelerator Science and Simulation (ComPASS) is developing a comprehensive set of interoperable components for beam dynamics, electromagnetics, electron cooling, and laser/plasma acceleration modelling. ComPASS is providing accelerator scientists the tools required to enable the necessary accelerator simulation paradigm shift from high-fidelity single physics process modeling (covered under SciDAC1) to high-fidelity multiphysics modeling. Our computational frameworks have been used to model the behavior of a large number of accelerators and accelerator R&D experiments, assisting both their design and performance optimization. As parallel computational applications, the ComPASS codes have been shown to make effective use of thousands of processors.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

A modular GUI-based program for genetic algorithm-based feedback-assisted wavefront shaping

Abstract We have developed a modular graphical user interface (GUI)-based program for use in genetic algorithm-based feedback-assisted wavefront shaping. The program uses a class-based structure to separate out the universal modules (e.g. GUI, multithreading, optimization algorithms) and hardware-specific modules (e.g. code for different SLMs and cameras). This modular design makes the program easily adaptable to a wide range of lab equipment, while providing easy access to a GUI, multithreading, and three optimization algorithms (phase-stepping, simple genetic, and microgenetic).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Optimization of a buncher for High Energy H- Beam Transport [Slides]

Accelerators and beam transport systems typically contain numerous beamline devices, diagnostics, system study, and upgrade. Routinely TRACE, MADX and TRANSPORT codes are used for the beam envelope prediction and optimization. Unknown problems of accelerators and beam-transport are addressed based on R&D, such as LFB voltage variation effect study. A nominal diagnostic have used in non-traditional beam transport path to address scientific question (such as use of a BPM in the air). High energy beam operation is a methodical process and routine R&D is required to maintain high standards to deliver beam.

43 PARTICLE ACCELERATORS↗
DOE OSTI2011

Method for Determining Optimal Residential Energy Efficiency Retrofit Packages

Businesses, government agencies, consumers, policy makers, and utilities currently have limited access to occupant-, building-, and location-specific recommendations for optimal energy retrofit packages, as defined by estimated costs and energy savings. This report describes an analysis method for determining optimal residential energy efficiency retrofit packages and, as an illustrative example,applies the analysis method to a 1960s-era home in eight U.S. cities covering a range of International Energy Conservation Code (IECC) climate regions. The method uses an optimization scheme that considers average energy use (determined from building energy simulations) and equivalent annual cost to recommend optimal retrofit packages specific to the building, occupants, and location. Energysavings and incremental costs are calculated relative to a minimum upgrade reference scenario, which accounts for efficiency upgrades that would occur in the absence of a retrofit because of equipment wear-out and replacement with current minimum standards.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗