Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “source code”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 937 records · Page 52

Pressure distribution based optimization of phase-coded acoustical vortices

Based on the acoustic radiation of point source, the physical mechanism of phase-coded acoustical vortices is investigated with formulae derivations of acoustic pressure and vibration velocity. Various factors that affect the optimization of acoustical vortices are analyzed. Numerical simulations of the axial, radial, and circular pressure distributions are performed with different source numbers, frequencies, and axial distances. The results prove that the acoustic pressure of acoustical vortices is linearly proportional to the source number, and lower fluctuations of circular pressure distributions can be produced for more sources. With the increase of source frequency, the acoustic pressure of acoustical vortices increases accordingly with decreased vortex radius. Meanwhile, increased vortex radius with reduced acoustic pressure is also achieved for longer axial distance. With the 6-source experimental system, circular and radial pressure distributions at various frequencies and axial distances have been measured, which have good agreements with the results of numerical simulations. The favorable results of acoustic pressure distributions provide theoretical basis for further studies of acoustical vortices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Education and training in the field of nuclear instrumentation and measurement: CEA/INSTN (National Institute for Nuclear Sciences and Technologies) strategy to improve and develop new pedagogical tools and methods

Part of the French Alternative Energies and Atomic Energy Commission (CEA), the National Institute for Nuclear Science and Technology (INSTN) is a higher education institution whose mission is to provide students and professionals a high level of scientific and technological qualification in all disciplines related to nuclear energy applications. In this frame, INSTN carries out education and training (E and T) programs in nuclear instrumentation and radioprotection. Its strategy has always been to complete theoretical courses by training courses and laboratory works carried out on an extensive range of training tools that includes a large panel of nuclear instrumentation as well as software applications. Since its creation in 1956, the INSTN has conducted both education and vocational programs on ionizing radiation detection. An extensive range of techniques have commonly been used during practical works with students and employees of companies who need to get the knowledge and specialization in this field. Today, the INSTN is mainly equipped with usual detectors and electronics in large numbers in order to be able to accommodate up to 48 trainees at the same time in two classrooms, with only two trainees for one workstation in order to optimize their learning. In the field of the neutron detection systems, the INSTN has strongly developed its offer taking advantage of the use of research reactors, such as ISIS reactor (700 kW) at Saclay. The implementation of neutron detection systems specific to the courses offers a unique way of observing and analysing the signal coming from neutron detectors, as well as learning how to set the parameters of the detection system in real conditions. Providing the trainees with an extensive overview of each part of the neutron monitoring instrumentation apply to a nuclear reactor, hands-on measurements on the ISIS reactor play a major role in ensuring a practical and comprehensive understanding of the neutron detection system and their integration in the safety system of nuclear reactors. It also gives a solid background for the follow up and the development of the neutron detection systems. Another field of activity of the INSTN is the development of new teaching tools using software applications. We describe hereafter two applications that have recently been developed by the institute, O.S.I.R.I.S and DOSIMEX. The INSTN and the company OREKA have developed an innovative teaching tool named O.S.I.R.I.S. The tool is built on a virtual 3D environment in which users operate in a totally free way. The action is located in a pressurised water reactor in the environment of a steam generator building. The users are students or professionals who want to learn and train on concrete situations on how to protect the workers against radiation. Through the use of this serious game, users must: establish a predictive dose evaluation, implement the principles of radiation protection, perform the dose result of the operation, analyse the gap between the predictive collective dose and the collective dose achieved, and decide about the lessons to provide feedback. Taking into account the increasing constraints related to radiation protection, it is essential to use computer codes to evaluate the dose rates generated by sources of ionizing radiation. Many powerful codes exist but are often relatively complex to implement and are dedicated to experts. In addition, these codes are often 'black boxes' that does not allow the user to understand the underlying physics. The objective of DOSIMEX codes developed by the INSTN is to give, by the means of a VBA code with a simple interface, the ability to easily conduct a wide range of radiological exposure situations: calculation of gamma, beta, neutron equivalent dose rate, calculation of internal contamination and atmospheric transfer.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Detailed Characterization of CZT Detector Response for Improved Coded-Aperture Imaging Performance

