Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reactor Physics Method”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Practical Insights on Applying Simulation-Based Control Methods in Experimental Studies

Advanced nuclear reactors are crucial to the future of energy both in the United States and around the globe. In contrast to the current operating fleet, they are characterized as being deployable in remote locations and able to operate in semi-autonomous or autonomous fashion. This leap forward necessitates a new reactor control paradigm. Because advanced nuclear reactors are still under development in the United States, the creation of new control methods to achieve autonomous operations has been based on systems modeling and simulation. However, an important factor in successfully deploying these new control methods is the ability to seamlessly transition from simulation environments to real-world settings. Control methods tested in both simulation and experimental settings need to be investigated in the context of advanced reactor applications. This work developed a series of simple controllers for Idaho National Laboratory (INL)’s Microreactor Applications Research Validation and Evaluation (MARVEL) microreactor operating in load-following scenarios. These controllers were tested in both simulation and experimental settings, and a comparative performance analysis was performed. The simulation tests leveraged the Control and Optimization Modular Modeling Application for Nuclear Deployment (COMMAND) software developed in a previous stage of the current effort, along with the MARVEL Reactor Excursion and Leak Analysis Program (RELAP5-3D) and Monte Carlo N-Particle (MCNP) models. The experimental tests leveraged the COMMAND software, MARVEL models, and the U.S. Department of Energy Microreactor Program’s Microreactor Automated Control System (MACS). MACS was developed to serve as a control method testbed. It was customized to mirror the MARVEL microreactor, and COMMAND enabled MACS to emulate the physics of MARVEL. The load-following controller was developed using the simulation platform, with efforts to emulate real systems by introducing actuator saturation and noise. These factors were incrementally accounted for in the controller design. After finalizing the controller design, it was implemented with the experimental setup. The experimental conditions tested included an initial test under conditions similar to the final simulation test, and two additional scenarios. The first scenario introduced additional actuator saturation to account for equipment aging over time, which was unknown to the controller. The second scenario introduced sensor delay, a phenomenon anticipated with the use of remote operations or wireless communication in advanced reactors. These tests revealed several notable differences. While the controller performed well in simulation, it exhibited several limitations when transitioning to hardware. The main challenges involved maintaining the steady-state target power, as evidenced by larger error values between the true reactor power and setpoint power, as well as persistent oscillations in controlled reactor power. These issues could lead to unacceptable transient conditions in real reactor testing. Introducing actuator aging and stochastic delays in the experimental setup significantly impacted controller performance, resulting in increased overshoot and undershoot, and exacerbated error and oscillations previously mentioned. These findings underscore the importance of experimental testbeds for testing and validating control methods, as controllers developed using only theory and/or simulation may perform unexpectedly when applied to actual hardware. This research emphasizes the need for an experimental testbed for achieving such validation.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗
DOE OSTI2016

Initial Comparison of Direct and Legacy Modeling Approaches for Radial Core Expansion Analysis

Radial core expansion in sodium-cooled fast reactors provides an important reactivity feedback effect. As the reactor power increases due to normal start up conditions or accident scenarios, the core and surrounding materials heat up, causing both grid plate expansion and bowing of the assembly ducts. When the core restraint system is designed correctly, the resulting structural deformations introduce negative reactivity which decreases the reactor power. Historically, an indirect procedure has been used to estimate the reactivity feedback due to structural deformation which relies upon perturbation theory and coupling legacy physics codes with limited geometry capabilities. With advancements in modeling and simulation, radial core expansion phenomena can now be modeled directly, providing an assessment of the accuracy of the reactivity feedback coefficients generated by indirect legacy methods. Recently a new capability was added to the PROTEUS-SN unstructured geometry neutron transport solver to analyze deformed meshes quickly and directly. By supplying the deformed mesh in addition to the base configuration input files, PROTEUS-SN automatically processes material adjustments including calculation of region densities to conserve mass, calculation of isotopic densities according to material models (for example, sodium density as a function of temperature), and subsequent re-homogenization of materials. To verify the new capability of directly simulating deformed meshes, PROTEUS-SN was used to compute reactivity feedback for a series of contrived yet representative deformed configurations for the Advanced Burner Test Reactor design. The indirect legacy procedure was also performed to generate reactivity feedback coefficients for the same deformed configurations. Interestingly, the legacy procedure consistently overestimated reactivity feedbacks by 35% compared to direct simulations by PROTEUS-SN. This overestimation indicates that the legacy procedures are in fact not conservative and could be overestimating reactivity feedback effects that are closely tied to reactor safety. We conclude that there is indeed value in performing direct simulation of deformed meshes despite the increased computational expense. PROTEUS-SN is already part of the SHARP multi-physics toolkit where both thermal hydraulics and structural mechanical feedback modeling can be applied but this is the first comparison of direct simulation to legacy techniques for radial core expansion.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2017

