Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “analysis codes”

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 307 records · Page 17

A Data Analysis Center for Electromagnetic and Hadronic Interaction. Products of the DAC members

The Data Analysis Center (DAC) of the Center for Nuclear Studies (CNS) at the George Washington University (GW) has made significant progress in its program to enhance and expand the partial-wave (and multipole) analyses of fundamental two- and three-body reactions (such as pion-nucleon, photon-nucleon, and nucleon-nucleon scattering) by maintaining and augmenting the analysis codes and databases associated with these reactions. These efforts provide guidance to experimental groups at the international level, forming an important link between theory and experiment. A renaissance in light hadron spectroscopy is underway as a continuous stream of polarization data issues from existing precision electromagnetic facilities and the coming Jefferson Lab 12 GeV Upgrade. Our principal goals have been focused on supporting the national N* resonance physics program. We have also continued to study topics more generally related to the problems associated with partial-wave analysis. On the Experimental side of the CNS DAC. Its primary goal is the enhancement of the body of data necessary for our analyses of fundamental γ - N reactions. We perform experiments that study the dynamics responsible for the internal structure of the nucleon and its excitations. Our principal focus is on the N* programs at JLab and MAMI. At JLab we study spin-polarization observables using polarized photons, protons and neutrons and yielding charged final states. Similarly at MAMI we study neutral meson photoproduction off polarized protons and neutrons. We use the Crystal Ball and TAPS spectrometers (CBT) to detect photons and neutrons to measure the photoproduction of π0, η, 2π0, π0η, and K0 off the neutron. The CBT program complements our program at JLab, which studies reactions resulting in charged final states. We are also involved in a renewed effort to make neutral pion photoproduction measurements close to threshold at Mainz. In addition to the programs underway, we are contributing to the future by participation in preparations for the coming JLab 12 GeV Upgrade. GW students are involved in tests of the detectors proposed to be used with CLAS12, i.e., for the CentralTime-of-Flight Barrel (CTOF). WJB is heavily involved in the MUSE quest at PSI to solve the Proton Radius Puzzle.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Collated experience of the use of multi-group analysis (MGA) in non-destructive waste assay systems - 15257

The quantification of plutonium is an important aspect of non-destructive assay of Special Nuclear Material (SNM). The calculation of {sup 240}Pu{sub eff} can be done through prior knowledge of the composition of the waste but more often is done through the analysis of a measured plutonium gamma-ray spectrum performed using an isotopic analysis code. In this work, the {sup 240}Pu{sub eff} results from a variety of non-destructive waste assay systems in the UK are presented and compared using MGA v9.63C and v10. The body of the results include real measurement results of a variety of waste materials containing plutonium, the results of reanalysis of these initial waste measurements with changes to certain setup parameters and the results of reanalysis with a later version of MGA (v10). Additionally, measurements of test matrices with known source isotopics are measured and the results analysed and presented here. The use of the 'waste' flag parameter introduced to allow measurement of unusual spectra in v10 is also evaluated. When using MGA as part of the waste assay characterisation process in automatic or batch mode, results have proven to be very dependent on the quality of the spectra measured and on the setup parameters that the code is reliant upon. With a modern version of the code, the likelihood of under-reporting the total plutonium content of a 200 litre/55 gallon drum using PNCC and MGA isotopics in combination is low. Additionally, the latest 'hardened' code will also reduce the number of cases where MGA fails to run and the use of default isotopics results in an over-pessimistic evaluation of the total plutonium mass. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Emerging Jets Search, Triton Server Deployment, and Track Quality Development: Machine Learning Applications in High Energy Physics

Machine learning is becoming prevalent in high energy physics, with numerous applications in physics analyses and event reconstruction showing great improvements compared to traditional computing methods. This thesis studies three projects which each propose new avenues for machine learning applications within the high energy physics CMS experiment located at CERN. In the first project, a search for a dark matter signal called “emerging jets” is performed, using graph neural networks to greatly increase sensitivity to the signal’s signature within the data. The result of this dark matter search sets the most stringent exclusion limits to date on theoretical emerging jet models. Motivated by inefficiencies encountered when processing the emerging jet graph neural network at Fermi National Accelerator Laboratory’s computing centers, the second project re-optimizes the computing centers for machine learning inference. This re-optimization uses NVIDIA Triton Inference Servers to process users’ analysis code heterogeneously, therefore achieving high processing throughput and decreasing user time-to-insight. The last project focuses on an upgrade to the CMS experiment’s real-time event selection system which improves physics object reconstruction under harsh processing conditions. A boosted decision tree is used to quickly and efficiently quantify a reconstructed particle’s “track quality” in order to remove particle tracks reconstructed erroneously. In summary, this thesis will not only present examples of how high energy physics can greatly benefit by leveraging machine learning techniques for physics analysis and reconstruction, but will also provide guidance on how the field can prepare for the inevitable increase in machine learning applications.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Measuring Degradation of Fielded Systems on an Ongoing Basis with RdTools

