Search NASA⌕ Search

SEARCH · Search NASA

Results for “C Code”

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,243 records · Page 69

Application of performance portability solutions for GPUs and many-core CPUs to track reconstruction kernels

Next generation High-Energy Physics (HEP) experiments are presented with significant computational challenges, both in terms of data volume and processing power. Using compute accelerators, such as GPUs, is one of the promising ways to provide the necessary computational power to meet the challenge. The current programming models for compute accelerators often involve using architecture-specific programming languages promoted by the hardware vendors and hence limit the set of platforms that the code can run on. Developing software with platform restrictions is especially unfeasible for HEP communities as it takes significant effort to convert typical HEP algorithms into ones that are efficient for compute accelerators. Multiple performance portability solutions have recently emerged and provide an alternative path for using compute accelerators, which allow the code to be executed on hardware from different vendors. We apply several portability solutions, such as Kokkos, SYCL, C++17 std::execution::par and Alpaka, on two mini-apps extracted from the mkFit project: p2z and p2r. These apps include basic kernels for a Kalman filter track fit, such as propagation and update of track parameters, for detectors at a fixed z or fixed r position, respectively. The two mini-apps explore different memory layout formats. We report on the development experience with different portability solutions, as well as their performance on GPUs and many-core CPUs, measured as the throughput of the kernels from different GPU and CPU vendors such as NVIDIA, AMD and Intel.

Kwok, Ka Martin↗

CCPDVANALYSIS Documentation

This document contains all relevant code and descriptions of working components of CCPDVANALYSIS software toolkit.

97 MATHEMATICS AND COMPUTING↗

Effects of Geometric Discontinuities on Creep Behavior of Alloy 617

Due to the excellent long-term high-temperature performance of Alloy 617, the Alloy has been identified as the primary structural material for the intermediate heat exchanger of the Very High Temperature Reactor (VHTR). Both the base and welded alloy has been qualified to be used up to 950 °C up to 100,000 hours by the American Society of Mechanical Engineer (ASME) Boiler and Pressure Vessel Code (BPVC). The qualification properties of the material were exclusively determined by uniaxial creep tests. The components may, however, experience multiaxial stress conditions in service due to the reactor geometry. This study used notched creep tests to investigate the effect of geometric discontinuities on creep behavior of both the base and welded Alloy 617. The V-notch tests reveal that the geometric discontinuities do not reduce the creep rupture live under high and intermediate stress conditions. The U-notch tests show that a larger radius notch behave in a similar way as a straight gauge section of a uniaxial test. The creep rupture live of the welded alloy was slightly longer than the base metal.

36 - MATERIALS SCIENCE↗

Models to Incorporate Reaction Mechanisms into DG-OSPREY: Fixed-Bed Simulations for Organic Iodides Capture Using Ag 0 Z