Power control system modification and simulation of CPR1000 - 181

Due to development of electric power system, the proportion of nuclear power and its needs for load follow operation have become large. Then there would be higher requirements to the power control ability of the nuclear power system. The axial power distribution may change a lot in the operation process and be affected by many factors. In particular, as the power change is performed within a load follow maneuver, several modifications occur in the core from a neutronic view point: the fuel and moderator temperature change, the xenon concentration and distribution are modified, the power distribution skewed axially, etc. These changes need to be adequately counterbalanced to keep both the core critical and the power distribution acceptable. The traditional approach in pressurized water reactors (PWRs) is compensating for the reactivity change due to the power variation by adjusting the soluble boron concentration and moving a limited number of control banks. However, dilution/boronation strategy becomes expensive for short time changes and leads to large volume of liquid effluent, so it is eager about a capability to perform load follow operation without changing soluble boron concentration. CPR1000, one of the main reactor types in China national development, is taken an example in this paper. It adopts a control method called Mode G, using G banks, N banks, R bank and soluble boron to perform the load follow control. According to the analysis of the physics and thermal principle, a two-node dynamic function model was built, putting the coupling coefficient and mutual influence into consideration. By the calculation of these transient parameter values which are acquired by the steady-state calculation of single channel, the core system model of CPR1000 is established. In order to carry out the load follow operation without the boron adjustment, named as BTP mode, the control banks were regrouped to control the average temperature by M banks and axial offset by AO bank. After that, the average temperature control system (M banks control system) and the axial offset control system (AO bank control system) are built. When the optimal control system is connected to the core model simulation platform, the BTP mode load follow operation of the modified CPR1000 is simulated in the end of the paper. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

Single-Bubble Flow Boiling Phenomena - Interface Tracking Simulation

