Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Static code analysis”

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 73 records · Page 4

Time Dependent Data Mining in RAVEN

RAVEN is a generic software framework to perform parametric and probabilistic analysis based on the response of complex system codes. The goal of this type of analyses is to understand the response of such systems in particular with respect their probabilistic behavior, to understand their predictability and drivers or lack of thereof. Data mining capabilities are the cornerstones to perform such deep learning of system responses. For this reason static data mining capabilities were added last fiscal year (FY 15). In real applications, when dealing with complex multi-scale, multi-physics systems it seems natural that, during transients, the relevance of the different scales, and physics, would evolve over time. For these reasons the data mining capabilities have been extended allowing their application over time. In this writing it is reported a description of the new RAVEN capabilities implemented with several simple analytical tests to explain their application and highlight the proper implementation. The report concludes with the application of those newly implemented capabilities to the analysis of a simulation performed with the Bison code.

97 MATHEMATICS AND COMPUTING↗

Transient convective delayed neutron precursor of {sup 235}U for the molten salt reactor experiment

This work investigates the delayed neutron precursors from the fission of {sup 235}U in molten salt reactors. The six delayed neutron groups predict the spatial distribution of the advective transport for radioactive mass transfer in the molten-salt reactor experiment (MSRE). The Mole code, using the framework of the Multiphysics Object-Oriented Simulation Environment (MOOSE), was used to compute the delayed neutron concentration of the whole system. This approach can be used for analysis of reactivity in circulating conditions. The effect of delayed neutron concentrations on the static or dynamic behavior of the system in MSRE was analyzed. Using this approach, the pattern and validation of delayed neutron for a range of parameters of interest are reported and discussed herein. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

RAFT (Response Amplitudes of Floating Turbines) [SWR-22-40]

RAFT is a Python code for frequency-domain analysis of floating wind turbines. RAFT works in the frequency domain to provide efficient computations of a floating wind system's response spectra accounting for platform hydrodynamics, quasi-static mooring reactions, rotor aerodynamics, and linearized control. It can be used to calculate response amplitude operators (RAO's), power spectral densities, mean/static properties, and response metrics based on mean and root-mean-square values of each output. RAFT serves as the Level 1 model in the WEIS (Wind Energy with Integrated Servo-control) toolkit, and features an OpenMDAO wrapper.

Hall, Matthew↗

Screening sensitivity analysis of a PWR fuel assembly FEM structural model

Within the present project, computational methods to describe the static mechanical response of PWR fuel assemblies (FA) in the reactor core are developed. For this purpose a structural model of a PWR fuel assembly has been implemented with the Finite Element Method (FEM) using the ANSYS APDL Mechanical code. The purpose of the model is to investigate the long-term creep deformation of the fuel assemblies in the reactor core, commonly known as fuel assembly bow. FA bow is considered to be a complex process with a large number of influencing mechanisms and many only roughly known boundary conditions. Uncertainty and sensitivity analyses are a common way to investigate the predictive power and reliability of such complex models. Hence, as a first approach to understand the importance of the many model parameters, a qualitative sensitivity analysis of the model was conducted using the screening method proposed by Morris. In the first part of the paper, the FA structural model and its principal parameters are presented. The second part presents first the different simulation experiments used for the screening analyses. Based on the model description, relevant input parameters are then identified and a typical range of values is determined. The sensitivity analysis is performed in two subsequent steps. First, the model is screened for influential stiffness parameters with the objective of stiffening the FA structure. Second, the influence of the uncertain boundary conditions and material models on the FA creep deformation is tested. For the considered FA design it is concluded that the influence of the involved uncertainties is comparable to the influence of the parameters which can increase the stiffness. An integral approach accounting for these uncertainties is hence necessary to predict FA bow reliably in order to take additional measures for a further reduction of FA bow amplitudes. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A New Light Weight Structural Material for Nuclear Structures

