Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo Methods”

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,801 records · Page 100
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Model Independent Approach to the Single Photoelectron Calibration of Photomultiplier Tubes

The accurate calibration of photomultiplier tubes is critical in a wide variety of applications for which it is necessary to know the absolute number of detected photons or precisely determine the resolution of the signal. Conventional calibration methods rely on fitting the photomultiplier response to a low intensity light source with analytical approximations to the single photoelectron distribution. We show that this approach often leads to biased estimates due to an inability to model the full distribution accurately, especially at low charge values. We present a simple statistical method to extract the relevant single photoelectron calibration parameters (first two central moments) without making any assumptions about the underlying single photoelectron distribution. We illustrate the use of this method through the calibration of a Hamamatsu R11410 photomultiplier tube and study the accuracy and precision of the method using Monte Carlo simulations. The method is found to have significantly reduced bias compared to conventional methods and works under a wide range of light intensities, making it suitable for the simultaneous calibration of large arrays of photomultiplier tubes where uniform illumination may not be possible.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2012

High-resolution Thomson scattering system on the COMPASS tokamak: Evaluation of plasma parameters and error analysis

The electron density and temperature profiles measured by the Thomson scattering diagnostic on the COMPASS tokamak are used for estimation of electron kinetic energy, energy confinement time, and effective charge number Z{sub eff}. Data are compared with the line-integrated electron density measured by a microwave interferometer in an ohmically heated plasma with a circular cross section. An error analysis of both electron temperature and density are performed by two methods-a constant chi-square boundaries method and a Monte Carlo simulation, determining asymmetrical error bars for the electron temperature.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

Stochastic lag time in nucleated linear self-assembly

Protein aggregation is of great importance in biology, e.g., in amyloid fibrillation. The aggregation processes that occur at the cellular scale must be highly stochastic in nature because of the statistical number fluctuations that arise on account of the small system size at the cellular scale. We study the nucleated reversible self-assembly of monomeric building blocks into polymer-like aggregates using the method of kinetic Monte Carlo. Kinetic Monte Carlo, being inherently stochastic, allows us to study the impact of fluctuations on the polymerization reactions. One of the most important characteristic features in this kind of problem is the existence of a lag phase before self-assembly takes off, which is what we focus attention on. We study the associated lag time as a function of system size and kinetic pathway. We find that the leading order stochastic contribution to the lag time before polymerization commences is inversely proportional to the system volume for large-enough system size for all nine reaction pathways tested. Finite-size corrections to this do depend on the kinetic pathway.

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

RAPID COSMIC-RAY ACCELERATION AT PERPENDICULAR SHOCKS IN SUPERNOVA REMNANTS

Perpendicular shocks are shown to be rapid particle accelerators that perform optimally when the ratio u{sub s} of the shock speed to the particle speed roughly equals the ratio 1/η of the scattering rate to the gyro frequency. We use analytical methods and Monte-Carlo simulations to solve the kinetic equation that governs the anisotropy generated at these shocks, and find, for ηu{sub s} ≈ 1, that the spectral index softens by unity and the acceleration time increases by a factor of two compared to the standard result of the diffusive shock acceleration theory. These results provide a theoretical basis for the 30 year old conjecture that a supernova exploding into the wind of a Wolf–Rayet star may accelerate protons to an energy exceeding 10{sup 15} eV.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

Dosimetry revisited for the HDR {sup 192}Ir brachytherapy source model mHDR-v2

