Search NASA⌕ Search

SEARCH · Search NASA

Results for “high performance computing systems”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,729 records · Page 96

DFT-FE 1.0: A massively parallel hybrid CPU-GPU density functional theory code using finite-element discretization

In this work, we present DFT-FE 1.0, building on DFT-FE 0.6 [Comput. Phys. Commun. 246, 106853 (2020)], to conduct fast and accurate large-scale density functional theory (DFT) calculations (reaching ~ 100,000 electrons) on both many-core CPU and hybrid CPU-GPU computing architectures. This work involves improvements in the real-space formulation—via an improved treatment of the electrostatic interactions that substantially enhances the computational efficiency—as well high-performance computing aspects, including the GPU acceleration of all the key compute kernels in DFT-FE. We demonstrate the accuracy by comparing the ground-state energies, ionic forces and cell stresses on a wide-range of benchmark systems against those obtained from widely used DFT codes. Further, we demonstrate the numerical efficiency of our implementation, which yields ~ 20× CPU-GPU speed-up by using GPU acceleration on hybrid CPU-GPU nodes. Notably, owing to the parallel-scaling of the GPU implementation, we obtain wall-times of 80–140 seconds for full ground-state calculations, with stringent accuracy, on benchmark systems containing ~ 6, 000 – 15,000 electrons.

pseudopotential↗

Toward a high-fidelity tritium transport modeling for retention and permeation experiments

We report Tritium Migration Analysis Program version 8 (TMAP8), the latest version of TMAP, was developed within the framework of the Multiphysics Object-Oriented Simulation Environment (MOOSE). Created at Idaho National Laboratory (INL), MOOSE is an open-source, dimension-agnostic, fully coupled, and fully implicit multiphysics platform featuring massively parallel computation capabilities. Using TMAP8, tritium transport in a divertor monoblock was analyzed to elucidate the effects of pulsed operation (up to fifty 1,600 s plasma discharge and cool-down cycles) on the tritium in-vessel inventory source term and ex-vessel release term (i.e., tritium retention and permeation) for safety analysis. With its built-in Message Passing Interface capability, TMAP8 can, in under 2 h, simulate tritium transport in three different layered materials (i.e., tungsten, copper, and copper-chromium-zirconium alloy) in 2D geometry, using a single device/computer with 10 cores. The MOOSE-based TMAP8 code can leverage other MOOSE tools developed under the Nuclear Energy Advanced Modeling and Simulation program to perform tritium and thermal transport in complex geometries and multiphysics environments. And via its massively parallel computation, MOOSE will enable the fusion pilot plant designers to conduct high-fidelity multiphysics modeling for the design of the divertor and blanket systems as well as for the safety analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Full Implementation of Matrix Approach to Biogeochemistry Module of CLM5

Abstract Earth system models (ESMs) have been rapidly developed in recent decades to advance our understanding of climate change‐carbon cycle feedback. However, those models are massive in coding, require expensive computational resources, and have difficulty in diagnosing their performance. It is highly desirable to develop ESMs with modularity and effective diagnostics. Toward these goals, we implemented a matrix approach to the Community Land Model version 5 (CLM5) to represent carbon and nitrogen cycles. Specifically, we reorganized 18 balance equations each for carbon and nitrogen cycles among the 18 vegetation pools in the original CLM5 into two matrix equations. Similarly, 140 balance equations each for carbon and nitrogen cycles among the 140 soil pools were reorganized into two additional matrix equations. The vegetation carbon and nitrogen matrix equations are connected to soil matrix equations via litterfall. The matrix equations fully reproduce simulations of carbon and nitrogen dynamics by the original model. The computational cost for forwarding simulation of the CLM5 matrix model was 26% more expensive than the original model, largely due to calculation of additional diagnostic variables, but the spin‐up computational cost was significantly saved. We showed a case study on modeled soil carbon storage under two forcing data sets to illustrate the diagnostic capability that the matrix approach uniquely offers to understand simulation results of global carbon and nitrogen dynamics. The successful implementation of the matrix approach to CLM5, one of the most complex land models, demonstrates that most, if not all, the biogeochemical models can be reorganized into the matrix form to gain high modularity, effective diagnostics, and accelerated spin‐up.

