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,351 records · Page 75

SCALE Enhancements for Detailed Cask Dose Rate Analysis - Paper 37

The Nuclear Fuels Storage and Transportation Planning Project and the Used Fuel Disposition Campaign, both under the US Department of Energy (DOE) Office of Nuclear Energy, require the ability to predict k{sub eff}, heat generation and temperature distributions, and external dose rates from any combination of commercial spent nuclear fuel assemblies loaded into any type of licensed cask. Work in the Reactor and Nuclear Systems Division at Oak Ridge National Laboratory has begun to develop the UNF-ST[and]DARDS system (Used Nuclear Fuel Storage, Transportation and Disposal Analysis Resource and Data System) that uses many of the nuclear analysis codes from the SCALE package, a comprehensive database of used nuclear fuel inventories from utilities, and a library of model template files to produce complete input files for the various types of UNF computational analyses. UNF-ST[and]DARDS uses the template files and input parameters that are applicable to a safety analysis to assemble and execute a complete input file, collect the output files and relevant information for subsequent calculations, and illustrate calculation results as a function of decay time and other relevant analysis parameters. This paper will focus on the shielding aspects and the ability to predict external dose rates from as-loaded casks. Realistic dose rate values based on actual canister fuel contents may be used in developing appropriate plans, procedures, and methods to keep personnel radiation exposures as low as reasonably achievable during cask repackaging and transportation operations, as well as in performing total system analyses and transportation risk evaluations. The MAVRIC shielding sequence of the SCALE package was created for the US Nuclear Regulatory Commission for analyzing dose rates outside of storage, transfer, or transportation casks loaded with spent nuclear fuel. MAVRIC contains automated variance reduction capabilities based on the CADIS and FW-CADIS methods that are effective for optimizing deep penetration Monte Carlo shielding problems for achieving low relative uncertainties in short computing times. MAVRIC uses the Denovo discrete ordinates code to compute importance maps and biased source terms for the Monaco fixed-source Monte Carlo code, which can use multi-group or continuous-energy cross section data. ORIGEN (Oak Ridge Isotope Generation and Depletion) is the SCALE computer code for calculating time-dependent isotopic concentrations during irradiation and decay for more than 2,200 nuclides that can be produced by nuclear fuel irradiation and activation. A series of SCALE enhancements was implemented to enable automatic decay heat and radiation source term generation and shielding safety analysis of as-loaded casks. The enhancements to SCALE made under this project will be a part of the SCALE 6.2 release expected in fall 2014. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Design and performance testing of a novel far field Gamma system to assay radioactive waste in 400 and 200 litre drums - 15339

A Far Field Gamma ray Measurement system has been developed for measuring and assaying both 400 and 200 litre waste drums with a wide range of activity at nuclear power stations. These include radioactive waste drums with surface dose-rates of less than 2 mSv/h (200 mrem/h) and, at the higher range of activity, 400 litre drums with surface dose-rates of up to 50 mSv/h (5 rem/h). The assay system incorporates an electro-mechanically cooled high purity Germanium profile detector and advanced digital spectroscopic analysis electronics. A novel optimised conical-trapezoidal collimator has been incorporated to reduce the effect of the radioactive background in addition to the conventional detector lead shield with tin and copper graded lining to reduce the effect of lead X-rays. Two features have been incorporated to reduce detector dead time and maintain detector resolution when measuring high activity drums. The first is the inclusion of two automated tungsten filters of different thickness as an integral part of the detector collimator. The filters are used in conjunction with an automated rail system on which the detector platform is mounted. The detector can be positioned at different distances from the drum surface in order to reduce detector dead time. The deployment of the filters and the adjustment of the drum - detector position are automatically controlled based on a user defined dose-rate limit table. The data is provided by two Geiger-Muller dose-rate detectors, which are used to measure the drum surface dose-rate. This data is subsequently used to adjust the detector position and in cases of higher dose-rate deploy one or other of the tungsten filters. The drum rotation platform incorporates a load cell to determine drum weight and hence the drum density. This information is used to calculate an attenuation correction based on the assumption of uniform drum density and uniform distribution of activity within the drum. Where the waste drums have a more complicated and known regular internal structure, such as a small radioactive region surrounded by an annulus of shielding material, the analysis algorithms incorporated into the far field geometry and attenuation correction spectroscopic data analysis code are able to make a more accurate determination of drum activity than would be possible assuming a completely uniform drum. The detector energy and efficiency calibration is achieved using a point source with multiple gamma ray energy peaks. The response of the high purity Germanium detector is also modelled using the Monte Carlo Neutron Photon code and the model is benchmarked using the measured point source calibration data. The performance of the instrument has been determined for a range of uniform density matrix drums containing a set of Eu-152 line or rod sources located in re-entrant tubes positioned on an equal volume basis. When rotated, the test drums with volume distributed line sources simulate waste drums of different density each with a uniformly distributed radioactive source. The measurements have been validated using Monte Carlo simulations. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

