Search NASA⌕ Search

SEARCH · Search NASA

Results for “core 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 685 records · Page 38

The HEPCloud Facility: elastic computing for High Energy Physics – The NOvA Use Case

The need for computing in the HEP community follows cycles of peaks and valleys mainly driven by conference dates, accelerator shutdown, holiday schedules, and other factors. Because of this, the classical method of provisioning these resources at providing facilities has drawbacks such as potential overprovisioning. As the appetite for computing increases, however, so does the need to maximize cost efficiency by developing a model for dynamically provisioning resources only when needed. To address this issue, the HEPCloud project was launched by the Fermilab Scientific Computing Division in June 2015. Its goal is to develop a facility that provides a common interface to a variety of resources, including local clusters, grids, high performance computers, and community and commercial Clouds. Initially targeted experiments include CMS and NOvA, as well as other Fermilab stakeholders. In its first phase, the project has demonstrated the use of the “elastic” provisioning model offered by commercial clouds, such as Amazon Web Services. In this model, resources are rented and provisioned automatically over the Internet upon request. In January 2016, the project demonstrated the ability to increase the total amount of global CMS resources by 58,000 cores from 150,000 cores - a 25 percent increase - in preparation for the Recontres de Moriond. In March 2016, the NOvA experiment has also demonstrated resource burst capabilities with an additional 7,300 cores, achieving a scale almost four times as large as the local allocated resources and utilizing the local AWS s3 storage to optimize data handling operations and costs. NOvA was using the same familiar services used for local computations, such as data handling and job submission, in preparation for the Neutrino 2016 conference. In both cases, the cost was contained by the use of the Amazon Spot Instance Market and the Decision Engine, a HEPCloud component that aims at minimizing cost and job interruption. This paper describes the Fermilab HEPCloud Facility and the challenges overcome for the CMS and NOvA communities.

Fuess, S.↗

Learning to Branch with Interpretable Machine Learning Models

Machine learning is being increasingly used in improving decisions made within branch-and-bound algorithms for solving mixed-integer programs (MIPs). Branching is a key component in branch-and-bound algorithms, this work presents IDAES-core project update on building simple and interpretable machine learning models for branching and improving decision-making tools applied for the optimization of advanced energy systems.

Bayramoglu, Selin↗

Mixed-precision iterative refinement using tensor cores on GPUs to accelerate solution of linear systems

Double-precision floating-point arithmetic (FP64) has been the de facto standard for engineering and scientific simulations for several decades. Problem complexity and the sheer volume of data coming from various instruments and sensors motivate researchers to mix and match various approaches to optimize compute resources, including different levels of floating-point precision. In recent years, machine learning has motivated hardware support for half-precision floating-point arithmetic. A primary challenge in high-performance computing is to leverage reduced-precision and mixed-precision hardware. We show how the FP16/FP32 Tensor Cores on NVIDIA GPUs can be exploited to accelerate the solution of linear systems of equations Ax = b without sacrificing numerical stability. The techniques we employ include multiprecision LU factorization, the preconditioned generalized minimal residual algorithm (GMRES), and scaling and auto-adaptive rounding to avoid overflow. We also show how to efficiently handle systems with multiple right-hand sides. On the NVIDIA Quadro GV100 (Volta) GPU, we achieve a 4×-5× performance increase and 5× better energy efficiency versus the standard FP64 implementation while maintaining an FP64 level of numerical stability.

GMRES↗

Optimized pulses for Raman excitation through the continuum: Verification using the multiconfigurational time-dependent Hartree-Fock method

Here, we have verified a mechanism for Raman excitation of atoms through continuum levels previously obtained by quantum optimal control using the multiconfigurational time-dependent Hartree-Fock (MCTDHF) method. For the optimal control, which requires running multiple propagations to determine the optimal pulse sequence, we used the computationally inexpensive time-dependent configuration interaction singles (TDCIS) method. TDCIS captures all of the necessary correlation of the desired processes but assumes that ionization pathways reached via double excitations are not present. MCTDHF includes these pathways and all multiparticle correlations in a set of time-dependent orbitals. The mechanism that was determined to be optimal in the Raman excitation of the Ne 1s 2 2s 2 2p 5 3p 1 valence state via the metastable 1s 2 2s 1 2p 6 3p 1 resonance state involves a sequential resonance-valence excitation. First, a long pump pulse excites the core-hole state, and then a shorter Stokes pulse transfers the population to the valence state. This process represents the first step in a multidimensional x-ray spectroscopy scheme that will provide a local probe of valence electronic correlations. Although at the optimal pulse intensities at the TDCIS level of theory the MCTDHF method predicts multiple ionization or excitation ionization of the atom, at slightly lower intensities (reduced by a factor of about 4) the TDCIS mechanism is shown to hold qualitatively. Quantitatively, the MCTDHF populations are reduced from the TDCIS calculations by a factor of 4.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Exploiting Geometric Partitioning in Task Mapping for Parallel Computes