58 GEOSCIENCES↗

Combining RAVEN, RELAP5-3D, and PHISICS for Fuel Cycle and Core Design Analysis for New Cladding Criteria

The Nuclear Regulatory Commission (NRC) has considered revision of 10 CFR 50.46C rule to account for the effects of burn-up rate in future analysis of reactor accident scenarios so that safety margins may evolve with the operation and reloading of the reactor as dynamic limits. In order to find these limiting conditions, both cladding oxidation and maximum temperature must be cast as functions of fuel exposure. It is computationally intensive to run a plant model through a long duration operational transient to the time of fuel reload, and this must be repeated for each reload until the time of the accident scenario. Moreover, fuel reload patterns may differ in many ways so that calculations with many different core configurations must be considered. To perform such new analyses in a reasonable amount of computational time with good accuracy, Idaho National Laboratory (INL) has developed new multi-physics tools by combining existing codes and adding new capabilities. The Parallel Highly Innovative Simulation INL Code System (PHISICS) toolkit for neutronic and reactor physics is coupled with the Reactor Excursion and Leak Analysis Program (RELAP5-3D) for the Loss of Coolant Accident (LOCA) analysis and RAVEN for the Probabilistic Risk Assessment (PRA) and margin characterization analysis. Expansion of code capability was also necessary. For RELAP5-3D to process a single sequence of cores in a continuous run required modification of RELAP5-3D to restart a sequence of input decks, each with a slightly modified neutronics model or thermal-hydraulic flow region and culminating in an accident scenario. A new multi-deck input processing capability was developed and verified for this analysis. In conclusion, the combined RAVEN/PHISICS/RELAP5-3D tool is used to analyze a typical Pressurized Water Reactor (PWR).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Development of CFD-Based Simulation Tools for In-Situ Thermal Processing of Oil Shale/Sands

In our research, we are taking the novel approach of developing and applying high performance computing, computational fluid dynamics (CFD)-based simulation tools to a modified in-situ process for production of oil from oil shale. The simulation tools being developed capture the relevant physical processes and data from a large-scale system. The modified in-situ application is a pilot-scale heat transfer process inside Red Leaf Resources EcoShale capsule. We demonstrate the need to understand fluid flow behavior in the convective channels of the rubblized shale bed as convective heating greatly decreases the time required to heat the oil shale to the production temperature when compared with conductive heating alone. We have developed and implemented a geometry creation strategy for a representative section of the EcoShale capsule, developed a meshing approach to deal with the complicated geometry and produce a well-behaved mesh, analyzed the effects of boundary conditions on the simulation results, and devised a new operator splitting solution algorithm that reduces computational costs by taking advantage of the differing convective and conductive time scales occurring in the simulation. These simulation tools can be applied to a wide range of processes involving convective fluid flow heating in rubblized beds.

04 OIL SHALES AND TAR SANDS↗

OpenSHMEM Implementation of IOR Benchmark

The benchmarking effort within the Computational Research & Development Programs at the Oak Ridge National Laboratory (ORNL) seeks to design and enable High Performance Computing (HPC) benchmarks and test suites. The work described in this report is a part of the effort focusing on the comparison and analysis of OpenSHMEM implementations using the Interleave Or Random (IOR) software for benchmarking parallel file system using POSIX, MPIIO, or HDF5 interfaces. We describe the effort to emulate the MPIIO parallel collective capabilities in the IOR benchmark using OpenSHMEM communication. One development effort was in emulating the MPI derived datatype used in the read/write operations and in setting the file view. Another effort was in implementing an internal cache in OpenSHMEM distributed shared memory to facilitate global collective I/O operations. Experiments comparing collective I/O in MPIIO implementations with the OpenSHMEM implementations were performed on the SGI Turing Cluster and the Cray XK7 Titan supercomputer at the Oak Ridge Leadership Computing Facility (OLCF). The preliminary results suggest that on the Cray XK7 Titan, the MPIIO implementations obtained higher write performance and the OpenSHMEM version obtained slightly higher read performance. On the SGI Turing Cluster, the MPIIO implementations obtained slightly higher performance over the OpenSHMEM implementations on large files.

