Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computerized 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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

SU-E-T-593: Outcomes and Toxicities From a Clinical Trial of APBI Using MERT+IMRT with the Same XMLC

Purpose: We present the results from a clinical trial of accelerated partial breast irradiation (APBI), using mixed modulated photon and electron beams (MERT+IMRT) with the same photon multileaf collimator (xMLC). Methods: Seven patients were enrolled in the first year of the APBI clinical trial. Patients were selected following the conditions included in the NSABP B-39/RTOG 0413 protocol. The targets and clinically relevant normal structures were contoured on the CT images following this protocol for APBI-EBRT. All treatments were delivered using combined modulated electron and photon beams by means of the same xMLC installed in a SIEMENS Primus linac, with a reduced SSD equal to 60 cm for electron beams. The plans were performed with a treatment planning system based on full Monte Carlo simulations, called CARMEN, developed by our group. Simultaneously, an alternative IMRT plan was calculated with the commercial TPS PINNACLE v8.0m (Philips), and both plans were compared. An ad-hoc breast phantom with semi-spherical geometry called NAOMI was designed for a specific QA protocol. Patients received a total dose of 38.5 Gy, delivered in 10 fractions over 5 consecutive days, with a twice-a-day hypofractionated schema.Follow-up visits during 2.5 years on average were repeated at 1 month post-treatment, every 3 months for the first year, and every 6 months for the second year. Toxicity was scored according to National Cancer Institute Common Terminology Criteria for Adverse Events (CTCAE version 3.0). Results: This APBI technique achieved high loco-regional control rates and showed low acute toxicity (grade 1 of CTCAE) and no toxicities from first month onwards. Photographic assessment of cosmesis showed skin excellent results. Conclusion: The clinical results achieved with MERT+IMRT by using the same xMLC are comparable or even better than those obtained with other APBI techniques, thanks to a software solution without any additional equipment or specific device.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

SU-E-T-285: Dose Variation at Bone in Small-Animal Irradiation: A Monte Carlo Study Using Monoenergetic Photon Beams

Purpose: The aim of this study is to investigate the variation of bone dose on photon beam energy (keV – MeV) in small-animal irradiation. Dosimetry of homogeneous and inhomogeneous phantoms as per the same mouse computed tomography image set were calculated using the DOSCTP and DOSXYZnrc based on the EGSnrc Monte Carlo code. Methods: Monte Carlo simulations for the homogeneous and inhomogeneous mouse phantom irradiated by a 360 degree photon arc were carried out. Mean doses of the bone tissue in the irradiated volumes were calculated at various photon beam energies, ranging from 50 keV to 1.25 MeV. The effect of bone inhomogeneity was examined through the Inhomogeneous Correction Factor (ICF), a dose ratio of the inhomogeneous to the homogeneous medium. Results: From our Monte Carlo results, higher mean bone dose and ICF were found when using kilovoltage photon beams compared to megavoltage. In beam energies ranging from 50 keV to 200 keV, the bone dose was found maximum at 50 keV, and decreased significantly from 2.6 Gy to 0.55 Gy, when 2 Gy was delivered at the center of the phantom (isocenter). Similarly, the ICF were found decreasing from 4.5 to 1 when the photon beam energy was increased from 50 keV to 200 keV. Both mean bone dose and ICF remained at about 0.5 Gy and 1 from 200 keV to 1.25 MeV with insignificant variation, respectively. Conclusion: It is concluded that to avoid high bone dose in the small-animal irradiation, photon beam energy higher than 200 keV should be used with the ICF close to one, and bone dose comparable to the megavoltage beam where photoelectric effect is not dominant.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

TU-CD-207-10: Dedicated Cone-Beam Breast CT: Design of a 3-D Beam-Shaping Filter

