Search NASA⌕ Search

SEARCH · Search NASA

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,441 records · Page 80
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

SU-F-I-53: Coded Aperture Coherent Scatter Spectral Imaging of the Breast: A Monte Carlo Evaluation of Absorbed Dose

Purpose: Coherent scatter based imaging has shown improved contrast and molecular specificity over conventional digital mammography however the biological risks have not been quantified due to a lack of accurate information on absorbed dose. This study intends to characterize the dose distribution and average glandular dose from coded aperture coherent scatter spectral imaging of the breast. The dose deposited in the breast from this new diagnostic imaging modality has not yet been quantitatively evaluated. Here, various digitized anthropomorphic phantoms are tested in a Monte Carlo simulation to evaluate the absorbed dose distribution and average glandular dose using clinically feasible scan protocols. Methods: Geant4 Monte Carlo radiation transport simulation software is used to replicate the coded aperture coherent scatter spectral imaging system. Energy sensitive, photon counting detectors are used to characterize the x-ray beam spectra for various imaging protocols. This input spectra is cross-validated with the results from XSPECT, a commercially available application that yields x-ray tube specific spectra for the operating parameters employed. XSPECT is also used to determine the appropriate number of photons emitted per mAs of tube current at a given kVp tube potential. With the implementation of the XCAT digital anthropomorphic breast phantom library, a variety of breast sizes with differing anatomical structure are evaluated. Simulations were performed with and without compression of the breast for dose comparison. Results: Through the Monte Carlo evaluation of a diverse population of breast types imaged under real-world scan conditions, a clinically relevant average glandular dose for this new imaging modality is extrapolated. Conclusion: With access to the physical coherent scatter imaging system used in the simulation, the results of this Monte Carlo study may be used to directly influence the future development of the modality to keep breast dose to a minimum while still maintaining clinically viable image quality.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Minibeam radiation therapy for the management of osteosarcomas: A Monte Carlo study

Purpose: Minibeam radiation therapy (MBRT) exploits the well-established tissue-sparing effect provided by the combination of submillimetric field sizes and a spatial fractionation of the dose. The aim of this work is to evaluate the feasibility and potential therapeutic gain of MBRT, in comparison with conventional radiotherapy, for osteosarcoma treatments. Methods: Monte Carlo simulations (PENELOPE/PENEASY code) were used as a method to study the dose distributions resulting from MBRT irradiations of a rat femur and a realistic human femur phantoms. As a figure of merit, peak and valley doses and peak-to-valley dose ratios (PVDR) were assessed. Conversion of absorbed dose to normalized total dose (NTD) was performed in the human case. Several field sizes and irradiation geometries were evaluated. Results: It is feasible to deliver a uniform dose distribution in the target while the healthy tissue benefits from a spatial fractionation of the dose. Very high PVDR values (⩾20) were achieved in the entrance beam path in the rat case. PVDR values ranged from 2 to 9 in the human phantom. NTD{sub 2.0} of 87 Gy might be reached in the tumor in the human femur while the healthy tissues might receive valley NTD{sub 2.0} lower than 20 Gy. The doses in the tumor and healthy tissues might be significantly higher and lower than the ones commonly delivered used in conventional radiotherapy. Conclusions: The obtained dose distributions indicate that a gain in normal tissue sparing might be expected. This would allow the use of higher (and potentially curative) doses in the tumor. Biological experiments are warranted.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2016

SU-F-T-159: Monte Carlo Simulation Studies of Three-Dimensional Dose Distribution for Polymer Gel Dosimeter and Radiochromic Gel Dosimeter in a Proton Beam

