Search NASA⌕ Search

SEARCH · Search NASA

Results for “Material point method”

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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Statistical Study of Primary Knock-on Atoms under Neutron Irradiation in RMC Code

Materials in high flux neutron condition will suffer from an evolution of damage which ultimately restricts the lifetime of the fuel clad and reactor pressure vessel. Atoms in irradiated material transform to primary knock-on atoms (PKAs) by absorbing the kinetic energy transferred from neutrons, and the energy released from nuclear reactions. The secondary collision cascades which are induced by PKAs create point defects within the material, with residual vacancies and interstitials eventually causing swelling and embrittlement that degrades mechanical behavior. For a long time, Displacements per Atom (DPA) index has been identified as a measure of radiation damage in material. Nowadays, several advanced methods, such as Molecular Dynamics (MD) and Kinetic Monte Carlo (KMC) simulations, have been developed for more accurately simulating the evolution of the damage. However, all of these techniques need more comprehensive information of PKA, such as the type, energy, direction and spatial distribution. The limitations of DPA emerge and a method should be provided to simulate the generative process of PKAs. Recently, a deterministic analysis code SPECTRA-PKA has been developed by Gilbert to calculate the energy spectra of almost all types PKAs by solving the group-to-group cross-section matrices equation which includes neutron scattering, light charged particles, and recoils of the residual nuclei. But this deterministic analysis code could not give the direction, spatial and time distribution of PKAs. Based on the Geant4 program, Cui Zhenguo has simulated the irradiation by neutrons with 1 MeV energy in zirconium, while only the two adjacent elastic collisions between injected neutron and target atoms are taken into consideration. In this paper, Monte Carlo method to calculate PKAs information during neutron transport in RMC code is introduced. RMC is a continuous-energy Reactor Monte Carlo neutron and photon transport code being developed by Department of Engineering Physics at Tsinghua University, Beijing. PKAs in zirconium are calculated in fixed neutron source mode. Statistical information of PKAs in PWR fuel cell and SFR fuel cell is shown and compared. The results show that Monte Carlo method can provide more comprehensive information of PKAs. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Phonon modal analysis of thermal transport in ThO 2 with point defects using equilibrium molecular dynamics

Defects can significantly degrade the thermal conductivity of ThO 2 , an advanced nuclear fuel material as well as a surrogate for other fluorite-structured materials. Here, we investigate how point defects in ThO 2 impact phonon mode-resolved thermal transport. By incorporating phonon modes from lattice dynamics, we decompose the trajectory and heat flux to phonon normal mode space and extract key phonon properties, including phonon relaxation times and their contributions to thermal conductivity. We implement two methods. The first method is based on the Green Kubo formalism to resolve the contribution of each phonon mode to thermal conductivity. The second resolves the lifetime of individual phonon modes and the thermal conductivity is calculated using the Boltzmann transport equation within relaxation time approximation. Notably, a lower contribution of acoustic modes is revealed compared to perturbative approaches considering only three-phonon scattering processes. The effects of four types of point defects are evaluated. The strongest impact on a reduction in thermal conductivity is from Th interstitials, followed by Th vacancies. O interstitials/vacancies have a similar impact, albeit smaller than defects on the thorium sublattice. These observations are consistent with previous studies.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

International interlaboratory compilation of trace element concentrations in the CUP-2 uranium ore concentrate standard

Herein, a nuclear forensics investigation involving a uranium ore concentrate relies on accurate and precise analysis of impurities. Analytical data defensibility requires the use of reference materials as part of quality control. This study presents a compilation of trace element concentration results of the CUP-2 Uranium Ore Concentrate Standard measured by 11 different laboratories. The laboratories employed various dissolution methods, analytical preparation methods, and instrumental platforms. The data presented here contain concentrations of 66 impurities with up to 138 individual data points for each impurity. Consensus values have been assigned to each impurity following a statistical analysis of the data set.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2010

Luminescence and photosensitization properties of ensembles of silicon nanocrystals in terms of an exciton migration model

The relaxation processes that occur in ensembles of coupled silicon nanocrystals are described by a quantitative model that takes into account the energy transfer between them during exciton migration. This model is used to study the formation of singlet oxygen during the photoexcitation of silicon nanocrystals in porous silicon layers under various external conditions. It is experimentally found that, upon fine milling of as-deposited porous silicon films, the photoluminescence decay time increases despite an increase in the concentration of point defects. The photosensitized activity of ensembles of silicon nanocrystals degrades monotonically during their photostimulated oxidation. These experimental results agree completely with the conclusions of the model and are explained by a decrease in the number of exciton migration ways between nanocrystals when the granule size of a porous silicon powder decreases and by an increase in the efficiency of nonradiative recombination during the photostimulated oxidation of the nanocrystals. Our data indicate that nanocrystalline silicon is a promising material for the methods of nontoxic photodynamic therapy of oncological diseases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2012