Gamma-ray imaging is a powerful method for locating and quantifying sources of radiation. The coded-aperture technique demonstrates superior angular resolution in comparison to other methods (e.g., Compton reconstruction). In this method, a mask constructed of highly attenuating material encodes the scene as a shadow pattern on a position-sensitive detector; this pattern can then be used to recreate the origin(s) of incident radiation. This is typically done through convolution of the mask and shadow patterns. Iterative methods which attempt to reconstruct the observed shadow pattern using a weighted combination of simulated patterns may also be employed. In either case, errors in event position reconstruction due to detector imperfections alter the shadow pattern and will therefore degrade system performance and may introduce imaging artifacts. These effects can be mitigated with a detailed understanding of such errors – allowing for the generation of representative simulations that include the errors and/or correction of raw imager data to remove the errors. We present a calibration process for a commercially available cadmium zinc telluride (CZT) gamma imager which provides a comprehensive characterization of the spatial and energy dependence of event reconstruction. By illuminating a mask featuring a regular grid of pinholes with a calibration source, the localized response of the detector can be measured with fine granularity. These local responses are combined to generate a full detector response map which can be used to distort simulations in a manner that is representative of the observed detector data. Details of the calibration procedure and an assessment of the impact of its end products on the performance of iterative imaging methods will be presented.

Ziock, Klaus-Peter↗

Verification of Upcoming MCNP Features For Estimating Nuclear Data Sensitivities in Fixed Source Simulations [Abstract]

Predictive simulation codes, like the Monte Carlo N-Particle (MCNP) transport code, are used throughout the nuclear community. These simulations are based on nuclear data. Maximizing the accuracy and precision of nuclear data maximizes the accuracy and precision of the overall simulation. This is imperative to applications that rely on simulations. For example, improving nuclear data for special nuclear material improves simulation accuracy in stockpile stewardship applications, which results in larger safety margins and decreased operational costs. The improvement and validation of nuclear data is completed through integral benchmark experiments. Past benchmarks have primarily been limited to focus on the effective multiplication factor ($\kappa$ eff ); broadening the purview of benchmarks beyond $\kappa$ eff -dependent nuclear data addresses nuclear data deficiencies. Different response types depend on different areas of nuclear data. This dependence is quantified as nuclear data sensitivity: the change in response due to perturbation of a contributing parameter. The larger the nuclear data sensitivity of a response, the more the experiment is influenced by the uncertainties of the nuclear data. The optimization of nuclear data sensitivities in future benchmarks would result in more detailed validation of lesser studied areas of nuclear data. Currently, direct sensitivity capabilities are not easily found for all experiment types and parameters. An MCNP tool to directly estimate the cross section sensitivities of tallied values is under development. Additionally, updates have been made to the perturbation feature of MCNP, which can be used in a less direct approach to estimating sensitivities. This work verifies these features to estimate nuclear data sensitivities in fixed source simulations of a 4.5-kg sphere of alpha- phase weapons-grade plutonium surrounded by differing amounts of copper and polyethylene. Integrated estimates made using MCNP’s tools were found to statistically agree with integrated estimates made from manual perturbation of nuclear data proving the validity of the MCNP tools.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Coupled Neutron and Gamma Heating Calculation Based on VARIANT Nodal Transport Method