MONTE CARLO ANALYSES OF THE YALINA THERMAL FACILITY WITH SERPENT STEREOLITHOGRAPHY GEOMETRY MODEL

This paper analyzes the YALINA Thermal subcritical assembly of Belarus using two different Monte Carlo transport programs, SERPENT and MCNP. The MCNP model is based on combinatorial geometry and universes hierarchy, while the SERPENT model is based on Stereolithography geometry. The latter consists of unstructured triangulated surfaces defined by the normal and vertices. This geometry format is used by 3D printers and it has been created by: the CUBIT software, MATLAB scripts, and C coding. All the Monte Carlo simulations have been performed using the ENDF/B-VII.0 nuclear data library. Both MCNP and SERPENT share the same geometry specifications, which describe the facility details without using any material homogenization. Three different configurations have been studied with different number of fuel rods. The three fuel configurations use 216, 245, or 280 fuel rods, respectively. The numerical simulations show that the agreement between SERPENT and MCNP results is within few tens of pcms.

CUBIT↗

The MCNP ® 6 code: A decade of progress

After several years of effort involved in merging the Los Alamos National Laboratory MCNP5 and MCNPX codes, in 2013 the first production release of version 6 of the Monte Carlo N-Particle ® , or MCNP ® , code MCNP6.1 was distributed publicly. Since then, three significant releases have been issued: MCNP6.1.1beta in 2014, MCNP6.2 in 2018, and MCNP6.3 in 2023. While each release always contains new features, code enhancements, and bug fixes, each version has had a different primary focus, ranging from improved calculational efficiency to new powerful utilities and tools, to software modernization of the code base. With all that has been learned over the first decade of the MCNP6 code, continuous progress is being made toward a modernized, general-purpose Monte Carlo radiation transport code that remains a trusted resource for the global community of practitioners. This paper describes these first 10+ years of the MCNP6 code and its continually improving data libraries, and gives some insight into how the next decade is expected to unfold.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Demonstration of the On-the-Fly Shielding Analysis Method: Spent Fuel and Waste Disposition

This report documents work performed supporting the US Department of Energy (DOE) Office of Nuclear Energy (NE) Spent Fuel and Waste Disposition (SFWD) Integrated Waste Management activities under work breakdown structure element 1.08.02.04.01, “Data and Tools Development, Validation, and Maintenance.” In particular, this report fulfills milestone M3SF-21OR020401016, “Implement on-the-fly dose analysis methodology in UNF-ST&DARDS” within work package SF-21OR02040101, “Commercial SNF Characterization - ORNL.” The Used Nuclear Fuel - Storage, Transportation & Disposal Analysis Resource and Data System (UNFST& DARDS) enables automated dose rate calculations for spent nuclear fuel (SNF) transportation packages and storage casks using a Monte Carlo radiation transport code. The explicit method uses a detailed model of the SNF system and its contents. Therefore, a dose rate calculation is required for each as-loaded transportation package or storage cask because the SNF assemblies within a canister typically have unique irradiation characteristics. An alternate method, referred to as the “on-the-fly” shielding analysis method, has been proposed that requires only a set of Monte Carlo dose rate calculations for each transportation packaging/storage cask design. The results of the Monte Carlo dose rate calculations are independent of the SNF assembly irradiation and decay characteristics. The dose rate values may then be combined with the radiation source strength of the SNF assemblies associated with a particular transportation packaging/storage cask design to determine actual dose rates. This report presents on-the-fly dose rate calculations for a representative SNF storage cask and verification of the on-the-fly dose rate calculation results by comparison with reference dose rate calculations using the explicit Monte Carlo dose rate calculation. The on-the-fly shielding analysis method was implemented in UNF-ST&DARDS. A Python program was developed to process the MAVRIC dose rate results obtained by source particle type, energy group, and fuel geometry region. A Python processor created binary files, which were saved as a special UNF-ST&DARDS library for on-the- fly shielding analyses. UNF-ST&DARDS uses the precalculated on-the-fly binary libraries generated by the Python data processor and directly executes the Python code for on-the-fly dose analysis. This Python code unzips the pre-generated binary files mentioned above, reads the data, and combines them with user-specified sources for dose and uncertainty calculations. The Python programs were verified using Excel calculations and by comparison with the values obtained with the MAVRIC post-processing utilities applied to the 3dmap files. This method can currently be used to determine dose rates for as-loaded HI-STORM FW storage casks. The UNF-ST&DARDS analysis wizard for on-the-fly shielding analysis is described in this report.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Determination of transit dose profile for a {sup 192}Ir HDR source