Purpose: The purpose of this simulation study is to evaluate the proton detectability of gel dosimeters, and estimate the three-dimensional dose distribution of protons in the radiochromic gel and polymer gel dosimeter compared with the dose distribution in water. Methods: The commercial composition ratios of normoxic polymer gel and LCV micelle radiochromic gel were included in this simulation study. The densities of polymer and radiochromic gel were 1.024 and 1.005 g/cm3, respectively. The 50, 80 and 140 MeV proton beam energies were selected. The dose distributions of protons in the polymer and radiochromic gel were simulated using Monte Carlo radiation transport code (MCNPX 2.7.0, Los Alamos Laboratory). The water equivalent depth profiles and the dose distributions of two gel dosimeters were compared for the water. Results: In case of irradiating 50, 80 and 140 MeV proton beam to water phantom, the reference Bragg-peak depths are represented at 2.22, 5.18 and 13.98 cm, respectively. The difference in the water equivalent depth is represented to about 0.17 and 0.37 cm in the radiochromic gel and polymer gel dosimeter, respectively. The proton absorbed doses in the radiochromic gel dosimeter are calculated to 2.41, 3.92 and 6.90 Gy with increment of incident proton energies. In the polymer gel dosimeter, the absorbed doses are calculated to 2.37, 3.85 and 6.78 Gy with increment of incident proton energies. The relative absorbed dose in radiochromic gel (about 0.47 %) is similar to that of water than the relative absorbed dose of polymer gel (about 2.26 %). In evaluating the proton dose distribution, we found that the dose distribution of both gel dosimeters matched that of water in most cases. Conclusion: As the dosimetry device, the radiochromic gel dosimeter has the potential particle detectability and is feasible to use for quality assurance of proton beam therapy beam.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2015

COMPUTATIONAL MODELING OF THE CLASS I LOW-MASS PROTOSTAR ELIAS 29 APPLYING OPTICAL CONSTANTS OF ICES PROCESSED BY HIGH ENERGY COSMIC RAY ANALOGS

We present a study of the effects of high energy cosmic rays (CRs) over the astrophysical ices, observed toward the embedded class I protostar Elias 29, by using computational modeling and laboratory data. Its spectrum was observed with the Infrared Space Observatory (ISO) covering 2.3–190 μm. The modeling employed the three-dimensional Monte Carlo radiative transfer code RADMC-3D and laboratory data of bombarded ice grains by CR analogs and unprocessed ices (not bombarded). We are assuming that Elias 29 has a self-irradiated disk with inclination i = 60.°0, surrounded by an envelope with a bipolar cavity. The results show that absorption features toward Elias 29 are better reproduced by assuming a combination between unprocessed astrophysical ices at low temperature (H{sub 2}O, CO, CO{sub 2}) and bombarded ices (H{sub 2}O:CO{sub 2}) by high energy CRs. Evidences of the ice processing around Elias 29 can be observed by the good fitting around 5.5–8.0 μm, by polar and apolar ice segregation in 15.15–15.25 μm, and by the presence of the CH{sub 4} and HCOOH ices. Given that non-nitrogen compounds were employed in this work, we assume that absorption around 5.5–8.0 μm should not be associated with the NH{sub 4}{sup +} ion (see the 2003 work of Shutte and Khanna ), but more probably with aliphatic ethers (e.g., R1-OCH{sub 2}-R2), CH{sub 3}CHO, and related species. The results obtained in this paper are important because they show that the environment around protostars is better modeled considering processed samples and, consequently, demonstrate the chemical evolution of the astrophysical ices.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

Determination of peripheral underdosage at the lung-tumor interface using Monte Carlo radiation transport calculations