Purpose: To design a 3 -D beam-shaping filter for cone-beam breast CT for equalizing x-ray photon fluence incident on the detector along both fan and cone angle directions. Methods: The 3-D beam-shaping filter was designed as the sum of two filters: a bow-tie filter assuming cylindrical breast and a 3D difference filter equivalent to the difference in projected thickness between the cylinder and the real breast. Both filters were designed with breast-equivalent material and converted to Al for the targeted x-ray spectrum. The bow-tie was designed for the largest diameter cylindrical breast by determining the fan-angle dependent path-length and the filter thickness needed to equalize the fluence. A total of 23,760 projections (180 projections of 132 binary breast CT volumes) were averaged, scaled for the largest breast, and subtracted from the projection of the largest diameter cylindrical breast to provide the 3D difference filter. The 3 -D beam shaping filter was obtained by summing the two filters. Numerical simulations with semi-ellipsoidal breasts of 10–18 cm diameter (chest-wall to nipple length=0.75 x diameter) were conducted to evaluate beam equalization. Results: The proposed 3-D beam-shaping filter showed a 140% -300% improvement in equalizing the photon fluence along the chest-wall to nipple (cone-angle) direction compared to a bow-tie filter. The improvement over bow-tie filter was larger for breasts with longer chest-wall to nipple length. Along the radial (fan-angle) direction, the performance of the 3-D beam shaping filter was marginally better than the bow-tie filter, with 4%-10% improvement in equalizing the photon fluence. For a ray traversing the chest-wall diameter of the breast, the filter transmission ratio was >0.95. Conclusion: The 3-D beam shaping filter provided substantial advantage over bow-tie filter in equalizing the photon fluence along the cone-angle direction. In conjunction with a 2-axis positioner, the filter can accommodate breasts of varying dimensions and chest-wall inclusion. Supported in part by NIH R01 CA128906 and R21 CA134128. The contents are solely the responsibility of the authors and do not reflect the official views of the NIH or NCI.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Synchronized multiartifact reduction with tomographic reconstruction (SMART-RECON): A statistical model based iterative image reconstruction method to eliminate limited-view artifacts and to mitigate the temporal-average artifacts in time-resolved CT

Purpose: In x-ray computed tomography (CT), a violation of the Tuy data sufficiency condition leads to limited-view artifacts. In some applications, it is desirable to use data corresponding to a narrow temporal window to reconstruct images with reduced temporal-average artifacts. However, the need to reduce temporal-average artifacts in practice may result in a violation of the Tuy condition and thus undesirable limited-view artifacts. In this paper, the authors present a new iterative reconstruction method, synchronized multiartifact reduction with tomographic reconstruction (SMART-RECON), to eliminate limited-view artifacts using data acquired within an ultranarrow temporal window that severely violates the Tuy condition. Methods: In time-resolved contrast enhanced CT acquisitions, image contrast dynamically changes during data acquisition. Each image reconstructed from data acquired in a given temporal window represents one time frame and can be denoted as an image vector. Conventionally, each individual time frame is reconstructed independently. In this paper, all image frames are grouped into a spatial–temporal image matrix and are reconstructed together. Rather than the spatial and/or temporal smoothing regularizers commonly used in iterative image reconstruction, the nuclear norm of the spatial–temporal image matrix is used in SMART-RECON to regularize the reconstruction of all image time frames. This regularizer exploits the low-dimensional structure of the spatial–temporal image matrix to mitigate limited-view artifacts when an ultranarrow temporal window is desired in some applications to reduce temporal-average artifacts. Both numerical simulations in two dimensional image slices with known ground truth and in vivo human subject data acquired in a contrast enhanced cone beam CT exam have been used to validate the proposed SMART-RECON algorithm and to demonstrate the initial performance of the algorithm. Reconstruction errors and temporal fidelity of the reconstructed images were quantified using the relative root mean square error (rRMSE) and the universal quality index (UQI) in numerical simulations. The performance of the SMART-RECON algorithm was compared with that of the prior image constrained compressed sensing (PICCS) reconstruction quantitatively in simulations and qualitatively in human subject exam. Results: In numerical simulations, the 240{sup ∘} short scan angular span was divided into four consecutive 60{sup ∘} angular subsectors. SMART-RECON enables four high temporal fidelity images without limited-view artifacts. The average rRMSE is 16% and UQIs are 0.96 and 0.95 for the two local regions of interest, respectively. In contrast, the corresponding average rRMSE and UQIs are 25%, 0.78, and 0.81, respectively, for the PICCS reconstruction. Note that only one filtered backprojection image can be reconstructed from the same data set with an average rRMSE and UQIs are 45%, 0.71, and 0.79, respectively, to benchmark reconstruction accuracies. For in vivo contrast enhanced cone beam CT data acquired from a short scan angular span of 200{sup ∘}, three 66{sup ∘} angular subsectors were used in SMART-RECON. The results demonstrated clear contrast difference in three SMART-RECON reconstructed image volumes without limited-view artifacts. In contrast, for the same angular sectors, PICCS cannot reconstruct images without limited-view artifacts and with clear contrast difference in three reconstructed image volumes. Conclusions: In time-resolved CT, the proposed SMART-RECON method provides a new method to eliminate limited-view artifacts using data acquired in an ultranarrow temporal window, which corresponds to approximately 60{sup ∘} angular subsectors.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Will weight loss cause significant dosimetric changes of target volumes and organs at risk in nasopharyngeal carcinoma treated with intensity-modulated radiation therapy?

