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

SU-E-T-524: In-Vivo Diode Dosimetry Proton Therapy Range Verification Validation Study for Pediatric CSI

Purpose: With the aim of reducing acute esophageal radiation toxicity in pediatric patients receiving craniospinal irradiation (CSI), we investigated the implementation of an in-vivo, adaptive proton therapy range verification methodology. Simulation experiments and in-phantom measurements were conducted to validate the range verification technique for this clinical application. Methods: A silicon diode array system has been developed and experimentally tested in phantom for passively scattered proton beam range verification for a prostate treatment case by correlating properties of the detector signal to the water equivalent path length (WEPL). We propose to extend the methodology to verify range distal to the vertebral body for pediatric CSI cases by placing this small volume dosimeter in the esophagus of the anesthetized patient immediately prior to treatment. A set of calibration measurements was performed to establish a time signal to WEPL fit for a “scout” beam in a solid water phantom. Measurements are compared against Monte Carlo simulation in GEANT4 using the Tool for Particle Simulation (TOPAS). Results: Measurements with the diode array in a spread out Bragg peak of 14 cm modulation width and 15 cm range (177 MeV passively scattered beam) in solid water were successfully validated against proton fluence rate simulations in TOPAS. The resulting calibration curve allows for a sensitivity analysis of detector system response with dose rate in simulation and with individual diode position through simulation on patient CT data. Conclusion: Feasibility has been shown for the application of this range verification methodology to pediatric CSI. An in-vivo measurement to determine the WEPL to the inner surface of the esophagus will allow for personalized adjustment of the treatment plan to ensure sparing of the esophagus while confirming target coverage. A Toltz acknowledges partial support by the CREATE Medical Physics Research Training Network grant of the Natural Sciences and Engineering Research Council (Grant number: 432290)

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

WE-AB-207A-09: Optimization of the Design of a Moving Blocker for Cone-Beam CT Scatter Correction: Experimental Evaluation

Purpose: A moving blocker based strategy has shown promising results for scatter correction in cone-beam computed tomography (CBCT). Different geometry designs and moving speeds of the blocker affect its performance in image reconstruction accuracy. The goal of this work is to optimize the geometric design and moving speed of the moving blocker system through experimental evaluations. Methods: An Elekta Synergy XVI system and an anthropomorphic pelvis phantom CIRS 801-P were used for our experiment. A blocker consisting of lead strips was inserted between the x-ray source and the phantom moving back and forth along rotation axis to measure the scatter signal. Accoriding to our Monte Carlo simulation results, three blockers were used, which have the same lead strip width 3.2mm and different gap between neighboring lead strips, 3.2, 6.4 and 9.6mm. For each blocker, three moving speeds were evaluated, 10, 20 and 30 pixels per projection (on the detector plane). Scatter signal in the unblocked region was estimated by cubic B-spline based interpolation from the blocked region. CBCT image was reconstructed by a total variation (TV) based algebraic iterative reconstruction (ART) algorithm from the partially blocked projection data. Reconstruction accuracy in each condition is quantified as CT number error of region of interest (ROI) by comparing to a CBCT reconstructed image from analytically simulated unblocked and scatter free projection data. Results: Highest reconstruction accuracy is achieved when the blocker width is 3.2 mm, the gap between neighboring lead strips is 9.6 mm and the moving speed is 20 pixels per projection. RMSE of the CT number of ROIs can be reduced from 436 to 27. Conclusions: Image reconstruction accuracy is greatly affected by the geometry design of the blocker. The moving speed does not have a very strong effect on reconstruction result if it is over 20 pixels per projection.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-I-07: Response Characteristics and Signal Conversion Modeling of KV Flat-Panel Detector in Cone Beam CT System