The solar investment community needs up-to-date and data-derived system degradation rates to use in financial models that calculate the risk and expectation of energy revenue. To date, the majority of financial models choose degradation rates from past studies that are limited and may not be representative of future solar assets (systems) they are financing. As the solar industry expands in all sectors (residential, commercial, and utility), there is a growing amount of energy production available to analyze. With the help of the open source degradation analysis code (RdTools) it is now possible to derive degradation on an ongoing basis and continuously provide up-to-date durability statistics to the investment community. Ongoing and accurate statistics of fielded photovoltaic systems allows financial stakeholders to constrain energy revenue projections, lower financial risk, and increase the bankability of solar. This poster will present considerations and software engineering best practices to encourage data owners to measure degradation of their own photovoltaic assets on an ongoing basis. This poster will also preview an upcoming DOE-funded kWh Analytics degradation project which involves applying RdTools analysis on kWh Analytics’ industry database, aggregating results across variables such as weather, hardware, system configurations etc., and sharing the aggregated results with the broader community.

14 SOLAR ENERGY↗

DIC Challenge Sample Images Update

After the meetings at SEM and ICEM this year, which were both well attended by participants, it was decided by the participants that a first round of scoring the codes would be done using the Sample 14 and Sample 15 images. There was plenty of discussion on how we (the DIC Challenge Board) were going to score the results. What is going to be the balance between noise and filtering? And so forth. So it was decided to use a sub-group of the participants to help figure out if the submission guidelines were working, and how we would score the results. An additional benefit of this is that we can fix any submission guideline issues before getting results more broadly, and begin writing automated analysis codes. I expect that there will be a discussion on both subjects after I create a draft document of the scoring. This document is a draft of that report.

97 MATHEMATICS AND COMPUTING↗

Code Validation of SAM Using Forced and Natural Circulation Data from NACIE-UP Benchmark

Heavy liquid metals (HLMs) are promising candidates as coolants of Generation IV fast reactors due to their thermophysical properties. In the last decade, experimental work has been proposed as part of research and development efforts to develop such systems. In this context, researchers from the Brasimone Research Center have conducted many experiments using the Natural Circulation Experiment Upgrade (NACIE-UP) facility to study the thermofluid dynamic behavior of HLMs in rod bundle configurations with or without wire wrappers. This facility consists of a rectangular loop operated with lead-bismuth eutectic. Sensors across the loop monitor relevant parameters, i.e. temperatures, heat transfer, and flow conditions. Here, in the present work, we carefully select published data from NACIE-UP to validate the System Analysis Module (SAM), a modern system analysis code developed at Argonne National Laboratory. We developed one SAM model using specifications of the facility geometry and materials existing in relevant papers and reports. On top of that, these references provided the boundary conditions for simulating natural circulation and forced convection experiments in either steady or transient conditions. The SAM model simulates five test cases with diverse operating conditions. Ultimately, the code is proven to predict temperatures and mass flow rates that closely match the experiments. The discrepancies between numerical predictions and diverse transients are limited to a few degrees Celsius, showcasing that SAM is well suited for analyzing nuclear systems relying on HLM coolants.

advanced reactors↗

A logical progression of steps for implementation and testing of the 7-equation, two-phase model into a computational framework

A seven-step progression is given for the implementation of the seven-equation model, employed in the new RELAP-7 thermal hydraulics system analysis code, within the MOOSE-compatible framework. The model is well-posed and treats both phases as compressible flow in a one-dimensional, variable cross-sectional area duct/pipe/component. Some appropriate tests are also outlined to benchmark the model/implementation for each step of the progression. (authors)

97 MATHEMATICS AND COMPUTING↗

Thermal-hydraulic validation tests using the SMART-ITL facility with a single train passive safety injection system for the smart design

