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,117 records · Page 62

The Effect of the Flux Separability Approximation on Multigroup Neutron Transport

The angular dependence of flux-weighted multigroup cross sections is commonly neglected when generating multigroup libraries. The error of this flux separability approximation is typically not isolated from other error sources due to a lack of availability of library generation and corresponding solvers that cannot relax this approximation. These errors can now be isolated and quantified with the availability of a multigroup Monte Carlo transport and multigroup library-generation capability in the OpenMC Monte Carlo transport code. This work will discuss relevant details of the OpenMC implementation, provide an example case useful for detailing the type of errors one can expect from making the flux separability approximation, and end with more realistic problems which show the impact of the approximation and highlight how it can strongly arise from an energy-dependent resonance absorption effect. Since the angle-dependence is intrinsically linked to the energy group structure, these examples also show that relaxing the flux separability approximation with angle-dependent cross sections could be used to reduce either the fine-tuning required to set a multigroup energy structure for a specific reactor type or the number of energy groups required to obtain a desired level of accuracy for a given problem. This trade-off could increase the costs of generating multigroup cross sections, and has the potential to require more memory for storing the multigroup library during the transport calculations, but it can significantly reduce the computational time required since the runtime of a discrete ordinates or method of characteristics neutron transport solver scales roughly linearly with the number of groups.

42 ENGINEERING↗

Improved mesh based photon sampling techniques for neutron activation analysis

The design of fusion power systems requires analysis of neutron activation of large, complex volumes, and the resulting particles emitted from these volumes. Structured mesh-based discretization of these problems allows for improved modeling in these activation analysis problems. Finer discretization of these problems results in large computational costs, which drives the investigation of more efficient methods. Within an ad hoc subroutine of the Monte Carlo transport code MCNP, we implement sampling of voxels and photon energies for volumetric sources using the alias method. The alias method enables efficient sampling of a discrete probability distribution, and operates in 0(1) time, whereas the simpler direct discrete method requires 0(log(n)) time. By using the alias method, voxel sampling becomes a viable alternative to sampling space with the 0(1) approach of uniformly sampling the problem volume. Additionally, with voxel sampling it is straightforward to introduce biasing of volumetric sources, and we implement this biasing of voxels as an additional variance reduction technique that can be applied. We verify our implementation and compare the alias method, with and without biasing, to direct discrete sampling of voxels, and to uniform sampling. We study the behavior of source biasing in a second set of tests and find trends between improvements and source shape, material, and material density. Overall, however, the magnitude of improvements from source biasing appears to be limited. Future work will benefit from the implementation of efficient voxel sampling - particularly with conformal unstructured meshes where the uniform sampling approach cannot be applied. (authors)

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Proton Lateral Broadening Distribution Comparisons Between GRNTRN, MCNPX, and Laboratory Beam Measurements

Recent developments in NASA s deterministic High charge (Z) and Energy TRaNsport (HZETRN) code have included lateral broadening of primary ion beams due to small-angle multiple Coulomb scattering, and coupling of the ion-nuclear scattering interactions with energy loss and straggling. This new version of HZETRN is based on Green function methods, called GRNTRN, and is suitable for modeling transport with both space environment and laboratory boundary conditions. Multiple scattering processes are a necessary extension to GRNTRN in order to accurately model ion beam experiments, to simulate the physical and biological-effective radiation dose, and to develop new methods and strategies for light ion radiation therapy. In this paper we compare GRNTRN simulations of proton lateral broadening distributions with beam measurements taken at Loma Linda University Proton Therapy Facility. The simulated and measured lateral broadening distributions are compared for a 250 MeV proton beam on aluminum, polyethylene, polystyrene, bone substitute, iron, and lead target materials. The GRNTRN results are also compared to simulations from the Monte Carlo MCNPX code for the same projectile-target combinations described above.

Mertens, Christopher J.↗

Inspection of Dry Storage Casks Using PGAA and NAA Techniques