97 MATHEMATICS AND COMPUTING↗

A scalable messaging system for accelerating discovery from large scale scientific simulations

Emerging scientific and engineering simulations running at scale on leadership-class High End Computing (HEC) environments are producing large volumes of data, which has to be transported and analyzed before any insights can result from these simulations. The complexity and cost (in terms of time and energy) associated with managing and analyzing this data have become significant challenges, and are limiting the impact of these simulations. Recently, data-staging approaches along with in-situ and in-transit analytics have been proposed to address these challenges by offloading I/O and/or moving data processing closer to the data. However, scientists continue to be overwhelmed by the large data volumes and data rates. In this paper we address this latter challenge. Specifically, we propose a highly scalable and low-overhead associative messaging framework that runs on the data staging resources within the HEC platform, and builds on the staging-based online in-situ/in- transit analytics to provide publish/subscribe/notification-type messaging patterns to the scientist. Rather than having to ingest and inspect the data volumes, this messaging system allows scientists to (1) dynamically subscribe to data events of interest, e.g., simple data values or a complex function or simple reduction (max()/min()/avg()) of the data values in a certain region of the application domain is greater/less than a threshold value, or certain spatial/temporal data features or data patterns are detected; (2) define customized in-situ/in-transit actions that are triggered based on the events, such as data visualization or transformation; and (3) get notified when these events occur. The key contribution of this paper is a design and implementation that can support such a messaging abstraction at scale on high- end computing (HEC) systems with minimal overheads. We have implemented and deployed the messaging system on the Jaguar Cray XK6 machines at Oak Ridge National Laboratory and the Lonestar system at the Texas Advanced Computing Center (TACC), and we present the experimental performance evaluation using these HEC platforms in the paper.

Jin, Tong↗

Parallel Monte Carlo Methods for Heterogeneous Hardware Computer Systems Using GPUs and Coprocessors: Recent Development of ARCHER Code - Paper 54

Traditional methods of accelerating Monte Carlo (MC) calculations typically use CPUs as the computing unit. The parallelism is based upon multithreading on a multi-core processor, and message passing interface (MPI) across many processors and nodes. In recent years, the hardware accelerators have been introduced by the high performance computing industry as an alternative approach to achieve massive parallelism. Hardware accelerators have the advantage of high Floating-point Operations Per Second (FLOPS) and high energy efficiency. Examples include graphics processing units (GPUs) and coprocessors. Compute-intensive programs can be off-loaded to them to be concurrently executed by hundreds or thousands of active threads. Existing production MC codes, however, cannot be directly compiled or executed on the new platforms. It is therefore necessary to rewrite or redevelop the code using the hardware accelerator specific programming models and APIs. Several studies have been conducted on using GPUs to accelerate MC calculations in reactor analysis and medical physics. In this paper, we present the recent progress at Rensselaer Polytechnic Institute (RPI) in developing a new parallel MC code called ARCHER (Accelerated Radiation-transport Computations in Heterogeneous Environments) for coupled photon-electron transport simulations. The code is designed as a test bed to evaluate the efficacy of heterogeneous computing systems for low-cost MC-based radiation dosimetric calculations. The code has three variants corresponding to three parallel hardware platforms: ARCHER-CPU (multi-core CPU), ARCHER-GPU (Nvidia GPU) and ARCHER-COP (Intel Xeon Phi coprocessor). The code has been used for several applications, including CT imaging dosimetry, radiotherapy dosimetry and radiation shielding design. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗

Effective Exploration of New 760°C-Capability Steels for Coal Energy

Cost effective and high performance alloys that are capable of operating at 760 °C or higher for extended periods of time under a very aggressive environment are critically required for the design and development of advanced ultrasupercritical (AUSC) boilers and steam turbines. Finely dispersed Laves phase precipitates have been shown by Takeyama and co-workers to be a viable strengthening mechanism in high temperature austenitic steels. There is currently no straightforward theory that can predict what other intermetallic phases can serve as potent precipitation-strengthening phases for steels; thus we employed a highly effective dual-anneal diffusion multiple (DADM) approach to screen for viable strengthening precipitates over a wide range of compositions. From the Fe-Co-Cr-Ni-Mo DADMs, the Fe-Cr-Mo based Chi phase was identified as a new strengthening phase for high temperature ferritic steels; and from the Fe-Mn-Cr-Nb-Ni-Mo-FeAl DADMs, the Laves phase was identified as a viable strengthening precipitate in Fe-Mn and Fe-Ni based austenitic steels. After identification of viable strengthening phases from the DADMs that covered compositions in the basic ternary and quaternary systems, we employed computation thermodynamics to perform multicomponent alloy design and optimization. For the new the Chi-phase strengthened steels, we performed thermodynamic calculations to vary the volume fraction of the Chi phase and introduced Nb and carbon to promote the formation of stable carbides for grain size control during solution heat treatment. For the Fe-Ni-Mn based austenitic steels, we performed extensive parametric optimization of compositions in order to reduce the expensive Ni content, add Cr and Al for oxidation resistance, and balance the alloying contents (Ni, Mn, Cr, Al, Mo) to suppress the ferritic phase and promote the austenitic matrix phase. Four steels (two ferritic + two austenitic) were designed and tested. The two Chi-phase strengthened ferritic steels exhibited excellent oxidation resistance and good creep-rupture strength at moderate temperatures, considering their ferritic matrix that usually results in lower creep resistance than austenitic steels. These steels showed brittleness and sample-to-sample variability in ductility. The low ductility might be due to the macro segregation during solidification or the significant grain growth during the solution heat treatments. We believe there is no inherent brittleness based on the chemistry of the steels. The creep-rupture performance of the steels is comparable to the 9Cr steels. Due to their ferritic matrix, the new Chi-phase strengthened ferritic steels may not be suited for the 760 °C AUSC applications, but they are very good candidates for intermediate temperature applications due to their outstanding oxidation resistance and high strength. Further study is required to find the source of low and highly variable ductility. We believe the compositions of the Chi-phase strengthened steels are not inherently brittle. The Chi-phase strengthened ferritic steels may also be excellent candidates for intermediate-temperature and room-temperature cast stainless steels, thus we highly recommend further investigations. The two Mn-containing austenitic steels based on the Laves phase showed good ductility, excellent oxidation resistance (slightly inferior to the two ferritic steels) at high temperatures and moderate creep strength. The creep-strength of the two austenitic steels based on the Larson-Miller parameters is higher than that of the traditional 316 stainless steels, but lower than the alumina-forming alloys (AFAs) developed at Oak Ridge National Laboratories. We do not recommend high priority in further studying these compositions unless higher Cr alloys are required for hot-corrosion resistance.

01 COAL, LIGNITE, AND PEAT↗

Region templates: Data representation and management for high-throughput image analysis