A test program was launched using a large-scale thermal-hydraulic integral effect test facility of SMART-ITL (or FESTA) to validate the performance of the SMART passive safety system. In this paper thermal-hydraulic validation tests using the SMART-ITL facility equipped with a single train of a passive safety injection system were performed and their major results were summarized. The effect of the flow distributor, break size and different types of Core Makeup Tank (CMT) and Safety Injection Tank (SIT) on the thermal-hydraulic behavior during a Small-Break Loss of Coolant Accident (SBLOCA) scenario was investigated for the SMART passive safety system development. During the SBLOCA scenario, the water layer was well stratified thermally in the core makeup tank, and the safety injection water was injected efficiently into the reactor pressure vessel from the initial period. Eventually, the reactor coolant system was cooled down properly. From the test results, it was estimated that the SMART passive safety system has sufficient cooling capability to deal with the SBLOCA scenario of SMART. The acquired data will be used to validate the safety analysis code and its related models, to evaluate the performance of SMART PSS, and to provide the base data during the application phase of the revision of standard design approval and construction licensing. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Diagnostics: Chapter 8 of the special issue: on the path to tokamak burning plasma operation

This chapter presents the activity conducted by the ITPA topical group (TG) on Diagnostics over about the last 15 years. Following a general introduction of the ITER Diagnostics led by their measurement roles, the document is organized in several subchapters detailing the design support, research and development activity conducted by each of the specialist working groups (WGs) of the TG. Please note that the magnetic diagnostics were supported at the TG without a specific WG. Their status is included in the general introduction. In the following some highlights of the subchapter’s contents are provided. Recent advances in ITER first wall (FW) diagnostics for the measurements of plasma-metallic wall interaction in support of the ITER research plan are reported. An InfraRed imaging Video Bolometer for ITER has been developed and tested on several tokamaks to measure the radiated power loss. A laser-induced breakdown spectroscopy (LIBS) technique which utilizes a pulsed laser beam to ablate locally by forming a crater, will measure local tritium inventory in the FW material. Real-time Residual Gas Analyzers will measure the neutral gas composition in a divertor port and an equatorial port during plasma operation. Due to the full metallic FW environment, the plasma-wall interaction in ITER will face several challenges such as the compromised radiated power and divertor heat flux measurements by reflection. Ray tracing and analysis codes have been developed to eliminate and correct the effects of reflection in the measurements. The characteristics of the reflecting surfaces depending on the roughness and angle of the incidence have been measured by dedicated experiments, and the results were applied to the reflection elimination. For the measurement of the metallic impurity radiation induced by eroded metallic atoms, a vacuum ultraviolet spectrometer has been developed and tested. An extensive thermonuclear diagnostic suite will be required to support the operation of ITER and the planned experimental program for future burning plasma experiments. Due to the harsh environmental conditions, the implementation of diagnostic systems in ITER is a major challenge. These conditions include high levels of neutron and gamma fluxes, neutron heating, particle bombardment. Therefore, the selection and design of diagnostic systems must take into account a number of phenomena previously unseen in diagnostic design. For this reason, the measurement of neutrons and confined or lost fast ions, with particular emphasis on alpha particles, is critical to ITER. The diagnostics associated with these measurements will be important for future plasma-burning experiments at ITER. The high neutron emission and very large plasma size in ITER make neutron diagnostics the main diagnostic method used to measure plasma parameters such as fusion power, fusion power density, ion temperature, energy of fast ions and their spatial distributions in the plasma core. Active spectroscopy techniques are methods where a neutral particle beam is injected into the plasma and information on plasma parameters is extracted from the measurement of line emission resulting from the beam-plasma interaction, either by plasma ions or by beam atoms. Spatial localization is achieved by crossing the beamline and multiple observation lines. The ITER plasma will be a high temperature, moderately dense, fully ionized collisional plasma. The plasma facing surfaces are principally metallic being fashioned from beryllium or tungsten but many other elements, arising from either structural or from operational needs, may enter this plasma. The energy range of the emitted photons range from meV (infra-red) to multi keV (x-rays) and originate from all areas of the plasma volume. The primary role of passive emission diagnostics is to identify what is in the plasma from spectral signatures. Extracting quantitative information from these measurements such as impurity content, ion temperature, rotation, degree of detachment and radiated power depends on calibrated instruments, a physics model of the atomic and molecular processes and plasma transport and an analysis workflow that takes into account environmental effects such as reflections. The particular needs for ITER have prompted a multi-machine, many-year effort to address all these aspects and this chapter reviews the work on diagnostic design, experiments and new analysis techniques. An overview of the laser diagnostics to be implemented on ITER is also provided in this paper. This includes descriptions of the Thomson scattering in the core, edge and divertor regions, polarimetry and interferometry diagnostics used for measuring plasma density and also measurements of helium density in the divertor using Laser Induced Flourescence. Techniques which can allow improvements on current measurements are also addressed in particular expanding poloidal polarimetry measurements to measure field fluctuations and proposed use of dispersion interferometery which has a number of advantages over existing methods. This paper identifies particular areas where further research and testing on existing tokamaks is useful even at this advanced stage to inform the design of diagnostics for ITER. Outstanding areas of concern for the implementation of laser diagnostics, in particular with a view to reliable operation are identified. An overview of the latest developments of microwave diagnostic systems and techniques is given. The primary focus is the contributions for ITER—the next step burning plasma experiment—which is supplemented by describing recent progress of techniques applicable for fusion experiments beyond ITER. The contributions are intentionally kept concise, and are being supplemented by a rich list of references for further studies. Radiation induced effects are receiving continuous and well-deserved attention of the ITER diagnostic community and they are in many cases one of the primary design drivers of the ITER diagnostic systems. The paper summarizes recent progress in this area focusing primarily on the ITER diagnostics but in some cases provides also outlook for the possible solutions for even more demanding radiation environment of fusion reactors beyond ITER. Despite advancements in the area of modeling and simulation of various radiation induced effects, experimental testing in a nuclear environment as close as possible to the target one is still seen as unavoidable for proper qualification of particular diagnostic functional elements. Recent advancement within three diagnostic areas: optical diagnostics, magnetics and bolometers is covered. Encouraging results on qualification of silica glass vacuum window assemblies are presented. In the area of magnetic sensors, progress of irradiation tests performed on ITER in-vessel LTCC inductive sensors is presented with outlook for novel technological approaches to inductive sensors utilizing thick printing and photolithography technologies being highlighted. Summary of advancements in the area of steady state magnetic field sensors based on Hall effect is given. New results of neutron irradiation test of the ITER borosilicate glass inserts for vacuum electrical feedthroughs are summarized finding negligible swelling at target level of neutron fluence. Off-line irradiation tests of fiber optic current sensors for plasma current measurement demonstrated that both for gamma doses up to 5 MGy and a total neutron fluence up to 10 15 cm −2 , radiation induced changes are still compatible with required measurement accuracy on ITER. The ITER bolometers are given as an example how considering radiation effects may influence the diagnostic design. Finally, outlook for future main R&D directions is outlined. All optical and laser-based diagnostics in ITER will be using mirrors to guide plasma radiation toward detectors, cameras and sensors. In the hostile plasma, radiation and particle environment the optical characteristics of diagnostic mirrors will degrade directly affecting the entire performance of involved diagnostic systems. An assessment of factors affecting mirror performance is provided. Among the prime adverse factors are deposition of plasma impurities, sputtering of mirror surface and steam ingress in the vicinity of mirrors. Within the International Tokamak Physics Activity with active support by ITER central team and domestic agencies, the structured research and development (R&D) program on mitigation of risks for diagnostic mirrors is underway. Within this program the mirror material development, the passive mitigation of mirror degradation by using diagnostic ducts and shutters along with an active mirror recovery program comprising the in-situ mirror cleaning and calibration is underway. Recent developments in diagnostic mirror R&D are described in this Chapter along with an example of their implementation of R&D solutions in ITER Infrared Thermography diagnostic. An assessment of still open engineering and physics questions, considerations on mirror risks during an early phase of ITER operation are given along with an overview of diagnostic mirror evolution in the late ITER operation stage toward the demonstration fusion power plant. Several crucial areas of diagnostic R&D outlined in ITER Research Plan are addressed. The basic control groups in a fusion reactor can be broken-down in five categories: (1) plasma position, magnetic configuration, and plasma current control, (2) profile control and confinement optimization, (3) MHD control and suppression, (4) edge dissipation control, radiation and plasma exhaust control and (5) break-down optimization. These categories are coupled via the physics (a control action in one domain will affect the other domains) and via shared actuators (e.g. ECRH for impurity accumulation avoidance, current density distribution control and MHD suppression). Consequently, a supervisory control system should determine the priority of the various control tasks, their couplings, and the interfaces with the safety and interlock system. For the systematic development of the various controllers taking the complexity of the plasma and the control system into account, a model-based approach is required. A short historical overview is given of the developments in systems and control theory and control engineering with special emphasis on those developments that are most relevant for Nuclear Fusion research and operation. An overview is given of the state of the field of fusion plasma control for the control categories. It will be shown how synthetic diagnostics are being developed in ITER and how they are used in diagnostic design and design validation and how they can be in model-based controller synthesis using relatively simple models. In modern control methods, multiple diagnostics are used to constrain relatively simple models. The constrained models provide an estimate for the state. This opens the route to state controllers, such as model predictive control. A major challenge in nuclear fusion research is the coherent combination of data from heterogeneous diagnostics and modeling codes for machine control and safety as well as physics studies. Measured data from different diagnostics often provide information about the same subset of physical parameters. Additionally, information provided by some diagnostics might be needed for the analysis of other diagnostics. A joint analysis of complementary and redundant data allows, e.g. to improve the reliability of parameter estimation, to increase the spatial and temporal resolution of profiles, to obtain synergistic effects, to consider diagnostics interdependencies and to find and resolve data inconsistencies. Physics-based modeling and parameter relationships provide additional information improving the treatment of ill-posed inversion problems. A coherent combination of all kind of available information within a probabilistic framework allows for improved data analysis results. The concept of integrated data analysis (IDA) in the framework of Bayesian probability theory is outlined and contrasted with conventional data analysis. Components of the probabilistic approach are summarized and specific ingredients beneficial for data analysis at fusion devices are discussed.