In this effort, an assessment of a chloride detection method is performed for inspecting the surfaces of dry cask storage system (DCSS) canisters, with focus on the implementation of a non-destructive nuclear analysis and methodology. Some DCSS canisters (especially those located in coastal environments) will be exposed to marine environments, which can cause atmospheric stress corrosion cracking (SCC). Chloride deposits on canister surfaces provide the corrosive environment necessary for atmospheric SCC with the residual stresses and materials properties conditions. The most susceptible regions on the canister surface are the weld areas, or where cold work has been performed, and may provide the necessary driving force for crack growth. The heat-affected zones (HAZs) near welds are also potentially susceptible to atmospheric SCC due to sensitization of the material caused by chromium depletion during the welding process. Information collected from the field and from laboratory studies has not been able to rule out the possibility of atmospheric chloride induced stress corrosion cracking (CISCC), although no occurrences of atmospheric SCC in DCSS canisters have been detected to-date. A number of researchers have demonstrated that localized corrosion can take place under such situations on materials commonly used for interim storage containers, including 304SS (Cook et al. 2010; Shirai et al. 2011) and 316SS (Tani et al. 2009). In this summary, the preliminary results of a neutron activation analysis (NAA) and prompt gamma activation analysis (PGAA) are performed using MCNP6, a highly sophisticated, multi-particle Monte Carlo simulation code (with the employment of the high performance computing system at the Idaho National Laboratory), and confirm the feasibility of using PGAA and NAA as an environmental inspection methodology for detection and monitoring of small quantities of chlorine on canister surfaces. The specific signature of gamma rays due reactions of sea salt constituents are documented and summarized in this paper. The environmental constituent for passive oxide film breaching is chlorine, which can cause pitting initiation and incubation time that may lead to cracking in stainless steel. Therefore, an assessment of chloride presence at a canister surface is essential before moving to other nondestructive analysis (NDA) methods for detection of pitting and cracking. Two nonlinear acoustic techniques will also be used in our integrated research project (IRP) research study that proves to be useful in the detection of SCC: nonlinear resonant ultrasound spectroscopy (NRUS), and time revisited elastic nonlinear diagnostic (TREND). These two techniques are complementary to marine atmospheric detection techniques. In this presentation only PGAA and NAA are described in detail. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Poster - 40: Treatment Verification of a 3D-printed Eye Phantom for Proton Therapy

Purpose: Ocular melanoma is a form of eye cancer which is often treated using proton therapy. The benefit of the steep proton dose gradient can only be leveraged for accurate patient eye alignment. A treatment-planning program was written to plan on a 3D-printed anatomical eye-phantom, which was then irradiated to demonstrate the feasibility of verifying in vivo dosimetry for proton therapy using PET imaging. Methods: A 3D CAD eye model with critical organs was designed and voxelized into the Monte-Carlo transport code FLUKA. Proton dose and PET isotope production were simulated for a treatment plan of a test tumour, generated by a 2D treatment-planning program developed using NumPy and proton range tables. Next, a plastic eye-phantom was 3D-printed from the CAD model, irradiated at the TRIUMF Proton Therapy facility, and imaged using a PET scanner. Results: The treatment-planning program prediction of the range setting and modulator wheel was verified in FLUKA to treat the tumour with at least 90% dose coverage for both tissue and plastic. An axial isotope distribution of the PET isotopes was simulated in FLUKA and converted to PET scan counts. Meanwhile, the 3D-printed eye-phantom successfully yielded a PET signal. Conclusions: The 2D treatment-planning program can predict required parameters to sufficiently treat an eye tumour, which was experimentally verified using commercial 3D-printing hardware to manufacture eye-phantoms. Comparison between the simulated and measured PET isotope distribution could provide a more realistic test of eye alignment, and a variation of the method using radiographic film is being developed.

60 APPLIED LIFE SCIENCES↗

Use of freeze-casting in advanced burner reactor fuel design