Prediction of dose distributions in close proximity to interfaces is difficult. In the context of radiotherapy of lung tumors, this may affect the minimum dose received by lesions and is particularly important when prescribing dose to covering isodoses. The objective of this work is to quantify underdosage in key regions around a hypothetical target using Monte Carlo dose calculation methods, and to develop a factor for clinical estimation of such underdosage. A systematic set of calculations are undertaken using 2 Monte Carlo radiation transport codes (EGSnrc and GEANT4). Discrepancies in dose are determined for a number of parameters, including beam energy, tumor size, field size, and distance from chest wall. Calculations were performed for 1-mm{sup 3} regions at proximal, distal, and lateral aspects of a spherical tumor, determined for a 6-MV and a 15-MV photon beam. The simulations indicate regions of tumor underdose at the tumor-lung interface. Results are presented as ratios of the dose at key peripheral regions to the dose at the center of the tumor, a point at which the treatment planning system (TPS) predicts the dose more reliably. Comparison with TPS data (pencil-beam convolution) indicates such underdosage would not have been predicted accurately in the clinic. We define a dose reduction factor (DRF) as the average of the dose in the periphery in the 6 cardinal directions divided by the central dose in the target, the mean of which is 0.97 and 0.95 for a 6-MV and 15-MV beam, respectively. The DRF can assist clinicians in the estimation of the magnitude of potential discrepancies between prescribed and delivered dose distributions as a function of tumor size and location. Calculation for a systematic set of 'generic' tumors allows application to many classes of patient case, and is particularly useful for interpreting clinical trial data.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2015

ULTRAVIOLET RADIATIVE TRANSFER MODELING OF NEARBY GALAXIES WITH EXTRAPLANAR DUSTS

In order to examine their relation to the host galaxy, the extraplanar dusts of six nearby galaxies are modeled, employing a three-dimensional Monte Carlo radiative transfer code. The targets are from the highly inclined galaxies that show dust-scattered ultraviolet halos, and the archival Galaxy Evolution Explorer FUV band images were fitted with the model. The observed images are generally well-reproduced by two dust layers and one light source layer, whose vertical and radial distributions have exponential profiles. We obtained several important physical parameters, such as star formation rate (SFR{sub UV}), face-on optical depth, and scale-heights. Three galaxies (NGC 891, NGC 3628, and UGC 11794) show clear evidence for the existence of an extraplanar dust layer. However, it is found that the remaining three targets (IC 5249, NGC 24, and NGC 4173) do not necessarily need a thick dust disk to model the ultraviolet (UV) halo, because its contribution is too small and the UV halo may be caused by the wing part of the GALEX point spread function. This indicates that the galaxy samples reported to have UV halos may be contaminated by galaxies with negligible extraplanar (halo) dust. The galaxies showing evidence of an extraplanar dust layer fall within a narrow range on the scatter plots between physical parameters such as SFR{sub UV} and extraplanar dust mass. Several mechanisms that could possibly produce the extraplanar dust are discussed. We also found a hint that the extraplanar dust scale-height might not be much different from the polycyclic aromatic hydrocarbon emission characteristic height.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Optimization of a dry, mixed nuclear fuel storage array for nuclear criticality safety - 14339

A dry storage array of used nuclear fuel at the Idaho National Laboratory contains a mixture of more than twenty different types of research and test reactor fuels. Significant neutronic interaction occurs between fuel canisters within this large storage arrangement of 636 fuel storage canisters. Nuclear criticality safety is maintained primarily by moderation control and through control of neutron interaction and leakage provided by the credited geometry. Over the nearly forty years of operation the collected fuels are now arranged in a configuration that does not minimize the system neutron multiplication factor (k{sub eff}). This creates a concern since the entire remaining facility storage capacity cannot be utilized without exceeding the system minimum subcritical state (k{sub safe}) for ensuring nuclear criticality safety. This work determines a more optimal arrangement for the stored fuels with a goal to minimize the system k{sub eff}, but also with a minimum of potential fuel canister relocation movements. The ultimate goal is to establish a solution for this multiple-objective optimization problem by finding the Pareto Front for maximum k{sub eff} reduction versus the number of canister relocations. The process first required additional analysis of the array by Monte-Carlo neutron transport codes (specifically, MCNP) to gain a greater understanding of the interaction characteristics of the disparate fuel types. Some stored fuels are composed of primarily metallic and oxide components that provide very little inherent moderation, while other fuel types contain graphite or zirconium-hydride which can contribute a significant amount of moderation. These factors cause a complex relationship within a fuel arrangement with respect to thermal utilization and axial leakage probability. The new MCNP analysis was utilized to develop a simplified spreadsheet-based empirical model which simulates the array k{sub eff} result. This empirical model was benchmarked against additional MCNP calculations for various arrangements and then employed in the optimization analysis. A Tabu metaheuristic and stochastic approximation optimization was applied to the empirical model with the optimization function to minimize the array k{sub eff}. The optimization algorithm interrogated the empirical model to determine the optimal rearrangement for up to sixty canister movements. The Pareto Front was found, which shows that a limited number of canister relocation movements could allow the existing thirty-one empty array positions to be filled with common research reactor fuels while also yielding a reduction in the overall array k{sub eff} relative to the current arrangement. This will impart an improvement in the facility utilization while also offering an enhancement in the safety margin. The reduction in array reactivity has been confirmed by a final set of MCNP calculations with the proposed new arrangements. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Development of a Griffin model of the advanced test reactor

