Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo Code”

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

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

At least 1,657 records · Page 92
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Improved Monte Carlo Glauber predictions at present and future nuclear colliders

We present the results of an improved Monte Carlo Glauber (MCG) model of relevance for collisions involving nuclei at center-of-mass energies of the BNL Relativistic Heavy Ion Collider ( s NN = 0.2 TeV), CERN Large Hadron Collider (LHC) ( s NN = 2.76 – 8.8 TeV ), and proposed future hadron colliders ( s NN ≈ 10 – 63 TeV). The inelastic p p cross sections as a function of s NN are obtained from a precise data-driven parametrization that exploits the many available measurements at LHC collision energies. We describe the nuclear density of a lead nucleus with two separated two-parameter Fermi distributions for protons and neutrons to account for their different densities close to the nuclear periphery. Furthermore, we model the nucleon degrees of freedom inside the nucleus through a lattice with a minimum nodal separation, combined with a “recentering and reweighting” procedure, that overcomes some limitations of previous MCG approaches. The nuclear overlap function, number of participant nucleons and binary nucleon-nucleon collisions, participant eccentricity and triangularity, overlap area, and average path length are presented in intervals of percentile centrality for lead-lead (PbPb) and proton-lead ( p Pb ) collisions at all collision energies. We demonstrate for collisions at s NN = 5.02 TeV that the central values of the Glauber quantities change by up to 7% in a few bins of reaction centrality, due to the improvements implemented, though typically they remain within the previously assigned systematic uncertainties, while their new associated uncertainties are generally smaller (mostly below 5%) at all centralities than for earlier calculations. Tables for all quantities versus centrality at present and foreseen collision energies involving Pb nuclei, as well as for collisions of XeXe at s NN = 5.44 TeV , and AuAu and CuCu at s NN = 0.2 TeV , are provided. The source code for the improved Monte Carlo Glauber model is made publicly available.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

SU-F-T-361: Dose Enhancement Due to Nanoparticle Addition in Skin Radiotherapy: A Monte Carlo Study Using Kilovoltage Photon Beams

Purpose: This study investigated the dose enhancement due to addition of nanoparticles with different types and concentrations in skin radiotherapy using kilovoltage photon beams. Methods: An inhomogeneous water phantom (15×15×10 cm{sup 3}) having the skin target layer (0.5–5 mm), added with different concentrations (3–40 mg/ml) of nanoparticles (Au, Pt, I, Ag and Fe{sub 2}O{sub 3}), was irradiated by the 105 and 220 kVp photon beams produced by a Gulmay D3225 Orthovoltage unit. The circular cone of 5-cm diameter and source-to-surface distance of 20 cm were used. Doses in the skin target layer with and without adding the nanoparticles were calculated using Monte Carlo simulation (the EGSnrc code) through the macroscopic approach. Dose enhancement ratio (DER), defined as the ratio of dose at the target with nanoparticle addition to the dose without addition, was calculated for each type and concentration of nanoparticle in different target thickness. Results: For Au nanoparticle, DER dependence on target thickness for the 220 kVp photon beams was not significant. However, DER for Au nanoparticle was found decreasing with an increase of target thickness when the nanoparticle concentration was increased from 18 to 40 mg/ml using the 105 kVp photon beams. For nanoparticle concentration of 40 mg/ml, DER variation with target thickness was not significant for the 220 kVp photon beams, but DEF was found decreasing with the target thickness when lower energy of photon beam (105 kVp) was used. DEF was found increasing with an increase of nanoparticle concentration. The higher the DEF increasing rate, the higher the atomic number of the nanoparticle except I and Ag for the same target thickness. Conclusion: It is concluded that nanoparticle addition can result in dose enhancement in kilovoltage skin radiotherapy. Moreover, the DER is related to the photon beam energy, target thickness, atomic number and concentration of nanoparticles.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

Comparison of MCNP6 and experimental results for neutron counts, Rossi-{alpha}, and Feynman-{alpha} distributions