Purpose: Recently, the manufacturer of the HDR {sup 192}Ir mHDR-v2 brachytherapy source reported small design changes (referred to herein as mHDR-v2r) that are within the manufacturing tolerances but may alter the existing dosimetric data for this source. This study aimed to (1) check whether these changes affect the existing dosimetric data published for this source; (2) obtain new dosimetric data in close proximity to the source, including the contributions from {sup 192}Ir electrons and considering the absence of electronic equilibrium; and (3) obtain scatter dose components for collapsed cone treatment planning system implementation. Methods: Three different Monte Carlo (MC) radiation transport codes were used: MCNP5, PENELOPE2008, and GEANT4. The source was centrally positioned in a 40 cm radius water phantom. Absorbed dose and collision kerma were obtained using 0.1 mm (0.5 mm) thick voxels to provide high-resolution dosimetry near (far from) the source. Dose-rate distributions obtained with the three MC codes were compared. Results: Simulations of mHDR-v2 and mHDR-v2r designs performed with three radiation transport codes showed agreement typically within 0.2% for r{>=}0.25 cm. Dosimetric contributions from source electrons were significant for r<0.25 cm. The dose-rate constant and radial dose function were similar to those from previous MC studies of the mHDR-v2 design. The 2D anisotropy function also coincided with that of the mHDR-v2 design for r{>=}0.25 cm. Detailed results of dose distributions and scatter components are presented for the modified source design. Conclusions: Comparison of these results to prior MC studies showed agreement typically within 0.5% for r{>=}0.25 cm. If dosimetric data for r<0.25 cm are not needed, dosimetric results from the prior MC studies will be adequate.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Patient-based estimation of organ dose for a population of 58 adult patients across 13 protocol categories

Purpose: This study aimed to provide a comprehensive patient-specific organ dose estimation across a multiplicity of computed tomography (CT) examination protocols. Methods: A validated Monte Carlo program was employed to model a common CT system (LightSpeed VCT, GE Healthcare). The organ and effective doses were estimated from 13 commonly used body and neurological CT examination. The dose estimation was performed on 58 adult computational extended cardiac-torso phantoms (35 male, 23 female, mean age 51.5 years, mean weight 80.2 kg). The organ dose normalized by CTDI{sub vol} (h factor) and effective dose normalized by the dose length product (DLP) (k factor) were calculated from the results. A mathematical model was derived for the correlation between the h and k factors with the patient size across the protocols. Based on this mathematical model, a dose estimation iPhone operating system application was designed and developed to be used as a tool to estimate dose to the patients for a variety of routinely used CT examinations. Results: The organ dose results across all the protocols showed an exponential decrease with patient body size. The correlation was generally strong for the organs which were fully or partially located inside the scan coverage (Pearson sample correlation coefficient (r) of 0.49). The correlation was weaker for organs outside the scan coverage for which distance between the organ and the irradiation area was a stronger predictor of dose to the organ. For body protocols, the effective dose before and after normalization by DLP decreased exponentially with increasing patient's body diameter (r > 0.85). The exponential relationship between effective dose and patient's body diameter was significantly weaker for neurological protocols (r < 0.41), where the trunk length was a slightly stronger predictor of effective dose (0.15 < r < 0.46). Conclusions: While the most accurate estimation of a patient dose requires specific modeling of the patient anatomy, a first order approximation of organ and effective doses from routine CT scan protocols can be reasonably estimated using size specific factors. Estimation accuracy is generally poor for organ outside the scan range and for neurological protocols. The dose calculator designed in this study can be used to conveniently estimate and report the dose values for a patient across a multiplicity of CT scan protocols.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-I-67: X-Ray Fluorescence for Energy Response Calibration of a Photon Counting Detector: A Simulation Study