This study aimed to quantify dosimetric effects of weight loss for nasopharyngeal carcinoma (NPC) treated with intensity-modulated radiation therapy (IMRT). Overall, 25 patients with NPC treated with IMRT were enrolled. We simulated weight loss during IMRT on the computer. Weight loss model was based on the planning computed tomography (CT) images. The original external contour of head and neck was labeled plan 0, and its volume was regarded as pretreatment normal weight. We shrank the external contour with different margins (2, 3, and 5 mm) and generated new external contours of head and neck. The volumes of reconstructed external contours were regarded as weight during radiotherapy. After recontouring outlines, the initial treatment plan was mapped to the redefined CT scans with the same beam configurations, yielding new plans. The computer model represented a theoretical proportional weight loss of 3.4% to 13.7% during the course of IMRT. The dose delivered to the planning target volume (PTV) of primary gross tumor volume and clinical target volume significantly increased by 1.9% to 2.9% and 1.8% to 2.9% because of weight loss, respectively. The dose to the PTV of gross tumor volume of lymph nodes fluctuated from −2.0% to 1.0%. The dose to the brain stem and the spinal cord was increased (p < 0.001), whereas the dose to the parotid gland was decreased (p < 0.001). Weight loss may lead to significant dosimetric change during IMRT. Repeated scanning and replanning for patients with NPC with an obvious weight loss may be necessary.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

SU-E-J-08: Comparison of Unintended Radiation Doses to Organs at Risk Resulting From the Out-Of-Field Therapeutic Beams and From Image-Guidance X-Ray Procedures

Purpose: The unintended radiation dose to organs at risk (OAR) can be contributed from imaging guidance procedures as well as from leakage and scatter of therapeutic beams. This study compares the imaging dose with the unintended out-of-field therapeutic dose to patient sensitive organs. Methods: The Monte Carlo EGSnrc user codes, BEAMnrc and DOSXYZnrc, were used to simulate kV X-ray sources from imaging devices as well as the therapeutic IMRT/VMAT beams and to calculate doses to target and OARs on patient treatment planning CT images. The accuracy of the Monte Carlo simulations was benchmarked against measurements in phantoms. The dose-volume histogram was utilized in analyzing the patient organ doses. Results: The dose resulting from Standard Head kV-CBCT scans to bone and soft tissues ranges from 0.7 to 1.1 cGy and from 0.03 to 0.3 cGy, respectively. The dose resulting from Thorax scans on the chest to bone and soft tissues ranges from 1.1 to 1.8 cGy and from 0.3 to 0.6 cGy, respectively. The dose resulting from Pelvis scans on the abdomen to bone and soft tissues range from 3.2 to 4.2 cGy and from 1.2 to 2.2 cGy, respectively. The out-of-field doses to OAR are sensitive to the distance between the treated target and the OAR. For a typical Head-and-Neck IMRT/VMAT treatment the out-of-field doses to eyes are 1–3% of the target dose, or 2–6 cGy per fraction. Conclusion: The imaging doses to OAR are predictable based on the imaging protocols used when OARs are within the imaged volume and can be estimated and accounted for by using tabulated values. The unintended out-of-field doses are proportional to the target dose, strongly depend on the distance between the treated target and OAR, and are generally higher comparing to the imaging dose. This work was partially supported by Varian research grant VUMC40590.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-T-679: Monte Carlo Simulation On Surface Dose in Preclinical Irradiation Using Monoenergetic Photon Beams

Purpose: This study investigated the surface dose variation in preclinical irradiation using small animal, when monoenergetic photon beams with energy range from 50 keV to 1.25 MeV were used. Methods: Inhomogeneous, homogeneous and bone-tissue homogeneous mouse phantom based on the same CT image set were used. The homogeneous and bone-tissue homogeneous phantom were created with the relative electron density of all and only bone voxels of the mouse overridden to one, respectively. Monte Carlo simulation based on the EGSnrc-based code was used to calculate the surface dose, when the phantoms were irradiated by a 360° photon arc with energies ranging from 50 keV to 1.25 MeV. The mean surface doses of the three phantoms were calculated. In addition, the surface doses from partial arcs, 45°–315°, 125°–225°, 45°–125° and 225°–315° covering the anterior, posterior, right lateral and left lateral region of the mouse were determined using different photon beam energies. Results: When the prescribed dose at the isocenter of the mouse was 2 Gy, the maximum mean surface doses, found at the 50-keV photon beams, were 0.358 Gy, 0.363 Gy and 0.350 Gy for the inhomogeneous, homogeneous and bone-tissue homogeneous mouse phantom, respectively. The mean surface dose of the mouse was found decreasing with an increase of the photon beam energy. For surface dose in different orientations, the lateral regions of the mouse were receiving lower dose than the anterior and posterior regions. This may be due to the increase of beam attenuation along the horizontal (left-right) axis than the vertical (anterior-posterior) in the mouse. Conclusion: It is concluded that consideration of phantom inhomogeneity in the dose calculation resulted in a lower mean surface dose of the mouse. The mean surface dose also decreased with an increase of photon beam energy in the kilovoltage range.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