Compartment fires exhibit unique features associated with smoke accumulation and restricted air ventilation. The objective of the present study is to evaluate the ability of CFD-based fire models to simulate these features and in particular the effects of flame extinction that results from under-ventilated fire conditions. The study is performed using FireFOAM; FireFOAM is an advanced fire modeling software developed by FM Global and is based on a general-purpose open-source software called OpenFOAM. A new flame extinction model based on the concept of a critical value of the flame Damköhler number is incorporated into FireFOAM. The performance of the extinction model is evaluated via comparisons with a previously developed experimental database corresponding to a reduced-scale, heptane-fueled, compartment fire configuration. The numerical simulations also include a description of fuel evaporation driven by the computed gas-to-liquid heat feedback. Comparisons between experimental data and numerical results provide a suitable test bed to evaluate the ability of CFD-based fire models to describe the transition from over- to under-ventilated fire conditions, as well as the transition from extinction-free conditions to conditions in which the flame experiences partial or total quenching.

Energy & Fuels↗

Radiation Belt Remediation: A Complex Engineered System (RBR-ACES) (Institutional Computing Report)

We have performed the simulation campaigns in support of LDRD-DR project Radiation Belt Remediation: A Complex Engineered System (RBR-ACES). This project involves modeling wave generation and wave propagation in plasma environment and studying the effect of waves on the distribution of high energy electrons with the focus on electrons appeared after a high-altitude nuclear explosion (HANE). One of the major goal of this research is strong validation of the numerical models with laboratory experiments and space-based experiments and observations.

42 ENGINEERING↗

The PetscSF Scalable Communication Layer

PetscSF, the communication component of the Portable, Extensible Toolkit for Scientific Computation (PETSc), is designed to provide PETSc's communication infrastructure suitable for exascale computers that utilize GPUs and other accelerators. PetscSF provides a simple application programming interface (API) for managing common communication patterns in scientific computations by using a star-forest graph representation. PetscSF supports several implementations based on MPI and NVSHMEM, whose selection is based on the characteristics of the application or the target architecture. An efficient and portable model for network and intra-node communication is essential for implementing large-scale applications. The Message Passing Interface, which has been the de facto standard for distributed memory systems, has developed into a large complex API that does not yet provide high performance on the emerging heterogeneous CPU-GPU-based exascale systems. Here, we discuss the design of PetscSF, how it can overcome some difficulties of working directly with MPI on GPUs, and we demonstrate its performance, scalability, and novel features.

97 MATHEMATICS AND COMPUTING↗

Status on the Component Models Developed in the Modelica Framework: High-Temperature Steam Electrolysis Plant & Gas Turbine Power Plant

This report has been prepared as part of an effort to design and build a modeling and simulation (M&S) framework to assess the economic viability of a nuclear-renewable hybrid energy system (N-R HES). In order to facilitate dynamic M&S of such an integrated system, research groups in multiple national laboratories have been developing various subsystems as dynamic physics-based components using the Modelica programming language. In fiscal year (FY) 2015, Idaho National Laboratory (INL) performed a dynamic analysis of two region-specific N-R HES configurations, including the gas-to-liquid (natural gas to Fischer-Tropsch synthetic fuel) and brackish water reverse osmosis desalination plants as industrial processes. In FY 2016, INL has developed two additional subsystems in the Modelica framework: a high-temperature steam electrolysis (HTSE) plant and a gas turbine power plant (GTPP). HTSE has been proposed as a high priority industrial process to be integrated with a light water reactor (LWR) in an N-R HES. This integrated energy system would be capable of dynamically apportioning thermal and electrical energy (1) to provide responsive generation to the power grid and (2) to produce alternative industrial products (i.e., hydrogen and oxygen) without generating any greenhouse gases. A dynamic performance analysis of the LWR/HTSE integration case was carried out to evaluate the technical feasibility (load-following capability) and safety of such a system operating under highly variable conditions requiring flexible output. To support the dynamic analysis, the detailed dynamic model and control design of the HTSE process, which employs solid oxide electrolysis cells, have been developed to predict the process behavior over a large range of operating conditions. As first-generation N-R HES technology will be based on LWRs, which provide thermal energy at a relatively low temperature, complementary temperature-boosting technology was suggested for integration with the HTSE process that requires higher temperature input. Simulation results involving several case studies show that the suggested control scheme could maintain the controlled variables (including the steam utilization factor, cathode stream inlet composition, and temperatures of the process streams at various locations) within desired limits under various plant operating conditions. The results also indicate that the proposed HTSE plant could provide operational flexibility to participate in energy management at the utility scale by dynamically optimizing the use of excess plant capacity within an N-R HES. A natural-gas fired GTPP has been proposed as a secondary energy supply to be included in an N-R HES. This auxiliary generator could be used to cover rapid dynamics in grid demand that cannot be met by the remainder of the N-R HES. To evaluate the operability and controllability of the proposed process during transients between load (demand) levels, the dynamic model and control design were developed. Special attention was given to the design of feedback controllers to regulate the power frequency, and exhaust gas and turbine firing temperatures. Several case studies were performed to investigate the system responses to the major disturbance (power load demand) in such a control system. The simulation results show that the performance of the proposed control strategies was satisfactory under each test when the GTPP experienced high rapid variations in the load.

