Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “predictive simulation”

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,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Fiber motion in highly confined flows of carbon fiber and non-Newtonian polymer [Advanced modeling of non-Newtonian polymer and carbon fibers]

Inks compounded of short carbon fibers suspended in polymer resin can be extruded to produce composite materials during additive manufacturing or 3D printing processes. The flow process induces anisotropic orientation of the fibers which is set into the matrix and significantly affects the mechanical and physical properties of the final product. Therefore, the flow of fiber suspensions needs to be understood in order to predict the orientation distribution of the fibers during such manufacturing processes. There is still a lack of knowledge for extrusion of the complex mixture of a non-Newtonian polymer containing a high-volume fraction of fibers with high fiber aspect ratio, where both inter-particle, fluid-particle and fiber-wall interactions are computationally evaluated. This paper presents predictive numerical simulations and experimental results of such confined flow in a concentrated regime. The code is based on Lagrange multiplier technique and resolves each particle and interaction between the fibers and surrounding fluid and nozzle walls. We investigate numerically how the fiber length impacts the fiber alignment during extrusion. We found that a fiber length above 67% of the nozzle's diameter induces dramatic change in the fiber flow, causing fibers to concentrate at the nozzle boundaries with very low concentration of fibers in the center of the nozzle. It was found that the best fiber alignment is reached for fiber lengths equal to 40–50% of the nozzle diameters. Numerical findings are supported by experimental results. In conclusion, this work improves understanding of fiber orientation during 3D printing and is an important milestone for the prediction of the complex mechanical properties of additively manufactured fiber composites.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Assessment the leaching characteristics and long-term leaching behavior of some radionuclides from synthesized zeolite cement matrix

Highlights: • Synthesis and characterization of zeolite cement synthesized from fly ash • Leaching tests were used to assess the release rate of {sup 85}Sr and {sup 60}Co radionuclides. • Simple mathematical models were used to detect leaching mechanisms. • Long-term leaching behavior of {sup 85}Sr and {sup 60}Co from cementitious matrix was assessed. The development of cementitious materials remains a vital goal to produce valuable products with good mechanical, physical and chemical properties suitable for the safe disposal of concentrated radionuclides resulted from the treatment of contaminated solutions. In this paper, zeolite cement with good properties was laboratory synthesized from industrial by-products fly ash and characterized using various characterization techniques. The leaching characteristics of strontium and cobalt radionuclides from the synthesized zeolite cement were scrutinized conforming with the standard leaching methodology of International Atomic Energy Authority (IAEA). A mechanical strength assessment was executed to characterize the extent of immobilization process of the solidified matrices. The cumulative leaching fraction (φ) of the two studied radionuclides was found to be less than 5% in all examined conditions, which implies the applicability of the IAEA-recommended methodology for estimating the diffusion coefficient. The experimental leaching data were regressed nonlinearly to various mathematical kinetic models to assess the controlling leaching mechanism and to determine the leaching parameters. The regression results indicated that strontium ({sup 85}Sr) and cobalt ({sup 60}Co) leaching resulted from two succeeding mechanisms: as first order kinetic reaction, and then diffusion. The calculated values of leachability indices signify that the performance of each of the studied matrices is within an acceptable range. A simplified mathematical model, rooted in the first order reaction and diffusion mechanisms, was simulated to predict the radionuclides leaching rates from zeolite cement matrices. By comparing the synthesized zeolite cement with other cementitious materials, it can be concluded that the synthesized material can be classified as an efficient material suitable to immobilize {sup 85}Sr and {sup 60}Co from radioactive wastes. The acquired findings demonstrated that the studied immobilized waste matrices have acceptable mechanical effectiveness.

36 MATERIALS SCIENCE↗
DOE OSTI2017

N-body simulations of gravitational redshifts and other relativistic distortions of galaxy clustering

Large redshift surveys of galaxies and clusters are providing the first opportunities to search for distortions in the observed pattern of large-scale structure due to such effects as gravitational redshift. We focus on non-linear scales and apply a quasi-Newtonian approach using N-body simulations to predict the small asymmetries in the cross-correlation function of two galaxy different populations. Following recent work by Bonvin et al., Zhao and Peacock and Kaiser on galaxy clusters, we include effects which enter at the same order as gravitational redshift: the transverse Doppler effect, light-cone effects, relativistic beaming, luminosity distance perturbation and wide-angle effects. We find that all these effects cause asymmetries in the cross-correlation functions. Quantifying these asymmetries, we find that the total effect is dominated by the gravitational redshift and luminosity distance perturbation at small and large scales, respectively. By adding additional subresolution modelling of galaxy structure to the large-scale structure information, we find that the signal is significantly increased, indicating that structure on the smallest scales is important and should be included. We report on comparison of our simulation results with measurements from the SDSS/BOSS galaxy redshift survey in a companion paper.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Effects of Mineral Impurities and Heat on Uranium(VI) Sorption onto Bentonite (Final Report)