As one of the most potent radioisotopes released during spent nuclear fuel reprocessing, 129 I is strictly regulated and must be removed before discharge. Organic iodides (primarily alkyl iodides with different chain lengths, i.e., CH 3 I, C 4 H 9 I, and C 12 H 25 I) comprise ~2% of the total iodine in the reprocessing off-gases and are primarily present in vessel off-gas (VOG). Reduced silver mordenite (Ag 0 Z) is predominantly considered for the removal of radioiodine; however, its capture performance and underlying interaction processes with long-chain organic iodides are not fully understood. Two major tasks were accomplished in this study. First, to improve upon the previous experimental studies where Ag 0 Z was used to capture CH 3 I, C 4 H 9 I, and C 12 H 25 I at different concentrations, we comprehensively investigated the corresponding capture mechanisms by characterizing fully loaded Ag 0 Z samples. Second, computational codes were implemented to perform fixed-bed simulations that account for transport and reaction mechanisms. Scanning electron microscopy with energy dispersive X-ray analysis (SEM-EDX), powder X-ray diffraction (PXRD), UV-visible diffuse reflectance spectroscopy (UV-vis DRS), and thermogravimetric analysis (TGA) were conducted on Ag 0 Z samples that are saturated with I 2 , CH 3 I, C 4 H 9 I, and C 12 H 25 I), respectively. Results indicate that AgI is the predominant adsorption product regardless of the adsorbed iodine species, yet alkyl iodides with different carbon chain lengths may have different compositions of α- and γ-AgI. Synchrotron pair distribution function (PDF) measurements and TGA coupled with a Fourier transformed infrared detector (TGA FTIR) have been performed, and experimental data are currently analyzed. Results are expected to provide further insights into the adsorption mechanisms. The fixed-bed capture performance for CH 3 I was successfully simulated by solving mathematical equations that describe the underlying transport processes and adsorption reactions. The computational framework, Catalytic After Treatment System (CATS), that was originated in our research group, was used to solve the governing equations. Kinetic parameters, including the pore diffusivity and reaction rate constant were obtained by optimization techniques using data from thin-bed experiments performed at Oak Ridge National Laboratory. The performance of a deep bed predicted using the optimized parameters showed promising agreement with the experimental data.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

FOCUS: A full-orbit CUDA solver for particle simulations in magnetized plasmas

A Full-Orbit CUDA Solver for particle simulations in plasmas, FOCUS, has been developed. The code follows exact particle trajectories by solving Newton’s equations with the Lorentz force. The code can use fields calculated analytically or numerically and it is also possible to use the information provided by equilibrium reconstruction and transport codes. FOCUS has an elastic collisions module which covers the whole particle energy range in magnetic fusion devices. Moreover, an atomic collision module was also included to simulate the interaction of the test particles with neutral or partially ionized species. Regarding performance, the main feature is that the code runs on GPUs, allowing the simulation of a large number of particles using moderate computational resources. Different versions of the code have been used in several papers. Here we present a complete description of its capabilities, the basis physics included and detailed information about the numerical algorithms employed. Finally, some examples are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An independent analysis of bias sources and variability in wind plant pre-construction energy yield estimation methods

The wind resource assessment community has long had the goal of reducing the bias between wind plant pre-construction energy yield assessment (EYA) and the observed annual energy production (AEP). This comparison is typically made between the 50% probability of exceedance (P50) value of the EYA and the long-term corrected operational AEP (hereafter OA P50), and is known as the P50 bias. The industry has critically lacked an independent analysis of bias reduction investigated across multiple consultants to identify the greatest sources of uncertainty and variance in the EYA process and the best opportunities for uncertainty reduction. The present study addresses this gap by benchmarking consultant methodologies against each other and against operational data at a scale not seen before in industry collaborations. We consider data from 10 wind plants and evaluate discrepancies between eight consultancies in the steps taken from estimates of gross to net energy. Consultants tend to overestimate the gross energy produced at the turbines and then compensate by further overestimating downstream losses, leading to a mean P50 bias near zero, still with significant variability among the individual wind plants. Within our data sample, we find that consultant estimates of all loss categories, except environmental losses, tend to reduce the project-to-project variability of the P50 bias. The disagreement between consultants, however, remains flat throughout the addition of losses. Finally, we find that differences in consultants’ estimates of project performance can lead to differences up to $10/MWh in the levelized cost of energy for a wind plant.

Todd, Austin C.↗

Fatigue and Creep-Fatigue Evaluation of Alloy 709 at 760 and 816°C