We present a new method for mapping applications' MPI tasks to cores of a parallel computer such that applications' communication time is reduced. We address the case of sparse node allocation, where the nodes assigned to a job are not necessarily located in a contiguous block nor within close proximity to each other in the network, although our methods generalize to contiguous allocations as well. The goal is to assign tasks to cores so that interdependent tasks are performed by "nearby' cores, thus lowering the distance messages must travel, the amount of congestion in the network, and the overall cost of communication. Our new method applies a geometric partitioning algorithm to both the tasks and the processors, and assigns task parts to the corresponding processor parts. We also present a number of algorithmic optimizations that exploit specific features of the network or application. We show that, for the structured finite difference mini-application MiniGhost, our mapping methods reduced communication time up to 75% relative to MiniGhost's default mapping on 128K cores of a Cray XK7 with sparse allocation. For the atmospheric modeling code E3SM/HOMME, our methods reduced communication time up to 31% on 32K cores of an IBM BlueGene/Q with contiguous allocation.

97 MATHEMATICS AND COMPUTING↗

SymProp: Scaling Sparse Symmetric Tucker Decomposition via Symmetry Propagation

Sparse symmetric tensors are an important class of tensors, and their decompositions serve as powerful tools for revealing low-rank structures. This paper introduces SymProp, a novel approach for scaling sparse symmetric Tucker decomposition by propagating symmetry through intermediate computations. SymProp optimizes two key computational kernels: Sparse Symmetric Tensor Times Same Matrix chain (S3 TTMc) for Higher-Order Orthogonal Iteration (HOOI) and Sparse Symmetric Tensor Times Same Matrix chain Times Core (S3 TTMcTC) for Higher-Order QR Iteration (HOQRI). Our method employs a metaprogramming-based index iteration approach to efficiently handle the upper triangular parts of intermediate dense symmetric tensors. SymProp achieves up to 50.9× speedup over SPLATT and up to 360.8× over Compressed Sparse Symmetric (CSS) format on the S3 TTMc operation. Moreover, our S3 TTMc and S3 TTMcTC implementations support tensor orders four levels higher than state-of-the-art methods. Our HOQRI demonstrates superior scalability and up to a 33.6× speedup over optimized HOOI. By enabling more scalable Tucker decompositions for higher orders, decomposition ranks, and dimension sizes, SymProp opens new possibilities for analyzing complex hypergraph structures in fields such as network science, data mining, and machine learning.

Li, Zecheng [North Carolina State University]↗

Multispecies density peaking in gyrokinetic turbulence simulations of low collisionality Alcator C-Mod plasmas

Peaked density profiles in low-collisionality AUG and JET H-mode plasmas are probably caused by a turbulently driven particle pinch, and Alcator C-Mod experiments confirmed that collisionality is a critical parameter. Density peaking in reactors could produce a number of important effects, some beneficial, such as enhanced fusion power and transport of fuel ions from the edge to the core, while others are undesirable, such as lower beta limits, reduced radiation from the plasma edge, and consequently higher divertor heat loads. Fundamental understanding of the pinch will enable planning to optimize these impacts. We show that density peaking is predicted by nonlinear gyrokinetic turbulence simulations based on measured profile data from low collisionality H-mode plasma in Alcator C-Mod. Multiple ion species are included to determine whether hydrogenic density peaking has an isotope dependence or is influenced by typical levels of low-Z impurities, and whether impurity density peaking depends on the species. We find that the deuterium density profile is slightly more peaked than that of hydrogen, and that experimentally relevant levels of boron have no appreciable effect on hydrogenic density peaking. The ratio of density at r/a = 0.44 to that at r/a = 0.74 is 1.2 for the majority D and minority H ions (and for electrons), and increases with impurity Z: 1.1 for helium, 1.15 for boron, 1.3 for neon, 1.4 for argon, and 1.5 for molybdenum. The ion temperature profile is varied to match better the predicted heat flux with the experimental transport analysis, but the resulting factor of two change in heat transport has only a weak effect on the predicted density peaking.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Transformational Challenge Reactor – On the Application of Design for Additive Manufacturing (DfAM) Techniques to the Conception of Nuclear Core