WE-EF-207-06: Dedicated Cone-Beam Breast CT with Laterally-Shifted Detector: Monte Carlo Evaluation of X-Ray Scatter Distribution and Scatter-To-Primary Ratio

Purpose: To determine the spatial distribution of x-ray scatter and scatter-to-primary ratio (SPR) in projections during cone-beam breast CT (CBBCT) with laterally-shifted detector that results in coronal (fan-angle) truncation. Methods: We hypothesized that CBBCT with coronal truncation would lower SPR due to reduction in irradiated breast volume, and that the location of maximum x-ray scatter fluence (scatter-peak) in the detector plane can be determined from the ratio of irradiated-to-total breast volume, breast dimensions and system geometry. Monte Carlo simulations (GEANT4) reflecting a prototype CBBCT system were used to record the position-dependent primary and scatter x-ray photon fluence incident on the detector without coronal truncation (full fan-angle, 2f=24-degrees) and with coronal truncation (fan-angle, f+ f=12+2.7-degrees). Semi-ellipsoidal breasts (10/14/18-cm diameter, chest-wall to nipple length: 0.75xdiameter, 2%/14%/100% fibroglandular content) aligned with the axis-of-rotation (AOR) were modeled. Mono-energy photons were simulated and weighted for 2 spectra (49kVp, 1.4-mm Al HVL; 60kVp, 3.76-mm Al HVL). In addition to SPR, the scatter maps were analyzed to identify the location of the scatter-peak. Results: For CBBCT without fan-angle truncation, the scatter-peaks were aligned with the projection of the AOR onto the detector for all breasts. With truncated fan-beam, the scatter-peaks were laterally-shifted from the projection of the AOR along the fan-angle direction by 14/38/70-pixels for 10/14/18-cm diameter breasts. The corresponding theoretical shifts were 14.8/39.7/68-pixels (p=0.47, 2-tailed paired-ratio t-test). Along the cone-angle, the shift in scatter-peaks between truncated and full-fan angle CBBCT were 2/2/4 -pixels for 10/14/18-cm diameter breasts. CBBCT with fan-angle truncation reduced SPR by 14/22/28% for 10/14/18-cm diameter breasts. 60kVp reduced SPR by 21–25% compared to 49kVp. Peak SPR for CBBCT with fan-angle truncation (60kVp) were 0.09/0.25/0.73 for 10/14/18-cm diameter breasts. Conclusion: CBBCT with laterally-shifted detector geometry and with appropriate kVp/beam quality reduces SPR. If residual scatter needs correction, the location corresponding to scatter-peak can be analytically computed. This work was supported in part by NIH R01 CA128906. The contents are solely the responsibility of the authors and do not reflect the official views of the NIH or NCI.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Resting myocardial blood flow quantification using contrast-enhanced magnetic resonance imaging in the presence of stenosis: A computational fluid dynamics study

Purpose: The extent to which atherosclerotic plaques affect contrast agent (CA) transport in the coronary arteries and, hence, quantification of myocardial blood flow (MBF) using magnetic resonance imaging (MRI) is unclear. The purpose of this work was to evaluate the influence of plaque induced stenosis both on CA transport and on the accuracy of MBF quantification. Methods: Computational fluid dynamics simulations in a high-detailed realistic vascular model were employed to investigate CA bolus transport in the coronary arteries. The impact of atherosclerosis was analyzed by inserting various medium- to high-grade stenoses in the vascular model. The influence of stenosis morphology was examined by varying the stenosis shapes but keeping the area reduction constant. Errors due to CA bolus transport were analyzed using the tracer-kinetic model MMID4. Results: Dispersion of the CA bolus was found in all models and for all outlets, but with a varying magnitude. The impact of stenosis was complex: while high-grade stenoses amplified dispersion, mild stenoses reduced the effect. Morphology was found to have a marked influence on dispersion for a small number of outlets in the post-stenotic region. Despite this marked influence on the concentration–time curves, MBF errors were less affected by stenosis. In total, MBF was underestimated by −7.9% to −44.9%. Conclusions: The presented results reveal that local hemodynamics in the coronary vasculature appears to have a direct impact on CA bolus dispersion. Inclusion of atherosclerotic plaques resulted in a complex alteration of this effect, with both degree of area reduction and stenosis morphology affecting the amount of dispersion. This strong influence of vascular transport effects impairs the accuracy of MRI-based MBF quantification techniques and, potentially, other bolus-based perfusion measurement techniques like computed tomography perfusion imaging.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Monte Carlo simulations of patient dose perturbations in rotational-type radiotherapy due to a transverse magnetic field: A tomotherapy investigation