MCNP6, the general-purpose Monte Carlo N-Particle code, has the capability to perform time-dependent calculations by tracking the time interval between successive events of the neutron random walk. In fixed-source calculations for a subcritical assembly, the zero time value is assigned at the moment the neutron is emitted by the external neutron source. The PTRAC and F8 cards of MCNP allow to tally the time when a neutron is captured by {sup 3}He(n, p) reactions in the neutron detector. From this information, it is possible to build three different time distributions: neutron counts, Rossi-{alpha}, and Feynman-{alpha}. The neutron counts time distribution represents the number of neutrons captured as a function of time. The Rossi-a distribution represents the number of neutron pairs captured as a function of the time interval between two capture events. The Feynman-a distribution represents the variance-to-mean ratio, minus one, of the neutron counts array as a function of a fixed time interval. The MCNP6 results for these three time distributions have been compared with the experimental data of the YALINA Thermal facility and have been found to be in quite good agreement. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Evaluation of Attila and MCNP computational methods for dose and exposure estimation

Radiation transport calculations are often used to estimate dose or exposure to components and personnel surrounding a radiation source. The sources for these calculations are decaying radionuclides within various nuclear materials. Historically, dose calculations use MCNP (Monte Carlo N-Particle) transport code as the primary particle transport tool without a secondary computational tool to validate the results from the MCNP simulations [1]. The goal of this study is to make an independent check of the Monte Carlo solution from MCNP6 Version 6.2.1 with the discrete ordinates solution from Attila 10.2.0 Beta 3. As an example problem for this study, water-filled, stainless-steel vessels, modeled with an unstructured mesh (UM) with both MCNP and Attila [2], are exposed to 252Cf and 60Co point sources. This report also includes a discussion of the limitations of unstructured mesh in a MCNP calculation.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2016

Fe kα profiles from simulations of accreting black holes

We present the first results from a new technique for the prediction of Fe Kα profiles directly from general relativistic magnetohydrodynamic (GRMHD) simulations. Data from a GRMHD simulation are processed by a Monte Carlo global radiation transport code, which determines the X-ray flux irradiating the disk surface and the coronal electron temperature self-consistently. With that irradiating flux and the disk’s density structure drawn from the simulation, we determine the reprocessed Fe Kα emission from photoionization equilibrium and solution of the radiation transfer equation. We produce maps of the surface brightness of Fe Kα emission over the disk surface, which—for our example of a 10M{sub ⊙} Schwarzschild black hole accreting at 1% the Eddington value—rises steeply one gravitational radius outside the radius of the innermost stable circular orbit and then falls ∝r {sup −2} at larger radii. We explain these features of the Fe Kα radial surface brightness profile as consequences of the disk’s ionization structure and an extended coronal geometry, respectively. We also present the corresponding Fe Kα line profiles as would be seen by distant observers at several inclinations. Both the shapes of the line profiles and the equivalent widths of our predicted Kα lines are qualitatively similar to those typically observed from accreting black holes. Most importantly, this work represents a direct link between theory and observation: in a fully self-consistent way, we produce observable results—iron fluorescence line profiles—from the theory of black hole accretion with almost no phenomenological assumptions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

SU-E-T-132: Assess the Shielding of Secondary Neutrons From Patient Collimator in Proton Therapy Considering Secondary Photons Generated in the Shielding Process with Monte Carlo Simulation

Purpose: In this study we present Monte Carlo based evaluation of the shielding effect for secondary neutrons from patient collimator, and secondary photons emitted in the process of neutron shielding by combination of moderator and boron-10 placed around patient collimator. Methods: The PHITS Monte Carlo Simulation radiation transport code was used to simulate the proton beam (Ep = 64 to 93 MeV) from a proton therapy facility. In this study, moderators (water, polyethylene and paraffin) and boron (pure {sup 10}B) were placed around patient collimator in this order. The rate of moderator and boron thicknesses was changed fixing the total thickness at 3cm. The secondary neutron and photons doses were evaluated as the ambient dose equivalent per absorbed dose [H*(10)/D]. Results: The secondary neutrons are shielded more effectively by combination moderators and boron. The most effective combination of shielding neutrons is the polyethylene of 2.4 cm thick and the boron of 0.6 cm thick and the maximum reduction rate is 47.3 %. The H*(10)/D of secondary photons in the control case is less than that of neutrons by two orders of magnitude and the maximum increase of secondary photons is 1.0 µSv/Gy with the polyethylene of 2.8 cm thick and the boron of 0.2 cm thick. Conclusion: The combination of moderators and boron is beneficial for shielding secondary neutrons. Both the secondary photons of control and those emitted in the shielding neutrons are very lower than the secondary neutrons and photon has low RBE in comparison with neutron. Therefore the secondary photons can be ignored in the shielding neutrons.This work was supported by JSPS Core-to-Core Program (No.23003). This work was supported by JSPS Core-to-Core Program (No.23003)

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