Accurate determination of temperature distributions in fuel pins and duct walls is essential for fuel-pin performance and assembly bowing analyses. Currently, the multi-assembly, sub-channel analysis code SE2-ANL is used in thermalhydraulic analysis of sodium cooled fast reactors (SFR). SE2- ANL is interfaced with the DIF3D/GAMSOR/DIF3D procedure for coupled neutron and gamma heating calculations, developed at Argonne National Laboratory (ANL). In this procedure, the neutron flux distribution is first calculated using the finite-difference diffusion theory option of the DIF3D code. Using the neutron flux solution, the GAMSOR code computes the gamma source distribution for the subsequent gamma flux calculation with DIF3D. However, in the current heating calculation scheme, the neutron and gamma flux calculations are limited to six triangular meshes per hexagonal assembly. In addition, the pin power distributions are estimated by assuming that the power distribution within an assembly is separable in the radial and axial directions. In each assembly, the axial power profile is obtained by averaging the power densities of six triangular meshes at each axial mesh interval. The radial power profile is evaluated at a single axial mesh (e.g., at the core mid-plane) by approximating the neutron and gamma flux distributions by a linear or quadratic polynomial of two variables x and y, depending on the polynomial order specified by the user. The coefficients of the linear and quadratic polynomials are determined by the least squares fit of six and 24 triangular mesh fluxes, respectively. Furthermore, the total heat generation rate in an assembly duct wall is simply determined by a user-specified fraction of the total gamma heating within that assembly. With the advancement of computer technology, the nodal transport theory code VARIANT is routinely used in core design analyses of SFR instead of DIF3D. It is therefore necessary to develop a new heating calculation procedure based on VARIANT neutron and gamma transport calculations. In this study, a new, coupled neutron and gamma heating calculation scheme has been developed to determine the pin power distributions more accurately by eliminating the aforementioned limitations of the current procedure. For this, the GAMSOR code has been extended to generate the intra-nodal gamma source distribution to be consistent with the trial functions of VARIANT. In the new heating calculation procedure, the neutron flux distribution is determined first by solving an eigenvalue transport problem using VARIANT. The obtained neutron flux is then used in GAMSOR to calculate the gamma source distribution. Using this gamma source distribution, the gamma flux distribution is determined by solving a fixed-source transport problem using VARIANT again. Using the computed neutron and gamma flux distributions, the power distribution in each fuel pin and duct wall is determined using a new utility code CURVE, which has been developed to evaluate the VARIANT flux solution at each fuel pin position and along the duct mid-wall. For each pin segment contained in a VARIANT node, the axial power shape is represented as a quadratic profile. The performance of the new heating calculation scheme based on VARIANT transport solutions has been tested against the existing method of SE2-ANL, using the ABR-1000 metal core design. The gamma source distribution and the assembly, fuel pin and duct wall power distributions were calculated and compared. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Validation of Light Water Reactor Ex-Core Calculations with VERA

The Consortium for Advanced Simulation of Light Water Reactors (CASL) Virtual Environment for Reactor Applications (VERA) is a reactor simulation software. It offers unique capabilities by combining high-fidelity in-core radiation transport with temperature feedback by using MPACT (a deterministic neutron transport code) and COBRA-TF (a thermal-hydraulic code) with follow-on, fixed-source transport calculations using the Shift Monte Carlo code to calculate ex-core quantities of interest. In these coupled calculations, MPACT provides Shift with the fission source for follow-on ex-core calculations. These ex-core simulations can be set up to calculate detector responses, as well as the flux and fluence in ex-core regions of interest, such as the reactor pressure vessel, nozzle, and irradiated capsules. A Watts Bar Nuclear Plant Unit 1 (WBN1) ex-core model was developed, as described in this paper, and this model was used to perform coupon calculations. The results for the coupon flux calculations show close agreement with the reference values for cycle 1 produced by the two-dimensional Discrete Ordinates Transport (DORT) code and presented in a BWXT Services Inc. report. However, differences in the results (10%) seen in cycles 2 and 3 and the reasons for these differences are discussed in this paper. The VERA WBN1 model was also used to perform a vessel fluence calculation for cycle 1. Additionally, a collaboration between CASL and Duke Energy led to the first code-to-code validation of VERA for reactor ex-core applications that used a model for the Shearon Harris reactor. Results from this collaboration show excellent agreement between VERA and the Monte Carlo N-Particle Transport Code for the detector response calculations. The work performed under this collaboration is also detailed in this paper.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Molecular dynamics simulation of thermal transport in UO 2 containing uranium, oxygen, and fission-product defects