Radiation shielding materials are commonly used in nuclear facilities to attenuate the background ionization radiations to a minimum level for creating a safer workplace, meeting regulatory requirements and maintaining high quality performance. The conventional radiation shielding materials have a number of drawbacks: heavy concrete contains a high amount of elements that are not desirable for an effective shielding such as oxygen, silicon, and calcium; a well known limitation of lead is its low machinability and toxicity, which is causing a major environmental concern. Therefore, an effective and environmentally friendly shielding material with increased attenuation and low mass density is desirable. Close-cell composite metal foams (CMFs) and open-cell Al foam with fillers are light-weight candidate materials that we have studied in this project. Close-cell CMFs possess several suitable properties that are unattainable by conventional radiation shielding materials such as low density and high strength for structural applications, high surface area to volume ratio for excellent thermal isolation with an extraordinary energy absorption capability. Open-cell foam is made up of a network of interconnected solid struts, which allows gas or fluid media to pass through it. This unique structure provided a further motive to investigate its application as radiation shields by infiltrating original empty pores with high hydrogen or boron compounds, which are well known for their excellent neutron shielding capability. The resulting open-cell foam with fillers will not only exhibit light weight and high specific surface area, but also possess excellent radiation shielding capability and good processability. In this study, all the foams were investigated for their radiation shielding efficiency in terms of X-ray, gamma ray and neutron. X-ray transmission measurements were carried out on a high-resolution microcomputed tomography (microCT) system. Gamma-emitting sources: 3.0mCi 60 Co, 1.8mCi 137 Cs, 13.5mCi 241 Am, and 5.0mCi 133 Ba were used for gamma-ray attenuation analysis. The evaluations of neutron transmission measurements were conducted at the Neutron Powder Diffractometer beam facility at North Carolina State University. The experimental results were verified theoretically through XCOM and Monte Carlo Z-particle Transport Code (MCNP). A mechanical investigation was performed by means of quasi-static compressive testing. Thermal characterizations were carried out through effective thermal conductivity and thermal expansion analyses in terms of high temperature guarded-comparative-longitudinal heat flow technique and thermomechanical analyzer (TMA), respectively. The experimental results were compared with analytical results obtained from, respectively, Brailsford and Major’s model and modified Turner’s model for verification. Flame test was performed in accordance with United States Nuclear Regulatory Commission (USNRC) standard. CMF sample and a 304L stainless steel control sample were subjected to a fully engulfing fire with an average flame temperature of 800°C for a period of 30 minutes. Finite Element Analysis was conducted to secure the credibility of the experimental results. This research indicates the potential of utilizing the light-weight close-cell CMFs and open-cell Al foam with fillers as shielding material replacing current heavy structures with additional advantage of high-energy absorption and excellent thermal characteristics.

36 MATERIALS SCIENCE↗

Fluoride salt static freezing study

Molten salt reactor designs such as the FHR must consider the safety implications of freezing, requiring the development of thermal hydraulics models that account for freezing behavior. Pipe models that couple heat transfer, fluid dynamics, and phase change physics will require verification and validation. A static freezing study is being conducted to provide an initial set of validation data with LiF-BeF{sub 2} (FLiBe). The study will also be performed with LiF-NaF-KF (FLiNaK) in order to quantify the distortion that arises from using FLiNaK as a non-toxic surrogate salt for thermal-hydraulic studies. In order to develop the numerical model, accurate inputs for the thermal hydraulic properties of the two salts of interest need to be measured near the freezing point. For this reason, a series of three experiments including Differential Scanning Calorimetry (DSC) measurements, the Density Measurement Experiment, and the Static Freezing Experiment are being conducted to determine these thermo-physical properties at temperatures near the freezing point. The Static Freezing Experiment is also designed to determine the degree of supercooling that FLiBe and FLiNaK experience. Supercooling is a temperature depression below the freezing point that the liquid experiences before freezing initiates. Increased supercooling in a liquid flowing through a pipe leads to increased dendritic freezing and an increased likelihood of pipe blockage and freeze-shut. Based on previous research on supercooling in water and non-fluoride salts, supercooling is likely affected most strongly by cooling rate and sample size (as well as salt purity). No data has been generated previously on supercooling in fluoride salts such as FLiBe and FLiNaK. The results of this set of DSC, Density Measurement and Static Freezing experiments will be used in an ongoing fluoride salt freezing study being conducted to develop pipe freezing models for molten salts to be used in safety analysis. The thermal hydraulic properties measured in this study will serve as inputs for all future models. With the supercooling behavior of both salts quantified, it can be assessed whether supercooling must be included in future thermal hydraulic models. The validated Static Freezing Experiment model can then serve as the framework for a pipe freezing model and eventually a 1D system code.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Trace-Driven Memory Access Pattern Recognition in Computational Kernels. In: WOSC '14 Proceedings of the Second Workshop on Optimizing Stencil Computations

Classifying memory access patterns is paramount to the selection of the right set of optimizations and determination of the parallelization strategy. Static analyses suffer from ambiguities present in source code, which modern compilation techniques, such as profile-guided optimization, alleviate by observing runtime behavior and feeding back into the compilation flow. This paper discusses a dynamic analysis technique for recognizing memory access patterns, with application to the stencils domain, and presents our design and C++ implementation using the memory-tracing tool Gleipnir. Finally, we evaluate and discuss the performance and matching capability of our classifiers in the context of the Polybench scientific benchmark suite, which includes both stencil and matrix computations.

