Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mcnp”

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

Differential cross sections of the 16 O (𝑛,𝛼) reaction at neutron energies from 3.8 to 15 MeV

The importance of studying the 16 O (𝑛,𝛼) reaction is motivated by multiple nuclear applications which rely on an accurate nuclear reaction data library for oxygen. Discrepancies between past experimental data on the 16 O (𝑛,𝛼) 13 C reaction and its time-reverse 13 C ⁡(𝛼,𝑛) 16 O reaction have led to various different nuclear data evaluations. Here, we have measured 16 O (𝑛,𝛼) reaction cross sections using the LENZ instrument with the unmoderated white neutron source at LANSCE. Results from 2016/2017 data are discussed and used to benchmark the mcnp and geant simulations of the LENZ experimental setup. We report partial differential cross sections of 16 O (𝑛,𝛼 0 ) at 𝐸 𝑛 = 3.8–15 MeV and 16 O ⁡(𝑛,𝛼 1 +𝛼 2 +𝛼 3 ) at 𝐸 𝑛 =9–15 MeV, based on the new measurement in 2021. The resonances that we observed are in good agreement with the levels in 17 O that were previously measured. The LENZ cross sections are in better agreement with the ENDF/B-VIII.0 evaluation than the with reduced cross section found in ENDF/B-VII.1 up to 6 MeV. However, the current results for 16 O ⁡(𝑛,𝛼 0 ) 13 C g.s. appear to be in the best agreement over the entire energy range with that of JENDL/AN-2005 (ENDF/B-VI.0).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Radioisotope power source dose estimation tool (RPS-DET)

The purpose of this work is to develop a radioisotope power source dose estimation tool (RPS-DET) capable of rapidly predicting accurate neutron and gamma-ray flux/dose implications from various radioisotope power sources (RPSs). RPS-DET can use this information to calculate the flux/dose incident on or absorbed by nearby humans, materials, and electronics over a complex mission lifetime. This simulation tool uses a Matlab user interface to combine common geometry scenarios for traditional and new RPS design configurations to write a particle-transport simulation using Monte Carlo N-Particle (MCNP) 6, which provides results that are automatically post-processed for interactive analysis. The user defines the average age of the PuO2fuel, which determines the appropriate time-dependent gamma-ray and neutron spectra, along with the desired RPS and local environment for the simulation. Results for a given simulation are gamma-ray, neutron, and combined flux spectra, whole-body effective dose for humans, and absorbed dose in silicon for electronics. This work has been benchmarked against measurements of existing RPS units and has been determined to predict neutron and gamma flux, absorbed dose, and dose equivalent within 50% of simulation uncertainties.

Smith, Michael↗
DOE OSTI2015

Neutron and gamma sensitivities of self-powered detectors: Monte Carlo modelling

This paper deals with the development of a detailed Monte Carlo approach for the calculation of the absolute neutron sensitivity of SPNDs, which makes use of the MCNP code. We will explain the calculation approach, including the activation and beta emission steps, the gamma-electron interactions, the charge deposition in the various detector parts and the effect of the space charge field in the insulator. The model can also be applied for the calculation of the gamma sensitivity of self-powered detectors and for the radiation-induced currents in signal cables. The model yields detailed information on the various contributions to the sensor currents, with distinct response times. Results for the neutron sensitivity of various types of SPNDs are in excellent agreement with experimental data obtained at the BR2 research reactor. For typical neutron to gamma flux ratios, the calculated gamma induced SPND currents are significantly lower than the neutron induced currents. The gamma sensitivity depends very strongly upon the immediate detector surroundings and on the gamma spectrum. Our calculation method opens the way to a reliable on-line determination of the absolute in-pile thermal neutron flux. (authors)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

Test of a prototype neutron spectrometer based on diamond detectors in a fast reactor