In the pursuit of a higher fidelity deterministic simulation capability of the Advanced Test Reactor, it is important to have a fast yet accurate deterministic neutronics model. Here, to achieve this, we employed an advanced two-step method. The first step involves generating homogenized cross sections using OpenMC, a cutting-edge Monte Carlo neutron transport code. OpenMC offers excellent modular capabilities, allowing for easy component integration and flexibility in incorporating new designs into the model. The second step involves deterministic transport calculations, which are performed using Griffin, a reactor physics application based on the Multiphysics Object-Oriented Simulation Environment (MOOSE). To ensure the accurate spatial resolution and assignment of material cross sections, a Cubit-generated mesh for the Advanced Test Reactor is utilized as an intermediate step between the OpenMC and Griffin models; Griffin utilizes the mesh for its finite element solution, while OpenMC material identifications are written to the mesh file to be used in Griffin material assignments. Additionally, a Python-based script converts the cross sections generated by OpenMC into the ISOXML format required by Griffin. Initial comparisons using the Griffin diffusion solver indicated good agreement between the neutron multiplication factors obtained from the standalone OpenMC model and the Griffin model, with differences of less than 10 pcm in the 2D geometry configuration; it was later determined that this agreement was likely due to compensating effect and was more likely on the order of –700 pcm relative to the OpenMC solution. However, in three-dimensional calculations, an unacceptably large error (almost 8,000 pcm) was found in the Griffin solution with the diffusion solver. Subsequent calculations using Griffin’s discrete ordinates solver demonstrated substantially improved agreement, within 116 pcm of the OpenMC solution used to generate the cross sections for Griffin. Building on this capability, future work will seek to perform more detailed validation calculations. The ultimate goal is to evaluate both transient and multiphysics simulations of the reactor.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2013

Obtaining identical results with double precision global accuracy on different numbers of processors in parallel particle Monte Carlo simulations

We describe and compare different approaches for achieving numerical reproducibility in photon Monte Carlo simulations. Reproducibility is desirable for code verification, testing, and debugging. Parallelism creates a unique problem for achieving reproducibility in Monte Carlo simulations because it changes the order in which values are summed. This is a numerical problem because double precision arithmetic is not associative. Parallel Monte Carlo, both domain replicated and decomposed simulations, will run their particles in a different order during different runs of the same simulation because the non-reproducibility of communication between processors. In addition, runs of the same simulation using different domain decompositions will also result in particles being simulated in a different order. In [1], a way of eliminating non-associative accumulations using integer tallies was described. This approach successfully achieves reproducibility at the cost of lost accuracy by rounding double precision numbers to fewer significant digits. This integer approach, and other extended and reduced precision reproducibility techniques, are described and compared in this work. Increased precision alone is not enough to ensure reproducibility of photon Monte Carlo simulations. Non-arbitrary precision approaches require a varying degree of rounding to achieve reproducibility. For the problems investigated in this work double precision global accuracy was achievable by using 100 bits of precision or greater on all unordered sums which where subsequently rounded to double precision at the end of every time-step.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

Development of array-type prompt gamma measurement system for in vivo range verification in proton therapy