MONTE CARLO NEUTRINO TRANSPORT THROUGH REMNANT DISKS FROM NEUTRON STAR MERGERS

We present Sedonu, a new open source, steady-state, special relativistic Monte Carlo (MC) neutrino transport code, available at bitbucket.org/srichers/sedonu. The code calculates the energy- and angle-dependent neutrino distribution function on fluid backgrounds of any number of spatial dimensions, calculates the rates of change of fluid internal energy and electron fraction, and solves for the equilibrium fluid temperature and electron fraction. We apply this method to snapshots from two-dimensional simulations of accretion disks left behind by binary neutron star mergers, varying the input physics and comparing to the results obtained with a leakage scheme for the cases of a central black hole and a central hypermassive neutron star. Neutrinos are guided away from the densest regions of the disk and escape preferentially around 45° from the equatorial plane. Neutrino heating is strengthened by MC transport a few scale heights above the disk midplane near the innermost stable circular orbit, potentially leading to a stronger neutrino-driven wind. Neutrino cooling in the dense midplane of the disk is stronger when using MC transport, leading to a globally higher cooling rate by a factor of a few and a larger leptonization rate by an order of magnitude. We calculate neutrino pair annihilation rates and estimate that an energy of 2.8 × 10{sup 46} erg is deposited within 45° of the symmetry axis over 300 ms when a central BH is present. Similarly, 1.9 × 10{sup 48} erg is deposited over 3 s when an HMNS sits at the center, but neither estimate is likely to be sufficient to drive a gamma-ray burst jet.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

Dosimetric Verification Using Monte Carlo Calculations for Tissue Heterogeneity-Corrected Conformal Treatment Plans Following RTOG 0813 Dosimetric Criteria for Lung Cancer Stereotactic Body Radiotherapy

Purpose: The recently activated Radiation Therapy Oncology Group (RTOG) studies of stereotactic body radiation therapy (SBRT) for non-small-cell lung cancer (NSCLC) require tissue density heterogeneity correction, where the high and intermediate dose compliance criteria were established based on superposition algorithm dose calculations. The study was aimed at comparing superposition algorithm dose calculations with Monte Carlo (MC) dose calculations for SBRT for NSCLC and to evaluate whether compliance criteria need to be adjusted for MC dose calculations. Methods and Materials: Fifteen RTOG 0236 study sets were used. The planning tumor volumes (PTV) ranged from 10.7 to 117.1 cm{sup 3}. SBRT conformal treatment plans were generated using XiO (CMS Inc.) treatment planning software with superposition algorithm to meet the dosimetric high and intermediate compliance criteria recommended by the RTOG 0813 protocol. Plans were recalculated using the MC algorithm of a Monaco (CMS, Inc.) treatment planning system. Tissue density heterogeneity correction was applied in both calculations. Results: Overall, the dosimetric quantities of the MC calculations have larger magnitudes than those of the superposition calculations. On average, R{sub 100%} (ratio of prescription isodose volume to PTV), R{sub 50%} (ratio of 50% prescription isodose volume to PTV), D{sub 2cm} (maximal dose 2 cm from PTV in any direction as a percentage of prescription dose), and V{sub 20} (percentage of lung receiving dose equal to or larger than 20 Gy) increased by 9%, 12%, 7%, and 18%, respectively. In the superposition plans, 3 cases did not meet criteria for R{sub 50%} or D{sub 2cm}. In the MC-recalculated plans, 8 cases did not meet criteria for R{sub 100%}, R{sub 50%}, or D{sub 2cm}. After reoptimization with MC calculations, 5 cases did not meet the criteria for R{sub 50%} or D{sub 2cm}. Conclusions: Results indicate that the dosimetric criteria, e.g., the criteria for R{sub 50%} recommended by RTOG 0813 protocol, may need to be adjusted when the MC dose calculation algorithm is used.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

SU-E-T-15: A Comparison of COMS and EP917 Eye Plaque Applicators Using Different Radionuclides