A prototype of neutron spectrometer based on diamond detectors has been developed. This prototype consists of a {sup 6}Li neutron converter sandwiched between two CVD diamond crystals. The radiation hardness of the diamond crystals makes it suitable for applications in low power research reactors, while a low sensitivity to gamma rays and low leakage current of the detector permit to reach good energy resolution. A fast coincidence between two crystals is used to reject background. The detector was read out using two different electronic chains connected to it by a few meters of cable. The first chain was based on conventional charge-sensitive amplifiers, the other used a custom fast charge amplifier developed for this purpose. The prototype has been tested at various neutron sources and showed its practicability. In particular, the detector was calibrated in a TRIGA thermal reactor (LENA laboratory, University of Pavia) with neutron fluxes of 10{sup 8} n/cm{sup 2}s and at the 3 MeV D-D monochromatic neutron source named FNG (ENEA, Rome) with neutron fluxes of 10{sup 6} n/cm{sup 2}s. The neutron spectrum measurement was performed at the TAPIRO fast research reactor (ENEA, Casaccia) with fluxes of 10{sup 9} n/cm{sup 2}s. The obtained spectra were compared to Monte Carlo simulations, modeling detector response with MCNP and Geant4. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2015

Neutron Measurement Instrumentation Development at KIT for the European ITER TBM

Fusion power reactors will rely on the internal production of the fuel tritium from lithium in the tritium breeding blanket. Test Blanket Modules (TBM) will be installed in ITER with the aim to investigate the nuclear performance of different breeding blanket designs. Currently there is no fully qualified nuclear instrumentation available for the measurement of neutron fluxes and tritium production rates which would be able to withstand the harsh environment conditions in the TBM such as high temperature (>400 deg. C) and, depending on the operation scenario, intense radiation levels. As partner of the European Consortium on Nuclear Data and Measurement Techniques in the framework of several F4E specific grants and contracts, KIT and ENEA have jointly studied the possibility to develop and test detectors suitable to operate in ITER-TBMs. Here we present an overview of ongoing work on three types of neutron flux monitors under development for the TBMs with focus on the KIT activities. A neutron activation system (NAS) with pneumatic sample transport could provide absolute neutron flux measurements in selected positions. A test system for investigating activation materials with short half-lives was constructed at the DT neutron generator laboratory of Technical University of Dresden to investigate the neutronics aspects. Several irradiations have been performed with focus on the simultaneous measurement of the extracted activated probes. An engineering assessment of a TBM NAS in the conceptual design phase has been done which considered issues of design requirements and integration. Last but not least, a mechanical test bench is under construction at KIT which will address issues of driving the activation probes, solutions for loading the system etc. experimentally. Self-powered neutron detectors (SPND) are widely applied in fission reactor monitoring, and the commercially available SPNDs are sensitive to thermal neutrons. We are investigating novel materials for SPND which would be sensitive also to the fast neutron flux expected in the TBMs. To this end simulations were done with the European Activation System EASY and neutron flux spectra which were calculated with MCNP for the HCPB TBM. Preliminary tests with commercial SPND in a fast reactor were performed. As a result of these activities, several materials have been found which may be suitable for the measurement of fast neutron fluxes in the TBM. Test detectors are under preparation for testing with DT neutron generators. Within the I{sub S}MART project, funded by KIC InnoEnergy, KIT is developing an online detector based on silicon carbide electronics for the TBMs. The operation of such detectors at TBM relevant temperatures is expected to incur lower accumulated radiation damage to them than at room temperature due to annealing effects. Detectors of several designs have been already irradiated with DT neutrons. Irradiation tests at elevated temperatures have been done and further tests are currently underway. This paper summarizes the status of the work for these three neutron flux monitor systems. (authors)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Development of nuclear data processing code FRENDY version 2

Nuclear data processing is an important interface between an evaluated nuclear data library and neutronics calculation codes. JAEA has been developed the new nuclear data processing code FRENDY from 2013. FRENDY version 1 generates the ACE files used for the continuous-energy Monte Carlo codes including PHITS, Solomon, Serpent,and MCNP and it was released as the open-source software under the 2-clause BSD license in 2019.After we released FRENDY version 1, many functions, e.g., the multi-group neutron cross-section library generation, the statistical uncertainty quantification of the probability tables for unresolved resonance cross-section, the 'perturbation' of the ACE file for uncertainty quantification, and the modification of the ENDF-6 formatted nuclear data file, have been developed. We released FRENDY version 2 including these functions. The present paper gives an overview of FRENDY and presents features of new functions implemented in FRENDY version 2. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2017