Purpose: In vivo range verification is one of the most important parts of proton therapy to fully utilize its benefits delivering high radiation dose to tumor, while sparing the normal tissue with the so-called Bragg peak. Currently, however, range verification method is not used in clinics. The purpose of the present study is to optimize and evaluate the configuration of an array-type prompt gamma measurement system on determining distal dose edge for in vivo range verification of proton therapy. Methods: To effectively measure the prompt gammas against the background gammas, the Monte Carlo simulations with the MCNPX code were employed in optimizing the configuration of the measurement system, and the Monte Carlo method was also used to understand the effect of the background gammas, mainly neutron capture gammas, in the measured gamma distribution. To reduce the effect of the background gammas, the optimized energy window of 4-10 MeV in measuring the prompt gammas was employed. A parameterized source was used to maximize computation speed in the optimization study. A simplified test measurement system, using only one detector moving from one measurement location to the next, was constructed and applied to therapeutic proton beams of 80-220 MeV. For accurate determination of the distal dose edge, the sigmoidal curve-fitting method was applied to the measured distributions of the prompt gammas, and then, the location of the half-value between the maximum and minimum value in the curve-fitting was determined as the distal dose edge and compared with the beam range assessed by the proton dose distribution. Results: The parameterized source term employed in optimization process improved the calculation speed by up to {approx}300 times. The optimization study indicates that an array-type measurement system with 3, 2, 2, and 150 mm for scintillator thickness, slit width, septal thickness, and slit length, respectively, can effectively measure the prompt gamma distributions minimizing the contribution of background gammas. The present results show that a few hundred counts of prompt gammas can be easily obtained by measuring 10 s at each measurement location for proton beams of {approx}4 nA. The distal dose edges determined by the prompt gamma distribution are 5.45, 14.73, and 27.74 cm for proton beams of 5.17 (80 MeV), 14.99 (150 MeV), and 27.38 (220 MeV) cm, respectively. Conclusions: The results show that the array-type measurement system can measure prompt gamma distributions from a therapeutic proton beam within a short measurement time, and that the distal dose edge can be determined within a few millimeters of error without using any sophisticated analysis.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2019

About the Impact of the Unresolved Resonance Region in Monte Carlo Simulations of Sodium Fast Reactors

In the last few years, and within the framework of different European projects, KENO-VI code from SCALE system has been employed to perform detailed continuous-energy Mon-te Carlo transport calculations for advanced fast reactors. The core characterization of both the sodium-cooled ASTRID and the lead-cooled ALFRED reactors was performed during the FP7 cross-cutting ESNII+ project; more recently, core calculations for the sodium-cooled Superphénix reactor and the improved European Sodium Fast Reactor design were performed within the HORIZON2020 ESFR-SMART project. In all cases, the effective multiplication factor predicted by KENO-VI was systematically higher (around 400-500 pcm) than the values computed by MCNP and Serpent Monte Carlo codes, using the same nu-clear data library.In order to provide insight into the origin of the observed discrepancies, a simplified 2D MOX-fueled SFR pin-cell benchmark has been launched. The multiplication factor, as well as 1-group and VITAMINJ 175-group cross-sections computed by KENO-VI, Serpent and MCNP codes employing ENDF/B-VII.1 data library, have been compared.Significant differences between KENO-VI and the other codes have been found in the un-resolved resonance regions of 239Pu and 241Pu capture and production cross sections, while negligible differences appeared outside those energy ranges. On the other hand, calculations without using probability tables have shown very good agreement. Quantita-tive comparison is presented and analyzed, along with a discussion of the impact of the probability-table treatment in the three codes for MOX-fueled systems with typical SFR spectrum.

Jimenez carrascosa, A↗
DOE OSTI2018

An angular biasing method using arbitrary convex polyhedra for Monte Carlo radiation transport calculations