The first 3D malonate bridged copper [Cu(O{sub 2}C-CH{sub 2}-CO{sub 2}H){sub 2}{center_dot}2H{sub 2}O]: Structure, properties and electronic structure

A new inorganic-organic compound [Cu(O{sub 2}C-CH{sub 2}-CO{sub 2}H){sub 2}{center_dot}2H{sub 2}O] ([Cumal]) was hydrothermally synthesized and characterized by IR spectroscopy, thermal analysis and single crystal X-ray diffraction. [Cumal] is the first three-dimensional compound existing in the system Cu(II)-malonic acid-H{sub 2}O. Its framework is built up through carboxyl bridged copper where CuO{sub 6} octahedra are elongated with an almost D{sub 4h} symmetry (4+2) due to the Jahn-Teller effect. The magnetic properties were studied by measuring its magnetic susceptibility in the temperature range of 2-300 K indicating the existence of weak ferromagnetic interactions. The electronic structure of [Cumal] was calculated within the density functional theory (DFT) framework. Structural features are well reproduced using DFT structural optimizations and the optical spectra, calculated within the dielectric formalism, explain very well the light blue colour of the compound. It is shown that a GGA+U approach with a U{sub eff} value of about 6 eV is necessary for a better correlation with the experiment. - Graphical abstract: [Cu(O{sub 2}C-CH{sub 2}-CO{sub 2}H){sub 2}{center_dot}2H{sub 2}O]: the first 3D hybrid organic-inorganic compound built up carboxyl groups. The network presents a diamond-like structure achieved via carboxyl. Highlights: Black-Right-Pointing-Pointer A new organic-inorganic material with an unprecedented topology is synthesized. Black-Right-Pointing-Pointer Crystallographic structure is determined using single crystal X-ray diffraction. Black-Right-Pointing-Pointer Electronic structure is obtained from DFT, GGA+U calculation. Black-Right-Pointing-Pointer Framework can be described as formed from CuC{sub 4} tetrahedron sharing four corners. Black-Right-Pointing-Pointer This structure can be classified as an extended diamond structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Anomalous properties and the liquid-liquid phase transition in gallium

A group of materials including water and silicon exhibit many anomalous behaviors, e.g., density anomaly and diffusivity anomaly (increase upon compression). These materials are hypothesized to have a liquid-liquid phase transition (LLPT) and the critical fluctuation in the vicinity of the liquid-liquid critical point is considered as the origin of different anomalies. Liquid gallium was also reported to have a LLPT, yet whether it shows similar water-like anomalies is not yet studied. Using molecular dynamics simulations on a modified embedded-atom model, we study the thermodynamic, dynamic, and structural properties of liquid gallium as well as its LLPT. We find that, similar to water-like materials predicted to have the LLPT, gallium also shows different anomalous behaviors (e.g., density anomaly, diffusivity anomaly, and structural anomaly). We also find that its thermodynamic and structural response functions are continuous and show maxima in the supercritical region, the loci of which asymptotically approach to the other and merge to the Widom line. These phenomena are consistent with the supercritical phenomenon in a category of materials with a liquid-liquid critical point, which could be common features in most materials with a LLPT.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Effect of annealing on the laser induced damage of polished and CO{sub 2} laser-processed fused silica surfaces

We investigate the effect of different heat treatments on the laser-induced damage probabilities of fused silica samples. Isothermal annealing in a furnace is applied, with different temperatures in the range 700–1100 °C and 12 h annealing time, to super-polished fused silica samples. The surface flatness and laser damage probabilities at 3 ns, 351 nm are measured before and after the different annealing procedures. We have found a significant improvement of the initial laser damage probabilities of the silica surface after annealing at 1050 °C for 12 h. A similar study has been conducted on CO{sub 2} laser-processed sites on the surface of the samples. Before and after annealing, we have studied the morphology of the sites, the evolution of residual stress, and the laser-induced damage threshold measured at 351 nm, 3 ns. In this case, we observe that the laser damage resistance of the laser created craters can reach the damage level of the bare fused silica surface after the annealing process, with a complete stress relieve. The obtained results are then compared to the case of local annealing process by CO{sub 2} laser irradiation during 1 s, and we found similar improvements in both cases. The different results obtained in the study are compared to numerical simulations made with a thermo-mechanical model based on finite-element method that allows the simulation of the isothermal or the local annealing process, the evolution of stress and fictive temperature. The simulation results were found to be very consistent with experimental observations for the stresses evolution after annealing and estimation of the heat affected area during laser-processing based on the density dependence with fictive temperature. Following this work, the temperature for local annealing should reach 1330–1470 °C for an optimized reduction of damage probability and be below the threshold for material removal, whereas furnace annealing should be kept below the annealing point to avoid sample deformation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Atomic defects in monolayer titanium carbide (Ti 3 C 2 T x ) MXene