Purpose: The flat-panel detector response characteristics are investigated to optimize the scanning parameter considering the image quality and less radiation dose. The signal conversion model is also established to predict the tumor shape and physical thickness changes. Methods: With the ELEKTA XVI system, the planar images of 10cm water phantom were obtained under different image acquisition conditions, including tube voltage, electric current, exposure time and frames. The averaged responses of square area in center were analyzed using Origin8.0. The response characteristics for each scanning parameter were depicted by different fitting types. The transmission measured for 10cm water was compared to Monte Carlo simulation. Using the quadratic calibration method, a series of variable-thickness water phantoms images were acquired to derive the signal conversion model. A 20cm wedge water phantom with 2cm step thickness was used to verify the model. At last, the stability and reproducibility of the model were explored during a four week period. Results: The gray values of image center all decreased with the increase of different image acquisition parameter presets. The fitting types adopted were linear fitting, quadratic polynomial fitting, Gauss fitting and logarithmic fitting with the fitting R-Square 0.992, 0.995, 0.997 and 0.996 respectively. For 10cm water phantom, the transmission measured showed better uniformity than Monte Carlo simulation. The wedge phantom experiment show that the radiological thickness changes prediction error was in the range of (-4mm, 5mm). The signal conversion model remained consistent over a period of four weeks. Conclusion: The flat-panel response decrease with the increase of different scanning parameters. The preferred scanning parameter combination was 100kV, 10mA, 10ms, 15frames. It is suggested that the signal conversion model could effectively be used for tumor shape change and radiological thickness prediction. Supported by National Natural Science Foundation of China (81101132, 11305203) and Natural Science Foundation of Anhui Province (11040606Q55, 1308085QH138)

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Tomography based numerical simulation of the demagnetizing field in soft magnetic composites

The magneto-static behaviour of soft magnetic composites (SMCs) is investigated using tomography based direct numerical simulation. The microgeometry crucially affects the magnetic properties of the composite since a geometry dependent demagnetizing field is established inside the composite, which lowers the magnetic permeability. We determine the magnetic field information inside the SMC using direct numerical simulation of the magnetic field based on high resolution micro-computed tomography data of the SMC's microstructure as well as artificially generated data made of statistically homogeneous systems of identical fully penetrable spheres and prolate spheroids. Quasi-static electromagnetic behaviour and linear material response are assumed. The 3D magnetostatic Maxwell equations are solved using Whitney finite elements. Simulations show that clustering and percolation behaviour determine the demagnetizing factor of SMCs rather than the particle shape. The demagnetizing factor correlates with the slope of a 2-point probability function at its origin, which is related to the specific surface area of the SMC. Comparison with experimental results indicates that the relatively low permeability of SMCs cannot be explained by demagnetizing effects alone and suggests that the permeability of SMC particles has to be orders of magnitude smaller than the bulk permeability of the particle material.

36 MATERIALS SCIENCE↗
DOE OSTI2014

WE-G-BRE-03: Dose Painting by Numbers Using Targeted Gold Nanoparticles

Purpose: Homogeneous dose enhancement in tumor cells of lung cancer patients treated with conventional dose of 60–66 Gy in five fractions is limited due to increased risk of toxicity to normal structures. Dose painting by numbers (DPBN) is the prescription of a non-uniform radiation dose distribution in the tumor for each voxel based on the intensity level of that voxel obtained from the tumor image. The purpose of this study is to show that DPBN using targeted gold nanoparticles (GNPs) could enhance conventional doses in the more resistant tumor areas. Methods: Cone beam computed tomography (CBCT) images of GNPs after intratumoral injection into human tumor were taken at 0, 48, 144 and 160 hours. The dose enhancement in the tumor voxels by secondary electrons from the GNPs was calculated based on analytical microdosimetry methods. The dose enhancement factor (DEF) is the ratio of the doses to the tumor with and without the presence of GNPs. The DEF was calculated for each voxel of the images based on the GNP concentration in the tumor sub-volumes using 6-MV photon spectra obtained using Monte Carlo simulations at 5 cm depth (10×10 cm2 field). Results: The results revealed DEF values of 1.05–2.38 for GNPs concentrations of 1–30 mg/g which corresponds to 12.60 – 28.56 Gy per fraction for delivering 12 Gy per fraction homogenously to lung tumor region. Conclusion: Our preliminary results verify that DPBN could be achieved using GNPs to enhance conventional doses to high risk tumor sub-volumes. In practice, DPBN using GNPs could be achieved due to diffusion of targeted GNPs sustainably released in-situ from radiotherapy biomaterials (e.g. fiducials) coated with polymer film containing the GNPs.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Numerical prediction of fiber orientation in injection-molded short-fiber/thermoplastic composite parts with experimental validation