Purpose: Several groups are exploring the integration of magnetic resonance (MR) image guidance with radiotherapy to reduce tumor position uncertainty during photon radiotherapy. The therapeutic gain from reducing tumor position uncertainty using intrafraction MR imaging during radiotherapy could be partially offset if the negative effects of magnetic field-induced dose perturbations are not appreciated or accounted for. The authors hypothesize that a more rotationally symmetric modality such as helical tomotherapy will permit a systematic mediation of these dose perturbations. This investigation offers a unique look at the dose perturbations due to homogeneous transverse magnetic field during the delivery of Tomotherapy{sup ®} Treatment System plans under varying degrees of rotational beamlet symmetry. Methods: The authors accurately reproduced treatment plan beamlet and patient configurations using the Monte Carlo code GEANT4. This code has a thoroughly benchmarked electromagnetic particle transport physics package well-suited for the radiotherapy energy regime. The three approved clinical treatment plans for this study were for a prostate, head and neck, and lung treatment. The dose heterogeneity index metric was used to quantify the effect of the dose perturbations to the target volumes. Results: The authors demonstrate the ability to reproduce the clinical dose–volume histograms (DVH) to within 4% dose agreement at each DVH point for the target volumes and most planning structures, and therefore, are able to confidently examine the effects of transverse magnetic fields on the plans. The authors investigated field strengths of 0.35, 0.7, 1, 1.5, and 3 T. Changes to the dose heterogeneity index of 0.1% were seen in the prostate and head and neck case, reflecting negligible dose perturbations to the target volumes, a change from 5.5% to 20.1% was observed with the lung case. Conclusions: This study demonstrated that the effect of external magnetic fields can be mitigated by exploiting a more rotationally symmetric treatment modality.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

Investigation of energy weighting using an energy discriminating photon counting detector for breast CT

Purpose: Breast CT is an emerging imaging technique that can portray the breast in 3D and improve visualization of important diagnostic features. Early clinical studies have suggested that breast CT has sufficient spatial and contrast resolution for accurate detection of masses and microcalcifications in the breast, reducing structural overlap that is often a limiting factor in reading mammographic images. For a number of reasons, image quality in breast CT may be improved by use of an energy resolving photon counting detector. In this study, the authors investigate the improvements in image quality obtained when using energy weighting with an energy resolving photon counting detector as compared to that with a conventional energy integrating detector.Methods: Using computer simulation, realistic CT images of multiple breast phantoms were generated. The simulation modeled a prototype breast CT system using an amorphous silicon (a-Si), CsI based energy integrating detector with different x-ray spectra, and a hypothetical, ideal CZT based photon counting detector with capability of energy discrimination. Three biological signals of interest were modeled as spherical lesions and inserted into breast phantoms; hydroxyapatite (HA) to represent microcalcification, infiltrating ductal carcinoma (IDC), and iodine enhanced infiltrating ductal carcinoma (IIDC). Signal-to-noise ratio (SNR) of these three lesions was measured from the CT reconstructions. In addition, a psychophysical study was conducted to evaluate observer performance in detecting microcalcifications embedded into a realistic anthropomorphic breast phantom.Results: In the energy range tested, improvements in SNR with a photon counting detector using energy weighting was higher (than the energy integrating detector method) by 30%–63% and 4%–34%, for HA and IDC lesions and 12%–30% (with Al filtration) and 32%–38% (with Ce filtration) for the IIDC lesion, respectively. The average area under the receiver operating characteristic curve (AUC) for detection of microcalcifications was higher by greater than 19% (for the different energy weighting methods tested) as compared to the AUC obtained with an energy integrating detector.Conclusions: This study showed that breast CT with a CZT photon counting detector using energy weighting can provide improvements in pixel SNR, and detectability of microcalcifications as compared to that with a conventional energy integrating detector. Since a number of degrading physical factors were not modeled into the photon counting detector, this improvement should be considered as an upper bound on achievable performance.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2016

SU-F-T-62: Three-Dimensional Dosimetric Gamma Analysis for Impacts of Tissue Inhomogeneity Using Monte Carlo Simulation in Intracavitary Brachytheray for Cervix Carcinoma