08 HYDROGEN↗

On the adequacy of lightweight thread approaches for high-level parallel programming models

High-level parallel programming models (PMs) are becoming crucial in order to extract the computational power of current on-node multi-threaded parallelism. The most popular PMs, such as OpenMP or OmpSs, are directive-based: the complexity of the hardware is hidden by the underlying runtime system, improving coding productivity. The implementations of OpenMP usually rely on POSIX threads (pthreads), offering excellent performance for coarse-grained parallelism and a perfect match with the current hardware. OmpSs is a task oriented PM based on an ad hoc runtime solution called Nanos++; it is the precursor of the tasking parallelism in the OpenMP tasking specification. A recent trend in runtimes and applications points to leveraging massive on-node parallelism in conjunction with fine-grained and dynamic scheduling paradigms. In this paper we analyze the behavior of the OpenMP and OmpSs PMs on top of the recently emerged Generic Lightweight Threads (GLT) API. GLT exposses a common API for lightweight thread (LWT) libraries that offers the possibility of running the same application over different native LWT solutions. We describe the design details of those high-level PMs implemented on top of GLT and analyze different scenarios in order to assess where the use of LWTs may benefit application performance. Our work reveals those scenarios where LWTs overperform pthread-based solutions and compares the performance between an ad hoc solution and a generic implementation.

97 MATHEMATICS AND COMPUTING↗

Computational Design and Experimental Validation of New Thermal Barrier Systems

The focus of this project is on the development of a reliable and efficient ab initio based computational high temperature material design method which can be used to assist the Thermal Barrier Coating (TBC) bond-coat and top-coat design. Experimental evaluations on the new TBCs are conducted to confirm the new TBCs’ properties. Southern University is the subcontractor on this project with a focus on the computational simulation method development. We have performed ab initio density functional theory (DFT) method and molecular dynamics simulation on screening the top coats and bond coats for gas turbine thermal barrier coating design and validation applications. For experimental validations, our focus is on the hot corrosion performance of different TBC systems. For example, for one of the top coatings studied, we examined the thermal stability of TaZr 2.75 O 8 and confirmed it’s hot corrosion performance.

36 MATERIALS SCIENCE↗

Global Sensitivity Analysis and Estimation of Model Error, Toward Uncertainty Quantification in Scramjet Computations

The development of scramjet engines is an important research area for advancing hypersonic and orbital flights. Progress toward optimal engine designs requires accurate flow simulations together with uncertainty quantification. However, performing uncertainty quantification for scramjet simulations is challenging due to the large number of uncertain parameters involved and the high computational cost of flow simulations. These difficulties are addressed in this paper by developing practical uncertainty quantification algorithms and computational methods, and deploying them in the current study to large-eddy simulations of a jet in crossflow inside a simplified HIFiRE Direct Connect Rig scramjet combustor. First, global sensitivity analysis is conducted to identify influential uncertain input parameters, which can help reduce the system’s stochastic dimension. Second, because models of different fidelity are used in the overall uncertainty quantification assessment, a framework for quantifying and propagating the uncertainty due to model error is presented. In conclusion, these methods are demonstrated on a nonreacting jet-in-crossflow test problem in a simplified scramjet geometry, with parameter space up to 24 dimensions, using static and dynamic treatments of the turbulence subgrid model, and with two-dimensional and three-dimensional geometries.

