Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “physical 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 271 records · Page 15

Increasing the hardness of posiform planting using random QUBOs for programmable quantum annealer benchmarking

Posiform planting is a method for constructing QUBO instances with a unique planted solution that can be tailored to arbitrary connectivity graphs. In this study we investigate making posiform planted QUBOs computationally harder by fusing many smaller random Ising models, whose global minimum is computed classically, with posiform planted QUBOs. The unique ground state of the resulting QUBO is the concatenation of (exactly one of) the ground states of each smaller problem. Our method generates QUBO instances that have a unique solution, are native to the hardware graph, and have tunable computational hardness. We use our QUBOs to benchmark three D-Wave quantum annealing processors (with 563–5627 qubits), and compare them against simulated annealing and Gurobi. Surprisingly, we find that the D-Wave ground state sampling success rate is not dependent on the glued random QUBO size, and that some QUBO classes are solved at high success rates at short annealing times on the Zephyr processors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Stiff neural ordinary differential equations

Neural Ordinary Differential Equations (ODEs) are a promising approach to learn dynamical models from time-series data in science and engineering applications. This work aims at learning neural ODEs for stiff systems, which are usually raised from chemical kinetic modeling in chemical and biological systems. We first show the challenges of learning neural ODEs in the classical stiff ODE systems of Robertson’s problem and propose techniques to mitigate the challenges associated with scale separations in stiff systems. We then present successful demonstrations in stiff systems of Robertson’s problem and an air pollution problem. The demonstrations show that the usage of deep networks with rectified activations, proper scaling of the network outputs as well as loss functions, and stabilized gradient calculations are the key techniques enabling the learning of stiff neural ODEs. The success of learning stiff neural ODEs opens up possibilities of using neural ODEs in applications with widely varying time-scales, such as chemical dynamics in energy conversion, environmental engineering, and life sciences.

97 MATHEMATICS AND COMPUTING↗

Fuel Behavior Implications of Reactor Design Choices in Pressurized Water SMRs

Small pressurized water reactors (PWRs) can feature boron free operation, natural circulation mode, reduced height assemblies and/or long refueling cycles. This paper attempts to explore core design optimization for each of these design evolutions. In consequence, five core design layouts are developed incorporating boron free operation with continuous control rods insertion, natural circulation with low burnup/low power density design, natural circulation with high burnup/low power density design, forced circulation with standard core power density design, and forced circulation with high power density design. These cores’ performance is compared to a standard 4-loop PWR. The design process aims to improve the fuel cycle cost under safety constraints through core design optimization using CASMO4E/SIMULATE3 reactor physics codes and FRAPCON4.1 fuel performance assessment tool. Core modeling assumes standard 17x17 PWR fuel assemblies loaded with low enriched uranium (LEU) up to 5wt% or LEU+ (i.e., below 10wt% enrichment) pellets with gadolinium oxide (Gd2O3) as the burnable poison. Satisfactory core and fuel performances are obtained for all the designed cores under steady state and considered overpower transients. For low power density operation, long cycle lengths are achieved reaching a 2.5- and a 5-year cycles and peak rod-average burnup is pushed to 83 MWd/kgU. Other cycle lengths are maintained at 18 months. Boron free operation exhibits the ability to achieve longer cycle lengths at the cost of higher peaking factors leading to high local power and fuel temperatures which prevents sizable power uprates and is deemed uneconomical. Fuel assembly height reduction allows coolant velocity retrofit which enables higher core power density without violating structural integrity of the fuel assembly. As a result, a core power density of 123 kW/l is reached where total cladding hoop strain becomes the limiting parameter.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

An finite element analysis surrogate model with boundary oriented graph embedding approach for rapid design