Here, the 2D transition metal carbides or nitrides, or MXenes, are emerging as a group of materials showing great promise in lithium ion batteries and supercapacitors. Until now, characterization and properties of single-layer MXenes have been scarcely reported. Here, using scanning transmission electron microscopy, we determined the atomic structure of freestanding monolayer Ti 3 C 2 T x flakes prepared via the minimally intensive layer delamination method and characterized different point defects that are prevalent in the monolayer flakes. We determine that the Ti vacancy concentration can be controlled by the etchant concentration during preparation. Density function theory-based calculations confirm the defect structures and predict that the defects can influence the surface morphology and termination groups, but do not strongly influence the metallic conductivity. Using devices fabricated from single- and few-layer Ti 3 C 2 T x MXene flakes, the effect of the number of layers in the flake on conductivity has been demonstrated.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Quantum effects on dislocation motion from ring-polymer molecular dynamics

Quantum motion of atoms known as zero-point vibration was recently proposed to explain a long-standing discrepancy between theoretically computed and experimentally measured low-temperature plastic strength of iron and possibly other metals with high atomic masses. This finding challenges the traditional notion that quantum motion of atoms is relatively unimportant in solids comprised of heavy atoms. Here we report quantum dynamic simulations of quantum effects on dislocation motion within the exact formalism of Ring-Polymer Molecular Dynamics (RPMD). To extend the reach of quantum atomistic simulations to length and time scales relevant for extended defects in materials, we implemented RPMD in the open-source code LAMMPS thus making the RPMD method widely available to the community. We use our RPMD/LAMMPS approach for direct calculations of dislocation mobility and its effects on the yield strength of α-iron. Our simulation results establish that quantum effects are noticeable at temperatures below 50 K but account for only a modest (≈13% at T = 0 K) overall reduction in the Peierls barrier, at variance with the factor of two reduction predicted earlier based on the more approximate framework of harmonic transition state theory. Our results confirm that zero-point vibrations provide ample additional agitation for atomic motion that increases with decreasing temperature, however its enhancing effect on dislocation mobility is largely offset by an increase in the effective atom size, an effect known as quantum dispersion that has not been accounted for in the previous calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Development of a phantom to validate high-dose-rate brachytherapy treatment planning systems with heterogeneous algorithms

Purpose: This work presents the development of a phantom to verify the treatment planning system (TPS) algorithms used for high-dose-rate (HDR) brachytherapy. It is designed to measure the relative dose in a heterogeneous media. The experimental details used, simulation methods, and comparisons with a commercial TPS are also provided. Methods: To simulate heterogeneous conditions, four materials were used: Virtual Water™ (VM), BR50/50™, cork, and aluminum. The materials were arranged in 11 heterogeneity configurations. Three dosimeters were used to measure the relative response from a HDR {sup 192}Ir source: TLD-100™, Gafchromic{sup ®} EBT3 film, and an Exradin™ A1SL ionization chamber. To compare the results from the experimental measurements, the various configurations were modeled in the PENELOPE/penEasy Monte Carlo code. Images of each setup geometry were acquired from a CT scanner and imported into BrachyVision™ TPS software, which includes a grid-based Boltzmann solver Acuros™. The results of the measurements performed in the heterogeneous setups were normalized to the dose values measured in the homogeneous Virtual Water™ setup and the respective differences due to the heterogeneities were considered. Additionally, dose values calculated based on the American Association of Physicists in Medicine-Task Group 43 formalism were compared to dose values calculated with the Acuros™ algorithm in the phantom. Calculated doses were compared at the same points, where measurements have been performed. Results: Differences in the relative response as high as 11.5% were found from the homogeneous setup when the heterogeneous materials were inserted into the experimental phantom. The aluminum and cork materials produced larger differences than the plastic materials, with the BR50/50™ material producing results similar to the Virtual Water™ results. Our experimental methods agree with the PENELOPE/penEasy simulations for most setups and dosimeters. The TPS relative differences with the Acuros™ algorithm were similar in both experimental and simulated setups. The discrepancy between the BrachyVision™, Acuros™, and TG-43 dose responses in the phantom described by this work exceeded 12% for certain setups. Conclusions: The results derived from the phantom measurements show good agreement with the simulations and TPS calculations, using Acuros™ algorithm. Differences in the dose responses were evident in the experimental results when heterogeneous materials were introduced. These measurements prove the usefulness of the heterogeneous phantom for verification of HDR treatment planning systems based on model-based dose calculation algorithms.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2012