Most nuclear waste disposal options include compacted bentonite, consisting of montmorillonite clay, as part of a barrier system to minimize contaminant mobility. Uranium (U) is the primary element in spent nuclear fuel, and from an environmental perspective, a potential contaminant of water resources. Furthermore, uranium may control nuclear fuel degradation rates and the consequent release of other radioactive contaminants based on its diffusive transport away from waste canisters. Uranium sorption onto clay and other mineral surfaces is expected to limit U(VI) mobility in these systems. However, at this point a prediction of U(VI) sorption and transport behavior in performance assessment (PA) models is complicated by a series of factors, such as: (1) the presence of various U(VI) solution species with different charges and sorption characteristics; (2) the complex microstructure of montmorillonite clay resulting in two types of clay porosities and multiple sorption sites; (3) the largely unknown effects of bentonite mineral impurities on pore water chemistry, U(VI) sorption and diffusion behavior; and (4) the potential impacts of heat, generated by the decay of spent fuel, on mineralogical and microstructural transformations, and any subsequent effects on radionuclide sorption and mobility. In this project, we investigated the effects of calcite impurities on U(VI) sorption and diffusion onto montmorillonite before and after mineral exposure to heat. Overall, our main goals were to: 1. Gain a mechanistic understanding of the effects of calcite impurities on U(VI) sorption and diffusion in montmorillonite/bentonite at variable chemical conditions, and before and after mineral exposure to heat; and to 2. Work towards the development of new surface complexation and diffusion models to predict U(VI) sorption and mobility in montmorillonite/bentonite while accounting for the impacts of calcite impurities and heat. In order to achieve these goals, we pursued the following research objectives and tasks: 1. Characterize the effects of calcite impurities on U(VI) sorption onto Na-montmorillonite, and evaluate if and how these effects change due to solids’ exposure to heat (Task 1). 2. Characterize U(VI) surface speciation in solid samples from Task 1, with a particular focus on U(VI) precipitates and co-precipitates (Task 2). 3. Test the relevance of observed changes in U(VI) sorption characteristics for U(VI) diffusive transport behavior in bentonite barriers (Task 3). 4. Evaluate the exclusion of anionic U(VI) solution complexes from clay nanopores (Task 4). 5. Building on Tasks 1-4, start the development of new surface complexation and reactive transport models that can simulate and predict U(VI) sorption and diffusion in complex bentonite systems (Task 5). A detailed description of all tasks, and their specific results and conclusions is provided in this report.

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

A coarse-grained model of clay colloidal aggregation and consolidation with explicit representation of the electrical double layer

The aggregation of clay minerals in liquid water exemplifies colloidal self-assembly in nature. These negatively charged aluminosilicate platelets interact through multiple mechanisms with different sensitivities to particle shape, surface charge, aqueous chemistry, and interparticle distance and exhibit complex aggregation structures. Experiments have difficulty resolving the associated colloidal assemblages at the scale of individual particles. Conversely, all-atom molecular dynamics (MD) simulations provide detailed insight on clay colloidal interaction mechanisms, but they are limited to systems containing a few particles. We develop a new coarse-grained (CG) model capable of representing assemblages of hundreds of clay particles with accuracy approaching that of MD simulations, at a fraction of the computational cost. Our CG model is parameterized based on MD simulations of a pair of smectite clay particles in liquid water. A distinctive feature of our model is that it explicitly represents the electrical double layer (EDL), i.e., the cloud of charge-compensating cations that surrounds the clay particles. Our model captures the simultaneous importance of long-range colloidal interactions (i.e., interactions consistent with simplified analytical models, already included in extant clay CG models) and short-range interactions such as ion correlation and surface and ion hydration effects. The resulting simulations correctly predict, at low solid-water ratios, the existence of ordered arrangements of parallel particles separated by water films with a thickness up to ~10 nm and, at high solid-water ratios, the coexistence of crystalline and osmotic swelling states, in agreement with experimental observations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Search for interacting galaxy clusters from SDSS DR-17 employing optimized friends-of-friends algorithm and multimessenger tracers