Purpose: To investigate the effect of plaque design and radionuclides on eye plaque dosimetry. Methods: The Monte Carlo N-particle Code version 6 (MCNP6) was used for radiation transport simulations. The 14 mm and 16 mm diameter COMS plaques and the model EP917 plaque were simulated using brachytherapy seeds containing I-125, Pd-103, and Cs-131 radionuclides. The origin was placed at the scleral inner surface. The central axis (CAX) doses of both COMS plaques at −1 mm, 0 mm, 1 mm, 2 mm, 5 mm, 10 mm, 15 mm, 20 mm, and 22.6 mm were compared to the model EP917 plaque. Dose volume histograms (DVHs) were also created for both COMS plaques for the tumor and outer sclera then compared to results for the model EP917 plaque. Results: For all radionuclides, the EP917 plaque delivered higher dose (max 343%) compared to the COMS plaques, except for the 14 mm COMS plaque with Cs-131 at 1 mm and 2 mm depths from outer sclera surface. This could be due to source design. For all radionuclides, the 14 mm COMS plaque delivered higher doses compared to the 16 mm COMS plaque for the depths up to 5 mm. Dose differences were not significant beyond depths of 10 mm due to ocular lateral scatter for the different plaque designs. Tumor DVHs for the 16 mm COMS plaque with Cs-131 provided better dose homogeneity and conformity compared to other COMS plaques with I-125 and Pd-103. Using Pd-103, DVHs for the 16 mm COMS plaque delivered less dose to outer sclera compared to other plaques. Conclusion: This study identified improved tumor homogeneity upon considering radionuclides and plaque designs, and found that scleral dose with the model EP917 plaque was higher than for the 16 mm COMS plaque for all the radionuclides studied.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Dosimetric characterization of two radium sources for retrospective dosimetry studies

Purpose: During the first part of the 20th century, {sup 226}Ra was the most used radionuclide for brachytherapy. Retrospective accurate dosimetry, coupled with patient follow up, is important for advancing knowledge on long-term radiation effects. The purpose of this work was to dosimetrically characterize two {sup 226}Ra sources, commonly used in Sweden during the first half of the 20th century, for retrospective dose–effect studies. Methods: An 8 mg {sup 226}Ra tube and a 10 mg {sup 226}Ra needle, used at Radiumhemmet (Karolinska University Hospital, Stockholm, Sweden), from 1925 to the 1960s, were modeled in two independent Monte Carlo (MC) radiation transport codes: GEANT4 and MCNP5. Absorbed dose and collision kerma around the two sources were obtained, from which the TG-43 parameters were derived for the secular equilibrium state. Furthermore, results from this dosimetric formalism were compared with results from a MC simulation with a superficial mould constituted by five needles inside a glass casing, placed over a water phantom, trying to mimic a typical clinical setup. Calculated absorbed doses using the TG-43 formalism were also compared with previously reported measurements and calculations based on the Sievert integral. Finally, the dose rate at large distances from a {sup 226}Ra point-like-source placed in the center of 1 m radius water sphere was calculated with GEANT4. Results: TG-43 parameters [including g{sub L}(r), F(r, θ), Λ, and s{sub K}] have been uploaded in spreadsheets as additional material, and the fitting parameters of a mathematical curve that provides the dose rate between 10 and 60 cm from the source have been provided. Results from TG-43 formalism are consistent within the treatment volume with those of a MC simulation of a typical clinical scenario. Comparisons with reported measurements made with thermoluminescent dosimeters show differences up to 13% along the transverse axis of the radium needle. It has been estimated that the uncertainty associated to the absorbed dose within the treatment volume is 10%–15%, whereas uncertainty of absorbed dose to distant organs is roughly 20%–25%. Conclusions: The results provided here facilitate retrospective dosimetry studies of {sup 226}Ra using modern treatment planning systems, which may be used to improve knowledge on long term radiation effects. It is surely important for the epidemiologic studies to be aware of the estimated uncertainty provided here before extracting their conclusions.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2013

Monte Carlo simulation based study of a proposed multileaf collimator for a telecobalt machine