Fresh Fuel Measurements With the Differential Die-Away Self-Interrogation Instrument

The purpose of the Next Generation Safeguards Initiative (NGSI)-Spent Fuel (SF) Project is to strengthen the technical toolkit of safeguards inspectors and/or other interested parties. Thus the NGSI-SF team is working to achieve the following technical goals more easily and efficiently than in the past using nondestructive assay measurements of spent fuel assemblies: 1) verify the initial enrichment, burnup, and cooling time of facility declaration; 2) detect the diversion or replacement of pins; 3) estimate the plutonium mass; 4) estimate decay heat; and 5) determine the reactivity of spent fuel assemblies. The differential die-away self-interrogation (DDSI) instrument is one instrument that was assessed for years regarding its feasibility for robust, timely verification of spent fuel assemblies. The instrument was recently built and was tested using fresh fuel assemblies in a variety of configurations, including varying enrichment, neutron absorber content, and symmetry. The early die-away method, a multiplication determination method developed in simulation space, was successfully tested on the fresh fuel assembly data and determined multiplication with a root-mean-square (RMS) error of 2.9%. The experimental results were compared with MCNP simulations of the instrument as well. Low multiplication assemblies had agreement with an average RMS error of 0.2% in the singles count rate (i.e., total neutrons detected per second) and 3.4% in the doubles count rates (i.e., neutrons detected in coincidence per second). High-multiplication assemblies had agreement with an average RMS error of 4.1% in the singles and 13.3% in the doubles count rates.

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

Neutron Yield of Thermo Scientific P385 D-T Neutron Generator vs . Current and Voltage

The Thermo Scientific P385 Neutron Generator is a compact neutron source, producing 14 MeV neutrons through the deuterium-tritium (DT) fusion reaction. It is important to measure and understand the dependence of the neutron production rate on the accelerator current and voltage. In this study we evaluated neutron production with an absolutely calibrated liquid scintillator neutron spectrometer (BTI N-Probe), an absolutely calibrated He-3 detector surrounded by HDPE shells (Detec Nested Neutron Spectrometer, NNS), and two uncalibrated ZnS fast neutron scintillators (EJ-410), for both A3082 and A3083 sealed tubes. Here we also modeled the neutron yield using the TRIM code, which calculates the trajectory and the energy loss of deuterons and tritons within the target. Experimental results showed an essentially linear dependence on beam current, as expected. A 3.59 ±0.08 power law dependence on the operating voltage was measured, in effective agreement with the modeled value of 3.5. A series of absolute NNS and N-Probe measurements, matched against MCNP calculations, showed that the A3083 and A3082 tubes provide a maximum neutron yield of 8.2 × 10 8 n/s and 4.7 × 10 8 n/s respectively, with estimated uncertainty of ±10%.We showed, through modeling, that tritium decay is not a significant consideration for tubes, such as these, with lifetimes of less than 10 years.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Optimization of an Energy Tuning Assembly for High Explosives Detection

The Portable Isotopic Neutron Spectroscopy (PINS) system, employs neutron-induced gamma-ray spectroscopy and provides a nondestructive method for high explosives detection. In standard operation it uses Californium-252 as a neutron source. Operating PINS with a deuterium-tritium (DT) neutron generator has some advantages over Cf-252, including lifetime and ability to produce high-energy inelastic scattering gamma rays. However, current systems using DT neutron generators suffer from a high environmental background and reduced ability to induce neutron capture, reducing spectral quality and limiting nitrogen sensitivity. Here, this study presents the development of an energy-tuning assembly (ETA) designed to optimize the DT neutron energy spectrum to increase nitrogen reaction rates in a target, thereby improving high explosive detection capabilities. A metaheuristic optimization framework, MultiGNOWEE, coupled with MCNP, was employed to generate two ETA configurations: a single-objective ETA, which maximizes nitrogen capture reactions, and a multi-objective ETA, which balances neutron capture and inelastic scattering. Simulations demonstrated the optimized configurations achieved up to a 10-fold improvement in nitrogen capture rates compared to the bare configuration. Experimental validation was conducted using a DT neutron generator and a high-purity germanium (HPGe) detector. Two prototype ETAs were constructed and assessed on a melamine simulant. Measurements demonstrated improved nitrogen detection for both prototype ETA configurations when compared to the standard system.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Numerical Simulation of Irradiation Induced Swelling for STAR 4.1 Blanket