Boiling, as one of the most efficient heat transfer mechanisms, is widely used in various engineering systems. Better understanding and modeling of this process remains a major challenge in multiphase flow research. The distribution of vapor in boiling system affects the heat transfer rate and may cause unfavorable conditions, such as heater burn-out. A number of flow regimes have been observed experimentally. In general, there are three regimes under boiling conditions: nucleate boiling, transition boiling, and film boiling. The nucleate boiling regime is categorized into partial boiling regime and fully developed nucleate boiling, according to behavior of bubble dynamics and heat transfer characteristics. In the first region, as the wall temperature increases, the fraction of the wall surface subject to nucleate boiling increases until bubble formation occupies the entire heated surface. Heat transfer in this region is a complex mixture of single-phase forced convection and nucleate boiling. As bubble density increases, bubbles turn to coalesce and form insulating vapor patches, which impedes liquid flow back to the surface. The heat transfer rate under this condition reduces dramatically. This point is called Critical Heat Flux (CHF), which leads to boiling crisis. The behavior of boiling crisis is dependent on local fluid conditions due to bubble growth and departure. Experiments are the ultimate tool for boiling model development and verification, but they still lack some capabilities needed for better understanding of the local mechanism of this phenomena. It's difficult to capture many important details at the pressure/temperature conditions of interest using experimental techniques, such as full three dimensional measurements of velocity and temperature distribution and estimates of heat flux partitioning during the boiling process. The interface tracking simulation (ITS) is one of the promising approaches to describe heat transfer of boiling phenomena and their underlying mechanisms. The advances of high performance computing (HPC) in recent years made it possible to apply DNS to a wide variety of adiabatic flows. Some work has been done to implement evaporation / condensation models as well. Lee and Nydahl performed a numerical simulation of bubble growth and departure assuming that the bubble has the shape of a truncated sphere and a micro-layer exists during bubble growth. A 2-D cylindrical bubbles simulation has been done by Son with a simplified model to predict the heat flux in a thin liquid micro-layer. A direct simulation of 2D film boiling is conducted by Juric and Tryggvason. The film boiling, below critical pressure with a modified Level Set Method, is studied by Son and Dhir. Ose and Kunugi have developed a boiling and condensation model on subcooled boiling phenomena. Fuchs performed numerical simulations of bubble growth and departure considering wall conduction effects. Fully transient heat and fluid flow is modeled with a free surface of a rising bubble and a periodic calculation of the whole cycling of a growing, detaching, and rising bubble. Mukherjee and Kandikar have done a numerical simulation of an evaporating meniscus on a moving heated surface with level set method. The heat transfer associated with the advancing and receding contact angle has been studied. To develop the capability of physics-based boiling process, in the presented work the energy and momentum equations are coupled to simulate single-bubble growth scenarios. For a single bubble simulation in infinite liquid, the analytical solution exists and can be compared with the numerical results for verification purposes. The presented capabilities will be utilized in the large scale DNS of boiling flows in the future. The finite-element based code, PHASTA is used for the simulation. It is a parallel, hierarchic, higher-order accurate, adaptive, stabilized, transient analysis flow solver, which has been shown to be an effective tool for a wide range of single and two-phase flow problems. It has been shown to be highly scalable on top supercomputers ({approx}85% scaling on up to 3 million mesh partitions using a 92 billion element mesh . PHASTA has been shown to reliably predict the details of adiabatic bubbly flows as well as single-phase flows with heat transfer. To develop the capability of physics-based boiling process, in the presented work the energy and momentum equations are coupled to simulate a single-bubble growth scenario. For a single bubble growth in superheated water, the analytical solution exists and can be compared with the numerical results for verification purposes. The second-order accurate spatial and time discretization is used for the flow solver in the presented work. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2017

Industry Application ECCS / LOCA Integrated Cladding/Emergency Core Cooling System Performance: Demonstration of LOTUS-Baseline Coupled Analysis of the South Texas Plant Model

Under the auspices of the DOE LWRS Program RISMC Industry Application ECCS/LOCA, INL has engaged staff from both South Texas Project (STP) and the Texas A&M University (TAMU) to produce a generic pressurized water reactor (PWR) model including reactor core, clad/fuel design and systems thermal hydraulics based on the South Texas Project (STP) nuclear power plant, a 4-Loop Westinghouse PWR. A RISMC toolkit, named LOCA Toolkit for the U.S. (LOTUS), has been developed for use in this generic PWR plant model to assess safety margins for the proposed NRC 10 CFR 50.46c rule, Emergency Core Cooling System (ECCS) performance during LOCA. This demonstration includes coupled analysis of core design, fuel design, thermalhydraulics and systems analysis, using advanced risk analysis tools and methods to investigate a wide range of results. Within this context, a multi-physics best estimate plus uncertainty (MPBEPU) methodology framework is proposed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2015

Year End Progress Report on Rattlesnake Improvements