A significant research and development effort is underway to support the qualification of Alloy 709 as a Class A construction material in the American Society of Mechanical Engineers (ASME) Boiler and Pressure Vessel Code, Section III, Division 5, High Temperature Reactors. This initiative includes a comprehensive Alloy 709 code qualification plan aimed at generating extensive material testing data crucial for compiling the code case data package. The data package is essential in establishing material-specific design parameters for Alloy 709 to be used as Section III, Division 5 Class A construction material for fast reactors, molten salt reactors and gas-cooled reactors. An ASME Section III, Division 5 material code case requires the evaluation of mechanical properties from a minimum of three commercial heats, covering anticipated compositional ranges. A key part of the data package involves fatigue and creep-fatigue testing at elevated temperatures, needed for developing the fatigue design curves and the damage envelope of the creep-fatigue interaction diagram (D-diagram). This paper summarizes the strain-controlled fatigue testing on three commercial heats of Alloy 709 at 760 and 816°C with strain ranges between 0.25% and 3%. The fatigue failure data are used to generate a preliminary fatigue design curve. Additionally, the creep-fatigue testing results at 816°C with tensile hold times of 10, 30, and 60 minutes are presented in support of developing the D-diagram for Alloy 709.

Wang, Yanli↗

Innovative Method for Welding in Generation 3 CSP to Enable Reliable Manufacturing of Solar Receivers to withstand Daily Cycling at Temperatures Above 700°C (Final Technical Report)

Inconel® Alloy 740H® (alloy 740H) was the first age-hardenable nickel-based alloy approved by the ASME Boiler & Pressure Vessel Code for use in pressure-boundary applications. Over the past ~20 years the alloy has been optimized for weldability and high-temperature stability, approved for use in different applications. Development of a supply chain combined with the advantageous properties of the alloy (high-temperature creep strength, oxidation and corrosion resistance, etc.) have resulted in the alloy being applied to new high-temperature power cycle demonstration projects, and of particular interest are applications to concentrating solar power (CSP) to enable higher-efficiency Generation 3 CSP systems and the corresponding supercritical CO 2 (sCO 2 ) power cycle components (heat exchangers, piping, etc.). The high allowable stresses of alloy 740H also make it a desirable material for current Generation 2 CSP solar power receivers to improve cyclic capability and/or reduce receiver height. Recent experiences in demonstration projects utilizing alloy 740H identified cracking issues during welding and fabrication. In this project, a detailed study was done to confirm and clarify the Stress Relaxation Cracking (SRxC) mechanism, also known as stress relief cracking or strain-age cracking (SAC), during post-weld heat-treatment (PWHT). This involved detailed microscopy and advanced characterization to understand the root cause(s) of three failures obtained from industry. Based in-part on these findings, a targeted laboratory based SRxC test method was utilized to evaluate variables such as heat-to-heat variations, strain level, PWHT temperature, and starting material condition on three heats of alloy 740H. Industrial shop welding of cold-worked plates was also conducted. The research showed the following: SRxC was confirmed as the cracking mechanism for all field failures; Stress state (from residual stresses, constraint, deformation, and local stress concentrations) was playing a significant role in field failures and laboratory testing confirmed increasing susceptibility for all heats with increasing strain levels. High levels of microstructural strain were identified at crack initiation locations, in some cases leading to local recrystallization; Precipitate free zones (PFZs) at grain boundaries were found at relaxation cracks and crack initiation locations uniquely associated with SRxC in alloy 740H. Laboratory testing reproduced this microstructural feature which had only previously been reported in long-term creep testing of weldments. Advanced nano-scale characterization confirmed the presence of a moving boundary leading to coarsening of precipitates and PFZs where damage accumulated; The research suggested heat-to-heat variations due to local chemistry and processing may influence SRxC susceptibility, but more work is needed to fully clarify these effects. To disseminate the key learnings from this research to the scientific and engineering communities and alloy 740H end users, multiple technical publications and presentations were made, an industrial alloy 740H users meeting was held, and a new industry guideline specification document which can be directly implemented by end-users of alloy 740H was produced.

14 SOLAR ENERGY↗

3-D Deformable Grid Electomagnetic PIC for Parallel Computers

We have developed a three-dimensional non-orthogonal grid electromagnetic PIC code for parallel supercomputers. The numerical formulation is based on hexahedral cells which are logically connected cubic cells but distorted to body-fit complex geometries in physical space. An outstanding area of research is to improve the accuracy of the electromagnetic field solve on distorted meshes.