Here, this article presents the preliminary results of FEM implementation of irradiation-induced swelling, hardening, and creep effect using ANSYS user programmable features (UPFs) applied to material properties. The first liquid metal breeding blanket model for STAR 4.1 tokamak as a part of the virtual prototyping system is under development at Princeton Plasma Physics Laboratory (PPPL). Intense neutron irradiation produces significant changes in the physical and mechanical properties of Fe-(8%–9%)Cr-based reduced activation ferritic martensitic (RAFM) steels. An ANSYS model of the breeding blanket was built to involve these irradiation-induced material property changes with neutron fluence distribution mapped from the MCNP model. Simulation results show that significant structure deformation forms from nuclear swelling, stress, and plastic strain arising from DPA gradient through wall thickness. The material degradation effect is not negligible. This modeling feature can either assist in analyzing the structural behaviors with the influence of nuclear swelling or provide guidance to design the structure to withstand irradiation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Investigation of the effects of cell model and subcellular location of gold nanoparticles on nuclear dose enhancement factors using Monte Carlo simulation

Purpose: The authors’ aims were to model how various factors influence radiation dose enhancement by gold nanoparticles (AuNPs) and to propose a new modeling approach to the dose enhancement factor (DEF).Methods: The authors used Monte Carlo N-particle (MCNP 5) computer code to simulate photon and electron transport in cells. The authors modeled human breast cancer cells as a single cell, a monolayer, or a cluster of cells. Different numbers of 5, 30, or 50 nm AuNPs were placed in the extracellular space, on the cell surface, in the cytoplasm, or in the nucleus. Photon sources examined in the simulation included nine monoenergetic x-rays (10–100 keV), an x-ray beam (100 kVp), and {sup 125}I and {sup 103}Pd brachytherapy seeds. Both nuclear and cellular dose enhancement factors (NDEFs, CDEFs) were calculated. The ability of these metrics to predict the experimental DEF based on the clonogenic survival of MDA-MB-361 human breast cancer cells exposed to AuNPs and x-rays were compared.Results: NDEFs show a strong dependence on photon energies with peaks at 15, 30/40, and 90 keV. Cell model and subcellular location of AuNPs influence the peak position and value of NDEF. NDEFs decrease in the order of AuNPs in the nucleus, cytoplasm, cell membrane, and extracellular space. NDEFs also decrease in the order of AuNPs in a cell cluster, monolayer, and single cell if the photon energy is larger than 20 keV. NDEFs depend linearly on the number of AuNPs per cell. Similar trends were observed for CDEFs. NDEFs using the monolayer cell model were more predictive than either single cell or cluster cell models of the DEFs experimentally derived from the clonogenic survival of cells cultured as a monolayer. The amount of AuNPs required to double the prescribed dose in terms of mg Au/g tissue decreases as the size of AuNPs increases, especially when AuNPs are in the nucleus and the cytoplasm. For 40 keV x-rays and a cluster of cells, to double the prescribed x-ray dose (NDEF = 2) using 30 nm AuNPs, would require 5.1 ± 0.2, 9 ± 1, 10 ± 1, 10 ± 1 mg Au/g tissue in the nucleus, in the cytoplasm, on the cell surface, or in the extracellular space, respectively. Using 50 nm AuNPs, the required amount decreases to 3.1 ± 0.3, 8 ± 1, 9 ± 1, 9 ± 1 mg Au/g tissue, respectively.Conclusions: NDEF is a new metric that can predict the radiation enhancement of AuNPs for various experimental conditions. Cell model, the subcellular location and size of AuNPs, and the number of AuNPs per cell, as well as the x-ray photon energy all have effects on NDEFs. Larger AuNPs in the nucleus of cluster cells exposed to x-rays of 15 or 40 keV maximize NDEFs.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2014