Rattlesnake is a MOOSE-based radiation transport application developed at INL to support modern multi-physics simulations. At the beginning of the last year, Rattlesnake was able to perform steady-state, transient and eigenvalue calculations for the multigroup radiation transport equations. Various discretization schemes, including continuous finite element method (FEM) with discrete ordinates method (SN) and spherical harmonics expansion method (PN) for the self-adjoint angular flux (SAAF) formulation, continuous FEM (CFEM) with SN for the least square (LS) formulation, diffusion approximation with CFEM and discontinuous FEM (DFEM), have been implemented. A separate toolkit, YAKXS, for multigroup cross section management was developed to support Rattlesnake calculations with feedback both from changes in the field variables, such as fuel temperature, coolant density, and etc., and in isotope inventory. The framework for doing nonlinear diffusion acceleration (NDA) within Rattlesnake has been set up, and both NDA calculations with SAAF-SN-CFEM scheme and Monte Carlo with OpenMC have been performed. It was also used for coupling BISON and RELAP-7 for the full-core multiphysics simulations. Within the last fiscal year, significant improvements have been made in Rattlesnake. Rattlesnake development was migrated into our internal GITLAB development environment at the end of year 2014. Since then total 369 merge requests has been accepted into Rattlesnake. It is noted that the MOOSE framework that Rattlesnake is based on is under continuous developments. Improvements made in MOOSE can improve the Rattlesnake. It is acknowledged that MOOSE developers spent efforts on patching Rattlesnake for the improvements made on the framework side. This report will not cover the code restructuring for better readability and modularity and documentation improvements, which we have spent tremendous effort on. It only details some of improvements in the following sections.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

Improved Few-group Homogenized Cross-section Sensitivity Analysis for Lattice Physics Calculation

In the conventional reactor physics 'two-step' calculation scheme, lattice-physics calculation is followed by core calculation using neutron diffusion theory. The few-group spatially homogenized cross-sections are the outputs of the lattice-physics calculation as well as the inputs of the core calculation. The uncertainties of fundamental nuclear data will be propagated to the core calculation results through few-group spatially homogenized cross-sections. Consequently, assessing the sensitivity and uncertainty of few-group homogenized cross-sections is necessary for the core-calculation sensitivity and uncertainty analysis. Conventional method for sensitivity analysis is usually performed with respect to effective self-shielding cross sections but not continuous-energy cross sections by neglecting resonance self-shielding effects. The purpose of this work is to improve the sensitivity analysis of few-group cross sections and assess the effects introduced by resonance calculation with the subgroup method being used for resonance calculation. This paper is organized as follows: Firstly, theoretical models of this work are described. Secondly, validation and analysis are performed based on a boiling-water reactor (BWR) fuel-pin cell under hot zero power (HZP) condition. Finally, summaries and conclusions are given. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

Effects of fission yield data in the calculation of antineutrino spectra for U 235 ( n , fission ) at thermal and fast neutron energies

Fission yields form an integral part of the prediction of antineutrino spectra generated by nuclear reactors, but little attention has been paid to the quality and reliability of the data used in current calculations. Following a critical review of the thermal and fast ENDF/B-VII.1 235 U fission yields, deficiencies are identified and improved yields are obtained, based on corrections of erroneous yields, consistency between decay and fission yield data, and updated isomeric ratios. These corrected yields are used to calculate antineutrino spectra using the summation method. An anomalous value for the thermal fission yield of 86 Ge generates an excess of antineutrinos at 5–7 MeV, a feature which is no longer present when the corrected yields are used. Thermal spectra calculated with two distinct fission yield libraries (corrected ENDF/B and JEFF) differ by up to 6% in the 0–7 MeV energy window, allowing for a basic estimate of the uncertainty involved in the fission yield component of summation calculations. Lastly, the fast neutron antineutrino spectrum is calculated, which at the moment can only be obtained with the summation method and may be relevant for short baseline reactor experiments using highly enriched uranium fuel.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Validation of k eff Calculations for Extended BWR Burnup Credit

The US Nuclear Regulatory Commission and consensus standards recommend validation of the numerical methods used in criticality safety analyses. This validation requires the comparison of computational results with measurements of physical systems which are neutronically similar to those used in the safety analysis being performed. To this end, this document examines the methods available to generate sensitivity data to help identify systems similar to spent boiling- water reactor (BWR) fuel in a flooded spent nuclear fuel transportation or storage cask. A large number of critical benchmark experiments are surveyed using sensitivity/uncertainty (S/U) techniques to assess their applicability to BWR burnup credit (BUC) beyond the burnup of peak reactivity. Multiple burnups of BWR assemblies are considered herein, as well as the actinide-only (AO) and actinide-and-major-fission-product (AFP) isotope sets. Sample validations are completed for representative application models to demonstrate that appropriate validation is possible and to indicate the bias and bias uncertainty values expected for related applications.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2016