This paper presents a new method for performing angular biasing in Monte Carlo radiation transport codes using arbitrary convex polyhedra to define regions of interest toward which to project particles (DXTRAN regions). The method is derived and is implemented using axis-aligned right parallelepipeds (AARPPs) and arbitrary convex polyhedra. Attention is also paid to possible numerical complications and areas for future refinement. A series of test problems are executed with void, purely absorbing, purely scattering, and 50% absorbing/50% scattering materials. For all test problems tally results using AARPP and polyhedral DXTRAN regions agree with analog and/or spherical DXTRAN results within statistical uncertainties. In cases with significant scattering the figure of merit (FOM) using AARPP or polyhedral DXTRAN regions is lower than with spherical regions despite the ability to closely fit the tally region. This is because spherical DXTRAN processing is computationally less expensive than AARPP or polyhedral DXTRAN processing. Thus, it is recommended that the speed of spherical regions be considered versus the ability to closely fit the tally region with an AARPP or arbitrary polyhedral region. It is also recommended that short calculations be made prior to final calculations to compare the FOM for the various DXTRAN geometries because of the influence of the scattering behavior.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Criticality concerns of a group actinide co-crystallization separations approach to used nuclear fuel recycling

We report the co-crystallization of hexavalent actinides with uranyl nitrate hexahydrate has been established as a viable approach to achieving group actinide separation in used nuclear fuel recycle while avoiding the difficulties associated with the solvent-extraction process altogether. A preliminary neutronics safety analysis has been performed to alleviate criticality concerns associated with this group hexavalent actinide co-crystallization approach. The basic uranyl nitrate hexahydrate material compositions used are representative of used fuel found at commercial pressurized water reactors all around the U.S. This material is subsequently doped with various percentages of either typical reactor grade plutonium or with the total transuranic vector, replacing the uranium atoms in the crystal lattice structure. Infinite and effective neutron multiplication factors were calculated using the state-of-the-art Monte Carlo neutron transport code MCNP in a variety of plausible material configurations, assuming an acceptable safety limit for criticality to be an effective neutron multiplication factor less than 0.95. The results indicate that the crystal structures form a very subcritical intrinsic geometry even with six water molecules being incorporated into the crystal structure. In conclusion, the results of this research suggest that significant opportunity exists for the safe transfer of the group actinide separation approach from a laboratory to an operational scale.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

Measurement of leakage neutron spectra for zirconium with D-T neutrons and validation of evaluated nuclear data

Integral neutronics experiments have been a powerful tool and method for examining the reliability of evaluated nuclear data. In the present paper, the accuracy of evaluated nuclear data for Zirconium has been validated by comparing measured leakage neutron spectra with calculated ones. Leakage neutron spectra from the irradiation of D-T neutrons on Zirconium samples were experimentally measured at 60° and 120° by using a time-of-flight method. Theoretical calculations are carried out by Monte Carlo neutron transport code MCNP-4C with evaluated nuclear data of the ENDF/B-VII.0, ENDF/B-VIII.0, JENDL-4.0 and CENDL-3.1 libraries. From the comparisons, it is found that the calculations with JENDL-4.0 gives better agreements with the experiments at 60° while a large discrepancies are observed in the 4–13.5 MeV energy range at 120°. However, the calculated spectra using the data from ENDF/B-VII.0, ENDF/B-VIII.0 and CENDL-3.1 libraries showed some discrepancies with the measured ones, especially caused by the difference of the elastic and inelastic angular distributions and the secondary neutron energy distributions in the continuum inelastic scattering.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Subject-specific modeling framework for particle deposition using computational fluid dynamics