Abstract In this work, we present a boundary oriented graph embedding (BOGE) approach for the graph neural network to assist in rapid design and digital prototyping. The cantilever beam problem has been solved as an example to validate its potential of providing physical field results and optimized designs using only 10 ms. Providing shortcuts for both boundary elements and local neighbor elements, the BOGE approach can embed unstructured mesh elements into the graph and performs an efficient regression on large-scale triangular-mesh-based finite element analysis (FEA) results, which cannot be realized by other machine-learning-based surrogate methods. It has the potential to serve as a surrogate model for other boundary value problems. Focusing on the cantilever beam problem, the BOGE approach with 3-layer DeepGCN model achieves the regression with mean square error (MSE) of 0.011 706 (2.41% mean absolute percentage error) for stress field prediction and 0.002 735 MSE (with 1.58% elements having error larger than 0.01) for topological optimization. The overall concept of the BOGE approach paves the way for a general and efficient deep-learning-based FEA simulator that will benefit both industry and Computer Aided Design (CAD) design-related areas.

42 ENGINEERING↗

Data driven drift correction for complex optical systems

To exploit the thousand-fold increase in spectral brightness of modern light sources, increasingly intricate experiments are being conducted that demand extremely precise beam trajectory. Maintaining the optimal trajectory over several hours of an experiment with the needed precision necessitates active drift control. Here, we outline time varying Bayesian optimization (TVBO) as a data driven approach for robust drift correction, and illustrate its application for a split and delay optical system composed of six crystals and twelve input dimensions. Using numerical simulations, we exhibit the application of TVBO for linear drift, non-smooth temporal drift as well as constrained TVBO for multi-objective control settings, representing real-life operating conditions. This approach can be easily adapted to other X-ray beam conditioning and guidance systems, including multi-crystal monochromators and grazing-incidence mirrors, to maintain sub-micrometer and nanoradian beam stability over the course of an experiment spanning several hours.

Bayesian optimization↗

Findings from Design and Operation of Connected Neighborhoods

Distributed energy resources and load management are an important emerging part of smart grids. Concurrent management of homeowner preferences and utility objectives requires an intelligent negotiation strategy supported by physical infrastructure that learns, optimizes, and controls the system. While there is a wide body of research related to modelling and simulation of distributed resources, relatively little is known about their practical feasibility. One of the first attempts to address this knowledge gap was through a 62-house connected neighborhood in Alabama. This study makes one more step in advancing this knowledge through a 46-townhome demonstration neighborhood located in Atlanta, GA. It reports hardware design, system architecture, and the results from the summer experimental work. The findings are discussed in the context of earlier experience, and analysis is expanded to areas which were not the focus of earlier research.

Tsybina, Eve↗

Super-resolution within a one-dimensional phononic crystal of resonators using time reversal in an equivalent circuit model

An equivalent circuit model has been developed to model a one-dimensional waveguide with many side-branch Helmholtz resonators. This waveguide constitutes a phononic crystal that has been shown to have decreased phase speed below the resonance frequency of an individual resonator. This decreased phase speed can be exploited to achieve super-resolution using broadband time reversal focusing techniques. It is shown that the equivalent circuit model is capable of quantifying this change in phase speed of the crystal and also the small-scale wave-resonator interactions within the crystal. The equivalent circuit model enables the parameterization of the physical variables and the optimization of the focusing bandwidth by balancing the combination of increasing resolution and decreasing amplitude near the resonance frequency. It is shown that the quality factor—in this case, the quality factor determined by the geometric shape of each resonator—controls the range of frequencies that are strongly affected by the Helmholtz resonators.

Kingsley, Adam D.↗

Leveraging Application Data Constraints to Optimize Database-Backed Web Applications

Exploiting the relationships among data is a classical query optimization technique. As persistent data is increasingly being created and maintained programmatically, prior work that infers data relationships from data statistics misses an important opportunity. We present Coco, the first tool that identifies data relationships by analyzing database-backed applications. Once identified, Coco leverages the constraints to optimize the application's physical design and query execution. Instead of developing a fixed set of predefined rewriting rules, Coco employs an enumerate-test-verify technique to automatically exploit the discovered data constraints to improve query execution. Each resulting rewrite is provably equivalent to the original query. Using 14 real-world web applications, our experiments show that Coco can discover numerous data constraints from code analysis and improve real-world application performance significantly.

Computer Science↗

Collaborative Research: Improved Efficiency and Coupling of the Radiation Code in the ACME Earth System Model (Final Report)

