Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fission Calculations”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS1978

Direct nuclear pumping by a volume source of fission fragments

A detailed kinetic model is presented for the analysis of nuclear pumped lasers when the pumping is a result of a volume source of fission fragments. The results of the model are employed to study a He-3 - Xe laser. For the range of pressures, neutron fluxes and mixtures considered, the gain and power calculations are in good agreement with experiment. Moreover, based on these calculations, it appears that the collisional recombination is the dominant pumping mechanism for 7p-7s transitions while direct excitation is the dominant pumping mechanism for the 5d-6p transitions.

Deese, J. E.↗
DOE OSTI2021

Which nuclear data can be validated with LLNL pulsed-sphere experiments?

It is shown here that 14-MeV D+T LLNL pulsed-sphere experiments bring complementary information into the process of validating nuclear data compared to experiments that are traditionally used for this purpose—such as critical assemblies. To be more specific, the 14-MeV D+T LLNL pulsed-sphere neutron-leakage spectra enable to validate scattering and fission nuclear data up to 15 MeV (compared to approximately up to 5 MeV when using criticality experiments) and employ to this end simple compound targets containing only few isotopes. In this work, sensitivity profiles of the spectra to nuclear data are calculated in order to understand in detail which isotopes, observables, and energy ranges of nuclear data contribute significantly to their simulation. These profiles are presented for a few selected spheres containing 16 O, 12 C, 56 Fe, and 239 Pu. It is shown that the neutron-leakage spectra of spheres containing light elements are mostly sensitive to elastic- and inelastic-scattering cross sections on discrete levels and corresponding angular distributions. Spheres of structural materials are sensitive to elastic- and inelastic-scattering cross sections, including scattering on discrete levels and the continuum, and double-differential cross sections. Actinide spheres are also strongly sensitive to the fission observables, in particular to the total-fission neutron spectrum. Thin spheres (in which neutrons experience on average less than one scatter) are mostly sensitive to data near the elastic peak, in the energy range from 12–15 MeV, while thicker ones can be sensitive to data at lower incident-neutron energies due to multiple-scattering effects. This information is brought together with simulations of 71 pulsed-sphere neutron-leakage spectra using the ENDF/B-VII.1 and ENDF/B-VIII.0 nuclear-data libraries. This analysis highlights ENDF/B-VIII.0 data that could be further investigated for potential shortcomings ( 6 Li, 12 C, 16 O, 24-26 Mg, 27 Al, 48 Ti, 56 Fe, and 208 Pb) or are likely reliable ( 1,2 H, 7 Li, 9 Be, 14 N, 235,238 U, and 239 Pu) as indicated by validating with LLNL pulsed-sphere experiments.

14-MeV D+T LLNL pulsed-sphere neutron-leakage spec↗
DOE OSTI2020

Performing k eff Validation of As-Loaded Criticality Safety Calculations Using UNF-ST&DARDS: Applicable Experiment Selection

The general method for performing validation of as loaded criticality safety calculations using UNF ST&DARDS is presented in a paper by Clarity, which includes a description of the UNF-ST&DARDS system. Proof-of-principle analyses were performed in the summer of 2019 for MPC-32 dual purpose canisters (DPCs) containing pressurized water reactor (PWR) fuel assemblies. Summaries of these results are presented in this and a companion paper for this conference. The current paper describes the TSUNAMI-IP calculations performed to select applicable experiments for validation of 11 MPC-32 DPCs. The companion paper discusses the TSUNAMI-3D calculations used to generate sensitivity data to support the experiment selections discussed here. Experiment selection is based on the sensitivity/uncertainty (S/U) methods used to validate criticality safety calculations of as-loaded DPCs containing pressurized water reactor (PWR) spent nuclear fuel (SNF). This process has been demonstrated and is summarized in this paper. The approach is similar to that used in NUREG/CR-7109, which provides an approach for validation of PWR burnup credit (BUC), including major and minor actinides and major fission products. The premise of S/U-based validation is that applicable experiments—those having a similar bias to a given application system—will have similar sensitivities for each isotope and reaction in the two systems. It is assumed that cross sections with larger uncertainties are more likely to contain data errors which contribute to the bias. The integral index c k thus propagates the system sensitivities with the nuclear covariance data to calculate a correlation coefficient representing the similarity of the two systems. In this work, a c k value of 0.8 or higher is interpreted as identifying an experiment with sufficient similarity for use in validation. This paper presents a brief summary of the characteristics of the 11 MPC-32 DPCs used in the proof of-principle analysis for as-loaded criticality safety calculation validation and an overview of the critical experiment suite with which each of these DPC models was compared. A summary and discussion of c k results is also presented, followed by conclusions and a discussion of future work to be performed for validation of UNF ST&DARDS as-loaded criticality safety calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