ITER↗

EVALUATION OF PRESSURE DROP CORRELATIONS FOR THE WIRE-WRAPPED ROD BUNDLES

The accurate prediction of pressure drop in the fuel rod bundle is of importance for the design and thermal-hydraulic analysis of nuclear reactor system. Due to very complex geometry of wire-wrapped fuel rod bundle geometry for the Sodium Fast Reactor (SFR), a simple method using hydraulic (equivalent) diameter is not sufficient to predict the pressure drop of rod bundle accurately. For several decades, many experimental works have been performed to develop an accurate correlation for wire-wrapped fuel bundle friction factor. In these works, the friction factor correlation of wire-wrapped fuel bundle has been developed as a function of Reynolds number of the flow and the design parameters such as fuel pin diameter, wire diameter, rod pitch, and wire lead length. The Cheng and Todreas correlations are well-known and used in many subchannel analysis code, such as CADET, TRIO-U, MATRA-LMW, SE2-ANL, etc. Novendstern model, Rehme model, Baxi and Dalle Donne model are also well-known and widely referred. Computational Fluid Dynamics is another typical approach to predict the pressure drop of wire-wrapped fuel rod bundle, but it requires the model validation because the accuracy of pressure drop prediction by the CFD simulation depends on what turbulence model is employed. The present work aims at evaluating the existing pressure drop correlations for the wire-wrapped fuel rod bundle and evaluating the pressure drop of 217-pin rod bundle by the CFD analysis.