Purpose: The objective of the present work was to propose a design of a secondary multileaf collimator (MLC) for a telecobalt machine and optimize its design features through Monte Carlo simulation. Methods: The proposed MLC design consists of 72 leaves (36 leaf pairs) with additional jaws perpendicular to leaf motion having the capability of shaping a maximum square field size of 35 Multiplication-Sign 35 cm{sup 2}. The projected widths at isocenter of each of the central 34 leaf pairs and 2 peripheral leaf pairs are 10 and 5 mm, respectively. The ends of the leaves and the x-jaws were optimized to obtain acceptable values of dosimetric and leakage parameters. Monte Carlo N-Particle code was used for generating beam profiles and depth dose curves and estimating the leakage radiation through the MLC. A water phantom of dimension 50 Multiplication-Sign 50 Multiplication-Sign 40 cm{sup 3} with an array of voxels (4 Multiplication-Sign 0.3 Multiplication-Sign 0.6 cm{sup 3}= 0.72 cm{sup 3}) was used for the study of dosimetric and leakage characteristics of the MLC. Output files generated for beam profiles were exported to the PTW radiation field analyzer software through locally developed software for analysis of beam profiles in order to evaluate radiation field width, beam flatness, symmetry, and beam penumbra. Results: The optimized version of the MLC can define radiation fields of up to 35 Multiplication-Sign 35 cm{sup 2} within the prescribed tolerance values of 2 mm. The flatness and symmetry were found to be well within the acceptable tolerance value of 3%. The penumbra for a 10 Multiplication-Sign 10 cm{sup 2} field size is 10.7 mm which is less than the generally acceptable value of 12 mm for a telecobalt machine. The maximum and average radiation leakage through the MLC were found to be 0.74% and 0.41% which are well below the International Electrotechnical Commission recommended tolerance values of 2% and 0.75%, respectively. The maximum leakage through the leaf ends in closed condition was observed to be 8.6% which is less than the values reported for other MLCs designed for medical linear accelerators. Conclusions: It is concluded that dosimetric parameters and the leakage radiation of the optimized secondary MLC design are well below their recommended tolerance values. The optimized design of the proposed MLC can be integrated into a telecobalt machine by replacing the existing adjustable secondary collimator for conformal radiotherapy treatment of cancer patients.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Design of an Yb-169 source optimized for gold nanoparticle-aided radiation therapy

Purpose: To find an optimum design of a new high-dose rate ytterbium (Yb)-169 brachytherapy source that would maximize the dose enhancement during gold nanoparticle-aided radiation therapy (GNRT), while meeting practical constraints for manufacturing a clinically relevant brachytherapy source. Methods: Four different Yb-169 source designs were considered in this investigation. The first three source models had a single encapsulation made of one of the following materials: aluminum, titanium, and stainless steel. The last source model adopted a dual encapsulation design with an inner aluminum capsule surrounding the Yb-core and an outer titanium capsule. Monte Carlo (MC) simulations using the Monte Carlo N-Particle code version 5 (MCNP5) were conducted initially to investigate the spectral changes caused by these four source designs and the associated variations in macroscopic dose enhancement across the tumor loaded with gold nanoparticles (GNPs) at 0.7% by weight. Subsequent MC simulations were performed using the EGSnrc and NOREC codes to determine the secondary electron spectra and microscopic dose enhancement as a result of irradiating the GNP-loaded tumor with the MCNP-calculated source spectra. Results: Effects of the source filter design were apparent in the current MC results. The intensity-weighted average energy of the Yb-169 source varied from 108.9 to 122.9 keV, as the source encapsulation material changed from aluminum to stainless steel. Accordingly, the macroscopic dose enhancement calculated at 1 cm away from the source changed from 51.0% to 45.3%. The sources encapsulated by titanium and aluminum/titanium combination showed similar levels of dose enhancement, 49.3% at 1 cm, and average energies of 113.0 and 112.3 keV, respectively. While the secondary electron spectra due to the investigated source designs appeared to look similar in general, some differences were noted especially in the low energy region (<50 keV) of the spectra suggesting the dependence of the photoelectron yield on the atomic number of source filter material, consistent with the macroscopic dose enhancement results. A similar trend was also shown in the so-called microscopic dose enhancement factor, for example, resulting in the maximum values of 138 and 119 for the titanium- and the stainless steel-encapsulated Yb-169 sources, respectively. Conclusions: The current results consistently show that the dose enhancement achievable from the Yb-169 source is closely related with the atomic number (Z) of source encapsulation material. While the observed range of improvement in the dose enhancement may be considered moderate after factoring all uncertainties in the MC results, the current study provides a reasonable support for the encapsulation of the Yb-core with lower-Z materials than stainless steel, for GNRT applications. Overall, the titanium capsule design can be favored over the aluminum or dual aluminum/titanium capsule designs, due to its superior structural integrity and improved safety during manufacturing and clinical use.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2015