Park, Eunjung↗

A Program Analysis Tool to Query Static and Dynamic Characteristics of HPC Applications

Emerging HPC platforms are becoming extraordinarily difficult to program as a result of complex, deep and heterogeneous memory hierarchies, heterogeneous cores, and the need to divide work and data among them. New programming models and libraries are being developed to aid in porting efforts of application developers, but substantial code restructuring is still necessary to fully make use of these new technologies. To do this effectively, these developers need information about their source code characteristics, including static and dynamic (e.g. performance) information to direct their optimisation efforts and make key decisions. On the other hand, system administrators need to understand how users are using the software stack and resources on their systems to understand what software they need to provide to improve the productivity of the users on a platform.In this paper we describe a tool that combines compiler and profiler information to query performance characteristics within an application or across multiple applications on a given platform. Static and dynamic data about applications are collected and stored together in an SQL database that can be queried by either a developer or a system administrator. We will demonstrate the capabilities of this tool via a real-world example from application-driven case studies that aims at understanding the use of scientific libraries on a routine from the molecular simulation application CP2K.

Welch, Aaron↗

Automated Hybrid Variance Reduction on Advanced Architectures in the Shift Monte Carlo Code

Monte Carlo transport methods are the most accurate schemes for solving problems with complex energy and spatial features, but they come with a high computational cost. Although hybrid methods have enabled the use of Monte Carlo transport for a large class of problems, they still require significant computing resources. Modern multicore CPUs with large numbers of compute cores and graphical processing units (GPUs) provide opportunities to optimize the memory and run-time costs of hybrid Monte Carlo methods. This paper documents the development and analysis of three Monte Carlo transport algorithms that support hybrid transport using the consistent adjoint-driven importance sampling (CADIS) and forward-weighted CADIS methods in the Shift Monte Carlo code: history-based transport using static and dynamic threading on multicore CPUs and event-based transport enabling weight window tracking on GPUs. The results are shown for two challenging hybrid problems on the Frontier supercomputer at the Oak Ridge Leadership Computing Facility. The results show that all three methods yield good performance and enable solutions of difficult fixed-source transport problems in less than 2 min on 20 nodes of Frontier. Dynamic threading was observed to give up to 20% better scaling behavior than static threading. Moreover, the AMD Instinct 250X GPU was found to give 9 to 11 times greater throughput per graphics compute die than the best CPU performance. In conclusion, additional opportunities for optimization of hybrid transport on GPUs are discussed.

Denovo↗

Evaluating the effectiveness of program data features for guiding memory management

Recent trends have led to the adoption of larger and more complex memory systems, often with multiple tiers of memory performance within the same platform. To utilize complex memory systems efficiently, current data management strategies must be altered to map usage demands to the underlying hardware. Applications, as the generators of memory accesses, are well-suited to guide this process, but building and maintaining separate source code versions for different memory systems is not feasible in most cases. One potential solution is to employ automated program profiling and analysis to facilitate the production of application-based guidance. By attaching memory usage information to static or lightweight program features, compilers and runtime systems can generate fine-grained guidance without additional efforts from users or developers. Recent works have employed this approach with some success, but it is not clear which program features are most useful for guiding data management.This work evaluates the effectiveness of using different program data features to predict memory usage and guide memory management. It employs a custom set of simulation tools, based in the Intel® Pin framework, to collect and analyze the usage characteristics of application data associated with common program data features, such as allocation sites, types, and context. The results show that even relatively simple features, such as object size, are useful for selecting data with similar usage properties, but finer-grained features, such as the instructions that access a particular object, are often much more effective. Additionally, this work evaluates the performance of using different data features and different program inputs to guide data placement on a heterogeneous memory platform with a limited amount of high performance memory.

Effler, T. Chad↗

ROSE

Developed at Lawrence Livermore National Laboratory (LLNL), ROSE is an open source compiler infrastructure to build source-to-source program transformation and analysis tools for large-scale C (C89 to C23), C++ (C++98 to C++23), UPC, Fortran (Fortran4, 66, 77, 95, 2003), OpenMP, Java, Python, and Binary applications. ROSE users range from experienced compiler researchers to library and tool developers who may have minimal compiler experience. ROSE is particularly well suited for building custom tools for static analysis, program optimization, arbitrary program transformation, domain-specific optimizations, complex loop optimizations, performance analysis, and cyber-security. ROSE is: A library (and set of associated tools) to quickly and easily apply compiler techniques to one's code in order to improve application performance and developer productivity. A research and development compiler infrastructure for for writing custom source-to-source translators to perform source code transformations, analysis, and optimizations. Is