Additive manufacturing (AM) technologies are radically changing the way objects are designed and manufactured. They allow building by deposition and solidification of material layer by layer, enabling the possibility to create simple and complex features alike, almost seamlessly. Generally, the design optimization process requires to define objectives, design variables and constraints. Additive manufacturing does not challenge this process per se but does allow designers to completely redefine the constraints space as the ones originating from fabrication can be considerably relaxed compared to more “traditional” manufacturing. Thus, design optimization becomes naturally far more responsive to the actual physics being solved and considerably less influenced by fabrication limitations, leading to dramatically different designs. To take advantage of these new opportunities, so-called Designing for Additive Manufacturing (DfAM) techniques are emerging. Development of design techniques specifically tailored for additive manufacturing is warranted because, considering AM, the design space is typically considerably larger than with traditional manufacturing. The ability to explore the design space efficiently is of paramount importance for designers. This study proposes to investigate and apply some of these DfAM techniques to the conception of nuclear core. The goal being to assess if these new methods can be applied to core design and if core design could benefits from additive manufacturing technologies. After a brief investigation on the pertinence of some DfAM techniques for core design, algorithms are proposed and a workflow is established to carry neutronics and steady-state thermal-hydraulics analyses. To diminish the work load, the workflow has been automated using python modules. These modules allow the rapid creation of input files, post-treatment of output files and visualization. To test the pertinence of the proposed workflow, three test cases have been investigated: a research and test reactor, a micro-reactor and a space propulsion reactor. These test cases offered a variety of objectives, constraints and operating conditions. It is observed that the proposed workflow is capable of converging quickly and efficiently to valid design solutions. It is then concluded that DfAM techniques can be applied to core design.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Reduced Order Models Generation for HTGRs Pebble Shuffling Procedure Optimization Studies

This report provides an initial study for producing reduced-order models (ROMs) of pebble-bed high temperature gas reactor (HTGR) models for the purposes of design optimization. As an initial study, this work is meant to be exploratory---identifying useful workflows and methods for ROM generation---and not meant to be a catch-all analysis of HTGR ROM generation and usage for optimization. This report summarizes three tasks performed in Fiscal Year 2022: 1) the creation of HTGR model, 2) the sensitivity analysis of model design parameters, and 3) an introduction to ROM generation techniques. The representative HTGR model created in this work is a multiphysics equilibrium-core using the BlueCRAB (comprehensive reactor analysis bundle) reactor analysis application, coupling four physical phenomena: neutronics, streamline depletion, porous flow thermal hydraulics, and pebble heat conduction. Part of the model creation was identifying some design parameters and quantities of interest that are relevant in an optimization analysis and adjustable in the model. The sensitivity analysis utilized a polynomial chaos expansion methodology to compute global sensitivity metrics. This analysis showed that thermal hydraulics parameters and quantities of interest had a relatively small impact on simulation results. Finally, the ROM generation work involved exploring three different ROM methodologies: polynomial regression, a Gaussian process, and artificial neural networks. Using a cross-validation technique to characterize ROM performance, the Gaussian process and single-layer artificial neural networks showed the most promising results. Overall, this study was insightful and the lessons learned will be invaluable for the eventual development of an HTGR design optimization workflow.

97 MATHEMATICS AND COMPUTING↗

A Dynamic Risk Framework for the Physical Security of Nuclear Power Plants

This paper describes ongoing work within the Light Water Reactor Sustainability pathway at Idaho National Laboratory (INL) to optimize the security and cost of nuclear power plants. It introduces the dynamic risk assessment tool developed at INL, Event Modeling Risk Assessment using Linked Diagrams (EMRALD). EMRALD is leveraged to optimize the security posture of a nuclear power plant by integrating force-on-force (FOF) simulations and operator mitigation actions, including dynamic and flexible coping strategies (FLEX). To illustrate the methodology, four attack scenarios are modeled in a commercially available FOF simulation tool using a hypothetical nuclear power plant facility. The simulation results provide valuable insights into possible attack outcomes, as well as the probabilistic risk of a core damage event given these outcomes. Safety mitigation procedures are modeled in EMRALD dependent on the attack outcomes by considering human operator uncertainties. The results demonstrate that the number of armed responders can be optimized, while still maintaining the same protection level as the initial security posture. The proposed modeling and simulation framework of integrating FLEX equipment with FOF models enables the nuclear power plants to credit FLEX portable equipment in the plant security posture, resulting in an efficient and optimized physical security system.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A Dynamic Risk Framework for the Optimization of Physical Security Posture of Nuclear Power Plants