The complexity of radiative transfer, its importance to the exchange of energy in the climate system, and its high computational cost establishes importance of an accurate and efficient radiative transfer parameterization for climate simulation. Previous work by the proposing team at AER led to the development of the radiation code, RRTMG, which has been widely accepted by the global modeling community as a fast and accurate advancement over the previous generation of radiation codes. It has been in use in the NCAR CESM for many years, and it has been implemented in the initial version the DOE E3SM model. However, its computational cost remains high relative to other components in part due to its complexity and to its inefficient use of modern optimization strategies, and this project helped address this limitation for the code’s application in E3SM. Under other funding, the Investigators of this project led an effort to develop a high-performance broadband radiation code, called RTE+RRTMGP, which is a completely restructured code that will take advantage of modern computational capabilities to enhance its performance while retaining the strengths and accuracy of the original code. Designed to perform over a range of computer architectures, RTE+RRTMGP makes extensive use of Fortran 2003 features to improve both its efficiency and its use of memory. RTE+RRTMGP is expected to be adopted in the next generation of RTE+RRTMGP. This project allowed for advancements to the code’s computational capabilities and adding new and enhanced features. This included revising RTE+RRTMGP to run on GPU processors, analysis of and improvements to the code’s timing, modifying the gas optics in RRTMGP to increase the code’s accuracy, extensive validation of computed fluxes, heating rates and forcings, generating a cloud optical property and vertical sampling capabilities for RTE+RRTMGP, implementing a fast and accurate longwave scattering capability, and groundwork for the inclusion of a capability to specify solar variability. Many of the accomplishment in this project necessitated significant collaboration with the E3SM development team. The result of this project was optimization of a key physical component (radiative transfer calculations) of E3SM, directly supporting E3SM’s overarching global modeling objectives. More broadly, this project provided overall advancements in the use of radiative transfer calculations in atmospheric modeling and simulation, particularly for climate.

58 GEOSCIENCES↗

Collaborative Research: Improved Efficiency and Coupling of the Radiation Code in the ACME Earth System Model. Final Report

This final report details all work performed on the project by both project partners. This project provided support to properly couple RTE+RRTMGP, a high-performance broadband radiation code, within DOE’s Energy Exascale Earth System Model (E3SM). RTE+RRTMGP is a successor to the RRTMG radiation code, which has been widely accepted for its speed and accuracy by the global modeling community, and has been in use in the NCAR CESM for many years and was implemented in the initial version of E3SM. However, the computational cost of RRTMG remains high relative to other components in part due to its complexity and to its inefficient use of modern optimization strategies, issues that were rectified by the development of RTE+RRTMGP. Many of the accomplishment in this project necessitated significant collaboration with the E3SM development team. One focus of the project was to enhance the code’s optimization on the limited number of emerging computing systems on which the model is expected be used, including Many Integrated Core (MIC) architectures and Graphics Processing Unit (GPU) hardware. We also developed additional capabilities for RTE+RRTMGP that E3SM scientists identified as important for the planned applications of the model. The result of our project was optimization of a key physical component (radiative transfer calculations) of E3SM, directly supporting E3SM’s overarching global modeling objectives. More broadly, this project provided overall advancements in the use of radiative transfer calculations in atmospheric modeling and simulation, particularly for climate.

54 ENVIRONMENTAL SCIENCES↗

LANL update to 239 Pu in the fast energy range [Slides]

Researchers are in the process of updating the Pu cross sections in the fast energy range (complementing the work of IAEA / INDEN / ORNL lower energy work in the resonance range). LANL is overhauling its evaluation tools (CoH, CGMF, DeCE, Kalman, NEXUS, PySOK, SOK) and focusing on consistency throughout evaluation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Spiner vs EOSPAC6: capabilities, performance, and accuracy considerations