ABSTRACT In the theoretical framework of hierarchical structure formation, galaxy clusters evolve through continuous accretion and mergers of substructures. Cosmological simulations have revealed the best picture of the universe as a 3D filamentary network of dark-matter distribution called the cosmic web. Galaxy clusters are found to form at the nodes of this network and are the regions of high merging activity. Such mergers being highly energetic, contain a wealth of information about the dynamical evolution of structures in the Universe. Observational validation of this scenario needs a colossal effort to identify numerous events from all-sky surveys. Therefore, such efforts are sparse in literature and tend to focus on individual systems. In this work, we present an improved search algorithm for identifying interacting galaxy clusters and have successfully produced a comprehensive list of systems from SDSS DR-17. By proposing a set of physically motivated criteria, we classified these interacting clusters into two broad classes, ‘merging’ and ‘pre-merging/postmerging’ systems. Interestingly, as predicted by simulations, we found that most cases show cluster interaction along the prominent cosmic filaments of galaxy distribution (i.e. the proxy for dark matter filaments), with the most violent ones at their nodes. Moreover, we traced the imprint of interactions through multiband signatures, such as diffuse cluster emissions in radio or X-rays. Although we could not find direct evidence of diffuse emission from connecting filaments and ridges; our catalogue of interacting clusters will ease locating such faintest emissions as data from sensitive telescopes such as eROSITA or SKA, becomes accessible.

Oak, Tejas↗
DOE OSTI2020

Faint and Fading Tails: The Fate of Stripped H i Gas in Virgo Cluster Galaxies

Although many galaxies in the Virgo cluster are known to have lost significant amounts of H i gas, only about a dozen features are known where the H i extends significantly outside its parent galaxy. Previous numerical simulations have predicted that H i removed by ram pressure stripping should have column densities far in excess of the sensitivity limits of observational surveys. We construct a simple model to try and quantify how many streams we might expect to detect. This accounts for the expected random orientation of the streams in position and velocity space as well as the expected stream length and mass of stripped H i. Using archival data from the Arecibo Galaxy Environment Survey, we search for any streams that might previously have been missed in earlier analyses. We report the confident detection of 10 streams as well as 16 other less-certain detections. We show that these well match our analytic predictions for which galaxies should be actively losing gas; however, the mass of the streams is typically far below the amount of missing H i in their parent galaxies, implying that a phase change and/or dispersal renders the gas undetectable. By estimating the orbital timescales, we estimate that dissolution rates of 1–10 M {sub ⊙} yr{sup −1} are able to explain both the presence of a few long, massive streams and the greater number of shorter, less-massive features.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Radiative Energy Balance in the Tropical Tropopause Layer: An Investigation with ARM Data

The overall objective of this project is to use the ARM observational data to improve our understanding of cloud-radiation effects in the tropical tropopause layer (TTL), which is crucial for improving the simulation and prediction of climate and climate change. In last four and half years, we have been concentrating on (i) performing the comparison of the ice cloud properties from the ground-based lidar observations with those from the satellite CALIPSO lidar observations at the ARM TWP sites; (ii) analyzing TTL cirrus and its relation to the tropical planetary waves; (iii) calculating the radiative heating rates using retrieved cloud microphysical properties by combining the ground-based lidar and radar observations at the ARM TWP sites and comparing the results with those using cloud properties retrieved from CloudSat and CALIPSO observations; (iv) comparing macrophysical properties of tropical cirrus clouds from the CALIPSO satellite and from ground-based micropulse and Raman lidar observations; (v) improving the parameterization of optical properties of cirrus clouds with small effective ice particle sizes; and (vi) evaluating the enhanced maximum warming in the tropical upper troposphere simulated by the GCMs. The main results of our research efforts are reported in the 12 referred journal publications that acknowledge the DOE Grant No. DE-FG02-09ER64769.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

PSA 2025 Presentation: "Modeling and Sensitivity Analysis of a Generation IV Pebble Bed Reactor Using MELCOR 2.2"

