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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

An investigation of kV CBCT image quality and dose reduction for volume-of-interest imaging using dynamic collimation

Purpose: The focus of this work was to investigate the improvements in image quality and dose reduction for volume-of-interest (VOI) kilovoltage-cone beam CT (CBCT) using dynamic collimation. Methods: A prototype iris aperture was used to track a VOI during a CBCT acquisition. The current aperture design is capable of 1D translation as a function of gantry angle and dynamic adjustment of the iris radius. The aperture occupies the location of the bow-tie filter on a Varian On-Board Imager system. CBCT and planar image quality were investigated as a function of aperture radius, while maintaining the same dose to the VOI, for a 20 cm diameter cylindrical water phantom with a 9 mm diameter bone insert centered on isocenter. Corresponding scatter-to-primary ratios (SPR) were determined at the detector plane with Monte Carlo simulation using EGSnrc. Dose distributions for various sizes VOI were modeled using a dynamic BEAMnrc library and DOSXYZnrc. The resulting VOI dose distributions were compared to full-field distributions. Results: SPR was reduced by a factor of 8.4 when decreasing iris diameter from 21.2 to 2.4 cm (at isocenter). Depending upon VOI location and size, dose was reduced to 16%–90% of the full-field value along the central axis plane and down to 4% along the axis of rotation, while maintaining the same dose to the VOI compared to full-field techniques. When maintaining constant dose to the VOI, this change in iris diameter corresponds to a factor increase of approximately 1.6 in image contrast and a factor decrease in image noise of approximately 1.2. This results in a measured gain in contrast-to-noise ratio by a factor of approximately 2.0. Conclusions: The presented VOI technique offers improved image quality for image-guided radiotherapy while sparing the surrounding volume of unnecessary dose compared to full-field techniques.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Spectra of clinical CT scanners using a portable Compton spectrometer

Purpose: Spectral information of the output of x-ray tubes in (dual source) computer tomography (CT) scanners can be used to improve the conversion of CT numbers to proton stopping power and can be used to advantage in CT scanner quality assurance. The purpose of this study is to design, validate, and apply a compact portable Compton spectrometer that was constructed to accurately measure x-ray spectra of CT scanners. Methods: In the design of the Compton spectrometer, the shielding materials were carefully chosen and positioned to reduce background by x-ray fluorescence from the materials used. The spectrum of Compton scattered x-rays alters from the original source spectrum due to various physical processes. Reconstruction of the original x-ray spectrum from the Compton scattered spectrum is based on Monte Carlo simulations of the processes involved. This reconstruction is validated by comparing directly and indirectly measured spectra of a mobile x-ray tube. The Compton spectrometer is assessed in a clinical setting by measuring x-ray spectra at various tube voltages of three different medical CT scanner x-ray tubes. Results: The directly and indirectly measured spectra are in good agreement (their ratio being 0.99) thereby validating the reconstruction method. The measured spectra of the medical CT scanners are consistent with theoretical spectra and spectra obtained from the x-ray tube manufacturer. Conclusions: A Compton spectrometer has been successfully designed, constructed, validated, and applied in the measurement of x-ray spectra of CT scanners. These measurements show that our compact Compton spectrometer can be rapidly set-up using the alignment lasers of the CT scanner, thereby enabling its use in commissioning, troubleshooting, and, e.g., annual performance check-ups of CT scanners.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2015

SU-E-T-335: Dosimetric Investigation of An Advanced Rotating Gamma Ray System for Imaged Guided Radiation Therapy

Purpose: Co-60 beams have unique dosimetric properties for cranial treatments and thoracic cancers. The conventional concern about the high surface dose is overcome by modern system designs with rotational treatment techniques. This work investigates a novel rotational Gamma ray system for image-guided, external beam radiotherapy. Methods: The CybeRT system (Cyber Medical Corp., China) consists of a ring gantry with either one or two treatment heads containing a Gamma source and a multileaf collimator (MLC). The MLC has 60 paired leaves, and the maximum field size is either 40cmx40cm (40 pairs of 0.5cm central leaves, 20 pairs of 1cm outer leaves), or 22cmx40cm (32 pairs of 0.25cm central leaves, 28 pairs of 0.5cm outer leaves). The treatment head(s) can swing 35° superiorly and 8° inferiorly, allowing a total of 43° non-coplanar beam incident. The treatment couch provides 6-degrees-of-freedom motion compensation and the kV cone-beam CT system has a spatial resolution of 0.4mm. Monte Carlo simulations were used to compute dose distributions and compare with measurements. A retrospective study of 98 previously treated patients was performed to compare CybeRT with existing RT systems. Results: Monte Carlo results confirmed the CybeRT design parameters including output factors and 3D dose distributions. Its beam penumbra/dose gradient was similar to or better than that of 6MV photon beams and its isocenter accuracy is 0.3mm. Co-60 beams produce lower-energy secondary electrons that exhibit better dose properties in low-density lung tissues. Because of their rapid depth dose falloff, Co-60 beams are favorable for peripheral lung tumors with half-arc arrangements to spare the opposite lung and critical structures. Superior dose distributions were obtained for head and neck, breast, spine and lung tumors. Conclusion: Because of its accurate dose delivery and unique dosimetric properties of C-60 sources, CybeRT is ideally suited for advanced SBRT as well as conventional RT. This work was partially supported by Cyber Medical Corp.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