TEX-HEU: Integral Experiments with Highly Enriched Uranium and Polyethylene (IER-297 CED-3b Report)

This report documents the experimental critical configurations and measurements for IER-297, Thermal/Epithermal eXperiments (TEX) using highly enriched uranium (HEU) fuel, moderated and reflected by polyethylene. These experiments establish the baseline TEX-HEU configurations spanning thermal, intermediate, and fast fission energy regimes. The experiment campaign was completed in FY2020 at the National Critical Experiments Research Center at the Nevada National Security Site. The campaign produced six documented critical configurations (0”, ¹/₈” A, ¹/₈” B, ¹/₄”, ¹/₂”, and 1¹/₂”) and one reproducibility measurement (¹/₈” R). Table 1 summarizes these seven experiments along with their period and calculated reactivity, in order of execution. The Comet General Purpose Critical Assembly Machine was used to conduct the TEX-HEU experiments. The HEU fuel consists of U(93+) plates, collectively known as the “Jemima” plates. Each experiment consisted of an alternating stack of these HEU plates with or without polyethylene moderator plates. By varying the thicknesses of the polyethylene plates, the neutron spectrum of the experiment could be fine-tuned to a specific fission energy regime. The experiments were reflected by 1” of polyethylene consisting of a bottom and top plate and outer reflector rings to surround the HEU and polyethylene plates. The fine reactivity control of the system was provided by the thickness of the top reflector plate, which could be varied in increments of ¹/₃₂”.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Extending the Nuclide Inventory Validation Basis for High-Burnup Fuel with New Radiochemical Assay Data

Efforts are underway at Oak Ridge National Laboratory to improve the nuclide inventory validation basis for spent nuclear fuel at high burnups. Recently conducted radiochemical assay experiments provided new measurement data for nine samples of fuel irradiated in a pressurized water reactor, with estimated sample burnups in the 30 to 70 GWd/t range. This type of destructive assay data is essential for validating computational methods, tools, and nuclear data applied in nuclear safety analyses and for improving our understanding of the bias and uncertainty in code predictions. The measurement data include key actinides and fission products that span a gamut of needs and interests for nuclear science and engineering applications in criticality safety, reactor physics, nuclide inventory, decay heat, and radiation shielding. The SCALE 6.3 code system with ENDF/B-VII.1 cross-section libraries was used to simulate the irradiation histories of the measured fuel samples. The calculated nuclide concentrations are compared to corresponding measurement data. The significance of the comparisons is discussed, emphasizing how the addition of the new measurement data fills gaps in the validation basis at high burnups and contributes to the decrease in bias and uncertainty for predicted nuclide concentrations. The discussion addresses the effect of the sample burnup used in the simulation—which is based on reactor operator records or on calibration to measured data for burnup indicator fission products—on the validation results.

Nuclide inventory↗
DOE OSTI2025

First-principles investigation of cerium and neodymium diffusion in BCC chromium and vanadium via vacancy-mediated transport