3D deformable non-orthogonal grid code PIC codes p↗

WHONDRS Summer 2019 Sampling Campaign: Global River Corridor Sediment FTICR-MS, Dissolved Organic Carbon, Aerobic Respiration, Elemental Composition, Grain Size, Total Nitrogen and Organic Carbon Content, Bacterial Abundance, and Stable Isotopes (v10)

The WHONDRS Summer 2019 Sampling (S19S) Campaign collected samples in 97 globally distributed river corridor systems between July and September 2019. Surficial streambed sediments were collected at three locations within each site (upstream, midstream, and downstream). Surface water was collected at the downstream site. This dataset includes a portion of the data types produced from the sediment samples and does not include any results from the surface water. Surface water data can be found at https://data.ess-dive.lbl.gov/view/doi:10.15485/1603775. The S19S campaign was designed with the science community to ask questions associated with links among core/transient metabolomes, microbial metabolism, biogeochemical function, and physical properties of watershed and river corridor systems. This data package was originally published in December 2020. It was updated in April 2021 (v2; modified files), May 2022 (v3; new and modified files), November 2022 (v4; new and modified files), January 2023 (v5 and v6; new and modified files), February 2023 (v7; new and modified files), August 2023 (v8; new and modified files), July 2025 (v9; new and modified files), and June 2026 (v10; new and modified files). Version 10 includes major restructuring of the files and their formats. See the change history in the readme for details. In addition to a readme, this data package also includes a file-level metadata (FLMD) file that describes each file and a data dictionary (DD) that describes all column/row headers and variable definitions. This dataset contains (1) readme, (2) data dictionary, (3) file level metadata, (4) sampling protocol, (5) a data checks report, (6) field metadata, (7) international generic sample number (IGSN) mapping file, (8) joint genome institute (JGI) mapping file, (9) a folder of site photos, and (10) a folder of sample data. The sample data folder contains (1) dissolved organic carbon (DOC, measured as non-purgeable organic carbon, NPOC) and FTICR methods; (2) respiration rates and normalized respiration rates calculated from laboratory sediment incubations, (3) X-ray fluorescence (XRF) elemental composition, (4) grain size distribution from the < 2 millimeter fraction of sediment, (5) total organic carbon and total nitrogen content, (6) bacterial abundance flow cytometry data, (7) sediment stable isotopes (C and N), (8) laboratory methods, (9) water and sediment mass and volume used in incubations, (10) a folder with code, inputs, and plots for the respiration rates, (11) a folder containing XRF data, (12) a folder containing unaligned, peak-picked high resolution characterization of dissolved organic matter from sediment via 12 Tesla Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS) through the Environmental Molecular Sciences Laboratory (EMSL; https://www.pnnl.gov/environmental-molecular-sciences-laboratory), (12) and a folder containing processed FTICR-MS data. All files are .msa, .txt, .spx, .xml, .csv, .pdf, .txt, .jpg, or .jpeg.

13C↗

Low level coupling scheme between neutronics and thermal-hydraulics based on Anderson acceleration

The simulation of nuclear reactors is a multiphysics problem mixing, amongst other fields, neutron transport and thermal-hydraulics. The simplest and most used approach in multiphysics simulation is based on the coupling of single-physics codes in a black-box fashion. However, in order to reduce the computational time needed for such simulations, case-dependent optimizations are often required. In this paper, we aim at reducing the computational time required to solve a coupled neutronic/thermal-hydraulic steady-state problem on a simplified Pressurized Water Reactor (PWR) core. The idea is to deal simultaneously with the coupling of the energy groups of the deterministic neutronic description of the core and its thermal-hydraulic description with the Anderson acceleration. By doing so, the fission source terms are directly accelerated instead of the power map as done in most cases. The power method used to solve the k-eigenvalue problem inside the neutronic solver is thus accelerated with the Anderson acceleration. The numerical experimentations conducted in this work are performed using APOLLO3 and THEDI, and indicate that such coupling strategy improves the convergence rates in terms of number of iterations required and the total computational time. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Fabrication and evaluation of cold/formed/weldbrazed beta-titanium skin-stiffened compression panels