Purpose: To investigate the feasibility of energy calibration and energy response characterization of a photon counting detector using x-ray fluorescence. Methods: A comprehensive Monte Carlo simulation study was done to investigate the influence of various geometric components on the x-ray fluorescence measurement. Different materials, sizes, and detection angles were simulated using Geant4 Application for Tomographic Emission (GATE) Monte Carlo package. Simulations were conducted using 100 kVp tungsten-anode spectra with 2 mm Al filter for a single pixel cadmium telluride (CdTe) detector with 3 × 3 mm2 in detection area. The fluorescence material was placed 300 mm away from both the x-ray source and the detector. For angular dependence measurement, the distance was decreased to 30 mm to reduce the simulation time. Compound materials, containing silver, barium, gadolinium, hafnium, and gold in cylindrical shape, were simulated. The object size varied from 5 to 100 mm in diameter. The angular dependence of fluorescence and scatter were simulated from 20° to 170° with an incremental step of 10° to optimize the fluorescence to scatter ratio. Furthermore, the angular dependence was also experimentally measured using a spectrometer (X-123CdTe, Amptek Inc., MA) to validate the simulation results. Results: The detection angle between 120° to 160° resulted in more optimal x-ray fluorescence to scatter ratio. At a detection angle of 120°, the object size did not have a significant effect on the fluorescence to scatter ratio. The experimental results of fluorescence angular dependence are in good agreement with the simulation results. The Kα and Kβ peaks of five materials could be identified. Conclusion: The simulation results show that the x-ray fluorescence procedure has the potential to be used for detector energy calibration and detector response characteristics by using the optimal system geometry.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

SU-E-J-36: Combining CBCT Dose Into IMRT Treatment Planning

Purpose: Cone beam CT (CBCT) is increasingly used in patient setup for IMRT. Daily CBCT may provide effective localization, however, it introduces concern over excessive imaging dose. Previous studies investigated the calculation of excess CBCT dose, however, no study has yet treated this dose as a source of therapeutic radiation, optimized in consideration of PTV and OARs constrains. Here we present a novel combined MV+kV inverse optimization engine to weave the CBCT and MV dose together such that CBCT dose is used for both imaging and therapeutic purposes. This may mitigate some of the excess imaging dose effects of daily CBCT and allow complete evaluation of the CBCT dose prior to treatment. Methods: The EGSnrc Monte Carlo system was used to model a Varian Trilogy CBCT system and 6 MV treatment beam. Using the model, the dose to patient from treatment beam and imaging beam was calculated for ten patients. The standard IMRT objective function was modified to include CBCT dose. Treatment plan optimization using the MOSEK optimization tool was performed retrospectively with and without assuming kV radiation dose from CBCT, assuming one CBCT per fraction. Results: Across ten patients, the CBCT delivered peaks of between .4% and 3.0% of the prescription dose to the PTV, with average CBCT dose to the PTV between .3% and .8%. By including CBCT dose to skin as a constraint during optimization, peak skin dose is reduced by between 1.9% and 7.4%, and average skin dose is reduced by .2% to 3.3%. Conclusions: Pre-treatment CBCT may deliver a substantial amount of radiation dose to the target volume. By considering CBCT dose to skin at the point of treatment planning, it is possible to reduce patient skin dose from current clinical levels, and to provide patient treatment with the improved accuracy that daily CBCT provides.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Compartment modeling of dynamic brain PET—The impact of scatter corrections on parameter errors

Purpose: The aim of this study was to investigate the effect of scatter and its correction on kinetic parameters in dynamic brain positron emission tomography (PET) tumor imaging. The 2-tissue compartment model was used, and two different reconstruction methods and two scatter correction (SC) schemes were investigated. Methods: The GATE Monte Carlo (MC) software was used to perform 2 × 15 full PET scan simulations of a voxelized head phantom with inserted tumor regions. The two sets of kinetic parameters of all tissues were chosen to represent the 2-tissue compartment model for the tracer 3′-deoxy-3′-({sup 18}F)fluorothymidine (FLT), and were denoted FLT{sub 1} and FLT{sub 2}. PET data were reconstructed with both 3D filtered back-projection with reprojection (3DRP) and 3D ordered-subset expectation maximization (OSEM). Images including true coincidences with attenuation correction (AC) and true+scattered coincidences with AC and with and without one of two applied SC schemes were reconstructed. Kinetic parameters were estimated by weighted nonlinear least squares fitting of image derived time–activity curves. Calculated parameters were compared to the true input to the MC simulations. Results: The relative parameter biases for scatter-eliminated data were 15%, 16%, 4%, 30%, 9%, and 7% (FLT{sub 1}) and 13%, 6%, 1%, 46%, 12%, and 8% (FLT{sub 2}) for K{sub 1}, k{sub 2}, k{sub 3}, k{sub 4}, V{sub a}, and K{sub i}, respectively. As expected, SC was essential for most parameters since omitting it increased biases by 10 percentage points on average. SC was not found necessary for the estimation of K{sub i} and k{sub 3}, however. There was no significant difference in parameter biases between the two investigated SC schemes or from parameter biases from scatter-eliminated PET data. Furthermore, neither 3DRP nor OSEM yielded the smallest parameter biases consistently although there was a slight favor for 3DRP which produced less biased k{sub 3} and K{sub i} estimates while OSEM resulted in a less biased V{sub a}. The uncertainty in OSEM parameters was about 26% (FLT{sub 1}) and 12% (FLT{sub 2}) larger than for 3DRP although identical postfilters were applied. Conclusions: SC was important for good parameter estimations. Both investigated SC schemes performed equally well on average and properly corrected for the scattered radiation, without introducing further bias. Furthermore, 3DRP was slightly favorable over OSEM in terms of kinetic parameter biases and SDs.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-J-06: Additional Imaging Guidance Dose to Patient Organs Resulting From X-Ray Tubes Used in CyberKnife Image Guidance System

