Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Calculations”

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,423 records · Page 79

SU-F-T-125: Radial Dose Distributions From Carbon Ions of Therapeutic Energies Calculated with Geant4-DNA

Purpose: Radial dose distribution D(r) is the dose as a function of lateral distance from the path of a heavy charged particle. Its main application is in modelling of biological effects of heavy ions, including applications to hadron therapy. It is the main physical parameter of a broad group of radiobiological models known as the amorphous track models. Our purpose was to calculate D(r) with Monte Carlo for carbon ions of therapeutic energies, find a simple formula for D(r) and fit it to the Monte Carlo data. Methods: All calculations were performed with Geant4-DNA code, for carbon ion energies from 10 to 400 MeV/u (ranges in water: ∼ 0.4 mm to 27 cm). The spatial resolution of dose distribution in the lateral direction was 1 nm. Electron tracking cut off energy was 11 eV (ionization threshold). The maximum lateral distance considered was 10 µm. Over this distance, D(r) decreases with distance by eight orders of magnitude. Results: All calculated radial dose distributions had a similar shape dominated by the well-known inverse square dependence on the distance. Deviations from the inverse square law were observed close to the beam path (r<10 nm) and at large distances (r >1 µm). At small and large distances D(r) decreased, respectively, slower and faster than the inverse square of distance. A formula for D(r) consistent with this behavior was found and fitted to the Monte Carlo data. The accuracy of the fit was better than 10% for all distances considered. Conclusion: We have generated a set of radial dose distributions for carbon ions that covers the entire range of therapeutic energies, for distances from the ion path of up to 10 µm. The latter distance is sufficient for most applications because dose beyond 10 µm is extremely low.

60 APPLIED LIFE SCIENCES↗

SU-F-T-152: Experimental Validation and Calculation Benchmark for a Commercial Monte Carlo Pencil BeamScanning Proton Therapy Treatment Planning System in Heterogeneous Media

Purpose: Eclipse AcurosPT 13.7, the first commercial Monte Carlo pencil beam scanning (PBS) proton therapy treatment planning system (TPS), was experimentally validated for an IBA dedicated PBS nozzle in the CIRS 002LFC thoracic phantom. Methods: A two-stage procedure involving the use of TOPAS 1.3 simulations was performed. First, Geant4-based TOPAS simulations in this phantom were experimentally validated for single and multi-spot profiles at several depths for 100, 115, 150, 180, 210 and 225 MeV proton beams, using the combination of a Lynx scintillation detector and a MatriXXPT ionization chamber array. Second, benchmark calculations were performed with both AcurosPT and TOPAS in a phantom identical to the CIRS 002LFC, with the exception that the CIRS bone/mediastinum/lung tissues were replaced with similar tissues that are predefined in AcurosPT (a limitation of this system which necessitates the two stage procedure). Results: Spot sigmas measured in tissue were in agreement within 0.2 mm of TOPAS simulation for all six energies, while AcurosPT was consistently found to have larger spot sigma (<0.7 mm) than TOPAS. Using absolute dose calibration by MatriXXPT, the agreements between profiles measurements and TOPAS simulation, and calculation benchmarks are over 97% except near the end of range using 2 mm/2% gamma criteria. Overdosing and underdosing were observed at the low and high density side of tissue interfaces, respectively, and these increased with increasing depth and decreasing energy. Near the mediastinum/lung interface, the magnitude can exceed 5 mm/10%. Furthermore, we observed >5% quenching effect in the conversion of Lynx measurements to dose. Conclusion: We recommend the use of an ionization chamber array in combination with the scintillation detector to measure absolute dose and relative PBS spot characteristics. We also recommend the use of an independent Monte Carlo calculation benchmark for the commissioning of a commercial TPS. Partially supported by Varian Medical System under the master agreement between Varian and University of pennsylvania.

60 APPLIED LIFE SCIENCES↗

SU-F-T-686: Considerations About Dose Protraction Factor in TCP Calculations for Prostate VMAT Treatments