The room temperature and elevated temperature buckling behavior of cold formed beta titanium hat shaped stiffeners joined by weld brazing to alpha-beta titanium skins was determined. A preliminary set of single stiffener compression panels were used to develop a data base for material and panel properties. These panels were tested at room temperature and 316 C (600 F). A final set of multistiffener compression panels were fabricated for room temperature tests by the process developed in making the single stiffener panels. The overall geometrical dimensions for the multistiffener panels were determined by the structural sizing computer code PASCO. The data presented from the panel tests include load shortening curves, local buckling strengths, and failure loads. Experimental buckling loads are compared with the buckling loads predicted by the PASCO code. Material property data obtained from tests of ASTM standard dogbone specimens are also presented.

Royster, D. M.↗

Funtools: Fits Users Need Tools for Quick, Quantitative Analysis

The Funtools project arose out of conversations with astronomers about the decline in their software development efforts over the past decade. A stated reason for this decline is that it takes too much effort to master one of the existing FITS libraries simply in order to write a few analysis programs. This problem is exacerbated by the fact that astronomers typically develop new programs only occasionally, and the long interval between coding efforts often necessitates re-learning the FITS interfaces. We therefore set ourselves the goal of developing a minimal buy-in FITS library for researchers who are occasional (but serious) coders. In this case, "minimal buy-in" meant "easy to learn, easy to use, and easy to re-learn next month". Based on conversations with astronomers interested in writing code, we concluded that this goal could be achieved by emphasizing two essential capabilities. The first was the ability to write FITS programs without knowing much about FITS, i.e., without having to deal with the arcane rules for generating a properly formatted FITS file. The second was to support the use of already-familiar C/Unix facilities, especially C structs and Unix stdio. Taken together, these two capabilities would allow researchers to leverage their existing programming expertise while minimizing the need to learn new and complex coding rules.

Mandel, Eric↗

DUKSUP: A Computer Program for High Thrust Launch Vehicle Trajectory Design and Optimization

From the late 1960's through 1997, the leadership of NASA's Intermediate and Large class unmanned expendable launch vehicle projects resided at the NASA Lewis (now Glenn) Research Center (LeRC). One of LeRC's primary responsibilities --- trajectory design and performance analysis --- was accomplished by an internally-developed analytic three dimensional computer program called DUKSUP. Because of its Calculus of Variations-based optimization routine, this code was generally more capable of finding optimal solutions than its contemporaries. A derivation of optimal control using the Calculus of Variations is summarized including transversality, intermediate, and final conditions. The two point boundary value problem is explained. A brief summary of the code's operation is provided, including iteration via the Newton-Raphson scheme and integration of variational and motion equations via a 4th order Runge-Kutta scheme. Main subroutines are discussed. The history of the LeRC trajectory design efforts in the early 1960's is explained within the context of supporting the Centaur upper stage program. How the code was constructed based on the operation of the Atlas/Centaur launch vehicle, the limits of the computers of that era, the limits of the computer programming languages, and the missions it supported are discussed. The vehicles DUKSUP supported (Atlas/Centaur, Titan/Centaur, and Shuttle/Centaur) are briefly described. The types of missions, including Earth orbital and interplanetary, are described. The roles of flight constraints and their impact on launch operations are detailed (such as jettisoning hardware on heating, Range Safety, ground station tracking, and elliptical parking orbits). The computer main frames on which the code was hosted are described. The applications of the code are detailed, including independent check of contractor analysis, benchmarking, leading edge analysis, and vehicle performance improvement assessments. Several of DUKSUP's many major impacts on launches are discussed including Intelsat, Voyager, Pioneer Venus, HEAO, Galileo, and Cassini.