TU-FG-BRB-09: Thermoacoustic Range Verification with Perfect Co-Registered Overlay of Bragg Peak onto Ultrasound Image

Purpose: The potential of particle therapy has not yet been fully realized due to inaccuracies in range verification. The purpose of this work was to correlate the Bragg peak location with target structure, by overlaying thermoacoustic localization of the Bragg peak onto an ultrasound image. Methods: Pulsed delivery of 50 MeV protons was accomplished by a fast chopper installed between the ion source and the inflector of the 88″ cyclotron at Lawrence Berkeley National Lab. 2 Gy were delivered in 2 µs by a beam with peak current of 2 µA. Thermoacoustic emissions were detected by a cardiac array and Verasonics V1 ultrasound system, which also generated a grayscale ultrasound image. 1024 thermoacoustic pulses were averaged before filtering and one-way beamforming focused signal onto the Bragg peak location with perfect co-registration to the ultrasound images. Data was collected in a room temperature water bath and gelatin phantom with a cavity designed to mimic the intestine, in which gas pockets can displace the Bragg peak. Experiments were performed with the cavity both empty and filled with olive oil. Results: In the waterbath overlays of the Bragg peak agreed with Monte Carlo simulations to within 800±170 µm. Agreement within 1.3 ± 0.2 mm was achieved in the gelatin phantom, although relative stopping powers were estimated only to first order from CT scans. Protoacoustic signals were detected after travel from the Bragg peak through 29 mm and 65 mm of phantom material when the cavity was empty and full of olive oil, respectively. Conclusion: Protoacoustic range verification is feasible with a commercial clinical ultrasound array, but at doses exceeding the clinical realm. Further optimization of both transducer array and injection line chopper is required to enable range verification within a 2 Gy dose limit, which would enable online adaptive treatment. This work was supported in part by a UWM Intramural Instrumentation Grant and by the Director, Office of Science, Office of Nuclear Physics, of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231. YMQ was supported by a UWM-OUR summer fellowship.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

TH-AB-207A-06: The Use of Realistic Phantoms to Predict CT Dose to Pediatric Patients

Purpose: To predict pediatric patient dose from diagnostic CT scans using Monte Carlo simulation of realistic reference phantoms of various ages, weights, and heights. Methods: A series of deformable pediatric reference phantoms using Non-Uniform Rational B-Splines (NURBS) was developed for a large range of ages, percentiles, and reference anatomy. Individual bones were modeled using age-dependent factors, and red marrow was modeled as functions of age and spatial distribution based on Cristy1. Organ and effective doses for the phantom series were calculated using Monte Carlo simulation of chest, abdominopelvic, and chest-abdomen-pelvis CT exams. Non-linear regression was performed to determine the relationship between dose-length-product (DLP)-normalized organ and effective doses and phantom diameter. Patient-specific voxel computational phantoms were also created by manual segmentation of previously acquired CT images for 40 pediatric patients (0.7 to 17 years). Organ and effective doses were determined by Monte Carlo simulation of these patient-specific phantoms. Each patient was matched to the closest pediatric reference phantom based primarily on age and diameter for all major organs within the torso. Results: A total of 80 NURBS phantoms were created ranging from newborn to 15 years with height/weight percentiles from 10 to 90%. Organ and effective dose normalized by DLP correlated strongly with exponentially decreasing average phantom diameter (R{sup 2} > 0.95 for most organs). A similar relationship was determined for the patient-specific voxel phantoms. Differences between patient-phantom matched organ-dose values ranged from 0.37 to 2.39 mGy (2.87% to 22.1%). Conclusion: Dose estimation using NURBS-based pediatric reference phantoms offers the ability to predict patient dose before and after CT examinations, and physicians and scientists can use this information in their analysis of dose prescriptions for particular subjects and study types. This may lead to practices that minimize radiation dose while still achieving high quality images and, ultimately, improved patient care. NIH/NCI 1 R01 CA155400-01A1.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