Purpose: The aim of this study was to investigate the impacts of tissue inhomogeneity on dose distributions using a three-dimensional (3D) gamma analysis in cervical intracavitary brachytherapy using Monte Carlo (MC) simulations. Methods: MC simulations for comparison of dose calculations were performed in a water phantom and a series of CT images of a cervical cancer patient (stage: Ib; age: 27) by employing a MC code, Particle and Heavy Ion Transport Code System (PHIT) version 2.73. The {sup 192}Ir source was set at fifteen dwell positions, according to clinical practice, in an applicator consisting of a tandem and two ovoids. Dosimetric comparisons were performed for the dose distributions in the water phantom and CT images by using gamma index image and gamma pass rate (%). The gamma index is the minimum Euclidean distance between two 3D spatial dose distributions of the water phantom and CT images in a same space. The gamma pass rates (%) indicate the percentage of agreement points, which mean that two dose distributions are similar, within an acceptance criteria (3 mm/3%). The volumes of physical and clinical interests for the gamma analysis were a whole calculated volume and a region larger than t% of a dose (close to a target), respectively. Results: The gamma pass rates were 77.1% for a whole calculated volume and 92.1% for a region within 1% dose region. The differences of 7.7% to 22.9 % between two dose distributions in the water phantom and CT images were found around the applicator region and near the target. Conclusion: This work revealed the large difference on the dose distributions near the target in the presence of the tissue inhomogeneity. Therefore, the tissue inhomogeneity should be corrected in the dose calculation for clinical treatment.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

WE-AB-BRB-05: Toward a 2D Water-Equivalent Dosimetry Panel Using KCl:Eu2+

Purpose: KCl:Eu2+ storage phosphor shows promise for radiation therapy dosimetry. The purpose of this work is to investigate several important aspects of this material for potential commercial use. Methods: KCl:Eu2+ chips were fabricated and a conformal coating using Parylene was applied. Material’s dose response in a 6 MV beam was investigated using Monte-Carlo simulations. We attempted to micronize the materials using a spiral jet mill. As we did not have a water-free glovebox, we used commercially available non-hygroscopic BaFBr0.85I0.15:Eu2+ computed radiography material to test if a homogeneous panel can be made using micron-sized phosphors. Results: Dosimeters remained intact and showed no change in PSL intensity after eight hrs of submersion in water. We then optically bleached the samples for reuse, irradiated and immersed for another 24 hrs. We observed marginal worsening of the PSL signal for both the soaked and un-soaked chips. By contrast, we were unable to measure PSL intensity of the un-coated pellets due to these pellets dissolving within minutes of being immersed in water. MC data indicate that the micron-sized KCl:Eu2+ is predicted to have a nearly water-equivalent response. KCl:Eu2+ particles with a median size of 3 microns can be produced using a jet mill, which could be reduced further if necessary. While the particles tend to agglomerate over time when stored in a desiccator, they still possess favorable d50’s and d99’s even after 100 minutes, providing an adequate time window for making a panel via tape casting. A panel cast using optimized methods exhibits nearly perfect particle arrangement. Conclusions: Data shown here support ongoing efforts in fabricating a reusable, high resolution dosimetry panel in a water-free glovebox using micron-sized KCl:Eu2+ particles separated by water-equivalent polymers. The conformal coating thereafter will provide good humidity resistance. HL is the founder of DoseImaging, LLC that is exclusively dedicated to commercializing this technology.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-I-38: Patient Organ Specific Dose Assessment in Coronary CT Angiograph Using Voxellaized Volume Dose Index in Monte Carlo Simulation

Purpose: Clinical use of multi detector computed tomography(MDCT) in diagnosis of diseases due to high speed in data acquisition and high spatial resolution is significantly increased. Regarding to the high radiation dose in CT and necessity of patient specific radiation risk assessment, the adoption of new method in the calculation of organ dose is completely required and necessary. In this study by introducing a conversion factor, patient organ dose in thorax region based on CT image data using MC system was calculated. Methods: The geometry of x-ray tube, inherent filter, bow tie filter and collimator were designed using EGSnrc/BEAMnrc MC-system component modules according to GE-Light-speed 64-slices CT-scanner geometry. CT-scan image of patient thorax as a specific phantom was voxellised with 6.25mm3 in voxel and 64×64×20 matrix size. Dose to thorax organ include esophagus, lung, heart, breast, ribs, muscle, spine, spinal cord with imaging technical condition of prospectively-gated-coronary CT-Angiography(PGT) as a step and shoot method, were calculated. Irradiation of patient specific phantom was performed using a dedicated MC-code as DOSXYZnrc with PGT-irradiation model. The ratio of organ dose value calculated in MC-method to the volume CT dose index(CTDIvol) reported by CT-scanner machine according to PGT radiation technique has been introduced as conversion factor. Results: In PGT method, CTDIvol was 10.6mGy and Organ Dose/CTDIvol conversion factor for esophagus, lung, heart, breast, ribs, muscle, spine and spinal cord were obtained as; 0.96, 1.46, 1.2, 3.28. 6.68. 1.35, 3.41 and 0.93 respectively. Conclusion: The results showed while, underestimation of patient dose was found in dose calculation based on CTDIvol, also dose to breast is higher than the other studies. Therefore, the method in this study can be used to provide the actual patient organ dose in CT imaging based on CTDIvol in order to calculation of real effective dose(ED) based on organ dose. This work has been supported by the research chancellor of tehran university of medical sciences(tums), school of medicine, Tehran, Iran.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