Dynamic Contrast-Enhanced MRI of Cervical Cancers: Temporal Percentile Screening of Contrast Enhancement Identifies Parameters for Prediction of Chemoradioresistance

Purpose: To systematically screen the tumor contrast enhancement of locally advanced cervical cancers to assess the prognostic value of two descriptive parameters derived from dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI). Methods and Materials: This study included a prospectively collected cohort of 81 patients who underwent DCE-MRI with gadopentetate dimeglumine before chemoradiotherapy. The following descriptive DCE-MRI parameters were extracted voxel by voxel and presented as histograms for each time point in the dynamic series: normalized relative signal increase (nRSI) and normalized area under the curve (nAUC). The first to 100th percentiles of the histograms were included in a log-rank survival test, resulting in p value and relative risk maps of all percentile-time intervals for each DCE-MRI parameter. The maps were used to evaluate the robustness of the individual percentile-time pairs and to construct prognostic parameters. Clinical endpoints were locoregional control and progression-free survival. The study was approved by the institutional ethics committee. Results: The p value maps of nRSI and nAUC showed a large continuous region of percentile-time pairs that were significantly associated with locoregional control (p < 0.05). These parameters had prognostic impact independent of tumor stage, volume, and lymph node status on multivariate analysis. Only a small percentile-time interval of nRSI was associated with progression-free survival. Conclusions: The percentile-time screening identified DCE-MRI parameters that predict long-term locoregional control after chemoradiotherapy of cervical cancer.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

Power optics