Measurements of SNAC2 area dosimeters placed in different configurations around the PROSPERO reactor and comparison with TRIPOLI-4 calculations

In the context of criticality accident alarm system tests, several experiments were carried out in 2013 on the PROSPERO reactor to study the response to neutron and gamma of different devices and dosimeters, particularly on the SNAC2 dosimeter. This article presents the results of this criticality dosimeter in different configurations, and compares the experimental measurements with the results of calculation performed with the TRIPOLI-4 Monte-Carlo Neutral Particles transport code. PROSPERO is a metallic critical assembly managed by the Criticality, Neutron Science and Measurement Department located at the French CEA Research Center of Valduc. The core, surrounded by a reflector of depleted uranium, is composed of 2 horizontal cylindrical blocks made of a highly enriched uranium alloy which can be placed in contact, and of 4 depleted uranium control rods which allow the reactor to be driven. This reactor, placed in a cell 10 m x 8 m x 6 m high, with 1.4-meter-thick concrete walls, is used as a fast neutron spectrum source and is operated at stable power level in delayed critical state, which can vary from 3 mW to 3 kW. PROSPERO is extensively used for electronic hardening or to study the effect of the neutrons on various materials. The SNAC2 criticality dosimeter is a zone dosimeter allowing the off line measurement of criticality accident neutron doses. This dosimeter consists of the pile up of seven activation foils embedded into a 23 mm diameter x 21 mm height cadmium container. The activation measurement of each foil, using a gamma spectroscopy technique, gives information about the neutron reaction rates. The SNAC2 software allows the spectrum unfolding from these values, taking into account the hypothesis of a particular spectrum shape, in three components: a Maxwell spectrum component for the thermal range, a 1/E component for the epithermal range, and a Watt spectrum component for the high energy range. Moreover, from the neutron spectrum, the SNAC software can calculate the neutron fluence integrated by the dosimeter and the neutron dose. During the 3 weeks measurement campaign many radioprotection devices were used. To modify the spectrum seen by these devices, several shields of various thicknesses made of concrete or polyethylene, with or without cadmium covers, were placed in the PROSPERO cell. These devices allow the study of criticality accident spectra in several environments: from metal to pseudo liquid. The fluxes measured by the SNAC2 devices were compared with TRIPOLI-4 calculations. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2011

3D Monte Carlo Modeling of the EEDF in Negative Hydrogen Ion Sources

For optimization and accurate prediction of the amount of H{sup -} ion production in negative ion sources, analysis of electron energy distribution function (EEDF) is necessary. We are developing a numerical code which analyzes EEDF in tandem-type arc-discharge sources. It is a three-dimensional Monte Carlo simulation code with the realistic geometry and magnetic configuration. Coulomb collision between electrons is treated with 'Binary Collision' model and collisions with hydrogen species are treated with 'Null-collision (NC)' method. We have applied this code to the analysis of the JAEA 10 ampere negative ion source. The numerical result shows that the obtained EEDFs reasonably agree with experimental results.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2015

penORNL: a parallel Monte Carlo photon and electron transport package using PENELOPE

The parallel Monte Carlo photon and electron transport code package penORNL was developed at Oak Ridge National Laboratory to enable advanced scanning electron microscope (SEM) simulations on high-performance computing systems. This paper discusses the implementations, capabilities and parallel performance of the new code package. penORNL uses PENELOPE for its physics calculations and provides all available PENELOPE features to the users, as well as some new features including source definitions specifically developed for SEM simulations, a pulse-height tally capability for detailed simulations of gamma and x-ray detectors, and a modified interaction forcing mechanism to enable accurate energy deposition calculations. The parallel performance of penORNL was extensively tested with several model problems, and very good linear parallel scaling was observed with up to 512 processors. penORNL, along with its new features, will be available for SEM simulations upon completion of the new pulse-height tally implementation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Design and validation of a gamma-ray scanning system for measuring irradiated nuclear fuel