This paper will detail the modeling of a fast reactor with fuel pins created using a freeze-casting process. Freeze-casting is a method of creating an inert scaffold within a fuel pin. The scaffold is created using a directional solidification process and results in open porosity for emplacement of fuel, with pores ranging in size from 300 microns to 500 microns in diameter. These pores allow multiple fuel types and enrichments to be loaded into one fuel pin. Also, each pore could be filled with varying amounts of fuel to allow for the specific volume of fission gases created by that fuel type. Currently fast reactors, including advanced burner reactors (ABR's), are not economically feasible due to the high cost of operating the reactors and of reprocessing the fuel. However, if the fuel could be very precisely placed, such as within a freeze-cast scaffold, this could increase fuel performance and result in a valid design with a much lower cost per megawatt. In addition to competitive costs, freeze-cast fuel would also allow for selective breeding or burning of actinides within specific locations in fast reactors. For example, fast flux peak locations could be utilized on a minute scale to target specific actinides for transmutation. Freeze-cast fuel is extremely flexible and has great potential in a variety of applications. This paper performs initial modeling of freeze-cast fuel, with the generic fast reactor parameters for this model based on EBR-II. The core has an assumed power of 62.5 MWt. The neutronics code used was Monte Carlo N-Particle (MCNP5) transport code. Uniform pore sizes were used in increments of 100 microns. Two different freeze-cast scaffold materials were used: ceramic (MgO-ZrO{sub 2}) and steel (SS316L). Separate models were needed for each material because the freeze-cast ceramic and metal scaffolds have different structural characteristics and overall porosities. Basic criticality results were compiled for the various models. Preliminary results show that criticality is achievable with freeze-cast fuel pins despite the significant amount of inert fuel matrix. Freeze casting is a promising method to achieve very precise fuel placement within fuel pins. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Calculations versus measurements for residual dose rates from SNS spent structures

Here, residual dose rate measurements were conducted on target vessel #13 and proton beam window #5 after extraction from their service locations. These measurements are used to verify calculation methods of radionuclide inventory assessment that are typically performed for nuclear waste characterization and transportation of these structures. Neutronics analyses for predicting residual dose rates are carried out using the transport code MCNPX and the transmutation code CINDER90. For transport analyses a complex and rigorous geometry model of the structures and their surroundings are applied. The neutronics analyses are carried out using the Bertini and CEM high energy physics models for simulating particles interactions above the table-based cross section range. Finally, obtained calculational results are analyzed and compared to the measured dose rates and overall show good agreement within 25%, which shows applicability of the methods used in analyses.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

ANALYSIS OF THE MSL/MEDLI ENTRY DATA WITH COUPLED CFD AND MATERIAL RESPONSE.

The Mars Science Laboratory (MSL) was protected during its atmospheric entry by an instrumented heat-shield using NASA's Phenolic Impregnated Carbon Ablator (PICA) material. PICA is a lightweight carbon fiber/polymeric resin material that offers out-standing performances for protecting probes during planetary entry. The Mars Entry Descent and Landing Instrument (MEDLI) suite on MSL offers unique in-flight validation data for models of material response and atmospheric entry. MEDLI recorded, among other things, time-resolved in-depth temperature data of PICA using thermocouple sensors assembled in the MEDLI Integrated Sensor Plugs (MISP). The objective of this work is to showcase and analyze the coupling between the material response and the aerothermal environment. As shown in Figure 1, the workflow is divided into the following steps. First, the aerothermal properties are computed in the Data Parallel Line Relaxation (DPLR) code [3] and used with the Nonequilibrium air radiation (NEQAIR) program [8] to compute radiative heating. Second, the thermal response inside the material is computed in the Porous material Analysis Toolbox based on Open-FOAM (PATO) using a fixed blowing correction parameter. Third, the pyrolysis gases computed in PATO are used as inputs to a blowing boundary condition within DPLR. Fourth, the new environment properties from DPLR are used in NEQAIR to provide an updated solution, then both the updated aerothermal environment and radiative heating are used in PATO without blowing correction. The third and fourth steps are then repeated until convergence in surface temperature is obtained. Convergence in the radiative heating is generally achieved before surface temperature, at which point the radiative heating is no longer updated. Char mass loss rates are forced to zero to produce a non-receding surface condition. For early time points in the trajectory, where flow around the MSL aeroshell is rarefied, the Direct Simulation Monte Carlo (DSMC) code, SPARTA, is used to compute the aerothermal environment. Iteration between PATO and SPARTA is not performed due to the computational cost of DSMC simulations. Preliminary results of the coupling between PATO and DPLR for the MSL heatshield atmospheric entry model are presented in Figures 2-4 at 65 seconds after entry interface. Figure 2 shows the surface temperature results from an uncoupled simulation in PATO with the blowing correction parameter applied (left) along with the coupled surface temperature after iteration (right). Figure 3 shows the surface temperature along the centerline from windward to leeward for easier comparison. Figure 4 shows the coupled and uncoupled pyrolysis gas blowing rate. Mars 2020 used a similar heatshield consisting of PICA for thermal protection during entry, descent, and landing. In preparation for Mars 2020 post-flight analysis, the predictive material response capability is benchmarked against flight data from MEDLI. This work represents an important milestone toward the development of validated predictive capabilities for designing thermal protection systems for planetary probes.

Mars Science Laboratory↗

Analysis of MSL/MEDLI Entry Data with Coupled CFD and Material Response

The Mars Science Laboratory (MSL) was protected during its atmospheric entry by an instrumented heatshield using NASA's Phenolic Impregnated Carbon Ablator (PICA) material [1]. PICA is a lightweight carbon fiber/polymeric resin material that offers outstanding performances for protecting probes during planetary entry. The Mars Entry Descent and Landing Instrument (MEDLI) suite on MSL offers unique in-flight validation data for models of material response and atmospheric entry. MEDLI recorded, among other things, time-resolved in-depth temperature data of PICA using thermocouple sensors assembled in the MEDLI Integrated Sensor Plugs (MISP) [2]. The objective of this work is to showcase and analyze the coupling between the material response and the aerothermal environment. As shown in Figure 1, the workflow is divided into the following steps. First, the aerothermal properties are computed in the Data Parallel Line Relaxation (DPLR) code [3] and used with the Nonequilibrium air radiation (NEQAIR) program [8] to compute radiative heating. Second, the thermal response inside the material is computed in the Porous material Analysis Toolbox based on OpenFOAM (PATO) [4,5,6] using a fixed blowing correction parameter. Third, the pyrolysis gases computed in PATO are used as inputs to a blowing boundary condition within DPLR. Fourth, the new environment properties from DPLR are used in NEQAIR to provide an updated solution, then both the updated aerothermal environment and radiative heating are used in PATO without blowing correction. The third and fourth steps are then repeated until convergence in surface temperature is obtained. Convergence in the radiative heating is generally achieved before surface temperature, at which point the radiative heating is no longer updated. Char mass loss rates are forced to zero to produce a non-receding surface condition. For early time points in the trajectory, where flow around the MSL aeroshell is rarefied, the Direct Simulation Monte Carlo (DSMC) code, SPARTA [7], is used to compute the aerothermal environment. Iteration between PATO and SPARTA is not performed due to the computational cost of DSMC simulations. Preliminary results of the coupling between PATO and DPLR for the MSL heatshield atmospheric entry model are presented in Figures 2-4 at 65 seconds after entry interface. Figure 2 shows the surface temperature results from an uncoupled simulation in PATO with the blowing correction parameter applied (left) along with the coupled surface temperature after iteration (right). Figure 3 shows the surface temperature along the centerline from windward to leeward for easier comparison. Figure 4 shows the coupled and uncoupled pyrolysis gas blowing rate. Mars 2020 used a similar heatshield consisting of PICA for thermal protection during entry, descent, and landing. In preparation for Mars 2020 post-flight analysis, the predictive material response capability is benchmarked against flight data from MEDLI. This work represents an important milestone toward the development of validated predictive capabilities for designing thermal protection systems for planetary probes.

Thermal Protection Systems↗

Effect of the size of nanoparticles on their dissolution within metal-glass nanocomposites under sustained irradiation

We investigate the dissolution law of metallic nanoparticles (NPs) under sustained irradiation. The system is composed of isolated spherical gold NPs (4–100 nm) embedded in an amorphous silica host matrix. Samples are irradiated at room temperature in the nuclear stopping power regime with 4 MeV Au ions for fluences up to 8 × 10{sup 16 }cm{sup −2}. Experimentally, the dependence of the dissolution kinetics on the irradiation fluence is linear for large NPs (45–100 nm) and exponential for small NPs (4–25 nm). A lattice-based kinetic Monte Carlo (KMC) code, which includes atomic diffusion and ballistic displacement events, is used to simulate the dynamical competition between irradiation effects and thermal healing. The KMC simulations allow for a qualitative description of the NP dissolution in two main stages, in good agreement with the experiment. Moreover, the perfect correlation obtained between the evolution of the simulated flux of ejected atoms and the dissolution rate in two stages implies that there exists an effect of the size of NPs on their dissolution and a critical size for the transition between the two stages. The Frost-Russell model providing an analytical solution for the dissolution rate, accounts well for the first dissolution stage but fails in reproducing the data for the second stage. An improved model obtained by including a size-dependent recoil generation rate permits fully describing the dissolution for any NP size. This proves, in particular, that the size effect on the generation rate is the principal reason for the existence of two regimes. Finally, our results also demonstrate that it is justified to use a unidirectional approximation to describe the dissolution of the NP under irradiation, because the solute concentration is particularly low in metal-glass nanocomposites.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

SU-G-JeP3-06: Lower KV Image Dose Are Expected From a Limited-Angle Intra-Fractional Verification (LIVE) System for SBRT Treatments

Purpose: In order to track the tumor movement for patient positioning verification during arc treatment delivery or in between 3D/IMRT beams for stereotactic body radiation therapy (SBRT), the limited-angle kV projections acquisition simultaneously during arc treatment delivery or in-between static treatment beams as the gantry moves to the next beam angle was proposed. The purpose of this study is to estimate additional imaging dose resulting from multiple tomosynthesis acquisitions in-between static treatment beams and to compare with that of a conventional kV-CBCT acquisition. Methods: kV imaging system integrated into Varian TrueBeam accelerators was modeled using EGSnrc Monte Carlo user code, BEAMnrc and DOSXYZnrc code was used in dose calculations. The simulated realistic kV beams from the Varian TrueBeam OBI 1.5 system were used to calculate dose to patient based on CT images. Organ doses were analyzed using DVHs. The imaging dose to patient resulting from realistic multiple tomosynthesis acquisitions with each 25–30 degree kV source rotation between 6 treatment beam gantry angles was studied. Results: For a typical lung SBRT treatment delivery much lower (20–50%) kV imaging doses from the sum of realistic six tomosynthesis acquisitions with each 25–30 degree x-ray source rotation between six treatment beam gantry angles were observed compared to that from a single CBCT image acquisition. Conclusion: This work indicates that the kV imaging in this proposed Limited-angle Intra-fractional Verification (LIVE) System for SBRT Treatments has a negligible imaging dose increase. It is worth to note that the MV imaging dose caused by MV projection acquisition in-between static beams in LIVE can be minimized by restricting the imaging to the target region and reducing the number of projections acquired. For arc treatments, MV imaging acquisition in LIVE does not add additional imaging dose as the MV images are acquired from treatment beams directly during the treatment.

60 APPLIED LIFE SCIENCES↗

TU-F-BRE-03: Application of a Novel Mass-Density Compensation Optimization Method to Improve the Response of a Liquid-Filled Ionization Chamber in Nonstandard Fields

Purpose: To improve the response of the microLion detector (PTW 31018) in small field conditions by optimizing the density of the detector's non-sensitive components using a novel method based on the detector's dose response function. Methods: The central values h(0,0) of the perturbation functions for the microLion detector and a volume of water equivalent to its sensitive volume were calculated using the Monte Carlo user code egs-chamber by scoring the dose absorbed by a 6 MV photon pencil beam incident on their centroids. Values of h(0,0) were plotted as a function of the density of the microLion's graphite electrode with the detector placed in the axial orientation. The optimized density was found by finding the minimal value of h(0,0). Results: A density of 1.37 g/cm3 was found to minimize the perturbation function of the microLion detector. The modified microLion's response was then evaluated in small square fields with sides in the range of 5 – 40 mm and found to be consistent with highly watere-quivalent detectors such as a scintillating detector (Exradin W1) and a generic alanine detector in both axial and radial orientations. Conclusion: This work illustrates a novel method which can used to optimize the design of radiation detectors in small fields. This method should also work with other optimization parameters (e.g. thickness of electrode). Density-compensated detectors have the potential to eliminate the need to evaluate nonstandard field correction factors as described by the IAEA-AAPM formalism (Alfonso et al.) and simplify future dosimetry protocols for SRS/SBRT modalities. Finally, we also expect an improvement in the response of density-compensated detectors for composite IMRT fields.

07 ISOTOPE AND RADIATION SOURCES↗

A Mercury Model of the Molly-G Fast Burst Reactor (FBR)

In order to design neutron irradiation experiments at a given nuclear reactor, one requires a representative neutron spectrum for the configuration of the reactor that will exist during the experiment. This usually implies modeling the operation of the reactor using a neutron transport code. Since the modeling of pulsed neutron reactors is quite complex, requiring coupled modeling of neutronics and thermo-structural response, the typical methodology employs analysis of the neutron spectrum which is obtained from a static keff calculation. While this time-independent Eigenfunction is not truly representative of that during a reactor pulse, we postulate that it is a reasonable representation for our purposes. The Nuclear Survivability (NS) Program at LLNL is considering fielding such experiments at the Molly-G (‘Molybdenum Godiva-II’) Fast Burst Reactor (FBR) at the White Sands Missile Range (WSMR). The goal of these experiments is to (a) test / proof / calibrate diagnostics, and (b) test the efficacy of neutron shield configurations, before fielding at other nuclear reactors, such as the Annular Core Research Reactor (ACRR) at SNL/NM. After attempting to secure either drawings or a copy of an existing MCNP model of the Molly-G FBR from various institutions, it was decided to develop a new model for use with the Mercury Monte Carlo transport code. The model described herein will be used to design these experiments.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Exploring {sup 57}Co as a new isotope for brachytherapy applications

Purpose: The characteristics of the radionuclide {sup 57}Co make it interesting for use as a brachytherapy source. {sup 57}Co combines a possible high specific activity with the emission of relatively low-energy photons and a half-life (272 days) suitable for regular source exchanges in an afterloader. {sup 57}Co decays by electron capture to the stable {sup 57}Fe with emission of 136 and 122 keV photons. Methods: A hypothetical {sup 57}Co source based on the Flexisource brachytherapy encapsulation with the active core set as a pure cobalt cylinder (length 3.5 mm and diameter 0.6 mm) covered with a cylindrical stainless-steel capsule (length 5 mm and thickness 0.125 mm) was simulated using Geant4 Monte Carlo (MC) code version 9.4. The radial dose function, g(r), and anisotropy function F(r,{theta}), for the line source approximation were calculated following the TG-43U1 formalism. The results were compared to well-known {sup 192}Ir and {sup 125}I radionuclides, representing the higher and the lower energy end of brachytherapy, respectively. Results: The mean energy of photons in water, after passing through the core and the encapsulation material was 123 keV. This hypothetical {sup 57}Co source has an increasing g(r) due to multiple scatter of low-energy photons, which results in a more uniform dose distribution than {sup 192}Ir. Conclusions: {sup 57}Co has many advantages compared to {sup 192}Ir due to its low-energy gamma emissions without any electron contamination. {sup 57}Co has an increasing g(r) that results in a more uniform dose distribution than {sup 192}Ir due to its multiple scattered photons. The anisotropy of the {sup 57}Co source is comparable to that of {sup 192}Ir. Furthermore, {sup 57}Co has lower shielding requirements than {sup 192}Ir.

07 ISOTOPE AND RADIATION SOURCES↗

Experiment Design of Secondary Neutron and Charged Particles Measurement with Stopping Targets Bombarded by 100- and 230-MeV/amu Helium Ions - Paper 124

Radiotherapy using ion beams has gained momentum in radiation oncology. While the number of proton therapy facilities that have been established or are under construction is rapidly growing, interest in heavy ion (whose atomic number Z is greater than one) radiotherapy has also been revitalized in the U.S. Helium ions, with mass and charge between protons and carbon ions, have advantages of deep penetration distances with significant peak-to-entrance dose ratios and limited dose contributions from secondary neutrons and other secondary charged fragments. In consideration of these advantages, the U.S. ion beam therapy community has expressed interests in helium ions as a potential candidate for ion beam therapy, especially pediatric cases. As such, characterizing the secondary nuclear fragment field created by primary helium ion interactions along the beam delivery devices and in human tissue is important for dose assessment and patient protection. Up to now the experimental data with respect to helium ions over the therapeutic energy range, 100 - 230 MeV/amu, is scarce compared to protons and carbon ions. Kurosawa et al. measured stopping target yields and angular spectra of neutrons, protons, deuterons and tritons produced by various target materials bombarded by 100- and 180- MeV/amu helium ions; Heilbronn et al. and Cecil et al. also measured the stopping target neutron yields by 155-, 160- and 177.5-MeV/amu helium ion bombardments. We have proposed to measure angular distributions of secondary neutrons and light charged ions produced by 100- and 230-MeV/amu helium ions bombarding stopping water and iron targets. The experiment was conducted at the Heavy Ion Medical Accelerator in Chiba (HIMAC) of the National Institute of Radiological Sciences (NIRS) in Japan this June. During the literature survey for experimental setup, we found that the ratio of the ion beam range (R) to the target depth (d) varies from 0.5 - 0.98 in the stopping target measurements either for helium ions incident on various targets or carbon ions bombarding targets as tissue surrogates. The difference in the values of ratio R/d makes comparison of experimental results complicated and require accurate corrections. Hence, how the R/d ratio and cross sectional area of the targets influence the angular distributions of secondary particles is investigated by using a Monte Carlo method in this study. The influence of the depth and the cross sectional area of stopping water and iron targets on the yield and the spectra of secondary particles induced by 100- and 230-MeV/amu He ions was studied by a Monte Carlo simulation code, PHITS. The secondary charged particles and secondary neutrons who come from breakup of projectiles are highly forward-focused, while the secondary particles produced by target breakup are nearly isotropic. As the energy of {sup 4}He ions increases from 100 MeV/amu to 230 MeV/amu, the normalized intensities of charged particles decrease more gently with increasing angles, whereas the neutron angular distribution falls quickly at angles smaller than 30 deg. and starts to exceed the normalized profile of 100 MeV/amu He from 60 deg. When the target depth increases, the impact on attenuation of secondary particles is more significant for lighter target mass and higher-energy projectiles at forward angles. Also, with deeper targets, more interactions occur between the secondary particles and the target elements, which results in more target-like fragments at large and backward angles. In respect of the cross sectional area of the stopping targets, the forward angular distributions are similar to the system with smaller cross sectional area of the targets; however, no general rule was found for how the large-angle distributions are affected by the target cross sectional area. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗

Accelerated Coupled Monte Carlo-Thermal Hydraulic Calculations using a Hybrid GTF-Diffusion-based Prediction Block: First Results

Accurate predictions of spatial power and temperature distributions require the coupling of a neutron transport solver with a thermal-hydraulic (TH) feedback. Nowadays, Monte Carlo (MC) codes are widely coupled to TH solvers, typically via a Picard iteration (PI) method, due to the higher fidelity that such frameworks can produce. To speed up a PI, a prediction step can produce an improved initial guess for a source distribution and feed it to the MC code. Recent work [1, 2] investigated a prediction step that uses generalized transfer functions (GTFs) to predict the macroscopic cross sections’ variations following a perturbation in TH properties, such as coolant density. The previous method also relied on first order perturbation (FOP) theory to predict perturbed power profiles, rather than using an expensive MC iterate. The implemented FOP method relied on generating a fission matrix from which the forward and adjoint eigenmodes were extracted and later used to for power calculations. The generation of the fission matrix can introduce a significant computational overhead, therefore undermining the performance of the proposed hybrid technique when applied to high-dimensional problems, e.g., full core calculations. This work attempts to improve the GTF-FOP prediction step by replacing the FOP solver with a nodal diffusion solver, thus eliminating the need to calculate a fission matrix. The GTF-diffusion step was tested for various moderator density perturbations. In each case, the predicted power distribution showed good agreement with the reference case. The latter is attributed to the generally good prediction of most spatially distributed macroscopic cross sections, except the transport cross section, which will become the focus of future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Prediction of neutron production and energy spectrum by the inverse kinematic reaction between an incident 7 Li 3+ beam and a proton target in PHITS

A neutron source using the inverse kinematic reaction between lithium and proton, p( 7 Li, n) 7 Be, achieves forward-directed neutrons, potentially enhancing neutron yield in the forward direction. Despite the advantage, no evaluated-cross-section data for this reaction can be used in Monte Carlo simulation codes, such as PHITS. To solve this problem, this study aims to evaluate the applicability of the user-defined cross-section data, Frag data, for p ( 7 Li, n) 7 Be in PHITS. The simulations reproduced collisions between 7 Li 3+ ions and polypropylene targets. The Frag data was edited based on the JENDL-5 by utilizing the two-body collision kinematics. The neutron yield and angular distribution were investigated in the simulation. As a result, the forward neutron convergence with a reasonable neutron yield and energy spectrum was observed. The expected neutron yield in the forward 1-steradian area is 2.46 × 10 10 n/s when lithium-ion energy and current are 16.45 MeV and 0.1 mA.

43 PARTICLE ACCELERATORS↗