Four-dimensional Monte Carlo simulations demonstrating how the extent of intensity-modulation impacts motion effects in proton therapy lung treatments

Purpose: To compare motion effects in intensity modulated proton therapy (IMPT) lung treatments with different levels of intensity modulation.Methods: Spot scanning IMPT treatment plans were generated for ten lung cancer patients for 2.5Gy(RBE) and 12Gy(RBE) fractions and two distinct energy-dependent spot sizes (σ∼8–17 mm and ∼2–4 mm). IMPT plans were generated with the target homogeneity of each individual field restricted to <20% (IMPT{sub 20%}). These plans were compared to full IMPT (IMPT{sub full}), which had no restriction on the single field homogeneity. 4D Monte Carlo simulations were performed upon the patient 4DCT geometry, including deformable image registration and incorporating the detailed timing structure of the proton delivery system. Motion effects were quantified via comparison of the results of the 4D simulations (4D-IMPT{sub 20%}, 4D-IMPT{sub full}) with those of a 3D Monte Carlo simulation (3D-IMPT{sub 20%}, 3D-IMPT{sub full}) upon the planning CT using the equivalent uniform dose (EUD), V{sub 95} and D{sub 1}-D{sub 99}. The effects in normal lung were quantified using mean lung dose (MLD) and V{sub 90%}.Results: For 2.5Gy(RBE), the mean EUD for the large spot size is 99.9%± 2.8% for 4D-IMPT{sub 20%} compared to 100.1%± 2.9% for 4D-IMPT{sub full}. The corresponding values are 88.6%± 8.7% (4D-IMPT{sub 20%}) and 91.0%± 9.3% (4D-IMPT{sub full}) for the smaller spot size. The EUD value is higher in 69.7% of the considered deliveries for 4D-IMPT{sub full}. The V{sub 95} is also higher in 74.7% of the plans for 4D-IMPT{sub full}, implying that IMPT{sub full} plans experience less underdose compared to IMPT{sub 20%}. However, the target dose homogeneity is improved in the majority (67.8%) of plans for 4D-IMPT{sub 20%}. The higher EUD and V{sub 95} suggests that the degraded homogeneity in IMPT{sub full} is actually due to the introduction of hot spots in the target volume, perhaps resulting from the sharper in-target dose gradients. The greatest variations between the IMPT{sub 20%} and IMPT{sub full} deliveries are observed for patients with the largest motion amplitudes. These patients would likely be treated using gating or another motion mitigation technique, which was not the focus of this study.Conclusions: For the treatment parameters considered in this study, the differences between IMPT{sub full} and IMPT{sub 20%} are only likely to be clinically significant for patients with large (>20 mm) motion amplitudes.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2015

Optimization of an on-board imaging system for extremely rapid radiation therapy

Purpose: Next-generation extremely rapid radiation therapy systems could mitigate the need for motion management, improve patient comfort during the treatment, and increase patient throughput for cost effectiveness. Such systems require an on-board imaging system that is competitively priced, fast, and of sufficiently high quality to allow good registration between the image taken on the day of treatment and the image taken the day of treatment planning. In this study, three different detectors for a custom on-board CT system were investigated to select the best design for integration with an extremely rapid radiation therapy system. Methods: Three different CT detectors are proposed: low-resolution (all 4 × 4 mm pixels), medium-resolution (a combination of 4 × 4 mm pixels and 2 × 2 mm pixels), and high-resolution (all 1 × 1 mm pixels). An in-house program was used to generate projection images of a numerical anthropomorphic phantom and to reconstruct the projections into CT datasets, henceforth called “realistic” images. Scatter was calculated using a separate Monte Carlo simulation, and the model included an antiscatter grid and bowtie filter. Diagnostic-quality images of the phantom were generated to represent the patient scan at the time of treatment planning. Commercial deformable registration software was used to register the diagnostic-quality scan to images produced by the various on-board detector configurations. The deformation fields were compared against a “gold standard” deformation field generated by registering initial and deformed images of the numerical phantoms that were used to make the diagnostic and treatment-day images. Registrations of on-board imaging system data were judged by the amount their deformation fields differed from the corresponding gold standard deformation fields—the smaller the difference, the better the system. To evaluate the registrations, the pointwise distance between gold standard and realistic registration deformation fields was computed. Results: By most global metrics (e.g., mean, median, and maximum pointwise distance), the high-resolution detector had the best performance but the medium-resolution detector was comparable. For all medium- and high-resolution detector registrations, mean error between the realistic and gold standard deformation fields was less than 4 mm. By pointwise metrics (e.g., tracking a small lesion), the high- and medium-resolution detectors performed similarly. For these detectors, the smallest error between the realistic and gold standard registrations was 0.6 mm and the largest error was 3.6 mm. Conclusions: The medium-resolution CT detector was selected as the best for an extremely rapid radiation therapy system. In essentially all test cases, data from this detector produced a significantly better registration than data from the low-resolution detector and a comparable registration to data from the high-resolution detector. The medium-resolution detector provides an appropriate compromise between registration accuracy and system cost.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