Implementation of the LAX-Wendroff Method in Cobra-TF for Solving Two-Phase Flow Transport Equations

COBRA-TF (Coolant Boiling in Rod Arrays Two Fluid), or CTF, is a subchannel code used to conduct the reactor core thermal hydraulic (T/H) solution in both standalone and coupled multi-physics applications. CTF applies the first-order upwind spatial discretization scheme for solving two-phase flow conservation equations. In this work, the second-order Lax-Wendroff (L-W) scheme has been implemented in CTF to solve the two-phase flow transport equations to improve numerical accuracy in both temporal and spatial discretization. To avoid the oscillation issue, a non-linear flux limiter VA (Van Albada) is employed for the convective terms in the transport equations. Assessments have been carried out to evaluate the performance and stability of the implemented second-order L-W scheme. It has been found that the L-W scheme performs better than the upwind scheme for the single-phase and two-phase flow problems in terms of numerical accuracy and computational efficiency.

CTF↗
DOE OSTI2017

Release of MCNP6.2 and WHISPER-1.1 - Guidance for NCS users

In 2017, a new version of the MCNP6 Monte Carlo code package was released: MCNP6.2. This is the first version of MCNP to include the WHISPER package to support sensitivity-uncertainty based nuclear criticality safety validation. The coding changes to MCNP6.2 physics are relatively insignificant. Corrections to the S(α,β) thermal scattering numerics are generally negligible relative to problem statistics (or, in rare cases, prevent aborts). Similarly, the changes to adjoint-weighting for computing kinetics parameters may result in small differences, generally negligible compared to problem statistics. The change to the MCNP6.2 geometry treatment to correctly handle coincident surfaces in problems with universe/fill features will produce different round-off in the geometry tracking. This will produce differences in results relative to previous versions, but those differences should be small relative to problem statistics, and are not a concern. An immediate benefit of using MCNP6.2 (rather than MCNP6.1) is that the new version is typically 1.5-2 times faster.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Dynamic Monte Carlo Transient Analysis with Thermal-Hydraulic Feedback by McCARD/CUPID

The reactor core transient analysis requires neutronics calculations with the thermal-hydraulic (TH) feedback. Recently the dynamic Monte Carlo (DMC) neutron transport algorithm controlling delayed neutron precursors was developed and successfully applied for core transient problems with help of a sub-channel code. In the Seoul National University Monte Carlo (MC) code, McCARD, a DMC scheme for the core transient analysis has been implemented by coupling with a computational fluid dynamics code, CUPID using the TCP/IP socket communication. In the McCARD/CUPID transient calculations, temperature-dependent cross sections are produced by a built-in on-the-fly Doppler broadening (OTF DB) module. This paper presents McCARD/CUPID transient analysis results for a triangular unit cell problem in the IAEA coordinated research project (CRP) benchmark on the high temperature gas cooled reactor (HTGR) uncertainty analysis. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Photoelectron spectroscopy in molecular physical chemistry

Photoelectron spectroscopy has long been a powerful method in the toolbox of experimental physical chemistry and molecular physics. Recent improvements in coincidence methods, charged-particle imaging, and electron energy resolution have greatly expanded the variety of environments in which photoelectron spectroscopy can be applied, as well as the range of questions that can now be addressed. In this Perspectives Article, we focus on selected recent studies that highlight these advances and research areas. Furthermore, the topics include reactive intermediates and new thermochemical data, high-resolution comparisons of experiment and theory using methods based on pulsed-field ionisation (PFI), and the application of photoelectron spectroscopy as an analytical tool to monitor chemical reactions in complex environments, like model flames, catalytic or high-temperature reactors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Research and Development to Reduce Impurity Production and Transport of the Impurities to the Target in Linear Plasma Devices Using Helicon Plasma Sources