A feasibility study of dynamic adaptive radiotherapy for nonsmall cell lung cancer

Purpose: The final state of the tumor at the end of a radiotherapy course is dependent on the doses given in each fraction during the treatment course. This study investigates the feasibility of using dynamic adaptive radiotherapy (DART) in treating lung cancers assuming CBCT is available to observe midtreatment tumor states. DART adapts treatment plans using a dynamic programming technique to consider the expected changes of the tumor in the optimization process. Methods: DART is constructed using a stochastic control formalism framework. It minimizes the total expected number of tumor cells at the end of a treatment course, which is equivalent to maximizing tumor control probability, subject to the uncertainty inherent in the tumor response. This formulation allows for nonstationary dose distributions as well as nonstationary fractional doses as needed to achieve a series of optimal plans that are conformal to the tumor over time, i.e., spatiotemporally optimal plans. Sixteen phantom cases with various sizes and locations of tumors and organs-at-risk (OAR) were generated using in-house software. Each case was planned with DART and conventional IMRT prescribing 60 Gy in 30 fractions. The observations of the change in the tumor volume over a treatment course were simulated using a two-level cell population model. Monte Carlo simulations of the treatment course for each case were run to account for uncertainty in the tumor response. The same OAR dose constraints were applied for both methods. The frequency of replanning was varied between 1, 2, 5 (weekly), and 29 times (daily). The final average tumor dose and OAR doses have been compared to quantify the potential dosimetric benefits of DART. Results: The average tumor max, min, mean, and D95 doses using DART relative to these using conventional IMRT were 124.0%–125.2%, 102.1%–114.7%, 113.7%–123.4%, and 102.0%–115.9% (range dependent on the frequency of replanning). The average relative maximum doses for the cord and esophagus, mean doses for the heart and lungs, and D05 for the unspecified tissue resulting 84%–102.4%, 99.8%–106.9%, 66.9%–85.6%, 58.2%–78.8%, and 85.2%–94.0%, respectively. Conclusions: It is feasible to apply DART to the treatment of NSCLC using CBCT to observe the midtreatment tumor state. Potential increases in the tumor dose and reductions in the OAR dose, particularly for parallel OARs with mean or dose–volume constraints, could be achieved using DART compared to nonadaptive IMRT.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

TU-AB-BRC-10: Modeling of Radiotherapy Linac Source Terms Using ARCHER Monte Carlo Code: Performance Comparison of GPU and MIC Computing Accelerators

Purpose: (1) To perform phase space (PS) based source modeling for Tomotherapy and Varian TrueBeam 6 MV Linacs, (2) to examine the accuracy and performance of the ARCHER Monte Carlo code on a heterogeneous computing platform with Many Integrated Core coprocessors (MIC, aka Xeon Phi) and GPUs, and (3) to explore the software micro-optimization methods. Methods: The patient-specific source of Tomotherapy and Varian TrueBeam Linacs was modeled using the PS approach. For the helical Tomotherapy case, the PS data were calculated in our previous study (Su et al. 2014 41(7) Medical Physics). For the single-view Varian TrueBeam case, we analytically derived them from the raw patient-independent PS data in IAEA’s database, partial geometry information of the jaw and MLC as well as the fluence map. The phantom was generated from DICOM images. The Monte Carlo simulation was performed by ARCHER-MIC and GPU codes, which were benchmarked against a modified parallel DPM code. Software micro-optimization was systematically conducted, and was focused on SIMD vectorization of tight for-loops and data prefetch, with the ultimate goal of increasing 512-bit register utilization and reducing memory access latency. Results: Dose calculation was performed for two clinical cases, a Tomotherapy-based prostate cancer treatment and a TrueBeam-based left breast treatment. ARCHER was verified against the DPM code. The statistical uncertainty of the dose to the PTV was less than 1%. Using double-precision, the total wall time of the multithreaded CPU code on a X5650 CPU was 339 seconds for the Tomotherapy case and 131 seconds for the TrueBeam, while on 3 5110P MICs it was reduced to 79 and 59 seconds, respectively. The single-precision GPU code on a K40 GPU took 45 seconds for the Tomotherapy dose calculation. Conclusion: We have extended ARCHER, the MIC and GPU-based Monte Carlo dose engine to Tomotherapy and Truebeam dose calculations.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Collimator design for a multipinhole brain SPECT insert for MRI