Accompanying the advancement of reactor technologies is the need for computational modeling and simulation to predict their behavior under normal operating conditions and accident scenarios. New Generation IV reactor designs which employ non-conventional fuel have a particular need for modeling the behavior and release of radionuclides and other material from the fuel. In this work, MELCOR version 2.2, a system-level safety and accident scenario code developed by Sandia National Laboratories, was used to model a 200-MWth pebble bed modular reactor and calculate the inventories of circulating and deposited graphite, metal dust, and elemental components released from the fuel elements. A base case modeling the reactor under standard operating conditions was calculated using MELCOR and the inventories were extrapolated to 30 years of operation time using a logarithmic regression fit. A sensitivity analysis was also performed in which several key parameters for the base case model were modified to explore the effect of these changes on the inventories calculated by MELCOR. A set of transient scenario simulations for a depressurized loss of forced cooling (DLOFC) accident were also performed. The results of the sensitivity analysis and transient simulations are reported and discussed in relation to the modeling techniques used for this study.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2013

Laser produced plasma sources for nanolithography—Recent integrated simulation and benchmarking

Photon sources for extreme ultraviolet lithography (EUVL) are still facing challenging problems to achieve high volume manufacturing in the semiconductor industry. The requirements for high EUV power, longer optical system and components lifetime, and efficient mechanisms for target delivery have narrowed investigators towards the development and optimization of dual-pulse laser sources with high repetition rate of small liquid tin droplets and the use of multi-layer mirror optical system for collecting EUV photons. We comprehensively simulated laser-produced plasma sources in full 3D configuration using 10–50 μm tin droplet targets as single droplets as well as, for the first time, distributed fragmented microdroplets with equivalent mass. The latter is to examine the effects of droplet fragmentation resulting from the first pulse and prior to the incident second main laser pulse. We studied the dependence of target mass and size, laser parameters, and dual pulse system configuration on EUV radiation output and on atomic and ionic debris generation. Our modeling and simulation included all phases of laser target evolution: from laser/droplet interaction, energy deposition, target vaporization, ionization, plasma hydrodynamic expansion, thermal and radiation energy redistribution, and EUV photons collection as well as detail mapping of photons source size and location. We also simulated and predicted the potential damage to the optical mirror collection system from plasma thermal and energetic debris and the requirements for mitigating systems to reduce debris fluence. The debris effect on mirror collection system is analyzed using our three-dimensional ITMC-DYN Monte Carlo package. Modeling results were benchmarked against our CMUXE laboratory experimental studies for the EUV photons production and for debris and ions generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Measurements of b -jet tagging efficiency with the ATLAS detector using $t\overline{t}$ events at $\sqrt{s}=13$ TeV

The efficiency to identify jets containing b-hadrons (b-jets) is measured using a high purity sample of dileptonic top quark-antiquark pairs ($t\bar{t}$) selected from the 36.1 fb -1 of data collected by the ATLAS detector in 2015 and 2016 from proton-proton collisions produced by the Large Hadron Collider at a centre-of-mass energy $\sqrt{s}=13$ TeV. Two methods are used to extract the efficiency from $t\bar{t}$ events, a combinatorial likelihood approach and a tag-and-probe method. A boosted decision tree, not using b-tagging information, is used to select events in which two b-jets are present, which reduces the dominant uncertainty in the modelling of the flavour of the jets. The efficiency is extracted for jets in a transverse momentum range from 20 to 300 GeV, with data-to-simulation scale factors calculated by comparing the efficiency measured using collision data to that predicted by the simulation. In conclusion, the two methods give compatible results, and achieve a similar level of precision, measuring data-to-simulation scale factors close to unity with uncertainties ranging from 2% to 12% depending on the jet transverse momentum.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Studies of the beam extraction system of the GTS-LHC electron cyclotron resonance ion source at CERN

The 14.5 GHz GTS-LHC Electron Cyclotron Resonance Ion Source (ECRIS) provides multiply charged heavy ion beams for the CERN experimental program. The GTS-LHC beam formation has been studied extensively with lead, argon, and xenon beams with varied beam extraction conditions using the ion optical code IBSimu. The simulation model predicts self-consistently the formation of triangular and hollow beam structures which are often associated with ECRIS ion beams, as well as beam loss patterns which match the observed beam induced markings in the extraction region. These studies provide a better understanding of the properties of the extracted beams and a way to diagnose the extraction system performance and limitations, which is otherwise challenging due to the lack of direct diagnostics in this region and the limited availability of the ion source for development work.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Mechanically and electrically tunable Rashba-Edelstein effect in ferroelectric semiconductors, CsGe⁢𝑋 3 (𝑋 = I, Br, Cl)