Purpose: Dose protraction factor should be considered in order to model the TCP calculations. Nevertheless, this study describes a brief discussion showing that the lack of its inclusion should not invalidate these calculations for prostate VMAT treatments. Methods: Dose protraction factor (G) modifies the quadratic term of the linear-quadratic expression in order to take into account the sublethal damage repair of protracting the dose delivery. If the delivery takes a short time (instantaneous), G = 1. For any other dose delivery pattern, G < 1. The Lea-Catcheside dose protraction factor for external beam radiotherapy contains terms depending of on the tissue specific repair parameter (λ) and the irradiation time (T). Expanding the exponential term using a Taylor’s series and neglecting terms of order (λT){sup 3}, the approximation leads to G = 1. The described situation occurs for 3DCRT techniques, where treatment times are about few minutes. For IMRT techniques, fraction times are prolonged compared to 3DCRT times. Wang et al. (2003) and Fowler et al. (2004) investigated the protraction effect with respect to IMRT treatments, reporting clinically significant loss in biological effect associated with IMRT delivery times. Results: Treatment times are noticeably reduced for prostate treatments using VMAT techniques. These times are comparable to 3DCRT times, leading to consider the previous approximation. Conclusion: Dose protraction factor can be approximated by G = 1 in TCP calculations for prostate treatments using VMAT techniques.

60 APPLIED LIFE SCIENCES↗

SU-G-IeP2-03: Comparison of Dose Calculation On MAR (metal Artifact Reduction) and Non-MAR Datasets for Pelvic Patients with Hip Prosthesis and Head and Neck Patients with Dental Filling

Purpose: Metal artifact reduction (MAR) software in computed tomography (CT) was previously evaluated with phantoms demonstrating the algorithm is capable of reducing metal artifacts without affecting the overall image quality. The goal of this study is to determine the dosimetric impact when calculating with CT datasets reconstructed with and without MAR software. Methods: Twelve head and neck cancer patients with dental fillings and four pelvic cancer patients with hip prosthesis were scanned with a GE Optima RT 580 CT scanner. Images were reconstructed with and without the MAR software. 6MV IMRT and VMAT plans were calculated with AAA on the MAR dataset until all constraints met our clinic’s guidelines. Contours from the MAR dataset were copied to the non-MAR dataset. Next, dose calculation on the non-MAR dataset was performed using the same field arrangements and fluence as the MAR plan. Conformality index, D99% and V100% to PTV were compared between MAR and non-MAR plans. Results: Differences between MAR and non-MAR plans were evaluated. For head and neck plans, the largest variations in conformality index, D99% and V100% were −3.8%, −0.9% and −2.1% respectively whereas for pelvic plans, the biggest discrepancies were −32.7%, −0.4% and -33.5% respectively. The dosimetric impact from hip prosthesis is greater because it produces more artifacts compared to dental fillings. Coverage to PTV can increase or decrease depending on the artifacts since dark streaks reduce the HU whereas bright streaks increase the HU. In the majority of the cases, PTV dose in the non-MAR plans is higher than MAR plans. Conclusion: With the presence of metals, MAR algorithm can allow more accurate delineation of targets and OARs. Dose difference between MAR and non-MAR plans depends on the proximity of the organ to the high density material, the streaking artifacts and the beam arrangements of the plan.

60 APPLIED LIFE SCIENCES↗

SU-G-IeP2-04: Dosimetric Accuracy of a Monte Carlo-Based Tool for Cone-Beam CT Organ Dose Calculation: Validation Against OSL and XRQA2 Film Measurements

Purpose: The intensive use of Cone-Beam Computed Tomography (CBCT) during radiotherapy treatments raise some questions about the dose to healthy tissues delivered during image acquisitions. We hence developed a Monte Carlo (MC)-based tool to predict doses to organs delivered by the Elekta XVI kV-CBCT. This work aims at assessing the dosimetric accuracy of the MC tool, in all tissue types. Methods: The kV-CBCT MC model was developed using the PENELOPE code. The beam properties were validated against measured lateral and depth dose profiles in water, and energy spectra measured with a CdTe detector. The CBCT simulator accuracy then required verification in clinical conditions. For this, we compared calculated and experimental dose values obtained with OSL nanoDots and XRQA2 films inserted in CIRS anthropomorphic phantoms (male, female, and 5-year old child). Measurements were performed at different locations, including bone and lung structures, and for several acquisition protocols: lung, head-and-neck, and pelvis. OSLs and film measurements were corrected when possible for energy dependence, by taking into account for spectral variations between calibration and measurement conditions. Results: Comparisons between measured and MC dose values are summarized in table 1. A mean difference of 8.6% was achieved for OSLs when the energy correction was applied, and 89.3% of the 84 dose points were within uncertainty intervals, including those in bones and lungs. Results with XRQA2 are not as good, because incomplete information about electronic equilibrium in film layers hampered the application of a simple energy correction procedure. Furthermore, measured and calculated doses (Fig.1) are in agreement with the literature. Conclusion: The MC-based tool developed was validated with an extensive set of measurements, and enables the organ dose calculation with accuracy. It can now be used to compute and report doses to organs for clinical cases, and also to drive strategies to optimize imaging protocols.