Purpose: Brain single photon emission computed tomography (SPECT) imaging is an important clinical tool, with unique tracers for studying neurological diseases. Nowadays, most commercial SPECT systems are combined with x-ray computed tomography (CT) in so-called SPECT/CT systems to obtain an anatomical background for the functional information. However, while CT images have a high spatial resolution, they have a low soft-tissue contrast, which is an important disadvantage for brain imaging. Magnetic resonance imaging (MRI), on the other hand, has a very high soft-tissue contrast and does not involve extra ionizing radiation. Therefore, the authors designed a brain SPECT insert that can operate inside a clinical MRI. Methods: The authors designed and simulated a compact stationary multipinhole SPECT insert based on digital silicon photomultiplier detector modules, which have shown to be MR-compatible and have an excellent intrinsic resolution (0.5 mm) when combined with a monolithic 2 mm thick LYSO crystal. First, the authors optimized the different parameters of the SPECT system to maximize sensitivity for a given target resolution of 7.2 mm in the center of the field-of-view, given the spatial constraints of the MR system. Second, the authors performed noiseless simulations of two multipinhole configurations to evaluate sampling and reconstructed resolution. Finally, the authors performed Monte Carlo simulations and compared the SPECT insert with a clinical system with ultrahigh-resolution (UHR) fan beam collimators, based on contrast-to-noise ratio and a visual comparison of a Hoffman phantom with a 9 mm cold lesion. Results: The optimization resulted in a stationary multipinhole system with a collimator radius of 150.2 mm and a detector radius of 172.67 mm, which corresponds to four rings of 34 diSPM detector modules. This allows the authors to include eight rings of 24 pinholes, which results in a system volume sensitivity of 395 cps/MBq. Noiseless simulations show sufficient axial sampling (in a Defrise phantom) and a reconstructed resolution of 5.0 mm (in a cold-rod phantom). The authors compared the 24-pinhole setup with a 34-pinhole system (with the same detector radius but a collimator radius of 156.63 mm) and found that 34 pinholes result in better uniformity but a worse reconstruction of the cold-rod phantom. The authors also compared the 24-pinhole system with a clinical triple-head UHR fan beam system based on contrast-to-noise ratio and found that the 24-pinhole setup performs better for the 6 mm hot and the 16 mm cold lesions and worse for the 8 and 10 mm hot lesions. Finally, the authors reconstructed noisy projection data of a Hoffman phantom with a 9 mm cold lesion and found that the lesion was slightly better visible on the multipinhole image compared to the fan beam image. Conclusions: The authors have optimized a stationary multipinhole SPECT insert for MRI and showed the feasibility of doing brain SPECT imaging inside a MRI with an image quality similar to the best clinical SPECT systems available.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-T-529: Dosimetric Investigation of a Rotating Gamma Ray System for ImagedGuided Modulated Arc Radiotherapy

Purpose: Because of their effectiveness and efficiency, rotational arc radiotherapy (MAT) techniques have been developed on both specialty machines such as Tomotherapy and conventional linear accelerators. This work investigates a new rotating Gamma therapy system for image guided MAT and SBRT of intra/extracranial tumors. Methods: The CyberMAT system (Cyber Medical Corp., China) consists of a ring gantry with a gamma source (effective source size 1cm and 1.7cm respectively), a 120-leaf MLC, a kV CBCT and an EPID. The treatment couch provides 6-degrees-of-freedom motion compensation and the kV CBCT system has a spatial resolution of 0.4mm for target localization. The maximum dose rate is >4.0 Gy/min and the maximum field size is 40cm × 40cm. Monte Carlo simulations were used to compute dose distributions and compare with measurements. A retrospective study of 125 previously treated SBRT patients was performed to evaluate the dosimetric characteristics of CyberMAT in comparison with existing VMAT systems. Results: Monte Carlo results confirmed the CyberMAT design parameters including output factors and 3D dose distributions. Its beam penumbra is 5mm to 10mm for field sizes 1cm to 10cm, respectively and its isocenter accuracy is <0.5mm. Compared to the 6 MV photons of Tomotherapy and conventional linacs, Cobalt beams produce lower-energy secondary electrons that exhibit better dose properties in low-density lung tissues. Cobalt beams are ideal for peripheral lung tumors with half-arc arrangements to spare the opposite lung and other critical structures. Superior dose distributions have been obtained for brain, head and neck, breast, spine and lung tumors with half/full arc arrangements. Conclusion: Because of the unique dosimetric properties of Cobalt sources and its accurate stereotaxy/dose delivery CyberMAT is ideally suited for image guided MAT and SBRT. Full-arc arrangements are superior for brain and H&N treatments while half-arc arrangements produce best dose distributions for thoracic tumors.

60 APPLIED LIFE SCIENCES↗
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 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↗