The ability of materials to convert charge current into spin current is fundamental to many spintronics applications. One means of realizing this conversion is via Rashba-Edelstein effect (REE). Using density functional theory simulations, we predict that REE can be induced in the recently discovered family of semiconducting ferroelectrics, CsGe⁢𝑋 3 (𝑋 = I, Br, Cl). The effect is quantified through Rashba-Edelstein coefficients, 𝜒 𝑥⁢𝑦 =−𝜒 𝑦⁢𝑥 , which are nonzero in valence and conduction bands. The largest values, obtained for CsGeI 3 , are 3.45 × 10 10 $\frac{ℏ}{Λ cm}$ and 0.97 × 10 10 $\frac{ℏ}{Λ cm}$ in the conduction and valence bands, respectively. The values are comparable to, and sometimes exceed, those of other inorganic materials, although the maximal values occur away from the band edges. The coefficients' sign couples to the direction of spontaneous polarization, offering opportunities for nonvolatile spin current manipulation via external electric field. Furthermore, these coefficients are highly tunable through strain engineering owing to strain-induced variations in spin textures and energy isosurfaces. Specifically, 𝜒 in the valence band of CsGeI 3 is enhanced to 3.61 × 10 10 $\frac{ℏ}{Λ cm}$ under 5% biaxial strain. Given this potential, we synthesized CsGeI 3 to validate our simulation structure and found excellent agreement between experiment and simulation, thereby allowing extrapolation of our findings to practically significant temperatures. Furthermore, our study identifies promising materials for nonvolatile, multifunctional spintronic applications.

Electronic structure↗
DOE OSTI2015

On the [C ii]–SFR relation in high redshift galaxies

After two Atacama Large Millimeter/submillimeter Array (ALMA) observing cycles, only a handful of [C ii] 158 μm emission line searches in z > 6 galaxies have reported a positive detection, questioning the applicability of the local [C ii]–star formation rate (SFR) relation to high-z systems. To investigate this issue we use the Vallini et al. (V13) model,based on high-resolution, radiative transfer cosmological simulations to predict the [C ii] emission from the interstellar medium of a z ≈ 7 (halo mass M{sub h} = 1.17 × 10{sup 11} M{sub ⊙}) galaxy. We improve the V13 model by including (a) a physically motivated metallicity (Z) distribution of the gas, (b) the contribution of photodissociation regions (PDRs), and (c) the effects of cosmic microwave background (CMB) on the [C ii] line luminosity. We study the relative contribution of diffuse neutral gas to the total [C ii] emission (F{sub diff}/F{sub tot}) for different SFR and Z values. We find that the [C ii] emission arises predominantly from PDRs: regardless of the galaxy properties, F{sub diff}/F{sub tot} ≤ 10%, since at these early epochs the CMB temperature approaches the spin temperature of the [C ii] transition in the cold neutral medium (T{sub CMB} ∼ T{sub s}{sup CNM} ∼ 20 K). Our model predicts a high-z [C ii]–SFR relation, consistent with observations of local dwarf galaxies (0.02 < Z/Z{sub ⊙} < 0.5). The [C ii] deficit suggested by actual data (L{sub Cii} < 2.0 × 10{sup 7} L{sub ⊙} in BDF3299 at z ≈ 7.1) if confirmed by deeper ALMA observations, can be ascribed to negative stellar feedback disrupting molecular clouds around star formation sites. The deviation from the local [C ii]–SFR would then imply a modified Kennicutt–Schmidt relation in z > 6 galaxies. Alternatively/in addition, the deficit might be explained by low gas metallicities (Z < 0.1 Z{sub ⊙}).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Stars get dizzy after lunch

Exoplanet searches have discovered a large number of {sup h}ot Jupiters{sup —}high-mass planets orbiting very close to their parent stars in nearly circular orbits. A number of these planets are sufficiently massive and close-in to be significantly affected by tidal dissipation in the parent star, to a degree parameterized by the tidal quality factor Q {sub *}. This process speeds up their star's rotation rate while reducing the planet's semimajor axis. In this paper, we investigate the tidal destruction of hot Jupiters. Because the orbital angular momenta of these planets are a significant fraction of their star's rotational angular momenta, they spin up their stars significantly while spiraling to their deaths. Using the Monte Carlo simulation, we predict that for Q {sub *} = 10{sup 6}, 3.9 × 10{sup –6} of stars with the Kepler Target Catalog's mass distribution should have a rotation period shorter than 1/3 day (8 hr) due to accreting a planet. Exoplanet surveys such as SuperWASP, HATnet, HATsouth, and KELT have already produced light curves of millions of stars. These two facts suggest that it may be possible to search for tidally destroyed planets by looking for stars with extremely short rotational periods, then looking for remnant planet cores around those candidates, anomalies in the metal distribution, or other signatures of the recent accretion of the planet.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Phase Switching as the Origin of Large Piezoelectric Response in Organic-Inorganic Perovskites: A First-Principles Study