We report a prototype gamma-ray scanning system was designed to perform post-irradiation measurements of nuclear fuel at Idaho National Laboratory (INL). The system is composed of a coaxial high-purity germanium (HPGe) detector, a collimator, and mechanical positioning stages that translate the fuel sample across the front of the collimator. This paper documents the design process, modeling and the laboratory characterization performed to validate the system. The modeling of the system was done using the Monte Carlo N-Particle code with vendor-supplied design specifications, dimensional measurements, and x-ray radiographs of the HPGe detector. The characterization included measurements of calibrated radiation sources in different geometries, and calculations for absolute efficiency and dead time. Benchmark measurements were performed scanning an irradiated fuel rodlet containing ten pellets of 4.9% enriched UO2 in zirconium alloy. These results were compared to simulated spectra to further help characterize the detector model. One goal of the system is to determine the number of fissions per gram of UO2 in the fuel; comparison of the results from this prototype system's assay of the irradiated fuel rodlet are in agreement with previously estimated results.

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

Software stewardship and advancement of a high-performance computing scientific application: QMCPACK

Here, we provide an overview of the software engineering efforts and their impact in QMCPACK, a production-level ab-initio Quantum Monte Carlo open-source code targeting high-performance computing (HPC) systems. Aspects included are: (i) strategic expansion of continuous integration (CI) targeting CPUs, using GitHub Actions own runners, and NVIDIA and AMD GPUs used in pre-exascale systems, (ii) incremental reduction of memory leaks using sanitizers, (iii) incorporation of Docker containers for CI and reproducibility, and (iv) refactoring efforts to improve maintainability, testing coverage, and memory lifetime management. We quantify the value of these improvements by providing metrics to illustrate the shift towards a predictive, rather than reactive, maintenance approach. Our goal, in documenting the impact of these efforts on QMCPACK, is to contribute to the body of knowledge on the importance of research software engineering (RSE) for the stewardship and advancement of community HPC codes to enable scientific discovery at scale.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Direction-Modulated Brachytherapy for High-Dose-Rate Treatment of Cervical Cancer. I: Theoretical Design

Purpose: To demonstrate that utilization of the direction-modulated brachytherapy (DMBT) concept can significantly improve treatment plan quality in the setting of high-dose-rate (HDR) brachytherapy for cervical cancer. Methods and Materials: The new, MRI-compatible, tandem design has 6 peripheral holes of 1.3-mm diameter, grooved along a nonmagnetic tungsten-alloy rod (ρ = 18.0 g/cm{sup 3}), enclosed in Delrin tubing (polyoxymethylene, ρ = 1.41 g/cm{sup 3}), with a total thickness of 6.4 mm. The Monte Carlo N-Particle code was used to calculate the anisotropic {sup 192}Ir dose distributions. An in-house-developed inverse planning platform, geared with simulated annealing and constrained-gradient optimization algorithms, was used to replan 15 patient cases (total 75 plans) treated with a conventional tandem and ovoids (T and O) applicator. Prescription dose was 6 Gy. For replanning, we replaced the conventional tandem with that of the new DMBT tandem for optimization but left the ovoids in place and kept the dwell positions as originally planned. All DMBT plans were normalized to match the high-risk clinical target volume V100 coverage of the T and O plans. Results: In general there were marked improvements in plan quality for the DMBT plans. On average, D2cc for the bladder, rectum, and sigmoid were reduced by 0.59 ± 0.87 Gy (8.5% ± 28.7%), 0.48 ± 0.55 Gy (21.1% ± 27.2%), and 0.10 ± 0.38 Gy (40.6% ± 214.9%) among the 75 plans, with best single-plan reductions of 3.20 Gy (40.8%), 2.38 Gy (40.07%), and 1.26 Gy (27.5%), respectively. The high-risk clinical target volume D90 was similar, with 6.55 ± 0.96 Gy and 6.59 ± 1.06 Gy for T and O and DMBT, respectively. Conclusions: Application of the DMBT concept to cervical cancer allowed for improved organ at risk sparing while achieving similar target coverage on a sizeable patient population, as intended, by maximally utilizing the anatomic information contained in 3-dimensional imaging. A series of mechanical and clinical validations are to be followed.

62 RADIOLOGY AND NUCLEAR MEDICINE↗