This paper describes an ongoing work within the Light Water Reactor Sustainability pathway at Idaho National Laboratory (INL) to optimize security and cost of nuclear power plants. It introduces the dynamic risk assessment tool developed at INL, Event Modeling Risk Assessment using Linked Diagrams (EMRALD). EMRALD was leveraged to optimize the security posture of a nuclear power plant by integrating force-on-force (FOF) simulations and operator mitigation actions including the dynamic and flexible coping strategies (FLEX). To illustrate the methodology, four attack scenarios were modeled in a commercially available FOF simulation tool using a hypothetical nuclear power plant facility. The simulation results provide valuable insights into possible attack outcomes, as well as the probabilistic risk of core damage event given these outcomes. Safety mitigation procedures were modeled in EMRALD dependent on the attack outcomes by considering human operator uncertainties. The results demonstrate that the number of armed responders can be optimized, while still maintaining the same protection level as the initial security posture. The proposed modeling and simulation framework of integrating FLEX equipment with FOF models enables the nuclear power plants to credit FLEX portable equipment in the plant security posture, resulting in an efficient and optimized physical security system.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A generalized approach for determining electric machine drive switching schemes using evolutionary optimization algorithms

The number of inverter switching events in one fundamental cycle of an electric machine can be limited by many factors such as inverter switching frequency, rotational speed, and computational time required by control algorithms. Even at low and moderate speeds, switching harmonics cause distortion in the phase current, which can lead to increased torque ripple and core losses when compared to operation with ideal waveforms. Furthermore, as rotational speed increases and voltage limitations are encountered, overmodulation schemes must be implemented to overcome increasing reactance and back-EMF (or pseudo back-EMF). The goal of this paper is to establish a framework and demonstrate methods for determining optimal switching times for a set number of switching events in a fundamental cycle using particle swarm optimization and multivariable fitness functions.

Burress, Timothy↗

Mesocrystallizing Nanograins for Enhanced Li + Storage

Morphology and crystallinity of electrode materials have a major effect on their charge carrier storage properties when applied in rechargeable batteries. While nanosizing electrode particles (with larger surface area) and maintaining electrode integrity are both good for performance enhancement, they seem to contradict each other and can be hardly well-balanced. Herein, we designed electrode particles consisting of numerous nanograins with uniform crystalline orientation to guarantee both high surface area and high structural integrity, allowing the significant improvement of Li + storage kinetics and performance. Applying this “mesocrystallizing” strategy to a NiCo 2 O 4 -based anode, resulting in various degrees of pseudocapacitance response, the long-term cyclability and rate performance of this material have been largely enhanced. Impressively, the mesocrystalline NiCo 2 O 4 electrode exhibits a high reversible capacity of 1403 mAh g -1 after 200 cycles at 1.6 A g -1 (a rate of 1.8 C). The growth mechanism of mesocrystalline materials with different morphologies is identified to be a topotactic strctural transition process featuring a gradual edge-to-core corrasion process. This work presents an important synthetic clue to balance the morpholgy and crystallinity of battery electrode materials for their perforamnce optimization and is expected to inspire future structural design for battery materials beyond the one prototyped here.

25 ENERGY STORAGE↗

Very High-Risk Localized Prostate Cancer: Outcomes Following Definitive Radiation

Purpose: Existing definitions of high-risk prostate cancer consist of men who experience significant heterogeneity in outcomes. As such, criteria that identify a subpopulation of National Comprehensive Cancer Network (NCCN) high-risk prostate cancer patients who are at very high risk (VHR) for poor survival outcomes following prostatectomy were recently developed at our institution and include the presence of any of the following disease characteristics: multiple NCCN high-risk factors, primary Gleason pattern 5 disease and/or ≥5 biopsy cores with Gleason sums of 8 to 10. Whether these criteria also apply to men undergoing definitive radiation is unclear, as is the optimal treatment regimen in these patients. Methods and Materials: All men consecutively treated with definitive radiation by a single provider from 1993 to 2006 and who fulfilled criteria for NCCN high-risk disease were identified (n=288), including 99 patients (34%) with VHR disease. Multivariate-adjusted competing risk regression models were constructed to assess associations between the VHR definition and biochemical failure (BF), distant metastasis (DM), and prostate cancer–specific mortality (PCSM). Multivariate-adjusted Cox regression analysis assessed the association of the VHR definition with overall mortality (OM). Cumulative incidences of failure endpoints were compared between VHR men and other NCCN high-risk men. Results: Men with VHR disease compared to other NCCN high-risk men experienced a higher 10-year incidence of BF (54.0% vs 35.4%, respectively, P<.001), DM (34.9% vs 13.4%, respectively, P<.001), PCSM (18.5% vs 5.9%, respectively, P<.001), and OM (36.4% vs 27.0%, respectively, P=.04). VHR men with a detectable prostate-specific antigen (PSA) concentration at the end of radiation (EOR) remained at high risk of 10-year PCSM compared to VHR men with an undetectable EOR PSA (31.0% vs 13.7%, respectively, P=.05). Conclusions: NCCN high-risk prostate cancer patients who meet VHR criteria experience distinctly worse outcomes following definitive radiation and long-term androgen deprivation therapy, particularly if an EOR PSA is detectable. Optimal use of local therapies for VHR patients should be explored further, as should novel agents.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Quantum Dynamics at Scale: Ultrafast Control of Emergent Functional Materials