99 GENERAL AND MISCELLANEOUS↗

Stabilizing effect of resistivity towards ELM-free H-mode discharge in lithium-conditioned NSTX

Linear stability analysis of the national spherical torus experiment (NSTX) Li-conditioned ELM-free H-mode equilibria is carried out in the context of the extended magneto-hydrodynamic (MHD) model in NIMROD. The purpose is to investigate the physical cause behind edge localized mode (ELM) suppression in experiment after the Li-coating of the divertor and the first wall of the NSTX tokamak. Besides ideal MHD modeling, including finite-Larmor radius effect and two-fluid Hall and electron diamagnetic drift contributions, a non-ideal resistivity model is employed, taking into account the increase of Z eff after Li-conditioning in ELM-free H-mode. Unlike an earlier conclusion from an eigenvalue code analysis of these equilibria, NIMROD results find that after reduced recycling from divertor plates, profile modification is necessary but insufficient to explain the mechanism behind complete ELMs suppression in ideal two-fluid MHD. After considering the higher plasma resistivity due to higher Z eff , the complete stabilization could be explained. A thorough analysis of both pre-lithium ELMy and with-lithium ELM-free cases using ideal and non-ideal MHD models is presented, after accurately including a vacuum-like cold halo region in NIMROD to investigate ELMs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Drift-flux Analysis in Horizontal Two-phase Flow