Purpose: Several studies have reported methodologies to calculate and correct the transit dose component of the moving radiation source for high dose rate (HDR) brachytherapy planning systems. However, most of these works employ the average source speed, which varies significantly with the measurement technique used, and does not represent a realistic speed profile, therefore, providing an inaccurate dose determination. In this work, the authors quantified the transit dose component of a HDR unit based on the measurement of the instantaneous source speed to produce more accurate dose values. Methods: The Nucletron microSelectron-HDR Ir-192 source was characterized considering the Task Group 43 (TG-43U1) specifications. The transit dose component was considered through the calculation of the dose distribution using a Monte Carlo particle transport code, MCNP5, for each source position and correcting it by the source speed. The instantaneous source speed measurements were performed in a previous work using two optical fibers connected to a photomultiplier and an oscilloscope. Calculated doses were validated by comparing relative dose profiles with those obtained experimentally using radiochromic films. Results: TG-43U1 source parameters were calculated to validate the Monte Carlo simulations. These agreed with the literature, with differences below 1% for the majority of the points. Calculated dose profiles without transit dose were also validated by comparison with ONCENTRA{sup Registered-Sign} Brachy v. 3.3 dose values, yielding differences within 1.5%. Dose profiles obtained with MCNP5 corrected using the instantaneous source speed profile showed differences near dwell positions of up to 800% in comparison to values corrected using the average source speed, but they are in good agreement with the experimental data, showing a maximum discrepancy of approximately 3% of the maximum dose. Near a dwell position the transit dose is about 22% of the dwell dose delivered by the source dwelling 1 s and reached 104.0 cGy per irradiation in a hypothetical clinical case studied in this work. Conclusions: The present work demonstrated that the transit dose correction based on average source speed fails to accurately correct the dose, indicating that the correct speed profile should be considered. The impact on total dose due to the transit dose correction near the dwell positions is significant and should be considered more carefully in treatments with high dose rate, several catheters, multiple dwell positions, small dwell times, and several fractions.

60 APPLIED LIFE SCIENCES↗

OBSERVATIONAL SIGNATURES OF CORONAL LOOP HEATING AND COOLING DRIVEN BY FOOTPOINT SHUFFLING

The evolution of a coronal loop is studied by means of numerical simulations of the fully compressible three-dimensional magnetohydrodynamic equations using the HYPERION code. The footpoints of the loop magnetic field are advected by random motions. As a consequence, the magnetic field in the loop is energized and develops turbulent nonlinear dynamics characterized by the continuous formation and dissipation of field-aligned current sheets: energy is deposited at small scales where heating occurs. Dissipation is nonuniformly distributed so that only a fraction of the coronal mass and volume gets heated at any time. Temperature and density are highly structured at scales that, in the solar corona, remain observationally unresolved: the plasma of our simulated loop is multithermal, where highly dynamical hotter and cooler plasma strands are scattered throughout the loop at sub-observational scales. Numerical simulations of coronal loops of 50,000 km length and axial magnetic field intensities ranging from 0.01 to 0.04 T are presented. To connect these simulations to observations, we use the computed number densities and temperatures to synthesize the intensities expected in emission lines typically observed with the Extreme Ultraviolet Imaging Spectrometer on Hinode. These intensities are used to compute differential emission measure distributions using the Monte Carlo Markov Chain code, which are very similar to those derived from observations of solar active regions. We conclude that coronal heating is found to be strongly intermittent in space and time, with only small portions of the coronal loop being heated: in fact, at any given time, most of the corona is cooling down.

79 ASTRONOMY AND ASTROPHYSICS↗

Advancement of Functional Expansion Tallies Capabilities in Serpent 2