Pinnow, NathanT [Lawrence Livermore National Labor↗

HPC-enabled computation of demand models at scale

The purpose of this project is to examine the energy impact of urban-scale traffic for the Los Angeles Basin by developing and implementing a scalable traffic assignment model. An energy optimization function will be posed and when integrated into the optimization code for travel assignment it can be mathematically proven to converge. The energy optimization function can then be compared to the typical travel time optimization that is traditionally used in traffic assignment models. The analysis will begin with static traffic assignment models with the routing for all origin and destinations computed in parallel on high performance computing facilities. Convergence of the numerical methods rely on the solution of convex programs (or extensions of these). This step will mostly consist of demonstrating the ability to parallelize the Frank Wolfe algorithm on various platforms.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

HPC4Mobilty w/ UCB

The purpose of this project is to examine the energy impact of urban-scale traffic for the Los Angeles Basin by developing and implementing a scalable traffic assignment model. An energy optimization function will be posed and when integrated into the optimization code for travel assignment it can be mathematically proven to converge. The energy optimization function can then be compared to the typical travel time optimization that is traditionally used in traffic assignment models. The analysis will begin with static traffic assignment models with the routing for all origin and destinations computed in parallel on high performance computing facilities. Convergence of the numerical methods rely on the solution of convex programs (or extensions of these). This step will mostly consist of demonstrating the ability to parallelize the Frank Wolfe algorithm on various platforms. This work will contribute to LBNL’s efforts to develop new processes, analytical tools, program designs, and business models to advance the state of the art in next-generation sustainable transportation solutions.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Expanded capabilities of the CarMa code in modeling resistive wall mode dynamics with 3-D conductors

In this work, an improved version of the CarMa code is presented, called CarMa-D, for the analysis of resistive wall modes (RWMs) in fusion devices, simultaneously considering the effects of volumetric three-dimensional conducting structures, plasma dynamics, toroidal rotation or drift-kinetic damping. Unlike static CarMa, CarMa-D does not rely on the simplifying assumptions such as neglecting the plasma mass, toroidal rotation and kinetic damping. The new coupling strategy is based on matrix-based Padé rational functions approximation of plasma a response. The arising mathematical model is formally equivalent to the original CarMa model, but with a higher number of degrees of freedom to model the dynamics of the plasma. CarMa-D overcomes the main limitations of the original CarMa, in particular: (i) the massless assumption for the plasma is removed, allowing modeling of global modes growing both on ideal kink time scales and in the typical RWM growth rate regime, with a suitable treatment of the model; (ii) the effects of toroidal plasma flow and drift kinetic damping can be included into the new model, providing a powerful tool to study macroscopic phenomena where both plasma dynamics and 3-D conducting structures play important roles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Creating an Interprocedural Analyst-Oriented Data Flow Representation for Binary Analysts (CIAO)

National security missions require understanding third-party software binaries, a key element of which is reasoning about how data flows through a program. However, vulnerability analysts protecting software lack adequate tools for understanding data flow in binaries. To reduce the human time burden for these analysts, we used human factors methods in a rolling discovery process to derive user-centric visual representation requirements. We encountered three main challenges: analysis projects span weeks, analysis goals significantly affect approaches and required knowledge, and analyst tools, techniques, conventions, and prioritization are based on personal preference. To address these challenges, we initially focused our human factors methods on an attack surface characterization task. We generalized our results using a two-stage modified sorting task, creating requirements for a data flow visualization. We implemented these requirements partially in manual static visualizations, which we informally evaluated, and partially in automatically generated interactive visualizations, which have yet to be integrated into workflows for evaluation. Our observations and results indicate that 1) this data flow visualization has the potential to enable novel code navigation, information presentation, and information sharing, and 2) it is an excellent time to pursue research applying human factors methods to binary analysis workflows.

42 ENGINEERING↗

Performance Improvements for the Griffin Transport Solvers