Lanthanide transport plays a crucial role in the performance and longevity of metallic nuclear fuels. This study examines the diffusion behavior of Ce and Nd—two major fission products—in body-centered cubic (BCC) Cr and V, which are potential liner or coating materials for mitigating fuel-cladding chemical interactions (FCCI). Using density functional theory (DFT) calculations and self-consistent mean-field (SCMF) analysis, the vacancy-mediated diffusion coefficients are evaluated. Our findings reveal that Ce and Nd act as oversized solutes and are strongly bound to vacancies in BCC Cr and V, with diffusivities in Cr significantly lower than in V and in hexagonal closed-packed (HCP) Zr, as investigated in our previous work. The activation energies for Ce and Nd diffusion are 3.39 and 3.32 eV, respectively, in BCC Cr, and 2.56 and 2.33 eV, respectively, in BCC V. Analysis of vacancy drag and partial diffusion coefficient ratios indicates a strong tendency for lanthanide enrichment at vacancy sinks in BCC Cr, and to a lesser extent in BCC V, with this effect persisting up to the melting point in Cr and remaining substantial for Nd in V at high temperatures. Under irradiation, the increase in vacancy concentration is expected to enhance lanthanide transport, potentially accelerating interactions at liner-cladding interfaces. Although BCC Cr exhibits relatively low lanthanide diffusivities under equilibrium conditions, the expected segregation tendencies under irradiation suggest that Zr liners may be a more favorable option. Further investigations using rate theory, cluster dynamics, and phase-field modeling are required to quantitatively assess the performance of these materials in reactor environments.

36 - MATERIALS SCIENCE↗
DOE OSTI2026

First Observation of Multiphonon 𝛾-Vibrations in an Odd-Odd Nuclear System

Here, the identification of the first multiphonon 𝛾-vibrational bands in an odd-odd neutron-rich nucleus of the nuclear chart is presented. These high-spin structures of hard to access $^{104}_{41}\textrm{Nb}_{63}$, produced in fission, were studied by combining a spectrometer with isotopic resolution coupled to a 𝛾-ray tracking array and independently from high-fold 𝛾 coincidence measurements. Triaxial projected shell model calculations for the high-spin states are in good agreement with the measured observables for the yrast, one-phonon, and two-phonon 𝛾 bands. The possibility of an oblate shape of an isomeric state and coexistence of triaxial and oblate configurations is investigated from the decay of the 141 keV isomer. The present Letter illustrates the robustness of vibration excitations in the presence of an odd valence proton and neutron as well as the possibly coexisting shapes beyond the 𝑁 = 60 transitional region.