The need for accident tolerant fuels (ATFs) became self-evident after the series of unfortunate nuclear events at Fukushima Daiichi in 2011. Restarting the Transient Reactor Test Facility (TREAT) at INL is one component in this concerted effort to establish the ATF development and testing program under the United States Department of Energy. Prior to the closure of TREAT in 1994, each experiment required multiple low-power calibration tests before and after each run. Since then, computing advancements have set computer-based modeling and simulation (M and S) as the preferred calibration method. The anticipated savings are both temporal and monetary. MAMMOTH, a collection of M and S tools built from the finite element-based MOOSE Framework, has been established to support the designed experiments intended for TREAT. This collection of tools includes: 1) Rattlesnake, a striking diffusion-based neutronics code, 2) BISON, an advanced materials and fuels performance code, and 3) RELAP-7, a powerful thermal-hydraulics solver. Each is capable of independently performing its function accurately and effectively; unified, they create a powerful multiphysics reactor core M and S application. One unfortunate weakness in this configuration lies in the unique geometry and design of TREAT's core. Although simple to model, this poses a challenge to the diffusion solver in Rattlesnake. One alternative is to replace Rattlesnake with Serpent, a fully-fledged Monte Carlo (MC) neutronics code. Coincidentally, Rattlesnake calculations are frequently validated against Serpent simulations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Examining Computational Assumptions For Godiva IV

Over the course of summer 2016, the effects of several computational modeling assumptions with respect to the Godiva IV reactor were examined. The majority of these assumptions pertained to modeling errors existing in the control rods and burst rod. The Monte Carlo neutron transport code, MCNP, was used to investigate these modeling changes, primarily by comparing them to that of the original input deck specifications.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Synchrotron polarization in blazars

We present a detailed analysis of time- and energy-dependent synchrotron polarization signatures in a shock-in-jet model for γ-ray blazars. Our calculations employ a full three-dimensional radiation transfer code, assuming a helical magnetic field throughout the jet. The code considers synchrotron emission from an ordered magnetic field, and takes into account all light-travel-time and other relevant geometric effects, while the relevant synchrotron self-Compton and external Compton effects are handled with the two-dimensional Monte-Carlo/Fokker-Planck (MCFP) code. We consider several possible mechanisms through which a relativistic shock propagating through the jet may affect the jet plasma to produce a synchrotron and high-energy flare. Most plausibly, the shock is expected to lead to a compression of the magnetic field, increasing the toroidal field component and thereby changing the direction of the magnetic field in the region affected by the shock. We find that such a scenario leads to correlated synchrotron + synchrotron-self-Compton flaring, associated with substantial variability in the synchrotron polarization percentage and position angle. Most importantly, this scenario naturally explains large polarization angle rotations by ≳ 180°, as observed in connection with γ-ray flares in several blazars, without the need for bent or helical jet trajectories or other nonaxisymmetric jet features.

79 ASTRONOMY AND ASTROPHYSICS↗

Investigation of mode activity in NBI-heated experiments of Wendelstein 7-X

The 2018 operation phase (OP 1.2b) of the stellarator Wendelstein 7-X (W7-X) included, for the first time, neutral beam injection (NBI) to heat the plasma. Since the injection geometry at W7-X is not parallel, this generates both passing and trapped fast particles. During longer phases of NBI injection, with the primary purpose to study the heating efficiency of this system, Alfvén eigenmodes (AEs) were observed by a number of diagnostics, including the phase contrast imaging (PCI) system, the magnetic pick-up coils (Mirnov coils), and the soft x-ray multi-camera tomography system (XMCTS). Alfvén eigenmodes are of great interest for future fusion reactors as it has been shown that the resonant interaction of fast ions with self-excited AEs can lead to enhanced transport of fast ions and potentially to energy losses. This is especially true for so-called gap-modes, Alfvén eigenmodes with frequencies in gaps of the continuous spectrum, since they lack continuum damping. These modes are commonly known to be excited by fast ions, but other destabilizing mechanisms, e.g. the electron-pressure gradient are also possible. In this article we present a first analysis of the experimentally observed frequencies from the theoretical side. The calculation of shear Alfvén wave continua for selected cases and the assignment of observed frequencies to the gaps of the continuous spectra are presented. Using the ideal-MHD code CKA (Könies A. 2007 10th IAEA TM on Energetic Particles in Magnetic Confinement System), we find gap modes that match the experimental measurements in terms of the observed frequencies. We emphasize the crucial roles played by the coupling of sound and Alfvén waves as well as of the Doppler shift arising as a consequence of the radial electric field in W7-X. We employ the perturbative gyrokinetic code CKA-EUTERPE (Feh´er 2013 Simulation of the interaction between Alfv´en waves and fast particles), using a slowing-down distribution function for the fast ions as calculated by the Monte-Carlo particle following code ASCOT (Hirvijoki et al 2014 Comput. Phys. Commun. 185 1310–21) to assess the fast-ion drive. Here we find that the fast-ion drive is insufficient to overcome the background-plasma damping. The fact that unstable modes were observed experimentally may point to problems with the modelling or indicate the existence of other destabilizing mechanisms, e.g. associated with the electron-pressure gradient (Windisch et al 2017 Plasma Phys. Control. Fusion 59 105002) that sensitively depend on the profiles of the background plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