Uranium dioxide (UO 2 ) is the most commonly used fuel in light-water nuclear reactors and thermal conductivity controls the removal of heat produced by fission, thereby governing fuel temperature during normal and accident conditions. The use of fuel performance codes by the industry to predict operational behavior is widespread. A primary source of uncertainty in these codes is thermal conductivity, and optimized fuel utilization may be possible if existing empirical models are replaced with models that incorporate explicit thermal-conductivity-degradation mechanisms during fuel burn up. This approach is able to represent the degradation of thermal conductivity due to each individual defect type, rather than the overall burn-up measure typically used, which is not an accurate representation of the chemical or microstructure state of the fuel that actually governs thermal conductivity and other properties. To generate a mechanistic thermal conductivity model, molecular dynamics (MD) simulations of UO 2 thermal conductivity including representative uranium and oxygen defects and fission products are carried out. These calculations employ a standard Buckingham-type interatomic potential and a potential that combines the many-body embedded-atom-method potential with Morse-Buckingham pair potentials. Potential parameters for UO 2+x and ZrO 2 are developed for the latter potential. Physical insights from the resonant phonon-spin-scattering mechanism due to spins on the magnetic uranium ions are introduced into the treatment of the MD results, with the corresponding relaxation time derived from existing experimental data. High defect scattering is predicted for Xe atoms compared to that of La and Zr ions. Uranium defects reduce the thermal conductivity more than oxygen defects. For each defect and fission product, scattering parameters are derived for application in both a Callaway model and the corresponding high-temperature model typically used in fuel-performance codes. The model is validated by comparison to low-temperature experimental measurements on single-crystal hyperstoichiometric UO 2+x samples and high-temperature literature data. Furthermore, this work will enable more accurate fuel-performance simulations and will extend to new fuel types and operating conditions, all of which improve the fuel economics of nuclear energy and maintain high fuel reliability and safety.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Summary report on UO 2 thermal conductivity model refinement and assessment studies

Uranium dioxide (UO 2 ) is the most commonly used fuel in light water nuclear reactors and thermal conductivity controls the removal of heat produced by fission, therefore, governing fuel temperature during normal and accident conditions. The use of fuel performance codes by the industry to predict operational behavior is widespread. A primary source of uncertainty in these codes is thermal conductivity, and optimized fuel utilization may be possible if existing empirical models were replaced with models that incorporate explicit thermal conductivity degradation mechanisms during fuel burn-up. This approach is able to represent the degradation of thermal conductivity due to each individual defect type, rather than the overall burn-up measure typically used which is not an accurate representation of the chemical or microstructure state of the fuel that actually governs thermal conductivity and other properties. To generate a mechanistic thermal conductivity model, molecular dynamics (MD) simulations of UO 2 thermal conductivity including representative uranium and oxygen defects and fission products are carried out. These calculations employ a standard Buckingham type interatomic potential and a potential that combines the many-body embedded atom method potential with Morse-Buckingham pair potentials. Potential parameters for UO 2+x and ZrO 2 are developed for the latter potential. Physical insights from the resonant phonon-spin scattering mechanism due to spins on the magnetic uranium ions have been introduced into the treatment of the MD results, with the corresponding relaxation time derived from existing experimental data. High defect scattering is predicted for Xe atoms compared to that of La and Zr ions. Uranium defects reduce the thermal conductivity more than oxygen defects. For each defect and fission product, scattering parameters are derived for application in both a Callaway model and the corresponding high-temperature model typically used in fuel performance codes. The model is validated by comparison to low-temperature experimental measurements on single crystal hyper-stoichiometric UO 2+x samples and high-temperature literature data. Ongoing works include investigation of the effect of phase separation to UO 2 +U 4 O 9 on the low temperature thermal conductivity of UO 2+x , and modeling of thermal conductivity using the Green-Kubo method. Ultimately, this work will enable more accurate fuel performance simulations as well as extension to new fuel types and operating conditions, all of which improve the fuel economics of nuclear energy and maintain high fuel reliability and safety.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Umpire: Integration with Other Software Technologies

Umpire is an open-source C++ library that provides memory and resource management functions for use with heterogeneous memory spaces. The Umpire library is being developed to simplify coordination and sharing of limited memory resources for complex integrated applications developed in the ECP/ATDM and ASC programs, and other multi-physics codes that combine many physics packages and libraries. The overarching goal of Umpire is to provide a single application-facing API for memory management on all platforms by decoupling memory operations from platform-specific memory spaces. It provides operations for: memory space discovery, memory provisioning, allocation and transfers, and introspection. Umpire development began in the Center for Applied Scientific Computing at Lawrence Livermore National Laboratory (LLNL). Umpire is funded primarily by LLNL ASC and ATMD portfolios and receives a small amount of support as part of ECP ST Project 2.3.1.06 STMP08-ISC. The source code is hosted and developed on GitHub. In this report, we explore possible integrations with other memory hierarchy management software. We examine both projects supported by ECP, as well as third-party and vendor provided tools.