Feasibility of using intermediate x-ray energies for highly conformal extracranial radiotherapy

Purpose: To investigate the feasibility of using intermediate energy 2 MV x-rays for extracranial robotic intensity modulated radiation therapy. Methods: Two megavolts flattening filter free x-rays were simulated using the Monte Carlo code MCNP (v4c). A convolution/superposition dose calculation program was tuned to match the Monte Carlo calculation. The modeled 2 MV x-rays and actual 6 MV flattened x-rays from existing Varian Linacs were used in integrated beam orientation and fluence optimization for a head and neck, a liver, a lung, and a partial breast treatment. A column generation algorithm was used for the intensity modulation and beam orientation optimization. Identical optimization parameters were applied in three different planning modes for each site: 2, 6 MV, and dual energy 2/6 MV. Results: Excellent agreement was observed between the convolution/superposition and the Monte Carlo calculated percent depth dose profiles. For the patient plans, overall, the 2/6 MV x-ray plans had the best dosimetry followed by 2 MV only and 6 MV only plans. Between the two single energy plans, the PTV coverage was equivalent but 2 MV x-rays improved organs-at-risk sparing. For the head and neck case, the 2MV plan reduced lips, mandible, tongue, oral cavity, brain, larynx, left and right parotid gland mean doses by 14%, 8%, 4%, 14%, 24%, 6%, 30% and 16%, respectively. For the liver case, the 2 MV plan reduced the liver and body mean doses by 17% and 18%, respectively. For the lung case, lung V20, V10, and V5 were reduced by 13%, 25%, and 30%, respectively. V10 of heart with 2 MV plan was reduced by 59%. For the partial breast treatment, the 2 MV plan reduced the mean dose to the ipsilateral and contralateral lungs by 27% and 47%, respectively. The mean body dose was reduced by 16%. Conclusions: The authors showed the feasibility of using flattening filter free 2 MV x-rays for extracranial treatments as evidenced by equivalent or superior dosimetry compared to 6 MV plans using the same inverse noncoplanar intensity modulated planning method.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

Absolute calibration of the Gamma Knife{sup ®} Perfexion™ and delivered dose verification using EPR/alanine dosimetry

Purpose: Elekta Leksell Gamma Knife{sup ®} (LGK) is a radiotherapy beam machine whose features are not compliant with the international calibration protocols for radiotherapy. In this scope, the Laboratoire National Henri Becquerel and the Pitié-Salpêtrière Hospital decided to conceive a new LKG dose calibration method and to compare it with the currently used one. Furthermore, the accuracy of the dose delivered by the LGK machine was checked using an “end-to-end” test. This study also aims to compare doses delivered by the two latest software versions of the Gammaplan treatment planning system (TPS). Methods: The dosimetric method chosen is the electron paramagnetic resonance (EPR) of alanine. Dose rate (calibration) verification was done without TPS using a spherical phantom. Absolute calibration was done with factors calculated by Monte Carlo simulation (MCNP-X). For “end-to-end” test, irradiations in an anthropomorphic head phantom, close to real treatment conditions, are done using the TPS in order to verify the delivered dose. Results: The comparison of the currently used calibration method with the new one revealed a deviation of +0.8% between the dose rates measured by ion chamber and EPR/alanine. For simple fields configuration (less than 16 mm diameter), the “end-to-end” tests showed out average deviations of −1.7% and −0.9% between the measured dose and the calculated dose by Gammaplan v9 and v10, respectively. Conclusions: This paper shows there is a good agreement between the new calibration method and the currently used one. There is also a good agreement between the calculated and delivered doses especially for Gammaplan v10.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