Numerical prediction of the fiber orientation in the short-glass fiber (GF) reinforced polyamide 6 (PA6) composites with the fiber weight concentration of 30%, 50%, and 70% manufactured by the injection molding process is presented. And the fiber orientation was also directly observed and measured through X-ray computed tomography. During the injection molding process of the short-fiber/thermoplastic composite, the fiber orientation is produced by the flow states and the fiber-fiber interaction. Folgar and Tucker equation is the well known for modeling the fiber orientation in a concentrated suspension. They included into Jeffrey’s equation a diffusive type of term by introducing a phenomenological coefficient to account for the fiber-fiber interaction. Our developed model for the fiber-fiber interaction was proposed by modifying the rotary diffusion term of the Folgar-Tucker equation. This model was presented in a conference paper of the 29{sup th} International Conference of the Polymer Processing Society published by AIP conference proceeding. For modeling fiber interaction, the fiber dynamic simulation was introduced in order to obtain a global fiber interaction coefficient, which is sum function of the fiber concentration, aspect ratio, and angular velocity. The fiber orientation is predicted by using the proposed fiber interaction model incorporated into a computer aided engineering simulation package C-Mold. An experimental program has been carried out in which the fiber orientation distribution has been measured in 100 x 100 x 2 mm injection-molded plate and 100 x 80 x 2 mm injection-molded weld by analyzed with a high resolution 3D X-ray computed tomography system XVA-160α, and calculated by X-ray computed tomography imaging. The numerical prediction shows a good agreement with experimental validation. And the complex fiber orientation in the injection-molded weld was investigated.

42 ENGINEERING↗
DOE OSTI2015

Degradation of proton depth dose distributions attributable to microstructures in lung-equivalent material

Purpose: The purpose of the work reported here was to investigate the influence of sub-millimeter size heterogeneities on the degradation of the distal edges of proton beams and to validate Monte Carlo (MC) methods’ ability to correctly predict such degradation. Methods: A custom-designed high-resolution plastic phantom approximating highly heterogeneous, lung-like structures was employed in measurements and in Monte Carlo simulations to evaluate the degradation of proton Bragg curves penetrating heterogeneous media. Results: Significant differences in distal falloff widths and in peak dose values were observed in the measured and the Monte Carlo simulated curves compared to pristine proton Bragg curves. Furthermore, differences between simulations of beams penetrating CT images of the phantom did not agree well with the corresponding experimental differences. The distal falloff widths in CT image-based geometries were underestimated by up to 0.2 cm in water (corresponding to 0.8–1.4 cm in lung tissue), and the peak dose values of pristine proton beams were overestimated by as much as ~35% compared to measured curves or depth-dose curves simulated on the basis of true geometry. The authors demonstrate that these discrepancies were caused by the limited spatial resolution of CT images that served as a basis for dose calculations and lead to underestimation of the impact of the fine structure of tissue heterogeneities. A convolution model was successfully applied to mitigate the underestimation. Conclusions: The results of this study justify further development of models to better represent heterogeneity effects in soft-tissue geometries, such as lung, and to correct systematic underestimation of the degradation of the distal edge of proton doses.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Effect of the mitral valve on diastolic flow patterns

The leaflets of the mitral valve interact with the mitral jet and significantly impact diastolic flow patterns, but the effect of mitral valve morphology and kinematics on diastolic flow and its implications for left ventricular function have not been clearly delineated. In the present study, we employ computational hemodynamic simulations to understand the effect of mitral valve leaflets on diastolic flow. A computational model of the left ventricle is constructed based on a high-resolution contrast computed-tomography scan, and a physiological inspired model of the mitral valve leaflets is synthesized from morphological and echocardiographic data. Simulations are performed with a diode type valve model as well as the physiological mitral valve model in order to delineate the effect of mitral-valve leaflets on the intraventricular flow. The study suggests that a normal physiological mitral valve promotes the formation of a circulatory (or “looped”) flow pattern in the ventricle. The mitral valve leaflets also increase the strength of the apical flow, thereby enhancing apical washout and mixing of ventricular blood. The implications of these findings on ventricular function as well as ventricular flow models are discussed.

42 ENGINEERING↗
DOE OSTI2014

SU-C-12A-03: The Impact of Contrast Medium On Radiation Dose in CT: A Systematic Evaluation Across 58 Patient Models