While the most detailed approach to model two-phase flows is the two-fluid model, the drift-flux model may also provide effective methods for estimating one-dimensional two-phase flow parameters. For example, in the current nuclear system analysis codes, the interfacial drag force is a function of the void fraction, drift velocity, and area-averaged relative velocity, which are estimated based on drift-flux approach. While correlations for these parameters exist for vertical two-phase flows, correlations for horizontal two-phase flows have not been developed. (authors)

42 ENGINEERING↗

Stabilizing effect of resistivity towards ELM-free H-mode discharge in lithium-conditioned NSTX

Linear stability analysis of the national spherical torus experiment (NSTX) Li-conditioned ELM-free H-mode equilibria is carried out in the context of the extended magneto-hydrodynamic (MHD) model in NIMROD. Our purpose is to investigate the physical cause behind edge localized mode (ELM) suppression in experiment after the Li-coating of the divertor and the first wall of the NSTX tokamak. Besides ideal MHD modeling, including finite-Larmor radius effect and two-fluid Hall and electron diamagnetic drift contributions, a non-ideal resistivity model is employed, taking into account the increase of Z eff after Li-conditioning in ELM-free H-mode. And unlike an earlier conclusion from an eigenvalue code analysis of these equilibria, NIMROD results find that after reduced recycling from divertor plates, profile modification is necessary but insufficient to explain the mechanism behind complete ELMs suppression in ideal two-fluid MHD. After considering the higher plasma resistivity due to higher Z eff, the complete stabilization could be explained. Furthermore, a thorough analysis of both pre-lithium ELMy and with-lithium ELM-free cases using ideal and non-ideal MHD models is presented, after accurately including a vacuum-like cold halo region in NIMROD to investigate ELMs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DYNAMIC PRA OF A MULTI-UNIT PLANT

Dynamic Probabilistic Risk Analysis (PRA) methods couple stochastic methods (e.g., RAVEN, ADAPT, ADS, MCDET) with safety analysis codes (RELAP5-3D, MELCOR, MAAP) to determine risk associated to complex systems such as nuclear plants. Compared to classical PRA methods, which are based on static logic structures (e.g., Event-Trees, Fault-Trees), they can evaluate with higher resolution the safety impact of timing and sequencing of events on the accident progression. Recently, special attention has been given to nuclear plants which consist of multiple units and, in particular, on the safety impact of system dependencies, shared systems and common resources on core damage frequencies. In the literature, while classical PRA methods have been employed to model multi-unit plants, Dynamic PRA methods have never been applied to analyze a full multi-unit model. This paper presents a PRA analysis of a multi-unit plant using Dynamic PRA methods. We employ RAVEN as stochastic method coupled with RELAP5-3D. The plant under consideration consists of the three units and their associated spent fuel pools (SFPs). The studied initiating event is a seismic induced station blackout event. We will describe in detail how the multi-unit plant has been constructed and, in particular, how unit dependencies and shared resources are modeled.

97 MATHEMATICS AND COMPUTING↗

Investigation of Abnormal Heat Transfer and Flow in a VHTR Reactor Core

The main objective of this project was to identify and characterize the conditions under which abnormal heat transfer phenomena would occur in a Very High Temperature Reactor (VHTR) with a prismatic core. High pressure/high temperature experiments have been conducted to obtain data that could be used for validation of VHTR design and safety analysis codes. The focus of these experiments was on the generation of benchmark data for design and off-design heat transfer for forced, mixed and natural circulation in a VHTR core. In particular, a flow laminarization phenomenon was intensely investigated since it could give rise to hot spots in the VHTR core.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Measurement Error Variance Estimation in Gamma-Spectroscopy Data using FRAM

The FRAM gamma-ray isotopic analysis code has recently been applied to four Pu data sets and four U data sets, each with relatively well-characterized standard items, in order to estimate the variances of random errors and of item-specific systematic errors. Assay results are available for Pu238, Pu239, Pu240, Pu241, Pu242 (operator entered, not independently measured), Am241, Specific Power, and Pu240eff. For the four U data sets, results are available for the measurands U234, U235, U236 (operator entered, not independently measured), and U238. For example, in Pu data set 1, there are 33 standard reference Pu items, which have well-known nominal/true values (except for a few items). FRAM was used to make approximately 15 repeat measurements for each of the 33 items, so empirical standard deviations can be calculated. In addition, FRAM internally estimates its random error variance using propagation of variance. This paper applies analysis of variance (ANOVA) to separately estimate the repeatability standard deviation and the standard deviation of the item-specific biases δ s in each of the 8 data sets.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Light Water Reactor Sustainability Program: Grizzly Year-End Progress Report