60 APPLIED LIFE SCIENCES↗

TH-AB-BRA-07: PENELOPE-Based GPU-Accelerated Dose Calculation System Applied to MRI-Guided Radiation Therapy

Purpose: The clinical commissioning of IMRT subject to a magnetic field is challenging. The purpose of this work is to develop a GPU-accelerated Monte Carlo dose calculation platform based on PENELOPE and then use the platform to validate a vendor-provided MRIdian head model toward quality assurance of clinical IMRT treatment plans subject to a 0.35 T magnetic field. Methods: We first translated PENELOPE from FORTRAN to C++ and validated that the translation produced equivalent results. Then we adapted the C++ code to CUDA in a workflow optimized for GPU architecture. We expanded upon the original code to include voxelized transport boosted by Woodcock tracking, faster electron/positron propagation in a magnetic field, and several features that make gPENELOPE highly user-friendly. Moreover, we incorporated the vendor-provided MRIdian head model into the code. We performed a set of experimental measurements on MRIdian to examine the accuracy of both the head model and gPENELOPE, and then applied gPENELOPE toward independent validation of patient doses calculated by MRIdian’s KMC. Results: We achieve an average acceleration factor of 152 compared to the original single-thread FORTRAN implementation with the original accuracy preserved. For 16 treatment plans including stomach (4), lung (2), liver (3), adrenal gland (2), pancreas (2), spleen (1), mediastinum (1) and breast (1), the MRIdian dose calculation engine agrees with gPENELOPE with a mean gamma passing rate of 99.1% ± 0.6% (2%/2 mm). Conclusions: We developed a Monte Carlo simulation platform based on a GPU-accelerated version of PENELOPE. We validated that both the vendor provided head model and fast Monte Carlo engine used by the MRIdian system are accurate in modeling radiation transport in a patient using 2%/2 mm gamma criteria. Future applications of this platform will include dose validation and accumulation, IMRT optimization, and dosimetry system modeling for next generation MR-IGRT systems.

60 APPLIED LIFE SCIENCES↗

Calculation of the Physicochemical Characteristics of a New Orthorhombic Form of Diamond

A geometrically optimized structure and some properties of a new diamond polymorph have been assessed using density-functional theory calculations in the local-density approximation. The new phase has a primitive orthorhombic lattice (Pmma). The calculation results indicate that the density, sublimation energy, bulk modulus, Vickers hardness, and band gap of this phase are lower than those of cubic diamond by 5–18%. The orthorhombic diamond polymorph can be obtained by compressing graphite in the [001] direction at a pressure of 59.4 GPa. For this phase transition, we calculated theoretical X-ray diffraction patterns of the new phase and intermediate states, which can be used for their identification in experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Calculation of Chemical Shifts of X-Ray-Emission Spectra of Niobium in Niobium(V) Oxides Relative to Metal

Chemical shifts of the K{sub α1} and K{sub β1} lines of X-ray-emission spectra of niobium in oxides (Nb{sub 2}O{sub 5}){sub n}, n = 1–4, relative to metal Nb have been calculated. Stoichiometric clusters (Nb2O5)n the electronic structure of which was calculated using accurate relativistic pseudopotentials and two-component version of the density functional theory are considered as prototypes for modeling different crystal forms of niobium(V) oxide. The chemical shifts were calculated using the method based on using the property of approximate proportionality of valence spinors in the core region of the heavy atom [11]. Corrections to values of chemical shifts have been determined with allowance for deviations from the abovementioned proportionality. Rapid convergence of results with respect to the size of the niobium oxide cluster has been demonstrated.

74 ATOMIC AND MOLECULAR PHYSICS↗

Precise calculations in simulations of the interaction of low energy neutrons with nano-dispersed media