MC21/CTF and VERA-CS Mutiphysics Solutions to VERA Core Physics Benchmark Problems 6 and 7

The continuous energy Monte Carlo neutron transport code, MC21, was coupled to the CTF subchannel thermal-hydraulics code using a combination of CASL tools and in-house Python scripts. An MC21/CTF solution for VERA Core Physics Benchmark Progression Problem 6 demonstrated excellent agreement with MC21/COBRA-IE and VERA solutions. The MC21/CTF solution for VERA Core Physics Benchmark Progression Problem 7, Watts Bar Unit 1 at beginning of cycle (BOC) hot full power (HFP) equilibrium xenon conditions, is the first published coupled Monte Carlo neutronics / subchannel T-H solution for this problem. MC21/CTF predicted a critical boron concentration of 854.5 ppm, yielding a critical eigenvalue of 0.99994 ± 6.8E-6 (95% CI). Excellent agreement with a VERA solution of Problem 7 was also demonstrated for integral and local power and temperature parameters.

Aviles, Brian↗

Reactivity worth measurements on fast burst reactor Caliban - description and interpretation of integral experiments for the validation of nuclear data

Reactivity perturbation experiments using various materials are being performed on the HEU fast core CALIBAN, an experimental device operated by the CEA VALDUC Criticality and Neutron Transport Research Laboratory. These experiments provide valuable information to contribute to the validation of nuclear data for the materials used in such measurements. This paper presents the results obtained in a first series of measurements performed with Au-197 samples. Experiments which have been conducted in order to improve the characterization of the core are also described and discussed. The experimental results have been compared to numerical calculation using both deterministic and Monte Carlo neutron transport codes with a simplified model of the reactor. This early work led to a methodology which will be applied to the future experiments which will concern other materials of interest. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Evaluating importance maps for TRIPOLI-4{sup R} using deterministic or on-line methods - 25446

Variance reduction is a key ingredient for solving radiation-protection problems with Monte-Carlo particle-transport codes. Many variance-reduction methods require the definition of an importance map and exhibit optimal performance if the importance map is given by the adjoint flux, the solution of the adjoint Boltzmann equation. This paper presents the implementation of the Consistent Adjoint-Driven Importance Sampling (CADIS) methodology via a coupling between the TRIPOLI-4 R Monte-Carlo particle-transport code and the IDT deterministic solver, Additionally, we describe the implementation of a new TRIPOLI-4 R score that makes it possible to estimate the adjoint flux during a direct Monte-Carlo calculation. These new features are expected to simplify the solution of difficult shielding problems. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗

GPU Acceleration of Kernel Density Estimators in Monte Carlo Neutron Transport Simulations

Kernel Density Estimators (KDEs) have recently been explored for use in Monte Carlo radiation transport simulations as an alternative to histogram tallies for capturing spatially-resolved quantities such as neutron flux and reaction rates. Histogram tallies suffer from large uncertainties when detailed spatial resolution of the quantity of interest is required; increasing the resolution of the underlying mesh increases the uncertainty in each histogram bin. KDEs obtain estimates of underlying densities on user-defined points, with the uncertainty at those points being independent of the spatial resolution desired since collisions and particle-track segments are allowed to contribute to the score at multiple tally points. Thus, KDEs show potential for obtaining estimates of a smoother spatially-resolved with reduced variance when compared to a histogram. However, this ability that enables KDEs to obtain lower uncertainties per particle history is more computationally expensive than using a traditional histogram bins on a structured mesh. While KDEs are capable of producing lower tally uncertainties compared to histogram binning for a given number of histories, their increased computational cost may make them less efficient. Even so, it is possible to reduce this additional cost by exporting the KDE tally routines onto Graphics Processing Units (GPUs). GPUs are commonly found on high-performance computing machines as a means of increasing performance without adding additional compute nodes. While it is difficult for existing Monte Carlo radiation transport codes to take advantage of GPUs without rewriting large portions of the code, it is possible to leverage GPUs through heterogeneous computing. This work uses a modified version of OpenMC and exports the KDE tally routines onto a GPU to obtain problem-dependent speedups ranging between 1.6 and 5.0 for the problem studied in this paper. While this methodology is applied to KDEs, the same ideas can be applied to other computationally expensive vectorizable portions of Monte Carlo algorithms. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Testing of Error-Correcting Sparse Permutation Channel Codes