Purpose: To assess the effect of contrast medium on radiation dose as a function of time via Monte Carlo simulation from the liver CT scan across a library of 5D XCAT models Methods: A validated Monte Carlo simulation package (PENELOPE) was employed to model a CT system (LightSpeed 64 VCT, GE Healthcare). The radiation dose was estimated from a common abdomen CT examination. The dose estimation was performed on a library of adult extended cardiac-torso (5D XCAT) phantoms (35 male, 23 female, mean age 51.5 years, mean weight 80.2 kg). The 5D XCAT models were created based on patient-specific iodine concentration-time results from our computational contrast medium propagation model for different intravenous injection protocols. To enable a dynamic estimation of radiation dose, each organ in the model was assigned to its own time-concentration curve via the PENELOPE package, material.exe. Using the Monte Carlo, for each scan time point after the injection, 80 million photons were initiated and tracked through the phantoms. Finally, the dose to the liver was tallied from the deposited energy. Results: Monte Carlo simulation results of radiation dose delivered to the liver from the XCAT models indicated up to 30% increase in dose for different time after the administration of contrast medium. Conclusion: The contrast enhancement is employed in over 60% of imaging modalities, which not only remarkably affects the CT image quality, but also increases the radiation dose by as much as 70%. The postinjection multiple acquisition in several enhanced CT protocols, makes the radiation dose increment through the use of contrast medium, an inevitable factor in optimization of these protocols. The relationship between radiation dose and injected contrast medium as a function of time studied in this work allows optimization of contrast administration for vulnerable individuals.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

TH-AB-207A-07: Radiation Dose Simulation for a Newly Proposed Dynamic Bowtie Filters for CT Using Fast Monte Carlo Methods

Purpose: Dynamic bowtie filter is an innovative design capable of modulating the X-ray and balancing the flux in the detectors, and it introduces a new way of patient-specific CT scan optimizations. This study demonstrates the feasibility of performing fast Monte Carlo dose calculation for a type of dynamic bowtie filter for cone-beam CT (Liu et al. 2014 9(7) PloS one) using MIC coprocessors. Methods: The dynamic bowtie filter in question consists of a highly attenuating bowtie component (HB) and a weakly attenuating bowtie (WB). The HB is filled with CeCl3 solution and its surface is defined by a transcendental equation. The WB is an elliptical cylinder filled with air and immersed in the HB. As the scanner rotates, the orientation of WB remains the same with the static patient. In our Monte Carlo simulation, the HB was approximated by 576 boxes. The phantom was a voxelized elliptical cylinder composed of PMMA and surrounded by air (44cm×44cm×40cm, 1000×1000×1 voxels). The dose to the PMMA phantom was tallied with 0.15% statistical uncertainty under 100 kVp source. Two Monte Carlo codes ARCHER and MCNP-6.1 were compared. Both used double-precision. Compiler flags that may trade accuracy for speed were avoided. Results: The wall time of the simulation was 25.4 seconds by ARCHER on a 5110P MIC, 40 seconds on a X5650 CPU, and 523 seconds by the multithreaded MCNP on the same CPU. The high performance of ARCHER is attributed to the parameterized geometry and vectorization of the program hotspots. Conclusion: The dynamic bowtie filter modeled in this study is able to effectively reduce the dynamic range of the detected signals for the photon-counting detectors. With appropriate software optimization methods, the accelerator-based (MIC and GPU) Monte Carlo dose engines have shown good performance and can contribute to patient-specific CT scan optimizations.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

MO-E-17A-11: In-Phantom Detector Response in a Kilovoltage X-Ray Beam

Purpose: To investigate changes in detector response in a kilovoltage x-ray beam as a function of depth in phantom. Methods: Detector response in a 120kVp x-ray beam as a function of depth in a PMMA phantom was determined for a reference class Exradin A12 ionization chamber, an Exradin A1SL ionization chamber, TLD-100 and EBT3 radiochromic film and normalized to a depth of 1cm. The output of each detector was referenced to the A12 ionization chamber reading to assess changes in detector response as a function of depth. The BEAMnrc, DOSxyz, and egs-chamber user codes of the EGSnrc Monte Carlo code were used to calculate the dose delivered to the detectors from an x-ray source with similar beam quality to the one used in the measurements. Results: Changes in detector response relative to the A12 ionization chamber of 2%, 7%, and 14% for the A1SL, TLDs, and film, respectively, were observed with the larger changes in response occurring at larger depths. The Monte Carlo calculated changes in dose response showed less deviation, with the largest difference being within 2%. Conclusion: The Monte Carlo simulations show that doses delivered to the different detectors are quite similar, while measurements show that the detector output varies drastically. This indicates that the detector output per unit dose is changing with depth, most likely due to changes in the photon spectrum that occur because of scatter and attenuation of the primary beam. Because detector calibration is often performed at a single point in the kilovoltage energy range, the change in detector response needs to be considered when the measurement conditions differ from those of the calibration. This is of particular concern in computed tomography dose measurement, where calibrations are performed at shallow depths or in air, and measurements are made at deeper locations.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2015