high thrust trajectory design↗

How initial conditions-, structural-, and parameter-based model uncertainty interact and influence predictions in permafrost ecosystems: Modeling Archive

This dataset contains model output and input data, as well as source code examples for the Terrestrial Ecosystem Model with the Dynamic Vegetation Model and Dynamic Organic Soil (DVM-DOS-TEM) for the field sites Imnavait creek and the Bonanza creek Long Term Ecological Research Network (LTER). The data covers simulations from the last glacial maximum (LGM) until 2100 for a selection of paleo scenarios, setting the mean temperature of the LGM up to 10°C lower than pre-industrial conditions. The model structure was modulated to represent various model versions, and this dataset contains the relevant changes in the source code. The raw output data, the processed statistical data, the setup and processing scripts as well as parameter value distribution files from a parameter sensitivity analysis are included as well. Model outputs include active layer depth, organic soil carbon, soil layer depths, gross primary productivity (GPP) with and without nitrogen limitation, net primary productivity (NPP), soil liquid water content, heterotrophic, maintenance, and growth respiration, soil temperature, and vegetation carbon (*.nc files). The Next-Generation Ecosystem Experiments in the Arctic (NGEE Arctic) project is a research effort to reduce uncertainty in the Department of Energy’s Energy Exascale Earth System Model (E3SM) by developing a predictive understanding of Arctic tundra ecosystems underlain by permafrost and to quantify feedbacks from the Arctic tundra to the Earth system. NGEE Arctic is supported by the Department of Energy's Office of Biological and Environmental Research.Over Phases 1–3, observations made by the NGEE Arctic team across a gradient of permafrost landscapes in Arctic Alaska improved the representation of tundra processes in the land surface component of E3SM (the E3SM Land Model, ELM). Model improvements emphasized unique aspects of permafrost environments and explored reductions in model complexity while retaining predictive power. The Arctic-informed ELM developed by NGEE Arctic has been used to make novel predictions on processes ranging from permafrost thaw to soil biogeochemical cycling to Earth system feedbacks associated with the unique characteristics of tundra plants. In Phase 4, the NGEE Arctic team is evaluating our new predictive understanding under novel conditions across the Arctic domain. In collaboration with partners at long-term pan-Arctic research sites we are examining whether an Arctic-informed ELM can faithfully simulate interactions among surface and subsurface processes at site, regional, and pan-Arctic scales. In turn, we are using variety of tools to dynamically extend and evaluate ELM inference, with an emphasis on data synthesis and pan-Arctic model evaluation, reintegration of code with an evolving E3SM, scaling across heterogeneous Arctic landscapes, and the appropriate representation of the impacts of increasingly frequent Arctic disturbances.

54 ENVIRONMENTAL SCIENCES↗

Helical core formation and evolution during current ramp-up in the high-field tokamak Alcator C-Mod

We report that large, spontaneous m/n = 1/1 helical cores are predicted in tokamaks with extended regions of low- or reversed-magnetic shear profiles in a region within the q = 1 surface and an onset condition determined by constant (dp/dρ)/B$2\atop{t}$ along the threshold. These 3D modes occurred frequently in Alcator C-Mod during ramp-up when slow current penetration results in a reversed shear q-profile. The onset and early development of a helical core in C-Mod were simulated using a new 3D time-dependent equilibrium reconstruction, based on the ideal MHD equilibrium code VMEC. The reconstruction used the experimental density, temperature, and soft-X-ray fluctuations. The pressure profile can become hollow due to an inverted, hollow electron temperature profile caused by molybdenum radiation in the plasma core during the current ramp-up phase before the onset of sawteeth, which may also occur in ITER with tungsten. Lastly, based on modeling, it is found that a reverse shear q-profile combined with a hollow pressure profile reduces the onset condition threshold, enabling helical core formation from an otherwise axisymmetric equilibrium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Final Report on CFD and Thermal-Mechanical Stress Analysis of PWR Surge Line under Transient Condition Thermal Stratification and an Evolutionary Cyclic Plasticity Based Transformative Fatigue Evaluation Approach without Using S~N Curve: Rev. 1