Purpose: The use of image-guided radiation therapy (IGRT) has become increasingly common, but the additional radiation exposure resulting from repeated image guidance procedures raises concerns. Although there are many studies reporting imaging dose from different image guidance devices, imaging dose for the CyberKnife Robotic Radiosurgery System is not available. This study provides estimated organ doses resulting from image guidance procedures on the CyberKnife system. Methods: Commercially available Monte Carlo software, PCXMC, was used to calculate average organ doses resulting from x-ray tubes used in the CyberKnife system. There are seven imaging protocols with kVp ranging from 60 – 120 kV and 15 mAs for treatment sites in the Cranium, Head and Neck, Thorax, and Abdomen. The output of each image protocol was measured at treatment isocenter. For each site and protocol, Adult body sizes ranging from anorexic to extremely obese were simulated since organ dose depends on patient size. Doses for all organs within the imaging field-of-view of each site were calculated for a single image acquisition from both of the orthogonal x-ray tubes. Results: Average organ doses were <1.0 mGy for every treatment site and imaging protocol. For a given organ, dose increases as kV increases or body size decreases. Higher doses are typically reported for skeletal components, such as the skull, ribs, or clavicles, than for softtissue organs. Typical organ doses due to a single exposure are estimated as 0.23 mGy to the brain, 0.29 mGy to the heart, 0.08 mGy to the kidneys, etc., depending on the imaging protocol and site. Conclusion: The organ doses vary with treatment site, imaging protocol and patient size. Although the organ dose from a single image acquisition resulting from two orthogonal beams is generally insignificant, the sum of repeated image acquisitions (>100) could reach 10–20 cGy for a typical treatment fraction.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-T-129: Impact of Radial Fluctuations in RBE for Therapeutic Proton Beams

Purpose: To evaluate the off axis relative biological effectiveness (RBE) for actively scanned proton beams and determine if a constant radial RBE can be assumed. Methods: The PHITS Monte Carlo code paired with a microscopic analytical function was used to determine probability distribution functions of the lineal energy in 0.3µm diameter spheres throughout a water phantom. Twenty million primary protons were simulated for a 0.6cm diameter pencil beam. Beam energies corresponding to Bragg Peak depths of 50, 100, 150, 200, 250, and 300mm were used and evaluated transversely every millimeter and radially for annuli of 1.0, 2.0, 3.0, 3.2, 3.4, 3.6, 4.0, 5.0, 10.0, 15.0, 20.0 and 25.0mm outer radius. The acquired probability distributions were reduced to dose-mean lineal energies and applied to the modified microdosimetric kinetic model, for human submandibular gland (HSG) cells, to calculate relative biological effectiveness (RBE) compared to 60Co beams at the 10% survival threshold. Results: RBE was generally seen to increase as distance from the central axis (CAX) increased. However, this increase was only seen in low dose regions and its overall effects on the transverse biological dose remains low. In the entrance region of the phantom (10mm depth), minimum and maximum calculated RBEs varied between 15.22 and 18.88% for different energies. At the Bragg peak, this difference ranged from 3.15 to 26.77%. Despite these rather large variations the dose-weighted RBE and the CAX RBE varied by less than 0.14% at 10mm depth and less than 0.16% at the Bragg peak. Similarly small variations were found at all depths proximal of the Bragg peak. Conclusion: Although proton RBE does vary radially, its overall effect on biological dose is minimal and the use of a radially constant RBE in treatment planning for scanned proton beams would not produce large errors.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-T-144: Analytical Closed Form Approximation for Carbon Ion Bragg Curves in Water