Evaluation of the OSC-TV iterative reconstruction algorithm for cone-beam optical CT

Purpose: The present work evaluates an iterative reconstruction approach, namely, the ordered subsets convex (OSC) algorithm with regularization via total variation (TV) minimization in the field of cone-beam optical computed tomography (optical CT). One of the uses of optical CT is gel-based 3D dosimetry for radiation therapy, where it is employed to map dose distributions in radiosensitive gels. Model-based iterative reconstruction may improve optical CT image quality and contribute to a wider use of optical CT in clinical gel dosimetry. Methods: This algorithm was evaluated using experimental data acquired by a cone-beam optical CT system, as well as complementary numerical simulations. A fast GPU implementation of OSC-TV was used to achieve reconstruction times comparable to those of conventional filtered backprojection. Images obtained via OSC-TV were compared with the corresponding filtered backprojections. Spatial resolution and uniformity phantoms were scanned and respective reconstructions were subject to evaluation of the modulation transfer function, image uniformity, and accuracy. The artifacts due to refraction and total signal loss from opaque objects were also studied. Results: The cone-beam optical CT data reconstructions showed that OSC-TV outperforms filtered backprojection in terms of image quality, thanks to a model-based simulation of the photon attenuation process. It was shown to significantly improve the image spatial resolution and reduce image noise. The accuracy of the estimation of linear attenuation coefficients remained similar to that obtained via filtered backprojection. Certain image artifacts due to opaque objects were reduced. Nevertheless, the common artifact due to the gel container walls could not be eliminated. Conclusions: The use of iterative reconstruction improves cone-beam optical CT image quality in many ways. The comparisons between OSC-TV and filtered backprojection presented in this paper demonstrate that OSC-TV can potentially improve the rendering of spatial features and reduce cone-beam optical CT artifacts.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

TH-CD-202-07: A Methodology for Generating Numerical Phantoms for Radiation Therapy Using Geometric Attribute Distribution Models

Purpose: A persistent challenge for the quality assessment of radiation therapy treatments (e.g. contouring accuracy) is the absence of the known, ground truth for patient data. Moreover, assessment results are often patient-dependent. Computer simulation studies utilizing numerical phantoms can be performed for quality assessment with a known ground truth. However, previously reported numerical phantoms do not include the statistical properties of inter-patient variations, as their models are based on only one patient. In addition, these models do not incorporate tumor data. In this study, a methodology was developed for generating numerical phantoms which encapsulate the statistical variations of patients within radiation therapy, including tumors. Methods: Based on previous work in contouring assessment, geometric attribute distribution (GAD) models were employed to model both the deterministic and stochastic properties of individual organs via principle component analysis. Using pre-existing radiation therapy contour data, the GAD models are trained to model the shape and centroid distributions of each organ. Then, organs with different shapes and positions can be generated by assigning statistically sound weights to the GAD model parameters. Organ contour data from 20 retrospective prostate patient cases were manually extracted and utilized to train the GAD models. As a demonstration, computer-simulated CT images of generated numerical phantoms were calculated and assessed subjectively and objectively for realism. Results: A cohort of numerical phantoms of the male human pelvis was generated. CT images were deemed realistic both subjectively and objectively in terms of image noise power spectrum. Conclusion: A methodology has been developed to generate realistic numerical anthropomorphic phantoms using pre-existing radiation therapy data. The GAD models guarantee that generated organs span the statistical distribution of observed radiation therapy patients, according to the training dataset. The methodology enables radiation therapy treatment assessment with multi-modality imaging and a known ground truth, and without patient-dependent bias.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-J-02: 4D Digital Tomosynthesis Based On Algebraic Image Reconstruction and Total-Variation Minimization for the Improvement of Image Quality