This report presents an evolutionary cyclic plasticity based transformative approach for fatigue evaluation of safety critical components. This approach is generic which can be used not only for the present generation light water reactor components (which is our focus) but also for advanced reactor components. The approach is based on the fundamental concept of time-evolution of progressive fatigue damage rather than the end-of-life data based conventional S~N curve approaches. Series of fatigue tests on 316 stainless steel uniaxial specimens were conducted under strain-controlled and stress-controlled cyclic loading in two different environment such as in-air at 300 °C and PWR coolant water condition at 300 °C. Material models are developed to capture the uniaxial test data and are subsequently programmed into commercially available ABAQUS code to transform the uniaxial material behavior to multiaxial loading domain. The developed evolutionary cyclic-plasticity approach is verified by both analytical and 3D-FE modeling of uniaxial fatigue test specimens. Results show that the developed material models can capture the time-dependence of material hardening and softening with great accuracy. In addition, results show that the material models not only can capture the material behavior under constant amplitude loading but also can capture the load-sequence effect under variable and random amplitude loading. The ANL developed approach is one of its kind, which shows first in the fatigue research community that fatigue evaluation of a component (in our case verified for laboratory scale fatigue test specimens with great accuracy) can be performed for its entire fatigue life, including the hardening and softening behavior, without using a conventional S~N curve based approach. Additionally results are presented related to fatigue evaluation of real reactor components such as of a pressurized water reactor (PWR) surge line (SL), under upset transient conditions. SLs of PWR experience complex thermal stratification due to the mixing of differential temperature flow from hot leg and pressurizer. A detailed parametric study of thermal stratification in SL is performed through both standalone computational fluid dynamics (CFD) analysis and fluid-structure-interaction based coupled CFD and heat transfer analysis. Results show that the mass flow rate in SL substantially affects the extent of its thermal stratification. The resulting thermal boundaries conditioned are used as inputs for thermal-mechanical stress analysis and component level fatigue evaluation. Fatigue evaluation of SL pipe under example loading conditions was performed using the above-mentioned evolutionary-plasticity based approach and ASME code (for in-air life) and NUREG-6909 (for PWR environment life) based approaches. Results show that the evolutionary plasticity based approach estimate substantially different life compared to the conventional ASME based approach. To show the entirety of the proposed evolutionary fatigue-modeling framework, this report includes the results presented in semi-annual March report (ANL/LWRS-17/01). Sections 2-5 of March report are directly included as Sections 2-5 in this report, without any significant changes. However, Sections 6-9 of this report present the work performed during March-September of 2017. Whereas, Sections 1 and 10 present the organization of this report and the summary and future work, respectively. This report (ANL/LWRS-17/03 Rev. 1) is a slightly revised version of the earlier published report (ANL/LWRS-17/03). In this revised version, some of the results and associated discussion have been corrected/updated.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Design and Evaluation of a Pixelated PSD-capable Scintillator Detector with SiPM Readout

We present the detector response comparison between a 10x10 pixellated array of scintillator read out with Anger logic using four 2" Hamamatsu R7724-100 super bialkali photomultiplier tubes (PMT) and a custom Silicon photomultiplier (SiPM) board consisting of 100 C-series 6x6 mm SiPMs from SensL. An array of these pixellated detectors are currently used in the Neutron Coded Aperture (NCA) imaging system. The energy, timing and pulse shape discrimination response using both readout schemes are presented, along with an analysis of multiple scatter events occurring within the block. An evaluation of the impact of photodetector readout on the overall detection efficiency and imaging accuracy is presented.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