By using the theory we developed in the early 1970s, a broad range of phenomena is considered for an optical surface of a solid body that is exposed to radiation arbitrarily varying in time and producing temperature fields, thermoelastic stresses and thermal deformations on the surface layer. The examination is based on the relations (which are similar to Duhamel's integral formula from the theory of heat conduction) between the quantities characterising the thermal stress state in any nonstationary regimes of energy input into a solid. A peculiar feature of the analysis of the thermal stress state in this case consists in the fact that this relation comprises time as a parameter, which in turn is a consequence of incoherence of the quasi-stationary problem of thermoelasticity. This phenomenon is particularly important for the optics of high-power, high-pulse repetition rate lasers, which are being actively developed. In the review, we have recently published in Laser Physics, the thermal stress state of a solid is analysed. In this state, time is treated as an independent variable used in differentiation. Such an approach greatly reduces the applicability of the method. The review published contains data on the use of capillary porous structures made of various materials with different degrees of the surface development. Moreover, such structures can be efficiently employed to increase the heat exchange at a temperature below the boiling point of the coolant. In the present review we discuss the dependences of the limiting laser intensities on the duration of a pulse or a pulse train, corresponding to the three stages of the state of the reflecting surface and leading to unacceptable elastic deformations of the surface, to the plastic yield of the material accompanied by the formation of residual stresses and to the melting of the surface layer. We also analyse the problem of heat exchange in the surface layer with a liquid metal coolant pumped through it. The theoretical estimates are compared with the experimental data. We discuss the issues related to the technology of fabrication of power optics elements based on materials with a porous structure, of lightweight highly stable large optics based on highly porous materials, multi-layer honeycomb structures and silicon carbide, as well as problems of application of physical and technical fundamentals of power optics in modern cutting-edge technology. (invited paper)

36 MATERIALS SCIENCE↗
DOE OSTI2024

Extraction of Vibration Data with Imaging

To date, the primary sensing technology used to measure the vibration response has been accelerometers and strain gages mounted directly to the structure and using either wired or, more recently, wireless telemetry. Cost issues with these sensors and the associated data acquisition systems typically limit the numbers that are deployed on in situ structures. Although there are a few structures with larger sensing counts that in some cases exceed over 1000 sensors, more typical numbers range from ten to one hundred sensors resulting in low spatial resolution when they are applied to physically large systems. When one considers that nuclear power plant structures usually have complex geometries, material properties, connectivity and boundary conditions, it is clear these current approaches to vibration measurements can only provide limited information about a system’s dynamics response characteristics. As an alternative, many non-contact measurement technologies have emerged, including point wise measurement methods such as Global Positioning System (GPS), microwave interferometry, and laser Doppler vibrometry (LDV), as well as simultaneous full-field measurement methods such as electronic speckle pattern interferometry, holography interferometry, and muon tomography, some of which can provide high spatial resolution measurements. Among these methods, digital video imaging techniques have emerged as a feasible solution for full-field vibration measurements that provide significantly more detailed dynamic response information because every pixel becomes a measurement point. Furthermore, recent advances in image processing and computer vision algorithms have been successfully used to process video data for experimental and operational modal analysis. Such full-field measurements have the potential to significantly improve many current structural assessment procedures including system identification (modal parameter estimation), structural health monitoring, load reconstruction, model validation, and model updating. Furthermore, more recent full-field imaging techniques can be accomplished with relatively low-cost, commercially-available off-the-shelf cameras. However, these measurement procedures have other limitations that must be considered such as the ability to only measure visibly accessible points on a structure and a more limited dynamic range and bandwidth than can be achieved with accelerometers or strain gages.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1978

Interpolation and extrapolation of creep rupture data by the minimum commitment method. III - Analysis of multiheats

An outline is presented of approaches for treating multiheats on the basis of the focal point convergence method. The method has been employed in the case of two highly characterized multiheats, including a 304 stainless steel studied in Japan and a low alloy carbon steel studied in England. The method makes use of the same functional form for all materials. Only the constants are varied for each multiheat. Completely computerized procedures are employed for the determination of the constants. Once the basic analysis has been performed, the representation of various members in the same system is achieved by adding linear expressions of log stress, changing only two constants in the equations to represent a selected heat.

Manson, S. S.↗
NASA NTRS2000

Microgravity Disturbance Characterization of the Quench Module Insert (QMI) Phase Change Device (PCD)

The Materials Science Research Facility (MSRF) is a multi-user, multi-purpose facility for materials science research. One experiment within the MSRF will be the Quench Module Insert (QMI), a high-temperature furnace with unique capabilities for processing different classes of materials. The primary functions of the QMI furnace are to melt, directionally solidify, and quench metallic samples, providing data to aid in understanding the effects of the microgravity environment on the characteristics of these processed metals. The QMI houses sealed individual sample ampoules containing material to be processed. Quenching of the samples in the QMI furnace is accomplished by releasing low-melting-point metallic shoes into contact with the outside of the sample ampoule, dissipating heat and cooling the sample inside. The impact from this method of quench will induce sample vibrations which could be large enough to adversely affect sample quality. Utilizing breadboard hardware, the sample quench sequence, releasing the shoes, was conducted. Data was collected from accelerometers located on the breadboard sample cartridge, indicating the maximum acceleration achieved by the sample. The primary objective of the test described in this presentation was to determine the acceleration imparted on the sample by the shoe contact. From this information, the science community can better assess whether this method of quench will allow them to obtain the data they need.

Gattis, Christy↗
NASA NTRS1999

Analytical Prediction of Damage Growth in Notched Composite Panels Loaded in Axial Compression

A progressive failure analysis method based on shell elements is developed for the computation of damage initiation and growth in stiffened thick-skin stitched graphite-epoxy panels loaded in axial compression. The analysis method involves a step-by-step simulation of material degradation based on ply-level failure mechanisms. High computational efficiency is derived from the use of superposed layers of shell elements to model each ply orientation in the laminate. Multiple integration points through the thickness are used to obtain the correct bending effects through the thickness without the need for ply-by-ply evaluations of the state of the material. The analysis results are compared with experimental results for three stiffened panels with notches oriented at 0, 15 and 30 degrees to the panel width dimension. A parametric study is performed to investigate the damage growth retardation characteristics of the Kevlar stitch lines in the pan

Ambur, Damodar R.↗
NASA NTRS2010

Corrosion-Resistant Container for Molten-Material Processing

In a carbothermal process, gaseous methane is passed over molten regolith, which is heated past its melting point to a temperature in excess of 1,625 C. At this temperature, materials in contact with the molten regolith (or regolith simulant) corrode and lose their structural properties. As a result, fabricating a crucible to hold the molten material and providing a method of contact heating have been problematic. Alternative containment approaches use a large crucible and limit the heat zone of the material being processed, which is inefficient because of volume and mass constraints. Alternative heating approaches use non-contact heating, such as by laser or concentrated solar energy, which can be inefficient in transferring heat and thus require higher power heat sources to accomplish processing. The innovation is a combination of materials, with a substrate material having high structural strength and stiffness and high-temperature capability, and a coating material with a high corrosion resistance and high-temperature capability. The material developed is a molybdenum substrate with an iridium coating. Creating the containment crucible or heater jacket using this material combination requires only that the molybdenum, which is easily processed by conventional methods such as milling, electric discharge machining, or forming and brazing, be fabricated into an appropriate shape, and that the iridium coating be applied to any surfaces that may come in contact with the corrosive molten material. In one engineering application, the molybdenum was fashioned into a container for a heat pipe. Since only the end of the heat pipe is used to heat the regolith, the container has a narrowing end with a nipple in which the heat pipe is snugly fit, and the external area of this nipple, which contacts the regolith to transfer heat into it, is coated with iridium. At the time of this reporting, no single material has been found that can perform the functions of this combination of materials, and other combinations of materials have not proven to be survivable to the corrosiveness of this environment. High-temperature processing of materials with similar constituencies as lunar regolith is fairly common. The carbo-thermal process is commonly used to make metallurgical-grade silicon for the semiconductor and solar-cell industries.

Stern, Theodore G.↗
DOE OSTI2016

Evaluation of a commercial MRI Linac based Monte Carlo dose calculation algorithm with GEANT 4

Purpose: This paper provides a comparison between a fast, commercial, in-patient Monte Carlo dose calculation algorithm (GPUMCD) and GEANT4. It also evaluates the dosimetric impact of the application of an external 1.5 T magnetic field. Methods: A stand-alone version of the Elekta™ GPUMCD algorithm, to be used within the Monaco treatment planning system to model dose for the Elekta™ magnetic resonance imaging (MRI) Linac, was compared against GEANT4 (v10.1). This was done in the presence or absence of a 1.5 T static magnetic field directed orthogonally to the radiation beam axis. Phantoms with material compositions of water, ICRU lung, ICRU compact-bone, and titanium were used for this purpose. Beams with 2 MeV monoenergetic photons as well as a 7 MV histogrammed spectrum representing the MRI Linac spectrum were emitted from a point source using a nominal source-to-surface distance of 142.5 cm. Field sizes ranged from 1.5 × 1.5 to 10 × 10 cm{sup 2}. Dose scoring was performed using a 3D grid comprising 1 mm{sup 3} voxels. The production thresholds were equivalent for both codes. Results were analyzed based upon a voxel by voxel dose difference between the two codes and also using a volumetric gamma analysis. Results: Comparisons were drawn from central axis depth doses, cross beam profiles, and isodose contours. Both in the presence and absence of a 1.5 T static magnetic field the relative differences in doses scored along the beam central axis were less than 1% for the homogeneous water phantom and all results matched within a maximum of ±2% for heterogeneous phantoms. Volumetric gamma analysis indicated that more than 99% of the examined volume passed gamma criteria of 2%—2 mm (dose difference and distance to agreement, respectively). These criteria were chosen because the minimum primary statistical uncertainty in dose scoring voxels was 0.5%. The presence of the magnetic field affects the dose at the interface depending upon the density of the material on either sides of the interface. This effect varies with the field size. For example, at the water-lung interface a 33.94% increase in dose was observed (relative to the D{sub max}), by both GPUMCD and GEANT4 for the field size of 2 × 2 cm{sup 2} (compared to no B-field case), which increased to 47.83% for the field size of 5 × 5 cm{sup 2} in the presence of the magnetic field. Similarly, at the lung-water interface, the dose decreased by 19.21% (relative to D{sub max}) for a field size of 2 × 2 cm{sup 2} and by 30.01% for 5 × 5 cm{sup 2} field size. For more complex combinations of materials the dose deposition also becomes more complex. Conclusions: The GPUMCD algorithm showed good agreement against GEANT4 both in the presence and absence of a 1.5 T external magnetic field. The application of 1.5 T magnetic field significantly alters the dose at the interfaces by either increasing or decreasing the dose depending upon the density of the material on either side of the interfaces.

61 RADIATION PROTECTION AND DOSIMETRY↗