Wang, E. H. [Shandong Univ., Jinan (China); Vander↗
DOE OSTI2022

Fitting $\overline{\nu}$ for minor Pu isotopes

After successful fitting of prompt $\overline{\nu}$ for 235 U, 238 U, and 239 Pu(n,f) using CGMF, we will move on to the minor plutonium isotopes. Minor isotopes pose a greater challenge both because there is less data available and we do not, by default, have parametrizations in CGMF already. To mitigate these challenges— and provide consistency within CGMF—we will take a stepped approach to the optimization. Continuing from the 239 Pu work, we will then fit 241 Pu(n,f) $\overline{\nu}$, where there is also a number of experimental measurements, keeping consistency between the parameters that are included in the calculation for 239 Pu and 241 Pu. With these parametrizations settled, we can consistently optimize 240 Pu(n,f) $\overline{\nu}$. Following that step, we will move on to 242 Pu(n,f) (and increasing neutron number) and 238 Pu(n,f) (and decreasing neutron number). By including the fissioning systems in this manner, we should be able to minimize the unknown parameters in CGMF. We will possibly also be able to develop systematics for the CGMF input parameters along the Pu isotopic chain. This work can serve as a guide to broadening the reactions available in CGMF. In this short report, we first give an example of how we have updated CGMF to include 240 Pu(n,f) and 242 Pu(n,f), keeping consistency with the current 239 Pu and 241 Pu calculations, but without rigorous optimization (Sec. 2). Then, we will shown in Section 3 what experimental data exist for the various isotopes to provide some insight into why 241 Pu and 239 Pu are used as anchor points.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2025

Analysis of 127 Xe tracer measurements using a net counts method

A suite of measurement systems were deployed as part of the Physical Experiment 1 series of experiments, which involved detonating chemical explosives along with radionuclide tracers in an underground cavity, at the Nevada National Security Site (NNSS) in the United States. One of the radionuclide tracers, 127 Xe was released from the containment following the explosion and detected on a SAUNA Q B sampler situated approximately 3.5 km away. The system uses a beta-gamma coincidence detector system to measure fission product radioisotopes of xenon relevant to nuclear explosion monitoring. In this work we use the coincidence measurement data to analyse and interpret the results from the SAUNA Q B system, to calculate the measured 127 Xe activity concentration(s).

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2021

Low-Energy Reactions of the n+ 233 U Nuclear Compound System and Its Initial Validation

The strong negative gradient as a function of the epithermal fission fraction (FEPIT) observed in both Evaluated Nuclear Data File (ENDF)/B-VII.1 and ENDF/B-VIII.0 nuclear data libraries poses new challenges for the verification of the n+ 233 U nuclear compound system, especially in the low energy range. The aim of this paper is to summarize the steps forward to consistently update the 233 U evaluated nuclear data to attain improved performance in benchmark calculations; this work is described in detail in a forthcoming paper by Pigni et al.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Influence of point defect accumulation on in-pile thermal conductivity degradation: Fuel rod defect distribution and deviation between in-pile and post irradiation thermal conductivity

As nuclear fuel burn-up increases, its thermal conductivity degrades due to the accumulation of defects that lead to increased phonon scattering rates. This results in a rise in the centerline temperature of the fuel rod, whereby heat generation must be decreased to avoid undesired behavior such as fuel melting and extensive fission gas release. Fuel performance codes are utilized to optimize the fuel's operational conditions; while they are based on established physical principles, their empirical nature limits their predictive capabilities. Recently, an effort has been made to develop predictive fuel performance codes for commonly used nuclear fuels, as well as accelerated qualification of advanced nuclear fuels. In this report, we elaborate on the importance of careful analysis of point defects’ impact on thermal conductivity in fuel performance analysis. Further, a model is presented where point defect concentration is estimated based on Rate Theory modeling and used as input to the Klemens-Callaway model to calculate their contribution to the degradation of thermal conductivity in UO 2 under prototypical irradiation conditions. This analysis suggests that point defect concentration is significant at the rim of the fuel pellet and neglecting this leads to underestimation of the centerline temperature, which may have consequences on fission gas behavior.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Determination of the 161 Tb half-life

There is significant interest in the use of terbium radioisotopes for applications in cancer therapy and diagnosis. Of these, 161 Tb, a medium energy beta- and prolific auger emitter, is being investigated as a potential alternative to 177 Lu. The higher proportion of conversion electron and Auger electron emissions, and variety of low-energy gammas make 161 Tb an attractive targeted theranostic. As a product of nuclear fission, 161 Tb is also of importance to nuclear forensics. The current evaluated half-life of 6.89(2) d has a standard uncertainty that contributes significantly to the standard uncertainty of decay calculated activity determination. Furthermore, the accuracy of this evaluated half-life has been called into question by measurements reported recently at the Institute of Radiation Physics (IRP), Switzerland. In this work, the half-life of 161 Tb was measured at three independent laboratories located in the United Kingdom and United States of America for a total of five determinations using three independent measurement techniques. In conclusion, the half-life determined for 161 Tb of 6.9637(29) d confirms the observed 1 % relative increase observed by IRP, though the reported half-lives in this work and at IRP are significantly different (ζ-score = 3.1).

half-life↗
DOE OSTI2020

In-pile thermal conductivity of uranium dioxide at low burnup

In-pile thermal conductivity of uranium dioxide (UO 2 ) was investigated at low burnup levels (<0.2 kWd/kg-UO 2 ) to elucidate the concurrent effects of the fission damage and the thermal recovery. In-pile experiments were performed in the Halden reactor. The fuel centerline and the cladding temperatures were measured during the experiments, and the power level of specially designed test rodlets were also monitored. The uncertainty of the pellet-cladding gap's thermal resistivity was minimized with a liquid-metal bond, while comparable measurements were made in reference test rodlets without liquid metal bonding. The experimental data were analyzed using the inverse heat transfer approach. The heat conduction equation was solved by applying the measured temperature of the cladding as a boundary condition, and the calculated fuel centerline temperature was compared with the measured fuel centerline temperature to determine the relative matrix resistivity. Results showed that the in-pile matrix resistivity, the phonon-lattice interaction terms of the thermal conductivity, was ~1.5–2.5 times higher than its unirradiated value due to fission-induced damage. Finally, the in-pile annealing caused a significant recovery of the matrix resistivity, and the amount of recovery increased from ~20 to ~33%, while the annealing temperature increased from 700 to 1,000 °C.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

A refined assessment of nuclear data target accuracy requirements for ALFRED

The interest in target accuracy assessment has been recently renewed within the OECD/NEA Working Party on International Nuclear Data Evaluation Cooperation with an update exercise, thanks to the invaluable initiative of Massimo Salvatores. This work, moving from evaluations on the main core integral parameters of the Advanced Lead-cooled Fast Reactor European Demonstrator (ALFRED), and the identification of the associated target accuracies needed for improving its design, approaches the retrieval of target accuracy requirements for nuclear data, as the inverse problem of uncertainty quantification. The preliminary results obtained in a previous work are here refined by explicitly taking into account the total contribution due to correlations among variables in the optimization constraint, while its various components are evaluated individually. Moreover, three different sets of cost parameters are considered to introduce a bias in the solution in order to take in due account the relative difficulties in the execution of new differential experiments as needed for achieving the proposed target accuracies. By this, the obtained results are deemed usefully informative, and proposed for the update of a High Priority List which identifies the isotopes-reaction couples of most interest for future refinement experiments, especially when integrated with further analogous information collected on other advanced nuclear systems. Among the cross-sections whose refinement contributes the most to the reduction of the multiplication factor, the fission of Pu{sup 239} was found, which needs to be fixed below 1% in the range of interest for ALFRED (2 keV-4 MeV). Another contribution for the amelioration of the multiplication factor calculation accuracy, is due to the capture channel of Pu{sup 239}, for which the initial high uncertainty (around 15%) can be significantly reduced resorting to easily practicable experiments. The same occurs for the inelastic of Pb{sup 207}, for which the initial uncertainty value of around 50% can be drastically reduced up to target values between 4% and 7%.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS2007

Testing in Support of Fission Surface Power System Qualification

The strategy for qualifying a FSP system could have a significant programmatic impact. The US has not qualified a space fission power system since launch of the SNAP-10A in 1965. This paper explores cost-effective options for obtaining data that would be needed for flight qualification of a fission system. Qualification data could be obtained from both nuclear and non-nuclear testing. The ability to perform highly realistic nonnuclear testing has advanced significantly throughout the past four decades. Instrumented thermal simulators were developed during the 1970s and 1980s to assist in the development, operation, and assessment of terrestrial fission systems. Instrumented thermal simulators optimized for assisting in the development, operation, and assessment of modern FSP systems have been under development (and utilized) since 1998. These thermal simulators enable heat from fission to be closely mimicked (axial power profile, radial power profile, temperature, heat flux, etc.) and extensive data to be taken from the core region. For transient testing, pin power during a transient is calculated based on the reactivity feedback that would occur given measured values of test article temperature and/or dimensional changes. The reactivity feedback coefficients needed for the test are either calculated or measured using cold/warm zero-power criticals. In this way non-nuclear testing can be used to provide very realistic information related to nuclear operation. Non-nuclear testing can be used at all levels, including component, subsystem, and integrated system testing. FSP fuels and materials are typically chosen to ensure very high confidence in operation at design burnups, fluences, and temperatures. However, facilities exist (e.g. ATR, HFIR) for affordably performing in-pile fuel and materials irradiations, if such testing is desired. Ex-core materials and components (such as alternator materials, control drum drives, etc.) could be irradiated in university or DOE reactors to ensure adequate radiation resistance. Facilities also exist for performing warm and cold zero-power criticals.

Houts, Mike↗
DOE OSTI2022

Report on Data Quality Required for Modeling an MSR

Calculations were performed to assess how uncertainties in molten salt property values impact the results of behavior models of molten salt reactor systems under steady-state conditions. The thermophysical and thermochemical properties of salt mixtures being considered for use in molten salt reactors are not yet well characterized and methods used to measure property values to molten salts are not yet well established. The limited property data for salts of interest that are available in the literature are inconsistent and unreliable. Many methods used to measure property values remain developmental in the sense that not all variables affecting measured values have been identified or controlled during the measurement and not all aspects of the measurements affecting data quality are calibrated. A prior report described the repeatability of measurements made at Argonne and the reproducibility of measurements made with the same salts at several institutions is being evaluated. There is significant uncertainty in measured property values due to instrumental limitations and the effects of contaminants in the salts used to make the measurements that impacts predictions of salt behavior in a reactor. Most salts of interest are multicomponent eutectic or near-eutectic mixtures and both uncertainty in the salt composition and uncertainty in measured values. Furthermore, the salt composition evolves over time due to the uptake of impurities from the atmosphere and corrosion of containment vessels. More significantly, fission products will accumulate during operation that change the salt compositions and property values. MSR developers need to understand how this uncertainty in the basic property data of the fuel will impact modeling predictions for normal reactor operation as well was in transient scenarios that will impact the safety basis. The calculations and analyses in this report address the impact of uncertainty in measured property values under steady-state conditions through a generalized parametric study. The pooled uncertainty from contributions by a number of factors were considered and quantified, including varying composition due to impurities and fission products, imprecision in the measured values, and uncontrolled aspects of the measurement. In this way, we can begin to gain an understanding of the magnitude of the impact uncertainty will have in modeling MSR systems.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Tracking Radioactive Isotopes in HVAC and Application for Hot Cell Analyses - 20155

In recent years, the surge in the number of isotope production facilities under design has increased the need to analyze isotope migration in facility Heating, Ventilation, and Air Conditioning (HVAC) during normal operations and accident scenarios. These facilities are used to produce isotopes for medical, security, and industrial applications and as such are subject to license and regulatory requirements 10CFR20, 10CFR30, 10CFR50, and 10CFR70. GOTHIC, a general-purpose thermal-hydraulics software package, includes the ability to model isotopic tracers, radioactive decay and isotope migration as well as HEPA and charcoal filters for isotope retention. A GOTHIC model was developed by Zachry Nuclear Engineering (ZNE) to examine the effect of negative room pressure, HEPA filtration, and HVAC fluctuations on radiation areas, hot cells and gloveboxes. Radioactive tracers were used to simulate the concentration of spills within contaminated areas, track the migration of isotopes of interest, and determine the isotopic retention and buildup on facility HEPA filters. A variety of isotopes with concern to dose (e.g., Kr-85, Sr-90, I-131, etc.), including their decay and progeny, are included in the analysis. Negative pressures are maintained in the regions of interest by a representative central HVAC system equipped with a volumetric fan that exhausts to the environment after a series of isolation valves and HEPA filters. GOTHIC is an industry trusted tool for providing engineering solutions for a variety of applications, including fission product tracking, aerosol and particulate transport and ventilation assessments. The software provides an integrated analysis environment that includes a graphical user interface (GUI) for constructing analysis models, a numerical solver that includes parallel processing capabilities and a post-processor for evaluating simulation results. It solves the conservation equations for mass, momentum and energy for multicomponent, multi-phase flow in lumped parameter and multi-dimensional geometries (1, 2, or full 3D), including the effects of turbulence, diffusion and buoyancy. It has been developed and maintained under a Quality Assurance program in compliance with the requirements of 10CFR50 Appendix B and applicable portions of ASME NQA-1 since 1995. GOTHIC has been used for assessing both forced and natural convection conditions for a wide range of applications, including: - Tracking concentration of hazardous gases and chemicals for habitability and safety assessments - Determining ventilation and filtration requirements and optimizing location and arrangement of these systems - Room heat-up, including diverse and Flexible coping strategies for Extended Loss of AC Power (FLEX/ELAP) - Equipment Qualification (EQ) A distinctive feature of GOTHIC is the ability to track many different fields/substances in a simulation, including user defined tracer elements, in the liquid, vapor and droplet fields as well as surfaces and filters. This capability allows GOTHIC to model fission product transport and release or the removal of particulates or harmful toxins from exhaust gases using a spray scrubber or other types of filtration systems. GOTHIC also includes models for engineered equipment, such as fans, filters, charcoal filters, dryers/demisters, dampers, etc. The aerosols and other filtered material are removed or accumulated in these components. The range of aerosol and radiological applications that GOTHIC has been used for includes: - Source Term: Primary Coolant (Equilibrium) Activity; Non-Water Coolant Source/Leakage. - Conditions for Iodine Re-evolution: Sump/Suppression Pool Conditions and pH; - RWST Conditions and pH. - Isotope Removal Mechanisms: Containment Sprayed and Unsprayed Region Mixing; Charcoal Filter Heating due to Iodine decay. - Radionuclide Transport and Decay: Post-LOCA Release in containment; Transport between connected Compartments and vent systems; Groundwater transport of radionuclides. - Non-Newtonian Fluid modeling for sludge, waste tanks, etc. ADAMS ML071581053 (titled 'Best Practice Guidelines for the use of CFD in Nuclear Reactor Safety Applications') poses guidelines for applying single phase CFD codes in nuclear reactor safety problems and GOTHIC is listed as a 'tool for 3D flows' and 'dispersal and deposition of radionuclides.' The Nuclear Quality Assurance (NQA) pedigree of GOTHIC is an important aspect for applications in the nuclear industry. The fundamental tracer models (convective transport, molecular and turbulent diffusion, removal mechanisms, etc.) have been verified using analytical solutions and validated against applicable separate effects tests. Also, GOTHIC has been benchmarked to many integrated effects tests, including Phebus FP (Fission Product). GOTHIC gives good agreement for the buildup and decay of fission products in Phebus Test 3. The model developed by ZNE demonstrates GOTHIC's applicability and acceptability for use in analyzing the migration and retention of radioactive isotopes and their progeny in normal operation and accident scenario analyses for isotope production facilities, hot cells, and gloveboxes. The tracer activities calculated by GOTHIC can then be used in downstream radiation transport and shielding codes like RADTRAD-NAI{sup C}, MCNP{sup R}, and MicroShield{sup R} to determine on-site and of-site doses. (authors)

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

On-the-fly response function generation method for composite coarse mesh

The hybrid stochastic deterministic transport code COMET, based on the incident response expansion theory, is used to model reactor cores with high fidelity and formidable computational speed. COMET models a reactor core using a library of incident flux response expansion coefficients that are pre computed for all the unique lattice cells (e.g., fuel assemblies, reflector blocks, etc.) in the core. In order to further improve its computational efficiency in pre-calculating the response library a new response function generation method is developed to compute the response functions for the composite coarse meshes made of a smaller set of unique lattices on the fly within the COMET's deterministic transport core sweep. The efficiency is achieved by eliminating a number of unique lattices that can be made up from the reduced set of unique meshes on the fly. The numerical process consists of the following steps. First, the boundary condition on composite coarse mesh boundaries is projected onto the expansion basis to compute the incident flux moments on external surfaces of all the basic (reduced set of unique) coarse meshes. Secondly, the deterministic sweeping solver in COMET is used to converge on the outgoing/incoming flux expansion moments crossing interfaces between the basic coarse meshes. Thirdly, the response functions for the composite coarse meshes are constructed as a superposition on the fly. The new response function generation method was tested on 88 composite coarse meshes consisting of CANDU fuel bundles and moderator blocks. It was found that response functions generated by the new method agree very well with those generated by direct Monte Carlo calculations. The average and maximum relative differences in the surface-to-surface response coefficients computed by the two methods are 0.10% and 0.20%, respectively. Similarly, the average and maximum relative differences in the response fission densities are 0.13% and 0.43%, respectively. These discrepancies are within one standard deviation of the stochastic uncertainties. The new method is five times faster than the original direct Monte Carlo method. The size of the response function library for the new method is five times smaller than that for the original method, leading to significantly less requirement for the computer hard drive space and memory. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