A computer program performs Monte Carlo direct numerical simulations for testing sparse permutation channel codes, which offer strong error-correction capabilities at high code rates and are considered especially suitable for storage of digital data in holographic and volume memories. A word in a code of this type is characterized by, among other things, a sparseness parameter (M) and a fixed number (K) of 1 or "on" bits in a channel block length of N.

Shcheglov, Kirill, V.↗

Shielding Design of Gantry 3 - Paper 104

The proton therapy complex at PSI will be extended with a third irradiation facility, Gantry 3. This new Gantry is being built in the south-west area of the experimental hall of the High Intensity Proton Accelerator (HIPA) complex at PSI. Originally, it was planned to restrict the operation of Gantry 3 to 180 deg. by hardware means, but due to new considerations, the Gantry will be able to be moved in the full 360 deg. field. This change requires the redesign and reinforcement of the too weak shielding in the north and east directions. Additionally, the isocenter has been fixed, thus allowing for the exact calculation of distances, angles and effective shielding thicknesses in order to determine the dose rates more accurately. For these calculations a fast and simple semi-empirical analytical formula and the more time consuming, but also more sophisticated best estimate code MCNPX 2.7.0 have been used. The two methods used for the calculations of the dose rates are in good agreement. Generally, the analytical method provides the results instantaneously, but predicts higher dose rate values of a factor 2-3 larger than the Monte-Carlo particle transport code MCNPX. Thus, MCNPX should be used to verify the final design of the shielding, including additional details in the simulation model (e.g. ventilation shafts) which have weakening effect on the shielding. Such influences are hard to quantify analytically. The analytical method fails when the geometry becomes complex, i.e. when it is difficult to determine the effective length of the shielding; in the case of position 10 the shielding effect of the inner wall is difficult to quantify. Due to the change of the specification of the Gantry, the shielding was reinforced at the north and east sides, and the HC-layer was extended in front of the control room. The dose rates are acceptable with the new shielding layout. A detailed dose rate mapping is planned in the commissioning phase of Gantry 3 to verify the simulation results and to identify any remaining weak spots of the shielding. (authors)

43 PARTICLE ACCELERATORS↗

Simulation of multi-photon emission isotopes using time-resolved SimSET multiple photon history generator

Abstract-Multiple-photon emitters, such as In-111 or Se-75, have enormous potential in the field of nuclear medicine imaging. For example, Se-75 can be used to investigate the bile acid malabsorption and measure the bile acid pool loss. The simulation system for emission tomography (SimSET) is a well-known Monte Carlo simulation (MCS) code in nuclear medicine for its high computational efficiency. However, current SimSET cannot simulate these isotopes due to the lack of modeling of complex decay scheme and the time-dependent decay process. To extend the versatility of SimSET for simulation of those multi-photon emission isotopes, a time-resolved multiple photon history generator based on SimSET codes is developed in present study. For developing the time-resolved SimSET (trSimSET) with radionuclide decay process, the new MCS model introduce new features, including decay time information and photon time-of-flight information, into this new code. The half-life of energy states were tabulated from the Evaluated Nuclear Structure Data File (ENSDF) database. The MCS results indicate that the overall percent difference is less than 8.5% for all simulation trials as compared to GATE. To sum up, we demonstrated that time-resolved SimSET multiple photon history generator can have comparable accuracy with GATE and keeping better computational efficiency. The new MCS code is very useful to study the multi-photon imaging of novel isotopes that needs the simulation of lifetime and the time-of-fight measurements. (authors)

07 ISOTOPE AND RADIATION SOURCES↗