We discuss properties of the interaction of slow neutrons with nano-dispersed media and their application for neutron reflectors. In order to increase the accuracy of model simulation of the interaction of neutrons with nanopowders, we perform precise quantum mechanical calculation of potential scattering of neutrons on single nanoparticles using the method of phase functions. We compare results of precise calculations with those performed within first Born approximation for nanodiamonds with the radius of 2–5 nm and for neutron energies 3 × 10{sup -7}–10{sup -3} eV. Born approximation overestimates the probability of scattering to large angles, while the accuracy of evaluation of integral characteristics (cross sections, albedo) is acceptable. Using Monte-Carlo method, we calculate albedo of neutrons from different layers of piled up diamond nanopowder.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

N-[4-(3-Methyl-3-phenyl-cyclobutyl)-thiazol-2-yl]-N′-pyridin-2-ylmethylene-chloro-acetic acid hydrazide: Synthesis and configurational assignment based on X-ray, {sup 1}H, and {sup 13}C NMR and theoretical calculations

In this study, quantum chemical calculations based on the density functional theory have been carried out to examine the effects of N-[4-(3-methyl-3-phenyl-cyclobutyl)-thiazol-2-yl]-N′-pyridin-2-ylmethylene- chloro-acetic acid hydrazide. The calculated values are compared with the experimental data available for these molecules as a mean of validation of our proposed chemistry model. Aided by normal coordinate analysis and potential energy distributions, a confident vibrational assignment of all fundamentals is proposed herein. Additional support is given by {sup 1}H and {sup 13}C NMR spectra recorded with the sample dissolved in CDCl{sub 3} and by predicted chemical shifts at the B3LYP/6-31G(d)/6-311G+(d) levels obtained using the gauge-invariant atomic orbital method. The calculated HOMO and LUMO energies also confirm that the charge transfer occurs within the molecule. Thiazole-based compounds are potential storehouse for exploiting CH···O and CH···N hydrogen bonding interactions for molecular self-assembly.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Structure of (2E)-Ethyl 2-((4-(cyanomethoxy)benzylidene)hydrazono)-3,4-dimethyl-2,3-dihydrothiazole-5-carboxylate Studied by X-ray and DFT Calculations

The molecular structure of (2E)-ethyl 2-((4-(cyanomethoxy)benzylidene)hydrazono)-3,4-dimethyl-2,3-dihydrothiazole-5-carboxylate was determined using X-ray diffraction. The crystals are triclinic: sp. gr. P1; Z = 2; the unit cell parameters a = 8.4747(7) Å, b = 8.9382(8) Å, c = 11.9913(10) Å. The title compound has two C–H···O type intramolecular hydrogen bonds, one C–H···O and one C–H···N type intermolecular hydrogen bonds. For theoretical calculations, the molecular structure was investigated by DFT/B3LYP method with 6-311G(d) and 6-311G(d, p) basis sets. The calculated and experumental results (bond lenghts, bond angles, and dihedral angles) were compared with each other. Total energy, dipole moment, and mulliken atomic charges were calculated using two different basis sets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ab initio calculations of magnetic properties of the interstitially doped YFe{sub 11}Mo compound

The recent increase in the number of studies of RFe{sub 11–x}M{sub x} compounds is related to their promising application as permanent magnets. However, the insufficiently high value of the Curie temperature T{sub C} of these compounds is a barrier to their widespread use. The increase in the Curie temperature of these compounds is achieved by doping with the light nonmetallic atoms such as hydrogen, nitrogen, and carbon. In this paper, it is shown numerically that this doping leads to drastic changes of the electronic band dispersions in a wide energy region around the Fermi level. This in turn changes values of the magnetic moments of ions and Heisenberg exchange interaction parameters. The values of ab initio calculated magnetic moments and direct exchange interaction parameters make it possible to calculate the Curie temperatures for both parent and nitrogen-doped compounds within the mean-field approach to the Heisenberg model in the sample of YFe{sub 11}Mo, a typical representative of the R(Fe,M){sub 12}L class. Theoretical values of T{sub C} obtained for YFe{sub 11}Mo and YFe{sub 11}MoN (514 and 723 K respectively) are consistent with experimental ones (472 and 664 K) with an accuracy of 10%. Also, the calculated increase in T{sub C} upon nitrogenization (about 200 K) is in good agreement with the experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Variational Monte Carlo calculations of lithium atom in strong magnetic field