33 ADVANCED PROPULSION SYSTEMS↗

Networked-based Cyber Analysis using Deep Packet Inspection (DPI) for High-Speed Networks

Today’s networked systems utilize advanced security components such as Next Generation Firewall (NGFW), Intrusion Detection Systems (IDS), Intrusion Prevention Systems (IPS), and methods for network traffic classification. A fundamental aspect of these security components and methods is network packet visibility and packet inspection. To achieve packet visibility, a compute mechanism used by these security components and methods is Deep Packet Inspection (DPI). DPI is used to obtain visibility into packet fields by looking deeper inside packets, beyond just IP address, port, and protocol. However, DPI is considered extremely expensive in terms of compute processing costs and very challenging to implement on high speed network systems. The fundamental scientific paradigm addressed in this research project is the application of greater network packet visibility and packet inspection at data rates greater than 40Gbps to secure computer network systems. The greater visibility and inspection will enable detection of advanced content-based threats that exploit application vulnerabilities and are designed to bypass traditional security approaches such as firewalls and antivirus scanners. Greater visibility and inspection are achieved through identification of the application protocol (e.g., HTTP, SMTP, Skype) and, in some cases, extraction and processing of the information contained in the packet payload. Analysis is then performed on the resulting DPI data to identify potentially malicious behavior. In order to obtain visibility and inspect the application protocol and contents at high speed data rates, advanced DPI technologies and implementations are developed.

97 MATHEMATICS AND COMPUTING↗

An Efficient Sampling-Based Method for Sensitivity and Uncertainty Analysis Through RAVEN

Sensitivity analysis (SA) and uncertainty quantification (UQ) are now widely recognized as essential parts of analyses for complex nuclear systems. Plenty of well-established methodologies are already available for performing SA and UQ, such as Monte Carlo analysis, adjoint-based method, and variance decomposition. The focus of this summary is on developing an efficient sampling-based approach to perform SA and UQ through RAVEN. RAVEN (Risk Analysis and Virtual control Environment) is a software framework able to perform SA and UQ, currently under development at the Idaho National Laboratory. In literature, a number of approaches have been developed to accelerate the process of SA and UQ for large complex systems, such as principal component analysis (PCA) of input space, input-output correlation analysis, and intersection subspace methods. The fundamental idea is to identify the dominant active directions in the parameter space prior to conducting SA and UQ. The premise is that the number of active directions is much less than the dimensionality of input space. In this work, the PCA or Karhunen-Loeve expansion or SVD method is implemented in RAVEN to identify the active directions. This approach is also combined with the forward sampling approaches, such as Monte Carlo, Grid, Latin Hypercube, and Response Surface Design sampling approaches that are currently available in RAVEN, to perform SA and UQ more efficiently. The summary is organized as follows. In the following section, we present the PCA-based forward sampling approach. The effectiveness and efficiency are discussed in the numerical results section. The paper then ends with a concluding summary. In this summary, we have presented PCA-based forward sampling approach to perform SA and UQ. The goal is to reduce the computational overhead for repetitive executions of complex nuclear systems that are generally required when performing SA and UQ. The control of correlations between input parameters is implemented in RAVEN to better represent the input relationships. In addition, we are currently examining our implementations with more realistic models, and we are also focusing on implementation of some new high-dimensional input reduction techniques in RAVEN. (authors)

97 MATHEMATICS AND COMPUTING↗