Purpose: Semi-empirical modeling is a powerful computational method in radiation dosimetry. A set of approximations exist for proton ion depth dose distribution (DDD) in water. However, the modeling is more complicated for carbon ions due to fragmentation. This study addresses this by providing and evaluating a new methodology for DDD modeling of carbon ions in water. Methods: The FLUKA, Monte Carlo (MC) general-purpose transport code was used for simulation of carbon DDDs for energies of 100–400 MeV in water as reference data model benchmarking. Based on Thomas Bortfeld’s closed form equation approximating proton Bragg Curves as a basis, we derived the critical constants for a beam of Carbon ions by applying models of radiation transport by Lee et. al. and Geiger to our simulated Carbon curves. We hypothesized that including a new exponential (κ) residual distance parameter to Bortfeld’s fluence reduction relation would improve DDD modeling for carbon ions. We are introducing an additional term to be added to Bortfeld’s equation to describe fragmentation tail. This term accounts for the pre-peak dose from nuclear fragments (NF). In the post peak region, the NF transport will be treated as new beams utilizing the Glauber model for interaction cross sections and the Abrasion- Ablation fragmentation model. Results: The carbon beam specific constants in the developed model were determined to be : p= 1.75, β=0.008 cm-1, γ=0.6, α=0.0007 cm MeV, σmono=0.08, and the new exponential parameter κ=0.55. This produced a close match for the plateau part of the curve (max deviation 6.37%). Conclusion: The derived semi-empirical model provides an accurate approximation of the MC simulated clinical carbon DDDs. This is the first direct semi-empirical simulation for the dosimetry of therapeutic carbon ions. The accurate modeling of the NF tail in the carbon DDD will provide key insight into distal edge dose deposition formation.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2016

SU-F-T-366: Dosimetric Parameters Enhancement of 120-Leaf Millennium MLC Using EGSnrc and IAEA Phase-Space Data

Purpose: Recently, the multileaf collimators (MLC) have become an important part of any LINAC collimation systems because they reduce the treatment planning time and improves the conformity. Important factors that affects the MLCs collimation performance are leaves material composition and their thickness. In this study, we investigate the main dosimetric parameters of 120-leaf Millennium MLC including dose in the buildup point, physical penumbra as well as average and end leaf leakages. Effects of the leaves geometry and density on these parameters are evaluated Methods: From EGSnrc Monte Carlo code, BEAMnrc and DOSXYZnrc modules are used to evaluate the dosimetric parameters of a water phantom exposed to a Varian xi for 100cm SSD. Using IAEA phasespace data just above MLC (Z=46cm) and BEAMnrc, for the modified 120-leaf Millennium MLC a new phase space data at Z=52cm is produces. The MLC is modified both in leaf thickness and material composition. EGSgui code generates 521ICRU library for tungsten alloys. DOSXYZnrc with the new phase space evaluates the dose distribution in a water phantom of 60×60×20 cm3 with voxel size of 4×4×2 mm3. Using DOSXYZnrc dose distributions for open beam and closed beam as well as the leakages definition, end leakage, average leakage and physical penumbra are evaluated. Results: A new MLC with improved dosimetric parameters is proposed. The physical penumbra for proposed MLC is 4.7mm compared to 5.16 mm for Millennium. Average leakage in our design is reduced to 1.16% compared to 1.73% for Millennium, the end leaf leakage suggested design is also reduced to 4.86% compared to 7.26% of Millennium. Conclusion: The results show that the proposed MLC with enhanced dosimetric parameters could improve the conformity of treatment planning.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-G-TeP3-13: The Role of Nanoscale Energy Deposition in the Development of Gold Nanoparticle-Enhanced Radiotherapy