Proton radiography in three dimensions: A proof of principle of a new technique

Purpose: Monte Carlo simulations were used to investigate a range of phantom configurations to establish enabling three-dimensional proton radiographic techniques.Methods: A large parameter space of stacked phantom geometries composed of tissue inhomogeneity materials such as lung, bone, and cartilage inserted within water background were simulated using a purposefully modified version of TOPAS, an application running on top of the GEANT4 Monte Carlo code. The phantoms were grouped in two classes, one with the inhomogeneity inserted only half-way in the lateral direction and another with complete inhomogeneity insertion. The former class was used to calculate the track count and the energy fluence of the protons as they exit the phantoms either having traversed the inhomogeneity or not. The latter class was used to calculate one yield value accounting for loss of protons due to physical processes only and another yield value accounting for deliberately discarded protons due to large scattering angles. A graphical fingerprinting method was developed to determine the inhomogeneity thickness and location within the phantom based on track count and energy fluence information. Two additional yield values extended this method to the general case which also determines the inhomogeneity material and the phantom thickness.Results: The graphical fingerprinting method was manually validated for two, and automatically tested for all, tissue materials using an exhaustive set of inhomogeneity geometries for 16 cm thick phantoms. Unique recognition of test phantom configurations was achieved in the large majority of cases. The method in the general case was further tested using an exhaustive set of inhomogeneity and phantom tissues and geometries where the phantom thicknesses ranged between 8 and 24 cm. Unique recognition of the test phantom configurations was achieved only for part of the phantom parameter space. The correlations between the remaining false positive recognitions were analyzed.Conclusions: The concept of 3D proton radiography for tissue inhomogeneities of simple geometries was established with the current work. In contrast to conventional 2D proton radiography, the main objective of the demonstrated 3D technique is not proton range. Rather, it is to measure the depth and thickness of an inhomogeneity located in an imaged geometry. Further work is needed to extend and apply the method to more complex geometries.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

Investigation of statistical iterative reconstruction for dedicated breast CT

Purpose: Dedicated breast CT has great potential for improving the detection and diagnosis of breast cancer. Statistical iterative reconstruction (SIR) in dedicated breast CT is a promising alternative to traditional filtered backprojection (FBP). One of the difficulties in using SIR is the presence of free parameters in the algorithm that control the appearance of the resulting image. These parameters require tuning in order to achieve high quality reconstructions. In this study, the authors investigated the penalized maximum likelihood (PML) method with two commonly used types of roughness penalty functions: hyperbolic potential and anisotropic total variation (TV) norm. Reconstructed images were compared with images obtained using standard FBP. Optimal parameters for PML with the hyperbolic prior are reported for the task of detecting microcalcifications embedded in breast tissue.Methods: Computer simulations were used to acquire projections in a half-cone beam geometry. The modeled setup describes a realistic breast CT benchtop system, with an x-ray spectra produced by a point source and an a-Si, CsI:Tl flat-panel detector. A voxelized anthropomorphic breast phantom with 280 μm microcalcification spheres embedded in it was used to model attenuation properties of the uncompressed woman's breast in a pendant position. The reconstruction of 3D images was performed using the separable paraboloidal surrogates algorithm with ordered subsets. Task performance was assessed with the ideal observer detectability index to determine optimal PML parameters.Results: The authors' findings suggest that there is a preferred range of values of the roughness penalty weight and the edge preservation threshold in the penalized objective function with the hyperbolic potential, which resulted in low noise images with high contrast microcalcifications preserved. In terms of numerical observer detectability index, the PML method with optimal parameters yielded substantially improved performance (by a factor of greater than 10) compared to FBP. The hyperbolic prior was also observed to be superior to the TV norm. A few of the best-performing parameter pairs for the PML method also demonstrated superior performance for various radiation doses. In fact, using PML with certain parameter values results in better images, acquired using 2 mGy dose, than FBP-reconstructed images acquired using 6 mGy dose.Conclusions: A range of optimal free parameters for the PML algorithm with hyperbolic and TV norm-based potentials is presented for the microcalcification detection task, in dedicated breast CT. The reported values can be used as starting values of the free parameters, when SIR techniques are used for image reconstruction. Significant improvement in image quality can be achieved by using PML with optimal combination of parameters, as compared to FBP. Importantly, these results suggest improved detection of microcalcifications can be obtained by using PML with lower radiation dose to the patient, than using FBP with higher dose.

62 RADIOLOGY AND NUCLEAR MEDICINE↗