Piezoelectrics are critical functional components of many practical applications such as sensors, ultrasonic transducers, actuators, medical imaging, telecommunications. So far, the best performing piezoelectrics are ferroelectric ceramics, many of which are toxic, heavy, hard, and cost-ineffective. Recently, a groundbreaking discovery of extraordinarily large piezoelectric coefficients in the family of organic-inorganic perovskites gave a hope for a cheaper, environmentally friendly, inexpensive, light-weight and flexible alternative. However, the origin of such response in organic-inorganic ferroelectrics whose spontaneous polarization is an order of magnitude smaller than for inorganic counterparts remains unclear. In this study we employ first-principles simulations to predict that the mechanism associated with large piezoelectric constants is of extrinsic origin and associated with switching between the stable phase and previously overlooked energetically competitive metastable phase that can be stabilized by the external stress. Here, the phase switching changes the polarization direction and, therefore, gives origin to large piezoelectric response, similar to PbZr$_{1-x}$Ti$_x$O$_3$ near morphotropic phase boundary. Existence of such metastable phases is likely to manifest as the dynamical molecular disorder above the Curie temperature and, therefore, could be intrinsic to the entire family of organic-inorganic ferroelectrics with such disorder.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

On the feasibility of polyurethane based 3D dosimeters with optical CT for dosimetric verification of low energy photon brachytherapy seeds

Purpose: To investigate the feasibility of and challenges yet to be addressed to measure dose from low energy (effective energy <50 keV) brachytherapy sources (Pd-103, Cs-131, and I-125) using polyurethane based 3D dosimeters with optical CT. Methods: The authors' evaluation used the following sources: models 200 (Pd-103), CS-1 Rev2 (Cs-131), and 6711 (I-125). The authors used the Monte Carlo radiation transport code MCNP5, simulations with the ScanSim optical tomography simulation software, and experimental measurements with PRESAGE{sup ®} dosimeters/optical CT to investigate the following: (1) the water equivalency of conventional (density = 1.065 g/cm{sup 3}) and deformable (density = 1.02 g/cm{sup 3}) formulations of polyurethane dosimeters, (2) the scatter conditions necessary to achieve accurate dosimetry for low energy photon seeds, (3) the change in photon energy spectrum within the dosimeter as a function of distance from the source in order to determine potential energy sensitivity effects, (4) the optimal delivered dose to balance optical transmission (per projection) with signal to noise ratio in the reconstructed dose distribution, and (5) the magnitude and characteristics of artifacts due to the presence of a channel in the dosimeter. Monte Carlo simulations were performed using both conventional and deformable dosimeter formulations. For verification, 2.8 Gy at 1 cm was delivered in 92 h using an I-125 source to a PRESAGE{sup ®} dosimeter with conventional formulation and a central channel with 0.0425 cm radius for source placement. The dose distribution was reconstructed with 0.02 and 0.04 cm{sup 3} voxel size using the Duke midsized optical CT scanner (DMOS). Results: While the conventional formulation overattenuates dose from all three sources compared to water, the current deformable formulation has nearly water equivalent attenuation properties for Cs-131 and I-125, while underattenuating for Pd-103. The energy spectrum of each source is relatively stable within the first 5 cm especially for I-125. The inherent assumption of radial symmetry in the TG43 geometry leads to a linear increase in sample points within the 3D dosimeter as a function of distance away from the source, which partially offsets the decreasing signal. Simulations of dose reconstruction using optical CT showed the feasibility of reconstructing dose out to a radius of 10 cm without saturating projection images using an optimal dose and high dynamic range scanning; the simulations also predicted that reconstruction artifacts at the channel surface due to a small discrepancy in refractive index should be negligible. Agreement of the measured with calculated radial dose function for I-125 was within 5% between 0.3 and 2.5 cm from the source, and the median difference of measured from calculated anisotropy function was within 5% between 0.3 and 2.0 cm from the source. Conclusions: 3D dosimetry using polyurethane dosimeters with optical CT looks to be a promising application to verify dosimetric distributions surrounding low energy brachytherapy sources.

07 ISOTOPE AND RADIATION SOURCES↗