Purpose: Gold nanoparticles (GNPs) can enhance radiotherapy effects. The high photoelectric cross section of gold relative to tissue, particularly at lower energies, leads to localized dose enhancement. However in a clinical context, photon energies must also be sufficient to reach a target volume at a given depth. These properties must be balanced to optimize such a therapy. Given that nanoscale energy deposition patterns around GNPs play a role in determining biological outcomes, in this work we seek to establish their role in this optimization process. Methods: The PENELOPE Monte Carlo code was used to generate spherical dose deposition kernels in 1000 nm diameter spheres around 50 nm diameter GNPs in response to monoenergetic photons incident on the GNP. Induced “lesions” were estimated by either a local effect model (LEM) or a mean dose model (MDM). The ratio of these estimates was examined for a range of photon energies (10 keV to 2 MeV), for three sets of linear-quadratic parameters. Results: The models produce distinct differences in expected lesion values, the lower the alpha-beta ratio, the greater the difference. The ratio of expected lesion values remained constant within 5% for energies of 40 keV and above across all parameter sets and rose to a difference of 35% for lower energies only for the lowest alpha-beta ratio. Conclusion: Consistent with other work, these calculations suggest nanoscale energy deposition patterns matter in predicting biological response to GNP-enhanced radiotherapy. However the ratio of expected lesions between the different models is largely independent of energy, indicating that GNP-enhanced radiotherapy scenarios can be optimized in photon energy without consideration of the nanoscale patterns. Special attention may be warranted for energies of 20 keV or below and low alpha-beta ratios.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

WE-H-BRA-09: Application of a Modified Microdosimetric-Kinetic Model to Analyze Relative Biological Effectiveness of Ions Relevant to Light Ion Therapy Using the Particle Heavy Ion Transport System

Purpose: To evaluate the dose and biological effectiveness of various ions that could potentially be used for actively scanned particle therapy. Methods: The PHITS Monte Carlo code paired with a microscopic analytical function was used to determine probability distribution functions of the lineal energy in 0.3µm diameter spheres throughout a water phantom. Twenty million primary particles for 1H beams and ten million particles for 4He, 7Li, 10B, 12C, 14N, 16O, and 20Ne were simulated for 0.6cm diameter pencil beams. Beam energies corresponding to Bragg peak depths of 50, 100, 150, 200, 250, and 300mm were used and evaluated transversely every millimeter and radially in annuli with outer radius of 1.0, 2.0, 3.0, 3.2, 3.4, 3.6, 4.0, 5.0, 10.0, 15.0, 20.0 and 25.0mm. The acquired probability distributions were reduced to dose-mean lineal energies and applied to the modified microdosimetric kinetic model for five different cell types to calculate relative biological effectiveness (RBE) compared to 60Co beams at the 10% survival threshold. The product of the calculated RBEs and the simulated physical dose was taken to create biological dose and comparisons were then made between the various ions. Results: Transversely, the 10B beam was seen to minimize relative biological dose in both the constant and accelerated dose change regions, proximal to the Bragg Peak, for all beams traveling greater than 50mm. For the 50mm beam, 7Li was seen to provide the most optimal biological dose profile. Radially small fluctuations (<4.2%) were seen in RBE while physical dose was greater than 1% for all beams. Conclusion: Even with the growing usage of 12C, it may not be the most optimal ion in all clinical situations. Boron was calculated to have slightly enhanced RBE characteristics, leading to lower relative biological doses.

60 APPLIED LIFE SCIENCES↗

Refine your search to explore more results.