97 MATHEMATICS AND COMPUTING↗

The HANDE-QMC Project: Open-Source Stochastic Quantum Chemistry from the Ground State Up

Building on the success of Quantum Monte Carlo techniques such as diffusion Monte Carlo, alternative stochastic approaches to solve electronic structure problems have emerged over the past decade. The full configuration interaction quantum Monte Carlo (FCIQMC) method allows one to systematically approach the exact solution of such problems, for cases where very high accuracy is desired. The introduction of FCIQMC has subsequently led to the development of coupled cluster Monte Carlo (CCMC) and density matrix quantum Monte Carlo (DMQMC), allowing stochastic sampling of the coupled cluster wave function and the exact thermal density matrix, respectively. In this Article, we describe the HANDE-QMC code, an open-source implementation of FCIQMC, CCMC and DMQMC, including initiator and semistochastic adaptations. We describe our code and demonstrate its use on three example systems; a molecule (nitric oxide), a model solid (the uniform electron gas), and a real solid (diamond). An illustrative tutorial is also included.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

SU-F-T-54: Determination of the AAPM TG-43 Brachytherapy Dosimetry Parameters for A New Titanium-Encapsulated Yb-169 Source by Monte Carlo Calculations

Purpose: To determine the AAPM TG-43 brachytherapy dosimetry parameters of a new titanium-encapsulated Yb-169 source designed to maximize the dose enhancement during gold nanoparticle-aided radiation therapy (GNRT). Methods: An existing Monte Carlo (MC) model of the titanium-encapsulated Yb-169 source, which was described in the current investigators’ published MC optimization study, was modified based on the source manufacturer’s detailed specifications, resulting in an accurate model of the titanium-encapsulated Yb-169 source that was actually manufactured. MC calculations were then performed using the MCNP5 code system and the modified source model, in order to obtain a complete set of the AAPM TG-43 parameters for the new Yb-169 source. Results: The MC-calculated dose rate constant for the new titanium-encapsulated Yb-169 source was 1.05 ± 0.03 cGy per hr U, indicating about 10% decrease from the values reported for the conventional stainless steel-encapsulated Yb-169 sources. The source anisotropy and radial dose function for the new source were found similar to those reported for the conventional Yb-169 sources. Conclusion: In this study, the AAPM TG-43 brachytherapy dosimetry parameters of a new titanium-encapsulated Yb-169 source were determined by MC calculations. The current results suggested that the use of titanium, instead of stainless steel, to encapsulate the Yb-169 core would not lead to any major change in the dosimetric characteristics of the Yb-169 source, while it would allow more low energy photons being transmitted through the source filter thereby leading to an increased dose enhancement during GNRT. Supported by DOD/PCRP grant W81XWH-12-1-0198 This investigation was supported by DOD/PCRP grant W81XWH-12-1- 0198.

60 APPLIED LIFE SCIENCES↗

MCNP6 Elevation Scaling of Cosmic Ray Backgrounds - Paper 29