Quantifying particle deposition and dose in the respiratory tract requires a physiologically realistic representation and reproducible computational workflows. However, existing modeling frameworks, such as the International Commission on Radiological Protection (ICRP) compartmental models and the Multiple Path Particle Dosimetry (MPPD) tool, lack detailed deposition profiles and subject-specific capabilities. The combination of advances in computer vision algorithms applied to the respiratory tract and Computational Fluid and Particle Dynamics (CFPD) allows high-fidelity simulations of particle behavior in anatomically accurate geometries derived from individual CT scans. The segmentation, preprocessing, and file preparation task for a CFPD simulation was often time-consuming, and no prior studies to-date have yet presented a fully automated framework. This work presents a fully automated workflow to obtain individualized particle deposition profiles in the human respiratory tract. The pipeline starts with segmenting upper and lower airway geometries using morphological and deep learning-based methods, generating three-dimensional (3D) models from CT imaging data. Next, a series of algorithms are presented to quality check and prepare the 3D geometry for a CFD or CFPD simulation. The preprocessing step includes correcting geometric artifacts, enforcing a physically consistent mesh, and automatically identifying and capping multiple outlets, which is required for CFD/CFPD simulations. These processed models are then input into open-source (OpenFOAM) or commercial (StarCCM+) CFD solvers, where flow and transient particle transport equations — including turbulence and particle–wall interactions are solved under realistic breathing conditions. Finally, the resulting particle deposition profiles can be integrated with Monte Carlo radiation transport codes and state-of-the-art computational phantoms to assess organ-specific absorbed doses in scenarios of radioactive aerosol inhalation. The presented work streamlines respiratory tract segmentation, preprocessing for CFD/CFPD simulations, and integration with dose assessment workflows, reducing manual intervention and improving access to high-fidelity, subject-specific modeling. The high precision in predicted particle deposition and dose distributions can improve personalized treatment strategies in respiratory medicine and refine dose estimates for radiation protection.

AI↗
DOE OSTI2013

Excited State Electronic Properties of Sodium Iodide and Cesium Iodide

We compute from first principles the dielectric function, loss function, lifetime and scattering rate of quasiparticles due to electronic losses, and secondary particle spectrum due to plasmon decay in two scintillating alkali halides, sodium iodide and cesium iodide. Particular emphasis is placed on quasiparticles within several multiples of the band gap from the band edges. A theory for the decay spectra of plasmons and other electronic excitations in crystals is presented. Applications to Monte Carlo radiation transport codes are discussed.

Campbell, Luke W↗
DOE OSTI2014

Reconstruction of NSTX midplane neutral density profiles from visible imaging data

In this report the experimental determination of neutral densities in tokamak plasmas from line radiation is only accurate in the narrow region in which both the excitation rate and neutral density are significant. We describe an alternative procedure using the DEGAS 2 Monte Carlo neutral transport code to invert light emission data obtained from a tangentially viewing camera, yielding absolute radial profiles of deuterium atoms and molecules at midplane. That the neutral source in these simulations can be adequately characterized as a uniform flux at the vacuum vessel wall is demonstrated by the similarity of the shapes of the simulated and observed brightness profiles. A second test is obtained by comparing the resulting neutral pressures at the vessel walls with data from midplane micro-ion gauges. We also show that the simulated camera image is insensitive to variations in the spatial distribution of the neutral source.

(PSI-21) NSTX DEGAS↗
DOE OSTI2026

Surrogate modeling of Monte Carlo radiation transport with convolutional neural networks for shielding optimization

Here, we present a machine learning (ML)-based surrogate model using convolutional neural networks (CNN) designed to emulate the attenuation of neutron fields as they pass through various shielding materials. This model can compute the outgoing neutron flux almost instantaneously and achieves reasonable accuracy compared to traditional Monte Carlo (MC)-based codes, which are computationally intensive. This emulator alleviates the complexity of neutron radiation transport through shielding materials by reducing the dimensionality and enables shielding optimization for a known radiation environment. This optimization process, which would have taken an unrealistic timeline due to several complex radiation transport simulations, can now be achieved in minutes, thus increasing computational capabilities in radiation shielding assessment. We demonstrate the applications of this emulator in computing effective dose rates and optimizing shielding solutions for a heavy-ion accelerator facility, such as the Facility for Rare Isotope Beams, where secondary neutrons produced via beam interactions dominate the radiation environment.

accelerator shielding↗