The variational Monte Carlo method is applied to investigate the ground state and some excited states of the lithium atom and its ions up to Z = 10 in the presence of an external magnetic field regime with γ = 0–100 arb. units. The effect of increasing field strength on the ground state energy is studied and precise values for the crossover field strengths were obtained. Our calculations are based on using accurate forms of trial wave functions, which were put forward in calculating energies in the absence of magnetic field. Furthermore, the value of Y at which ground-state energy of the lithium atom approaches to zero was calculated. The obtained results are in good agreement with the most recent values and also with the exact values.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Comparison of Different Methods of Atomistic Simulation To Calculate the Temperature of Phase Transition Using the Example of Zirconium

Comparison of various methods for calculating the phase transition temperature based on atomistic simulation such as two-phase simulation, single-phase simulation, calculations of thermodynamic potentials using phonon spectra is performed. Based on the simulation of zirconium, melting and transition between different crystalline phases are studied. It is shown that the Lindemann criterion is quite rigorous when describing the melting curve in the phase diagram. However, the threshold values of the Lindemann parameter corresponding to the phase transition may differ in different models. In addition, the results of the study allow us to estimate the accuracy of calculating the phase transition temperature based on the thermodynamic approach using the concept of entropy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Features of the application of the Monte Carlo method to calculations for large RBMK reactors and to model correction on the basis of data from in-core detectors

Statistical errors in sampling neutron fields in physically large systems like an RBMK are analyzed both qualitatively and quantitatively. Recommendations concerning the choice of parameters for calculations are given. A new procedure for Monte Carlo RBMK calculations with model corrections on the basis of data from in-core detectors is proposed. Dedicated software based on the CUDA software and hardware platform is developed for computational research. Results of testing the procedure and software in question via calculations for real RBMK reactors are discussed.

74 ATOMIC AND MOLECULAR PHYSICS↗

(γ, 2n)-Reaction cross-section calculations of several even-even lanthanide nuclei using different level density models

There are several level density models that can be used to predict photo-neutron cross sections. Some of them are Constant Temperature + Fermi Gas Model (CTFGM), Back-Shifted Fermi Gas Model (BSFM), Generalized Superfluid Model (GSM), Hartree-Fock-Bogoliubov microscopic Model (HFBM). In this study, the theoretical photo-neutron cross sections produced by (γ, 2n) reactions for several eveneven lanthanide nuclei such as {sup 140,142}Ce, {sup 142,144,146,148,150}Nd, {sup 144,148,150,152,154}Sm, and {sup 160}Gd have been calculated on the different level density models as mentioned above by using TALYS 1.6 and EMPIRE 3.1 computer codes for incident photon energies up to 30 MeV. The obtained results have been compared with each other and available experimental data existing in the EXFOR database. Generally, at least one level density model cross-section calculations are in agreement with the experimental results for all reactions except {sup 144}Sm(γ, 2n){sup 142}Sm along the incident photon energy, TALYS 1.6 BSFM option for the level density model cross-section calculations can be chosen if the experimental data are not available or are improbable to be produced due to the experimental difficulty.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Development of methods for calculating basic features of the nuclear contribution to single event upsets under the effect of protons of moderately high energy

As a continuation and a development of previous studies of our group that were devoted to the investigation of nuclear reactions induced by protons of moderately high energy (between 10 and 400 MeV) in silicon, aluminum, and tungsten atoms, the results obtained by exploring nuclear reactions on atoms of copper, which is among the most important components in materials for contact pads and pathways in modern and future ultralarge-scale integration circuits, especially in three-dimensional topology, are reported in the present article. The nuclear reactions in question lead to the formation of the mass and charge spectra of recoil nuclei ranging fromheavy target nuclei down to helium and hydrogen. The kineticenergy spectra of reaction products are calculated. The results of the calculations based on the procedure developed by our group are compared with the results of calculations and experiments performed by other authors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Substantiation of parameters of the geometric model of the research reactor core for the calculation using the Monte Carlo method

The choice of the spatial nodalization for the calculation of the power density and burnup distribution in a research reactor core with fuel assemblies of the IRT-3M and VVR-KN type using the program based on the Monte Carlo code is described. The influence of the spatial nodalization on the results of calculating basic neutronic characteristics and calculation time is investigated.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