SU-E-I-42: Normalized Embryo/fetus Doses for Fluoroscopically Guided Pacemaker Implantation Procedures Calculated Using a Monte Carlo Technique

Purpose: It is well known that pacemaker implantation is sometimes needed in pregnant patients with symptomatic bradycardia. To our knowledge, there is no reported experience regarding radiation doses to the unborn child resulting from fluoroscopy during pacemaker implantation. The purpose of the current study was to develop a method for estimating embryo/fetus dose from fluoroscopically guided pacemaker implantation procedures performed on pregnant patients during all trimesters of gestation. Methods: The Monte Carlo N-Particle (MCNP) radiation transport code was employed in this study. Three mathematical anthropomorphic phantoms representing the average pregnant patient at the first, second and third trimesters of gestation were generated using Bodybuilder software (White Rock science, White Rock, NM). The normalized embryo/fetus dose from the posteroanterior (PA), the 30° left-anterior oblique (LAO) and the 30° right-anterior oblique (RAO) projections were calculated for a wide range of kVp (50–120 kVp) and total filtration values (2.5–9.0 mm Al). Results: The results consist of radiation doses normalized to a) entrance skin dose (ESD) and b) dose area product (DAP) so that the dose to the unborn child from any fluoroscopic technique and x-ray device used can be calculated. ESD normalized doses ranged from 0.008 (PA, first trimester) to 2.519 μGy/mGy (RAO, third trimester). DAP normalized doses ranged from 0.051 (PA, first trimester) to 12.852 μGy/Gycm2 (RAO, third trimester). Conclusion: Embryo/fetus doses from fluoroscopically guided pacemaker implantation procedures performed on pregnant patients during all stages of gestation can be estimated using the method developed in this study. This study was supported by the Greek Ministry of Education and Religious Affairs, General Secretariat for Research and Technology, Operational Program ‘Education and Lifelong Learning’, ARISTIA (Research project: CONCERT)

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-I-76: Matching Primary and Scattered X-Ray Spectra for Use in Calculating the Diagnostic Radiation Index of Protection

Purpose: Lightweight lead-free or lead-composite protective garments exploit k-edge interactions to attenuate scattered X-rays. Manufacturers specify the protective value of garments in terms of lead equivalence at a single kVp. This is inadequate, as the protection provided by such garments varies with radiation quality in different use conditions. We present a method for matching scattered X-ray spectra to primary X-ray spectra. The resulting primary spectra can be used to measure penetration through protective garments, and such measurements can be weighted and summed to determine a Diagnostic Radiation Index for Protection (DRIP). Methods: Scattered X-ray spectra from fluoroscopic procedures were modeled using Monte Carlo techniques in MCNP-X 2.7. Data on imaging geometry, operator position, patient size, and primary beam spectra were gathered from clinical fluoroscopy procedures. These data were used to generate scattered X-ray spectra resulting from procedural conditions. Technical factors, including kV and added filtration, that yielded primary X-ray spectra that optimally matched the generated scattered X-ray spectra were identified through numerical optimization using a sequential quadratic programming (SQP) algorithm. Results: The primary spectra generated with shape functions matched the relative flux in each bin of the scattered spectra within 5%, and half and quarter-value layers matched within 0.1%. The DRIP for protective garments can be determined by measuring the penetration through protective garments using the matched primary spectra, then calculating a weighted average according to the expected clinical use of the garment. The matched primary spectra are specified in terms of first and second half-value layers in aluminum and acrylic. Conclusion: Lead equivalence is inadequate for completely specifying the protective value of garments. Measuring penetration through a garment using full scatter conditions is very difficult. The primary spectra determined in this work allow for practical primary penetration measurements to be made with equipment readily available to clinical medical physicists.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2014

SU-E-T-90: Concrete Forward-Scatter Fractions for Radiotherapy Shielding Applications