The Grizzly software application is being developed under the Light Water Reactor Sustainability (LWRS) program to address aging and material degradation issues that could potentially become an obstacle to life extension of nuclear power plants beyond 60 years of operation. Grizzly is based on INL’s MOOSE multiphysics simulation environment, and can simultaneously solve a variety of tightly coupled physics equations, and is thus a very powerful and flexible tool with a wide range of potential applications. Grizzly, the development of which was begun during fiscal year (FY) 2012, is intended to address degradation in a variety of critical structures. The reactor pressure vessel (RPV) was chosen for an initial application of this software. Because it fulfills the critical roles of housing the reactor core and providing a barrier to the release of coolant, the RPV is clearly one of the most safety-critical components of a nuclear power plant. In addition, because of its cost, size and location in the plant, replacement of this component would be prohibitively expensive, so failure of the RPV to meet acceptance criteria would likely result in the shutting down of a nuclear power plant. The current practice used to perform engineering evaluations of the susceptibility of RPVs to fracture is to use the ASME Master Fracture Toughness Curve (ASME Code Case N-631 Section III). This is used in conjunction with empirically based models that describe the evolution of this curve due to embrittlement in terms of a transition temperature shift. These models are based on an extensive database of surveillance coupons that have been irradiated in operating nuclear power plants, but this data is limited to the lifetime of the current reactor fleet. This is an important limitation when considering life extension beyond 60 years. The currently available data cannot be extrapolated with confidence further out in time because there is a potential for additional damage mechanisms (i.e. late blooming phases) to become active later in life beyond the current operational experience. To develop a tool that can eventually serve a role in decision-making, it is clear that research and development must be perfomed at multiple scales. At the engineering scale, a multiphysics analysis code that can capture the thermomechanical response of the RPV under accident conditions, including detailed fracture mechanics evaluations of flaws with arbitrary geometry and orientation, is needed to assess whether the fracture toughness, as defined by the master curve, including the effects of embrittlement, is exceeded. At the atomistic scale, the fundamental mechanisms of degradation need to be understood, including the effects of that degradation on the relevant material properties. In addition, there is a need to better understand the mechanisms leading to the transition from ductile to brittle fracture through improved continuum mechanics modeling at the fracture coupon scale. Work is currently being conducted at all of these levels with the goal of creating a usable engineering tool informed by lower length-scale modeling. This report summarizes progress made in these efforts during FY 2013.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

HYSPLIT/MACCS Atmospheric Dispersion Model Technical Documentation and Benchmark Analysis

The nuclear accident consequence analysis code MACCS has traditionally modeled dispersion during downwind transport using a Gaussian plume segment model. MACCS is designed to estimate consequence measures such as air concentrations and ground depositions, radiological doses, and health and economic impacts on a statistical basis over the course of a year to produce annualaveraged output measures. The objective of this work is to supplement the Gaussian atmospheric transport and diffusion (ATD) model currently in MACCS with a new option using the HYSPLIT model. HYSPLIT/MACCS coupling has been implemented, with HYSPLIT as an alternative ATD option. The subsequent calculations in MACCS use the HYSPLIT-generated air concentration, and ground deposition values to calculate the same range of output quantities (dose, health effects, risks, etc.) that can be generated when using the MACCS Gaussian ATD model. Based on the results from the verification test cases, the implementation of the HYSPLIT/MACCS coupling is confirmed. This report contains technical details of the HYSPLIT/MACCS coupling and presents a benchmark analysis using the HYSPLIT/MACCS coupling system. The benchmark analysis, which involves running specific scenarios and sensitivity studies designed to examine how the results generated by the traditional MACCS Gaussian plume segment model compare to the new, higher fidelity HYSPLIT/MACCS modeling option, demonstrates the modeling results that can be obtained by using this new option. The comparisons provided herein can also help decision-makers evaluate the potential benefit of using results based on higher fidelity modeling with the additional computational burden needed to perform the calculations. Three sensitivity studies to investigate the potential impact of alternative modeling options, regarding 1) input meteorological data set, 2) method to estimate stability class, and 3) plume dispersion model for larger distances, on consequence results were also performed. The results of these analyses are provided and discussed in this report.

97 MATHEMATICS AND COMPUTING↗