The Cross-Cutting Capabilities Project (XCAP) has significant interest in collecting and implementing a common set of libraries to be used across the various production hydro-codes at LANL. A common set of libraries will help facilitate comparisons between codes. In addition to minimizing variables for comparisons and physics validation, efforts toward optimization and porting of packages for future architectures will be more efficient. As part of this push Singularity is being considered as an inclusive materials interface library. One crucial part of this library is the importing, inverting, and interpolation of equation of state (EOS) data. Currently the Lagrangian Applications Project (LAP) and Safety Applications Project (SAP) are using EOSPAC6 for delivery EOS data. The Eulerian Applications Project (EAP) has historically used TEOS and more recently an implementation of Singularly via a package called Spiner for EOS interpolation. Singularity currently has the option of using several analytical EOS models, directly employing SESAME via Spiner, or directly using SESAME via EOSPAC6 (not optimized yet). Assuming Singularity moves forward as a common platform for the implementation of materials models, a decision will need to be made on which EOS interpolation package XCAP should move forward with, given resources and people are finite. We will attempt to address pros and cons of Spiner and EOSPAC6 and the trade-offs that should be considered during the decision making process. This report is intended as an ASC-PEM-EOS perspective on what is needed in an EOS interpolation package. Performance and Accuracy sections will mostly address data from comparison studies in LA-UR-22-22699.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗

Novel Materials R&D for Next-Generation Accelerator Target Facilities

High-Entropy Alloys and Electrospun Nanofiber materials are two novel classes of materials that can offer improved resistance to beam-induced radiation damage and thermal shock. Research to develop these new materials specifically for multi-megawatt accelerator target applications, such as beam windows and particle-production targets, has recently begun. The research program will combine in-beam experiments with complementary simulations to tailor the microstructures of these novel materials for use in next-generation accelerator target facilities. Iterative simulations to optimize the material composition, physics performance and beam-induced thermomechanical response will guide the material design and fabrication processes based on established figures of merit. Ensuing material irradiation experiments using low-energy ions and prototypic high-energy protons, followed by extensive post-irradiation material characterization, will then assess and qualify the selected novel materials. This talk will provide an overview of the novel materials development research program initiated at Fermilab through my DOE Early Career Research Program award.

43 PARTICLE ACCELERATORS↗

Novel Materials R&D for Next-Generation Accelerator Target Facilities

High-Entropy Alloys and Electrospun Nanofiber materials are two classes of novel materials that can offer improved resistance to beam-induced radiation damage and thermal shock. Research to develop these new materials specifically for multi-megawatt accelerator target applications, such as beam windows and particle-production targets, are ongoing at Fermilab within the scope of a DOE Early Career Research Program. The research program combines in-beam experiments with complementary simulations to tailor the microstructures of these novel materials for use in next-generation accelerator target facilities. Iterative simulations to optimize the material composition, physics performance, and beam-induced thermomechanical response will guide the material design and fabrication processes based on established figures of merit. This will be followed by material irradiation experiments using low-energy ions and prototypic high-energy protons with extensive post-irradiation material characterization to assess and qualify the selected novel materials. This talk will describe the alloy design and synthesis, microstructural pre-characterization of the alloys, and plans for the eventual down selection following low-energy ion irradiation studies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Multibody for Everybody (M4E): A Symbolic Dynamics Modeling Tool with Applications in Simulation, Control, and Optimization

Developing the analytical model of a multibody system is often the initial step in control and optimization. The analytical model (equations of motion) describes a system’s time evolution under specified forcing conditions. Although developing these equations is easy for simple systems, this process becomes more complex for systems composed of multiple bodies. Deriving equations of motion for complex multibody systems requires specialized expertise in multibody dynamics, is time-consuming, and is susceptible to error. To address this issue, this paper presents an open-source, easy-to-use, systematic framework to derive symbolic equations of motion in both Python and MATLAB using the joint coordinate formulation. This formulation results in a set of ordinary differential equations that use the minimum set of coordinates needed to model a system. The symbolic representation provides better insight into the influence of design parameters on system performance, facilitates sensitivity analysis and parameter studies, and supports direct implementation of control and optimization routines. The tool enables numerical simulation for specified parameter sets, is modular for straightforward integration with other tools and libraries, and allows incorporation of hydrodynamics, mooring, and other external forces. The result is a reproducible, extensible pipeline for modeling, simulation, and design of complex multibody systems. The proposed tool is versatile and can be applied to domains such as robotics, control, and design. In addition, we integrated external libraries that provide capabilities for modeling offshore systems such as underwater robots and marine energy converters.

16 TIDAL AND WAVE POWER↗