The galactic cosmic-ray (GCR) source option was implemented in the all-energy, all-particle transport code MCNP62 in 2010. This source option along with other newly added features to MCNP6 have been used to produce increasingly accurate global cosmic and terrestrial background flux data files (background.dat).3 This data file is produced on a grid and is read and sampled by MCNP6 whenever a user invokes the background source option. The most recent adaptation of this data file, Version 4, will be completed in 2014. This paper reports on current progress toward enhancements that have been made to the modeling and simulation of these background spectra, in particular the elevation scaling used to produce more accurate spectra at off-grid locations. Cosmic radiation continuously bombards Earth with various particles, such as protons, alpha particles, and heavier nuclei, some of which are deflected by the Earth's shielding magnetic field. Particles that carry sufficient momentum can overcome the deflection and penetrate into the atmosphere. The sufficient momentum is dependent on the terrestrial coordinates due to the shape of the Earth's magnetic field and the Lorentz force's proportionality to the sine of the angle between the velocity vector of the incoming particle and the magnetic field direction. As the particles propagate through the atmosphere, collisions with atmospheric molecules generate secondary particles such as neutrons, protons, photons, muons, pions, and other exotic particles. These secondary particles often have sufficient energy to undergo additional nuclear interactions, and so on, forming what is known as a cascade shower. The particles in the cascade shower are constantly undergoing energy losses due to interactions, especially charged particles whose paths are bent into helical trajectories around magnetic field lines. These energy losses make the integral flux of the cosmic background at the ground level highly dependent on the atmospheric depth traveled by the incident particle and its progeny, and therefore dependent on the elevation of the ground where the cascade deposits its energy. For relatively small depth changes of approximately 100 g/cm{sup 2} or less, the ground-level cosmic flux follows an exponential model. The tabulation of background particle fluxes on the surface of the earth is important for a variety of reasons, one of which is the design of nuclear material detection systems. This method is also useful in producing user-specified effective attenuation lengths for particles. Attenuation lengths are an essential parameter used to compare particle spectra from different locations, elevations, and energy ranges. Traditionally these attenuation lengths are inferred from measured data. Now with the ability to simulate an accurate galactic cosmic ray source and transport it and its shower through a modeled atmosphere, these attenuation lengths can be calculated for any set of conditions. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗

Development of a BlueCRAB/MELCOR Framework for Supporting Realistic Mechanistic Source Term Calculations in Microreactors

Efforts are currently underway to deploy microreactor modeling and simulation tools to better support vendors and regulatory authorities in submitting and reviewing licensing applications. In particular, the Nuclear Regulatory Commission is expected to rely on the Comprehensive Reactor Analysis Bundle (BlueCRAB) multiphysics toolset in performing design- and beyond-design-basis accident analyses. In addition, the Nuclear Regulatory Commission has been using the MELCOR code to estimate mechanistic source terms during accidents. As MELCOR relies on isotopic inventory and reactor temperature/power evolution profiles during accident conditions—all of which can theoretically be obtained from BlueCRAB—the ultimate goal of this activity is to establish a common BlueCRAB-MELCOR framework. However, prior to the present research, BlueCRAB had never been used to calculate such quantities of interest at the full-core level. While there are many Monte Carlo (MC) codes capable of computing such quantities of interest, they are unable to readily account for multiphysics feedback. BlueCRAB allows for the coupling of different physics codes together to perform multiphysics-informed calculations. Therefore, the purpose of this fiscal year 2023 work is to investigate the feasibility and challenges of performing such calculations within BlueCRAB so as to generate the data that MELCOR relies on. To demonstrate the methodology, the proposed workflow was applied to a prototypical heat pipe-cooled microreactor model. To predict isotopic concentrations (taking into account the ac- cumulation of fission products during operation), the necessary microscopic cross sections were generated via OpenMC and tabulated with respect to temperature and burnup. Next, a recently developed capability in Griffin (the reactor physics application in BlueCRAB) was used to convert the OpenMC output format into the ISOXML format used by Griffin. A multiphysics microscopic depletion calculation that involved performing a coupled full-core, heterogeneous neutron trans- port and thermal calculation at each depletion step was conducted to deplete the core to end of life (EOL) conditions so as to provide both isotopics and the initial condition for the transient calculation. Following a brief null-transient to verify that the initial condition had been properly restarted and was indeed in thermal equilibrium, a heat pipe failure transient was simulated. Thus, the entire workflow of using BlueCRAB to generate MELCOR inputs, from cross-section generation to producing isotopic inventory and power/temperature evolution profiles during transients, is demonstrated. This report also details the identified gaps in the workflow and how they were (for the most part) addressed. Future work should focus on directly including MEL- COR into the workflow by performing a MELCOR calculation using the BlueCRAB-generated input data. In addition, the heat pipe reactor design should be improved so as to reflect more prototypical burnup characteristics at EOL.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

SU-E-T-380: Evaluation of BEBIG HDR 60Co System for AccuBoost Therapy