Purpose: There is little instruction within the primary shielding guidance document NCRP 151 for vault designs where the primary beam intercepts the maze. We have conducted a Monte-Carlo study to characterize forward-scattered radiation from concrete barriers with the intent of quantifying what amount of additional shielding outside the primary beam is needed in this situation. Methods: We reproduced our vault in MCNP 5 and simulated spectra obtained from the literature and from our treatment planning system for 10 and 18 MV beams. Neutron and gamma-capture contributions were not simulated. Energy deposited was scored at isocenter in a water phantom, within various cells that comprised the maze, and within cells that comprised the vault door. Tracks were flagged that scattered from within the maze to the door and their contributions were tallied separately. Three different concrete mixtures found in the literature were simulated. An empirically derived analytic equation was used for comparison, utilizing patient scatter fractions to approximate the scatter from concrete. Results: Our simulated data confirms that maze-scattered radiation is a significant contribution to total photon dose at the door. It contributes between 20-35% of the photon shielding workload. Forward-scatter fractions for concrete were somewhat dependent on concrete composition and the relative abundance of higher-Z elements. Scatter fractions were relatively insensitive to changes in the primary photon spectrum. Analytic results were of the same magnitude as simulated results. Conclusions: Forward-scattered radiation from the maze barrier needs to be included in the photon workload for shielding calculations in non-standard vault designs. Scatter fractions will vary with concrete composition, but should be insensitive to spectral changes between machine manufacturers. Further plans for investigation include refined scatter fractions for various concrete compositions, scatter fraction measurement, and simulation of photo-nuclear interaction within the maze barrier for high-energy beams.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2014

Independent dosimetric assessment of the model EP917 episcleral brachytherapy plaque

Purpose: To investigate the influence of slot design on dose distributions and dose-volume histograms (DVHs) for the model EP917 plaque for episcleral brachytherapy. Methods: Dimensions and orientations of the slots were measured for three model EP917 plaques and compared to data in the Plaque Simulator (PS) treatment planning software (version 5.7.6). These independently determined coordinates were incorporated into the MCNP Monte Carlo simulation environment to obtain dose from the plaques in a water environment and in a clinical environment with ocular structures. A tumor volume was simulated as 5 mm in apical height and 11 mm in basal diameter. Variations in plaque mass density and composition; slot length, width, and depth; seed positioning; and Ag-marker rod positioning were simulated to examine their influence on plaque central axis (CAX) and planar dose distributions, and DVHs. Results: Seed shifts in a single slot toward the eye and shifts of the{sup 125}I-coated Ag rod within the capsule had the greatest impact on CAX dose distribution. A shift of 0.0994 mm toward the eye increased dose by 14%, 9%, 4.3%, and 2.7% at 1, 2, 5, and 10 mm, respectively, from the inner sclera. When examining the fully-modeled plaque in the ocular geometry, the largest dose variations were caused by shifting the Ag rods toward the sclera and shifting the seeds away from the globe when the slots were made 0.51 mm deeper, causing +34.3% and −69.4% dose changes to the outer sclera, respectively. At points along the CAX, dose from the full plaque geometry using the measured slot design was 2.4% ± 1.1% higher than the manufacturer-provided slot design and 2.2% ± 2.3% higher than the homogeneous calculation of PS treatment planning results. The ratio of D{sub 10} values for the measured slot design to the D{sub 10} values for the manufacturer-provided slot design was higher by 9%, 10%, and 19% for the tumor, inner sclera, and outer sclera, respectively. In comparison to the measured slot design, a theoretical plaque having narrower and deeper slots delivered 30%, 37%, and 62% lower D{sub 10} doses to the tumor, inner sclera, and outer sclera, respectively. Conclusions: While the measured positions of the slots on the model EP917 plaque were in close agreement (<0.7 mm) with the PS values, small differences in the slot shape caused substantial differences in dose distributions and DVH metrics. Increasing slot depth by 0.1 mm decreased outer scleral dose by 20%, yet shifting the Ag rods in the seeds toward the globe by 0.1 mm increased outer scleral dose by 35%. The clinical medical physicist is advised to measure these types of plaques upon acceptance testing before clinical use to inspect slot shape and position for comparison with data used for treatment planning purposes.

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↗