Confluence of extreme-scale quantum dynamics simulations (i.e. quantum@scale) and cutting-edge x-ray free-electron laser experiments are revolutionizing materials science. An archetypal example is the exciting concept of using picosecond light pulses to control emergent material properties on demand in atomically-thin layered materials. This paper describes efforts to scale our quantum molecular dynamics engine toward the United States' first exaflop/s computer, under an Aurora Early Science Program project named "Metascalable layered material genome". Key algorithmic and computing techniques incorporated are: (1) globally-scalable and locally-fast solvers within a linear-scaling divide-conquer-recombine algorithmic framework; (2) algebraic `BLASification' of computational kernels; and (3) data alignment and loop restructuring, along with register and cache blocking, for enhanced vectorization and efficient memory access. The resulting weak-scaling parallel efficiency was 0.93 on 131,072 Intel Xeon Phi cores for a 56.6 million atom (or 169 million valence-electron) system, whereas the various code transformations achieved 5-fold speedup. The optimized simulation engine allowed us for the first time to establish a significant effect of substrate on the dynamics of layered material upon electronic excitation.

Algebraic BLASification↗

Selective Equatorial Sclera Crosslinking in the Orbit Using a Metal-Coated Polymer Waveguide

Purpose: Photochemical crosslinking of the sclera is an emerging technique that may prevent excessive eye elongation in pathologic myopia by stiffening the scleral tissue. To overcome the challenge of uniform light delivery in an anatomically restricted space, we previously introduced the use of flexible polymer waveguides. We presently demonstrate advanced waveguides that are optimized to deliver light selectively to equatorial sclera in the intact orbit. Methods: Our waveguides consist of a polydimethylsiloxane cladding and a polyurethane core, coupled to an optical fiber. A reflective silver coating deposited on the top and side surfaces of the waveguide prevents light leakage to nontarget, periorbital tissue. Postmortem rabbits were used to test the feasibility of in situ equatorial sclera crosslinking. Tensometry measurements were performed on ex vivo rabbit eyes to confirm a biomechanical stiffening effect. Results: Metal-coated waveguides enabled efficient light delivery to the entire circumference of the equatorial sclera with minimal light leakage to the periorbital tissues. Blue light was delivered to the intact orbit with a coefficient of variation in intensity of 22%, resulting in a 45 ± 11% bleaching of riboflavin fluorescence. A 2-fold increase in the Young's modulus at 5% strain (increase of 92% P < 0.05, at 25 J/cm 2 ) was achieved for ex vivo crosslinked eyes. Conclusions: Flexible polymer waveguides with reflective, biocompatible surfaces are useful for sclera crosslinking to achieve targeted light delivery. We anticipate that our demonstrated procedure will be applicable to sclera crosslinking in live animal models and, potentially, humans in vivo.

59 BASIC BIOLOGICAL SCIENCES↗

Efficient (~10%) generation of vacuum ultraviolet femtosecond pulses via four-wave mixing in hollow-core fibers

Here we report the generation of the fifth harmonic of Ti:sapphire, at 160 nm, with more than 4 µJ of pulse energy and a pulse length of 37 fs with a 1 kHz repetition rate. The vacuum ultraviolet pulses are produced using four-wave difference frequency mixing in a He-filled stretched hollow-core fiber, driven by a pump at 267 nm and seeded at 800 nm. Guided by simulations using Luna.jl, we are able to optimize the process carefully. The result is a conversion efficiency of ~10% from the 267 nm pump beam.

47 OTHER INSTRUMENTATION↗

Passive seismic investigation of Harrat Rahat

Ambient noise correlation was applied to 18 months of continuous seismic data from 14 stations. The procedure of Bensen et al [2007] was followed with some changes to optimize signal-to-noise of the results. The 18 months of correlations (representing about 1 week of CPU time on a 12 core machine) were stacked and manually inspected to yield about 40 cross-correlations. These cross-correlations represent the Green’s function between the station pairs and will be analyzed in part two of this project to yield velocity structure.

15 GEOTHERMAL ENERGY↗