Purpose: In this paper, we implemented the four-dimensional (4D) digital tomosynthesis (DTS) imaging based on algebraic image reconstruction technique and total-variation minimization method in order to compensate the undersampled projection data and improve the image quality. Methods: The projection data were acquired as supposed the cone-beam computed tomography system in linear accelerator by the Monte Carlo simulation and the in-house 4D digital phantom generation program. We performed 4D DTS based upon simultaneous algebraic reconstruction technique (SART) among the iterative image reconstruction technique and total-variation minimization method (TVMM). To verify the effectiveness of this reconstruction algorithm, we performed systematic simulation studies to investigate the imaging performance. Results: The 4D DTS algorithm based upon the SART and TVMM seems to give better results than that based upon the existing method, or filtered-backprojection. Conclusion: The advanced image reconstruction algorithm for the 4D DTS would be useful to validate each intra-fraction motion during radiation therapy. In addition, it will be possible to give advantage to real-time imaging for the adaptive radiation therapy. This research was supported by Leading Foreign Research Institute Recruitment Program (Grant No.2009-00420) and Basic Atomic Energy Research Institute (BAERI); (Grant No. 2009-0078390) through the National Research Foundation of Korea(NRF) funded by the Ministry of Science, ICT and Future Planning (MSIP)

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

TH-A-18C-06: A Scatter Elimination Scheme for Cone Beam CT Using An Oscillating Narrow Beam

Purpose: While cone beam CT (CBCT) has been widely used in image guided radiation therapy, its low image quality, primarily caused by scattered x-rays, hinders advanced clinical applications, e.g., CBCT based on-line adaptive re-planning. We propose in this abstract a new scheme called oscillating narrow beam CBCT (ONB-CBCT) to eliminate scatter signals. Methods: ONB-CBCT consists of two major components. 1) Oscillating narrow beam (ONB) scan and 2) partitioned flat panel containing multiple individual detector strips and their own readouts. Both the beam oscillation and detector partition are along the superior-inferior (SI) direction. During data acquisition, at a given projection, the narrow beam sweep through the detector region, and different portions of the detector acquires projection data in synchrony with the narrow beam. ONB can be generated by a rotating slit collimator design with conventional tube with single focal spot, or by directly using a new source with multiple focal spots. A proof-of-principle study via Monte Carlo simulation is conducted to demonstrate the feasibility of ONB-CBCT. Results: As the beam becomes narrower, more and more scatter signals are eliminated. For the case with a bowtie filter and using 15 ONBs, the maximum and the average intensity error due to scatter are below 20 and 10 HU, respectively. Conclusion: ONB yields a narrowed exposure field at each snapshot and hence an inherently negligible scatter effect. Meanwhile, the individualized detector units guarantee high frame rate detection and hence a same large volume coverage as that in conventional CBCT. In summary, ONB-CBCT is a promising design to achieve high-quality CBCT imaging. This study is supported in part by NIH (1R01CA154747-01)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

Radiation doses in volume-of-interest breast computed tomography—A Monte Carlo simulation study

Purpose: Cone beam breast computed tomography (breast CT) with true three-dimensional, nearly isotropic spatial resolution has been developed and investigated over the past decade to overcome the problem of lesions overlapping with breast anatomical structures on two-dimensional mammographic images. However, the ability of breast CT to detect small objects, such as tissue structure edges and small calcifications, is limited. To resolve this problem, the authors proposed and developed a volume-of-interest (VOI) breast CT technique to image a small VOI using a higher radiation dose to improve that region’s visibility. In this study, the authors performed Monte Carlo simulations to estimate average breast dose and average glandular dose (AGD) for the VOI breast CT technique. Methods: Electron–Gamma-Shower system code-based Monte Carlo codes were used to simulate breast CT. The Monte Carlo codes estimated were validated using physical measurements of air kerma ratios and point doses in phantoms with an ion chamber and optically stimulated luminescence dosimeters. The validated full cone x-ray source was then collimated to simulate half cone beam x-rays to image digital pendant-geometry, hemi-ellipsoidal, homogeneous breast phantoms and to estimate breast doses with full field scans. 13-cm in diameter, 10-cm long hemi-ellipsoidal homogeneous phantoms were used to simulate median breasts. Breast compositions of 25% and 50% volumetric glandular fractions (VGFs) were used to investigate the influence on breast dose. The simulated half cone beam x-rays were then collimated to a narrow x-ray beam with an area of 2.5 × 2.5 cm{sup 2} field of view at the isocenter plane and to perform VOI field scans. The Monte Carlo results for the full field scans and the VOI field scans were then used to estimate the AGD for the VOI breast CT technique. Results: The ratios of air kerma ratios and dose measurement results from the Monte Carlo simulation to those from the physical measurements were 0.97 ± 0.03 and 1.10 ± 0.13, respectively, indicating that the accuracy of the Monte Carlo simulation was adequate. The normalized AGD with VOI field scans was substantially reduced by a factor of about 2 over the VOI region and by a factor of 18 over the entire breast for both 25% and 50% VGF simulated breasts compared with the normalized AGD with full field scans. The normalized AGD for the VOI breast CT technique can be kept the same as or lower than that for a full field scan with the exposure level for the VOI field scan increased by a factor of as much as 12. Conclusions: The authors’ Monte Carlo estimates of normalized AGDs for the VOI breast CT technique show that this technique can be used to markedly increase the dose to the breast and thus the visibility of the VOI region without increasing the dose to the breast. The results of this investigation should be helpful for those interested in using VOI breast CT technique to image small calcifications with dose concern.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-T-671: Range-Modulation Effects of Carbon Ion Beams in Lung Tissue