Purpose: In this project, the possibility of utilizing the BEBIG 60Co HDR system for AccuBoostTM treatment has been evaluated. Methods: Dose distributions in various breast sizes have been calculated for both Co-60 and Ir-192 sources using the MCNP5 code. These calculations were performed in breast tissues with thicknesses of 4cm, 6cm, and 8cm. The initial calculations were performed with the same applicator dimensions as the existing applicators used with the HDR Ir-192 system. The activity of the Co-60 source was selected such that the dose at the breast center was the same as the values from 192Ir. Then, the applicator thicknesses were increased to twice of those used with HDR Ir-192 system, for reducing skin and chest doses by Co-60 system. Dose to breast skin and chest wall were compared for both applicators types, with and without inclusion of a focusing cone at the applicator center. Results: The results showed that loading HDR Co-60 source inside the thin applicators impose higher doses to breast skin and chest wall compared to the 192Ir source. The area of the chest wall covered by 10Gy when treated by Co-60 with the thin and thick applicators, or treated by Ir-192 with thin applicator are 79cm2, 39cm2, and 3.8cm2, respectively. These values are reduced to 34cm2, 0cm2, and 0cm2 by using the focusing cone. It is worth noting that the breast skin areas covered by the 60Gy isodose line are 9.9cm2 and 7.8cm2 for Co-60 with the thin and thick applicators, respectively, while it is 20cm2 for Ir-192 when no focusing cone is present. These values are 0cm2, 0cm2, and 11cm2 in the presence of the focusing cone. Conclusion: The results indicate that using Co-60 with the thicker applicators is beneficial because of the higher half-life of Co-60, and the reduced maximum skin dose when compared with Ir-192.

60 APPLIED LIFE SCIENCES↗

Current State of Neutron Source Term Calculations

SOURCES-4C has some significant deficiencies that are mostly due to problems with the default data libraries that it uses. Some research groups have modified these libraries to fix most of the problems, but not all of them. Adding modern data libraries to SOURCES-4C is one of the most promising approaches to resolving these issues. Other promising options include using either NEDIS or a GEANT4 tool. Both of these codes approximately match the accuracy of the modified SOURCES-4C code. The NEDIS code is Russian. The GEANT4 tool might be computationally slow and it might also no longer be actively maintained. There are a handful of other (α,n) codes, but they are not accurate.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

MetaboDirect: an analytical pipeline for the processing of FT-ICR MS-based metabolomic data

Background: Microbiomes are now recognized as the main drivers of ecosystem function ranging from the oceans and soils to humans and bioreactors. However, a grand challenge in microbiome science is to characterize and quantify the chemical currencies of organic matter (i.e., metabolites) that microbes respond to and alter. Critical to this has been the development of Fourier transform ion cyclotron resonance mass spectrometry (FT-ICR MS), which has drastically increased molecular characterization of complex organic matter samples, but challenges users with hundreds of millions of data points where readily available, user-friendly, and customizable software tools are lacking. Results: Here, we build on years of analytical experience with diverse sample types to develop MetaboDirect, an open-source, command-line-based pipeline for the analysis (e.g., chemodiversity analysis, multivariate statistics), visualization (e.g., Van Krevelen diagrams, elemental and molecular class composition plots), and presentation of direct injection high-resolution FT-ICR MS data sets after molecular formula assignment has been performed. When compared to other available FT-ICR MS software, MetaboDirect is superior in that it requires a single line of code to launch a fully automated framework for the generation and visualization of a wide range of plots, with minimal coding experience required. Among the tools evaluated, MetaboDirect is also uniquely able to automatically generate biochemical transformation networks (ab initio) based on mass differences (mass difference network-based approach) that provide an experimental assessment of metabolite connections within a given sample or a complex metabolic system, thereby providing important information about the nature of the samples and the set of microbial reactions or pathways that gave rise to them. Finally, for more experienced users, MetaboDirect allows users to customize plots, outputs, and analyses. Conclusion: Application of MetaboDirect to FT-ICR MS-based metabolomic data sets from a marine phage-bacterial infection experiment and a Sphagnum leachate microbiome incubation experiment showcase the exploration capabilities of the pipeline that will enable the research community to evaluate and interpret their data in greater depth and in less time. It will further advance our knowledge of how microbial communities influence and are influenced by the chemical makeup of the surrounding system. The source code and User’s guide of MetaboDirect are freely available through (https://github.com/Coayala/MetaboDirect) and (https://metabodirect.readthedocs.io/en/latest/), respectively.

54 ENVIRONMENTAL SCIENCES↗