Linear plasma devices used to test plasma facing materials (PFMs) and components for fusion reactors are often suffering under the production of intrinsic impurities from the plasma source system. Most linear plasma devices use internal electrodes (hollow cathodes, reflex arc, or cascading arc), which typically are the source of the main impurities. The next-generation plasma generators use radio frequency (RF) plasma sources like helicons to avoid internal electrodes. However, high power operation of helicons has proven to result in impurity production due to the high rectified sheath voltages created. Depending on the plasma parameters and magnetic configuration, these impurities can be transported to the target and deposited there to unacceptable high levels. Here, in this contribution, the experimental results from Proto-MPEX are summarized, and the conclusions of the impurity source physics are given. Methods to reduce the impurity production, the impurity transport, and the net deposition on the target are presented. These methods to reduce the impurity production include Faraday screens to reduce the sheath voltage drop, high-Z refractory coatings to reduce the erosion yield, and wall conditioning methods. Methods to reduce the impurity transport include changes in the magnetic configuration as well as electron heating to change axial and radial temperature profiles. Preliminary results on the effectiveness of some of these methods are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Present status of sensitive detector of reactor’s antineutrinos using scintillating detectors

In 2011, the reanalysis of the reactor antineutrinos spectra led to the formulation of the Reactor Antineutrino Anomaly (RAA) [1], which indicates the discrepancy between measured and expected antineutrino fluxes on short baselines. This discrepancy appears to favor the existence of the fourth “sterile” neutrino with |Δm{sup 2}|>1 eV{sup 2}. To confirm or reject this hypothesis a high sensitive antineutrino detector located close to the reactor is required. In addition to that such a detector could be used to online monitor the isotopic composition of the reactor core and to prevent illegal production and removal of{sup 239}Pu, which is the essential part of nuclear weapons. Detector DANSSino [2] already proved that even a compact antineutrino detector (∼ 1 m{sup 3}) based on polystyrene is capable of antineutrino detection in the close vicinity of a reactor core (∼ 10 m) with signal to background ratio about one. As a common activity between JINR Dubna and IEAP CTU a new prototype of detector (called S{sup 3}) has been proposed and is under construction. The construction design, selected results of Monte Carlo simulations and results of benchmark tests are presented.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

Heat loads to divertor nearby components from secondary radiation evolved during plasma instabilities

A fundamental issue in tokamak operation related to power exhaust during plasma instabilities is the understanding of heat and particle transport from the core plasma into the scrape-off layer and to plasma-facing materials. During abnormal and disruptive operation in tokamaks, radiation transport processes play a critical role in divertor/edge-generated plasma dynamics and are very important in determining overall lifetimes of the divertor and nearby components. This is equivalent to or greater than the effect of the direct impact of escaped core plasma on the divertor plate. We have developed and implemented comprehensive enhanced physical and numerical models in the upgraded HEIGHTS package for simulating detailed photon and particle transport in the evolved edge plasma during various instabilities. The paper describes details of a newly developed 3D Monte Carlo radiation transport model, including optimization methods of generated plasma opacities in the full range of expected photon spectra. Response of the ITER divertor's nearby surfaces due to radiation from the divertor-developed plasma was simulated by using actual full 3D reactor design and magnetic configurations. We analyzed in detail the radiation emission spectra and compared the emission of both carbon and tungsten as divertor plate materials. The integrated 3D simulation predicted unexpectedly high damage risk to the open stainless steel legs of the dome structure in the current ITER design from the intense radiation during a disruption on the tungsten divertor plate.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Seventy-Five Years of Nuclear Criticality Safety Documents - A Bibliography

A collection of nuclear criticality documents was begun in the 1970s at the University of California’s Lawrence Livermore Laboratory (now Lawrence Livermore National Laboratory). Documents in this bibliography are dated as early as 1940. The early documents described experimental work that was of vital importance to the safe use of fissile materials at nuclear reactors, the facilities for preparing or transporting nuclear fuel, or its eventual storage. The experiments were of direct value to support safe operations with fissile material, but they were also needed for validating early physics data and calculational methods and then for validating computational methods used in criticality safety evaluations. Bibliographic information from the past seventy-five years have been entered into a computerized bibliographic database. The database continues to expand. Publications were issued several decades ago based on the information collected at that time. Our new publication is based on a current “snapshot” of the present status of the database. Contributions of information or corrections to this bibliography will be appreciated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