Griffin is a Multiphysics Object-Oriented Simulation Environment based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory. Griffin includes a variety of deterministic radiation transport solvers for fixed source, k-eigenvalue, adjoint, and subcritical multiplication, as well as transient solvers for point-kinetics, improved quasi-static, and spatial dynamics. A code assessment performed in FY-20 identified two significant issues with the transport solvers in Griffin: first, the primary heterogeneous SN (discrete ordinates) transport solver based on continuous finite element methods required significant mesh refinement and higher memory usage compared to solvers based on the method of characteristic for equivalent accuracy. Second, the homogeneous PN (spherical harmonics expansion) transport solver did not adequately support polynomial refinement, which is a feature usually required for problems with spatial homogenization and pronounced streaming, typical in fast or gas-cooled reactor systems. To address the first issue, the development effort focused on the more promising discontinuous finite element method (DFEM)-based SN transport solver in Griffin. The addition of an asynchronous parallel transport sweeper and coarse mesh finite difference (CMFD) acceleration have rendered a superior heterogeneous SN transport capability for multiphysics problems that requires far less computing resources in terms of both CPU time and memory usage. This is demonstrated with typical thermal- and fast-spectrum reactor benchmark problems, including 2D Transient Reactor Test, 3D Advanced Burner Test Reactor (ABTR), and 2D and 3D Empire microreactor. For the second issue, the development effort focused on a new transport solver based on the hybrid finite element PN method (HFEM-PN), equivalent to the variational nodal method, as well as a new diffusion solver based on HFEM-Diffusion. This solver is intended for homogenized domains with multiphysics coupling (i.e., supports mesh displacement, seamless temperature feedback, etc.). Initial calculations with the HFEM-Diffusion implementation show very good parallel efficiency for the residual evaluations with the 2D ABTR benchmark. A future development effort will be centered on further improvements to the CMFD, HFEM-PN, and DFEM diffusion solvers to ensure Griffin meets performance and software quality assurance requirements for advanced reactor design and analysis.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A Research Roadmap for Computation-Based Human Reliability Analysis

The United States (U.S.) Department of Energy (DOE) is sponsoring research through the Light Water Reactor Sustainability (LWRS) program to extend the life of the currently operating fleet of commercial nuclear power plants. The Risk Informed Safety Margin Characterization (RISMC) research pathway within LWRS looks at ways to maintain and improve the safety margins of these plants. The RISMC pathway includes significant developments in the area of thermalhydraulics code modeling and the development of tools to facilitate dynamic probabilistic risk assessment (PRA). PRA is primarily concerned with the risk of hardware systems at the plant; yet, hardware reliability is often secondary in overall risk significance to human errors that can trigger or compound undesirable events at the plant. This report highlights ongoing efforts to develop a computation-based approach to human reliability analysis (HRA). This computation-based approach differs from existing static and dynamic HRA approaches in that it: (i) interfaces with a dynamic computation engine that includes a full scope plant model, and (ii) interfaces with a PRA software toolset. The computation-based HRA approach presented in this report is called the Human Unimodels for Nuclear Technology to Enhance Reliability (HUNTER) and incorporates in a hybrid fashion elements of existing HRA methods to interface with new computational tools developed under the RISMC pathway. The goal of this research effort is to model human performance more accurately than existing approaches, thereby minimizing modeling uncertainty found in current plant risk models.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

The RISMC Approach to Perform Advanced PRA Analyses

The existing fleet of nuclear power plants is in the process of extending its lifetime and increasing the power generated from these plants via power uprates. In order to evaluate the impact of these two factors on the safety of the plant, the RISMC Pathway aims to develop simulation-based tools and methods to assess risks for existing nuclear power plants in order to optimize safety. This pathway, by developing new methods, is extending the state-of-the-practice methods that have been traditionally based on logic structures such as Event-Trees and Fault-Trees. These static types of models mimic system response in an inductive and deductive way respectively, yet are restrictive in the ways they can represent spatial and temporal constructs. RISMC analyses are performed by using a combination of thermal-hydraulic codes and a stochastic analysis tool currently under development at the Idaho National Laboratory, i.e. RAVEN. This paper presents a case study in order to show the capabilities of the RISMC methodology to assess impact of power uprate of a boiling water reactor system during a station blackout accident scenario. We employ the system simulator code, RELAP5-3D, coupled with RAVEN which perform the stochastic analysis. Our analysis is in fact performed by: 1) sampling values of a set of parameters from the uncertainty space of interest, 2) simulating the system behavior for that specific set of parameter values and 3) analyzing the set of simulation runs. Results obtained give a detailed investigation of the issues associated with a plant power uprate including the effects of station blackout accident scenarios. We are able to quantify how the timing of specific events was impacted by a higher nominal reactor core power. Such safety insights can provide useful information to the decision makers to perform risk informed margins management.

97 MATHEMATICS AND COMPUTING↗