Purpose: When particles traversing inhomogeneous materials like lung they show a characteristic range modulation which cannot be observed in homogeneous materials. It is possible to describe the range modulation by a convolution of an unperturbed Bragg-Curve and a normal distribution. The sigma of the normal distribution is a parameter for the strength of the modulation effect. A new material parameter (modulation power, P-mod) is introduced which is independent of the material thickness. It is defined as the square of sigma divided by the mean water equivalent thickness of the target (µ). Methods: The modulation power of lung tissue was determined by actual Bragg-peak measurements after traversing an ex-vivo porcine lung and by Monte-Carlo simulations with micro-CT data of human lung tissue. The determined modulation powers were used to show the effect of range modulation effects in a simplified treatment situation. A four centimeter spread-out Bragg-peak after traversing eight centimeter of lung tissue was simulated in FLUKA. The SOBP with and without consideration of range modulation effects were compared. Results: As well in the measurements as in the MC simulations range modulation effects of lung tissue were observed. The determined modulation powers showed a great range from 0.05 mm, in the micro-CT data, to 0.7 mm in the lung measurements. The SOBP comparison showed that range modulation effects Result in over- and underdosages at the distal and proximal edge of the SOBP. In the investigated case, the last 0.5 cm of the SOBP showed an underdosage of up to 50% at the distal edge, while 0.5 cm distal to the SOBP an overdosage of up to 50% was observed. Conclusion: Range modulation effects occur in inhomogeneous materials like lung. These modulation effects may Result in clinically relevant over- and underdosages but are currently not considered in commercially available treatment planning systems.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

TU-CD-207-05: A Novel Digital Tomosynthesis System Using Orthogonal Scanning Technique: A Feasibility Study

Purpose: As an alternative to full tomographic imaging technique such as cone-beam computed tomography (CBCT), there is growing interest to adopt digital tomosynthesis (DTS) for the use of diagnostic as well as therapeutic applications. The aim of this study is to propose a new DTS system using novel orthogonal scanning technique, which can provide superior image quality DTS images compared to the conventional DTS scanning system. Methods: Unlike conventional DTS scanning system, the proposed DTS is reconstructed with two sets of orthogonal patient scans. 1) X-ray projections that are acquired along transverse trajectory and 2) an additional sets of X-ray projections acquired along the vertical direction at the mid angle of the previous transverse scan. To reconstruct DTS, we have used modified filtered backprojection technique to account for the different scanning directions of each projection set. We have evaluated the performance of our method using numerical planning CT data of liver cancer patient and a physical pelvis phantom experiment. The results were compared with conventional DTS techniques with single transverse and vertical scanning. Results: The experiments on both numerical simulation as well as physical experiment showed that the resolution as well as contrast of anatomical structures was much clearer using our method. Specifically, the image quality comparing with transversely scanned DTS showed that the edge and contrast of anatomical structures along Left-Right (LR) directions was comparable however, considerable discrepancy and enhancement could be observed along Superior-Inferior (SI) direction using our method. The opposite was observed when vertically scanned DTS was compared. Conclusion: In this study, we propose a novel DTS system using orthogonal scanning technique. The results indicated that the image quality of our novel DTS system was superior compared to conventional DTS system. This makes our DTS system potentially useful in various on-line clinical applications.

60 APPLIED LIFE SCIENCES↗