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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
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 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↗
DOE OSTI2016

SU-F-J-178: A Computer Simulation Model Observer for Task-Based Image Quality Assessment in Radiation Therapy

Purpose: Traditionally, image quality in radiation therapy is assessed subjectively or by utilizing physically-based metrics. Some model observers exist for task-based medical image quality assessment, but almost exclusively for diagnostic imaging tasks. As opposed to disease diagnosis, the task for image observers in radiation therapy is to utilize the available images to design and deliver a radiation dose which maximizes patient disease control while minimizing normal tissue damage. The purpose of this study was to design and implement a new computer simulation model observer to enable task-based image quality assessment in radiation therapy. Methods: A modular computer simulation framework was developed to resemble the radiotherapy observer by simulating an end-to-end radiation therapy treatment. Given images and the ground-truth organ boundaries from a numerical phantom as inputs, the framework simulates an external beam radiation therapy treatment and quantifies patient treatment outcomes using the previously defined therapeutic operating characteristic (TOC) curve. As a preliminary demonstration, TOC curves were calculated for various CT acquisition and reconstruction parameters, with the goal of assessing and optimizing simulation CT image quality for radiation therapy. Sources of randomness and bias within the system were analyzed. Results: The relationship between CT imaging dose and patient treatment outcome was objectively quantified in terms of a singular value, the area under the TOC (AUTOC) curve. The AUTOC decreases more rapidly for low-dose imaging protocols. AUTOC variation introduced by the dose optimization algorithm was approximately 0.02%, at the 95% confidence interval. Conclusion: A model observer has been developed and implemented to assess image quality based on radiation therapy treatment efficacy. It enables objective determination of appropriate imaging parameter values (e.g. imaging dose). Framework flexibility allows for incorporation of additional modules to include any aspect of the treatment process, and therefore has great potential for both assessment and optimization within radiation therapy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Individual pore and interconnection size analysis of macroporous ceramic scaffolds using high-resolution X-ray tomography

The pore interconnection size of β-tricalcium phosphate scaffolds plays an essential role in the bone repair process. Although, the μCT technique is widely used in the biomaterial community, it is rarely used to measure the interconnection size because of the lack of algorithms. In addition, discrete nature of the μCT introduces large systematic errors due to the convex geometry of interconnections. We proposed, verified and validated a novel pore-level algorithm to accurately characterize the individual pores and interconnections. Specifically, pores and interconnections were isolated, labeled, and individually analyzed with high accuracy. The technique was verified thoroughly by visually inspecting and verifying over 3474 properties of randomly selected pores. This extensive verification process has passed a one-percent accuracy criterion. Scanning errors inherent in the discretization, which lead to both dummy and significantly overestimated interconnections, have been examined using computer-based simulations and additional high-resolution scanning. Then accurate correction charts were developed and used to reduce the scanning errors. Only after the corrections, both the μCT and SEM-based results converged, and the novel algorithm was validated. Material scientists with access to all geometrical properties of individual pores and interconnections, using the novel algorithm, will have a more-detailed and accurate description of the substitute architecture and a potentially deeper understanding of the link between the geometric and biological interaction. - Highlights: •An algorithm is developed to analyze individually all pores and interconnections. •After pore isolating, the discretization errors in interconnections were corrected. •Dummy interconnections and overestimated sizes were due to thin material walls. •The isolating algorithm was verified through visual inspection (99% accurate). •After correcting for the systematic errors, algorithm was validated successfully.

36 MATERIALS SCIENCE↗
DOE OSTI2015

GPU-accelerated regularized iterative reconstruction for few-view cone beam CT

Purpose: The present work proposes an iterative reconstruction technique designed for x-ray transmission computed tomography (CT). The main objective is to provide a model-based solution to the cone-beam CT reconstruction problem, yielding accurate low-dose images via few-views acquisitions in clinically acceptable time frames. Methods: The proposed technique combines a modified ordered subsets convex (OSC) algorithm and the total variation minimization (TV) regularization technique and is called OSC-TV. The number of subsets of each OSC iteration follows a reduction pattern in order to ensure the best performance of the regularization method. Considering the high computational cost of the algorithm, it is implemented on a graphics processing unit, using parallelization to accelerate computations. Results: The reconstructions were performed on computer-simulated as well as human pelvic cone-beam CT projection data and image quality was assessed. In terms of convergence and image quality, OSC-TV performs well in reconstruction of low-dose cone-beam CT data obtained via a few-view acquisition protocol. It compares favorably to the few-view TV-regularized projections onto convex sets (POCS-TV) algorithm. It also appears to be a viable alternative to full-dataset filtered backprojection. Execution times are of 1–2 min and are compatible with the typical clinical workflow for nonreal-time applications. Conclusions: Considering the image quality and execution times, this method may be useful for reconstruction of low-dose clinical acquisitions. It may be of particular benefit to patients who undergo multiple acquisitions by reducing the overall imaging radiation dose and associated risks.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Numerical simulation of fiber interaction in short-fiber injection-molded composite using different cavity geometries

The theoretical fiber-interaction model for calculating the fiber orientation in the injection molded short fiber/thermoplastic composite parts was proposed. The proposed model included the fiber dynamics simulation in order to obtain an equation of the global interaction coefficient and accurate estimate of the fiber interacts at all orientation states. The steps to derive the equation for this coefficient in short fiber suspension as a function of the fiber aspect ratio, volume fraction and general shear rate are delineated. Simultaneously, the high-resolution 3D X-ray computed tomography system XVA-160α was used to observe fiber distribution of short-glass-fiber-reinforced polyamide specimens using different cavity geometries. The fiber orientation tensor components are then calculated. Experimental orientation measurements of short-glass-fiber-reinforced polyamide is used to check the ability of present theory for predicting orientation. The experiments and predictions show a quantitative agreement and confirm the basic understanding of fiber orientation in injection-molded composites.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Experimental demonstration of passive acoustic imaging in the human skull cavity using CT-based aberration corrections

Purpose: Experimentally verify a previously described technique for performing passive acoustic imaging through an intact human skull using noninvasive, computed tomography (CT)-based aberration corrections Jones et al. [Phys. Med. Biol. 58, 4981–5005 (2013)]. Methods: A sparse hemispherical receiver array (30 cm diameter) consisting of 128 piezoceramic discs (2.5 mm diameter, 612 kHz center frequency) was used to passively listen through ex vivo human skullcaps (n = 4) to acoustic emissions from a narrow-band fixed source (1 mm diameter, 516 kHz center frequency) and from ultrasound-stimulated (5 cycle bursts, 1 Hz pulse repetition frequency, estimated in situ peak negative pressure 0.11–0.33 MPa, 306 kHz driving frequency) Definity™ microbubbles flowing through a thin-walled tube phantom. Initial in vivo feasibility testing of the method was performed. The performance of the method was assessed through comparisons to images generated without skull corrections, with invasive source-based corrections, and with water-path control images. Results: For source locations at least 25 mm from the inner skull surface, the modified reconstruction algorithm successfully restored a single focus within the skull cavity at a location within 1.25 mm from the true position of the narrow-band source. The results obtained from imaging single bubbles are in good agreement with numerical simulations of point source emitters and the authors’ previous experimental measurements using source-based skull corrections O’Reilly et al. [IEEE Trans. Biomed. Eng. 61, 1285–1294 (2014)]. In a rat model, microbubble activity was mapped through an intact human skull at pressure levels below and above the threshold for focused ultrasound-induced blood–brain barrier opening. During bursts that led to coherent bubble activity, the location of maximum intensity in images generated with CT-based skull corrections was found to deviate by less than 1 mm, on average, from the position obtained using source-based corrections. Conclusions: Taken together, these results demonstrate the feasibility of using the method to guide bubble-mediated ultrasound therapies in the brain. The technique may also have application in ultrasound-based cerebral angiography.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-D-209-06: Study On the Dose Conversion Coefficients in Pediatric Radiography with the Development of Children Voxel Phantoms

Purpose: Conversion coefficients of organ dose normalized to entrance skin dose (ESD) are widely used to evaluate the organ doses directly using ESD without time-consuming dose measurement, this work aims to investigate the dose conversion coefficients in pediatric chest and abdomen radiography with the development of 5 years and 10 years old children voxel phantoms. Methods: After segmentation of organs and tissues from CT slice images of ATOM tissue-equivalent phantoms, a 5-year-old and a 10-year-old children computational voxel phantoms were developed for Monte Carlo simulation. The organ doses and the entrance skin dose for pediatric chest postero-anterior projection and abdominal antero-posterior projection were simulated at the same time, and then the organ dose conversion coefficients were calculated.To verify the simulated results, dose measurement was carried out with ATOM tissue-equivalent phantoms for 5 year chest radiography. Results: Simulated results and experimental results matched very well with each other, the result differences of all the organs covered in radiation field were below 16% for 5-year-old child in chest projection. I showed that the conversion coefficients of organs covered in the radiation field were much larger than organs out of the field for all the study cases, for example, the conversion coefficients of stomach, liver intestines, and pancreas are larger for abdomen radiography while conversion coefficients of lungs are larger for chest radiography. Conclusion: The voxel children phantoms were helpful to evaluate the radiation doses more accurately and efficiently. Radiation field was the essential factor that affects the organ dose, use reasonably small field should be encouraged for radiation protection. This work was supported by the National Natural Science Foundation of China(11475047)

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Monte Carlo simulation of an x-ray luminescence optical tomography scanner prototype

In this work we report the calculation of the deposited energy distribution produced by an x-ray luminescence optical tomography (XLOT) system in a phantom containing different concentrations of Gd{sub 2}O{sub 2}S:Eu nanoparticles. The calculations were performed via Monte Carlo simulation considering spectra from a W target x-ray tube operating between 30 and 90 kVp, with 1.0 mm Al added filtration. CT and XLOT tomographic images were reconstructed from the same data. The results show that XLOT has better detectability than CT alone, that the dose scales linearly with kVp for a fixed concentration of Gd{sub 2}O{sub 2}S:Eu and air-kerma rate, the scattered radiation contribution to the total dose and signal is about 20% and that the dose ratio for a 3 mm diameter insert containing 10 mg/ml Gd{sub 2}O{sub 2}S embedded in a 30 mm diameter water phantom is 6:1. This ratio drops to less than 2:1 for a 1 mg/ml concentration. Finally we show that the method of conjugate images can be used to correct for artifacts due to attenuation effects in XLOT images.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

TU-AB-303-01: A Feasibility Study for Dynamic Adaptive Therapy of Non-Small Cell Lung Cancer

Purpose: To compare plans for NSCLC optimized using Dynamic Adaptive Therapy (DAT) with conventional IMRT optimization. DAT adapts plans based on changes in the target volume by using dynamic programing techniques to consider expected changes into the optimization process. Information gathered during treatment, e.g. from CBCT, is incorporated into the optimization. Methods and materials: DAT is formulated using stochastic control formalism, which minimizes the total expected number of tumor cells at the end of a treatment course subject to uncertainty inherent in the tumor response and organs-at-risk (OAR) dose constraints. This formulation allows for non-stationary dose distribution as well as non-stationary fractional dose as needed to achieve a series of optimal plans that are conformal to tumor over time. Sixteen phantom cases with various sizes and locations of tumors, and OAR geometries were generated. Each case was planned with DAT and conventional IMRT (60Gy/30fx). Tumor volume change over time was obtained by using, daily MVCT-based, two-level cell population model. Monte Carlo simulations have been performed for each treatment course to account for uncertainty in tumor response. Same OAR dose constraints were applied for both methods. The frequency of plan modification was varied to 1, 2, 5 (weekly), and 29 (daily). The final average tumor dose and OAR doses have been compared to quantify the potential benefit of DAT. Results: The average tumor max, min, mean, and D95 resulted from DAT were 124.0–125.2%, 102.1–114.7%, 113.7–123.4%, and 102.0–115.9% (range dependent on the frequency of plan modification) of those from conventional IMRT. Cord max, esophagus max, lung mean, heart mean, and unspecified tissue D05 resulted from AT were 84–102.4%, 99.8–106.9%, 66.9–85.6%, 58.2–78.8%, and 85.2–94.0% of those from conventional IMRT. Conclusions: Significant tumor dose increase and OAR dose reduction, especially with parallel OAR with mean or dose-volume constraints, can be achieved using DAT.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Radiation Exposure of Abdominal Cone Beam Computed Tomography

PurposeTo evaluate patients radiation exposure of abdominal C-arm cone beam computed tomography (CBCT).MethodsThis prospective study was approved by the institutional review board; written, informed consent was waived. Radiation exposure of abdominal CBCT was evaluated in 40 patients who underwent CBCT during endovascular interventions. Dose area product (DAP) of CBCT was documented and effective dose (ED) was estimated based on organ doses using dedicated Monte Carlo simulation software with consideration of X-ray field location and patients’ individual body weight and height. Weight-dependent ED per DAP conversion factors were calculated. CBCT radiation dose was compared to radiation dose of procedural fluoroscopy. CBCT dose-related risk for cancer was assessed.ResultsMean ED of abdominal CBCT was 4.3 mSv (95 % confidence interval [CI] 3.9; 4.8 mSv, range 1.1–7.4 mSv). ED was significantly higher in the upper than in the lower abdomen (p = 0.003) and increased with patients’ weight (r = 0.55, slope = 0.045 mSv/kg, p < 0.001). Radiation exposure of CBCT corresponded to the radiation exposure of on average 7.2 fluoroscopy minutes (95 % CI 5.5; 8.8 min) in the same region of interest. Lifetime risk of exposure related cancer death was 0.033 % or less depending on age and weight.ConclusionsMean ED of abdominal CBCT was 4.3 mSv depending on X-ray field location and body weight.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2016

SU-D-206-07: CBCT Scatter Correction Based On Rotating Collimator

Purpose: Scatter correction in cone-beam computed tomography (CBCT) has obvious effect on the removal of image noise, the cup artifact and the increase of image contrast. Several methods using a beam blocker for the estimation and subtraction of scatter have been proposed. However, the inconvenience of mechanics and propensity to residual artifacts limited the further evolution of basic and clinical research. Here, we propose a rotating collimator-based approach, in conjunction with reconstruction based on a discrete Radon transform and Tchebichef moments algorithm, to correct scatter-induced artifacts. Methods: A rotating-collimator, comprising round tungsten alloy strips, was mounted on a linear actuator. The rotating-collimator is divided into 6 portions equally. The round strips space is evenly spaced on each portion but staggered between different portions. A step motor connected to the rotating collimator drove the blocker to around x-ray source during the CBCT acquisition. The CBCT reconstruction based on a discrete Radon transform and Tchebichef moments algorithm is performed. Experimental studies using water phantom and Catphan504 were carried out to evaluate the performance of the proposed scheme. Results: The proposed algorithm was tested on both the Monte Carlo simulation and actual experiments with the Catphan504 phantom. From the simulation result, the mean square error of the reconstruction error decreases from 16% to 1.18%, the cupping (τcup) from 14.005% to 0.66%, and the peak signal-to-noise ratio increase from 16.9594 to 31.45. From the actual experiments, the induced visual artifacts are significantly reduced. Conclusion: We conducted an experiment on CBCT imaging system with a rotating collimator to develop and optimize x-ray scatter control and reduction technique. The proposed method is attractive in applications where a high CBCT image quality is critical, for example, dose calculation in adaptive radiation therapy. We want to thank Dr. Lei Xing and Dr. Yong Yang in the Stanford University School of Medicine for this work. This work was jointly supported by NSFC (61471226), Natural Science Foundation for Distinguished Young Scholars of Shandong Province (JQ201516), and China Postdoctoral Science Foundation (2015T80739, 2014M551949).

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

TH-EF-207A-05: Feasibility of Applying SMEIR Method On Small Animal 4D Cone Beam CT Imaging

Purpose: Small animal cone beam CT imaging has been widely used in preclinical research. Due to the higher respiratory rate and heat beats of small animals, motion blurring is inevitable and needs to be corrected in the reconstruction. Simultaneous motion estimation and image reconstruction (SMEIR) method, which uses projection images of all phases, proved to be effective in motion model estimation and able to reconstruct motion-compensated images. We demonstrate the application of SMEIR for small animal 4D cone beam CT imaging by computer simulations on a digital rat model. Methods: The small animal CBCT imaging system was simulated with the source-to-detector distance of 300 mm and the source-to-object distance of 200 mm. A sequence of rat phantom were generated with 0.4 mm{sup 3} voxel size. The respiratory cycle was taken as 1.0 second and the motions were simulated with a diaphragm motion of 2.4mm and an anterior-posterior expansion of 1.6 mm. The projection images were calculated using a ray-tracing method, and 4D-CBCT were reconstructed using SMEIR and FDK methods. The SMEIR method iterates over two alternating steps: 1) motion-compensated iterative image reconstruction by using projections from all respiration phases and 2) motion model estimation from projections directly through a 2D-3D deformable registration of the image obtained in the first step to projection images of other phases. Results: The images reconstructed using SMEIR method reproduced the features in the original phantom. Projections from the same phase were also reconstructed using FDK method. Compared with the FDK results, the images from SMEIR method substantially improve the image quality with minimum artifacts. Conclusion: We demonstrate that it is viable to apply SMEIR method to reconstruct small animal 4D-CBCT images.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Ultra Low Level Environmental Neutron Measurements Using Superheated Droplet Detectors

Through the application of superheated droplet detectors (SDDs), the SIMPLE project for the direct search for dark matter (DM) reached the most restrictive limits on the spin-dependent sector to date. The experiment is based on the detection of recoils following WIMP-nuclei interaction, mimicking those from neutron scattering. The thermodynamic operation conditions yield the SDDs intrinsically insensitive to radiations with linear energy transfer below ∼150 keVμm{sup -1} such as photons, electrons, muons and neutrons with energies below ∼40 keV. Underground facilities are increasingly employed for measurements in a low-level radiation background (DM search, gamma-spectroscopy, intrinsic soft-error rate measurements, etc.), where the rock overburden shields against cosmic radiation. In this environment the SDDs are sensitive only to α-particles and neutrons naturally emitted from the surrounding materials. Recently developed signal analysis techniques allow discrimination between neutron and α-induced signals. SDDs are therefore a promising instrument for low-level neutron and α measurements, namely environmental neutron measurements and α-contamination assays. In this work neutron measurements performed in the challenging conditions of the latest SIMPLE experiment (1500 mwe depth with 50-75 cm water shield) are reported. The results are compared with those obtained by detailed Monte Carlo simulations of the neutron background induced by {sup 238}U and {sup 232}Th traces in the facility, shielding and detector materials. Calculations of the neutron energy distribution yield the following neutron fluence rates (in 10{sup -8} cm{sup -2}s{sup -1}): thermal (<0.5 eV): 2.5; epithermal (0.5 eV-100 keV): 2.2; fast (>1 MeV): 3.9. Signal rates were derived using standard cross sections and codes routinely employed in reactor dosimetry. The measured and calculated neutron count rates per unit of active mass were 0.15 ct/kgd and 0.33 ct/kg-d respectively. As the major signal contribution (98%) originates from radio-impurities in the detector container, alternative materials will be employed in future devices. Latest results regarding the improvement of the detector characterization accuracy towards its application in environmental neutron detection are in progress and will be described. (authors)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

SU-D-207-07: Implementation of Full/half Bowtie Filter Model in a Commercial Treatment Planning System for Kilovoltage X-Ray Imaging Dose Estimation

Purpose: To implement full/half bowtie filter models in a commercial treatment planning system (TPS) to calculate kilovoltage (kV) x-ray imaging dose of Varian On-Board Imager (OBI) cone beam CT (CBCT) system. Methods: Full/half bowtie filters of Varian OBI were created as compensator models in Pinnacle TPS (version 9.6) using Matlab software (version 2011a). The profiles of both bowtie filters were acquired from the manufacturer, imported into the Matlab system and hard coded in binary file format. A Pinnacle script was written to import each bowtie filter data into a Pinnacle treatment plan as a compensator. A kV x-ray beam model without including the compensator model was commissioned per each bowtie filter setting based on percent depth dose and lateral profile data acquired from Monte Carlo simulations. To validate the bowtie filter models, a rectangular water phantom was generated in the planning system and an anterior/posterior beam with each bowtie filter was created. Using the Pinnacle script, each bowtie filter compensator was added to the treatment plan. Lateral profile at the depth of 3cm and percent depth dose were measured using an ion chamber and compared with the data extracted from the treatment plans. Results: The kV x-ray beams for both full and half bowtie filter have been modeled in a commercial TPS. The difference of lateral and depth dose profiles between dose calculations and ion chamber measurements were within 6%. Conclusion: Both full/half bowtie filter models provide reasonable results in kV x-ray dose calculations in the water phantom. This study demonstrates the possibility of using a model-based treatment planning system to calculate the kV imaging dose for both full and half bowtie filter modes. Further study is to be performed to evaluate the models in clinical situations.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

TU-CD-207-11: Patient-Driven Automatic Exposure Control for Dedicated Breast CT

Purpose: To implement automatic exposure control (AEC) in dedicated breast CT (bCT) on a patient-specific basis using only the pre-scan scout views. Methods: Using a large cohort (N=153) of bCT data sets, the breast effective diameter (D) and width in orthogonal planes (Wa,Wb) were calculated from the reconstructed bCT image and pre-scan scout views, respectively. D, Wa, and Wb were measured at the breast center-of-mass (COM), making use of the known geometry of our bCT system. These data were then fit to a second-order polynomial “D=F(Wa,Wb)” in a least squares sense in order to provide a functional form for determining the breast diameter. The coefficient of determination (R{sup 2}) and mean percent error between the measured breast diameter and fit breast diameter were used to evaluate the overall robustness of the polynomial fit. Lastly, previously-reported bCT technique factors derived from Monte Carlo simulations were used to determine the tube current required for each breast diameter in order to match two-view mammographic dose levels. Results: F(Wa,Wb) provided fitted breast diameters in agreement with the measured breast diameters resulting in R{sup 2} values ranging from 0.908 to 0.929 and mean percent errors ranging from 3.2% to 3.7%. For all 153 bCT data sets used in this study, the fitted breast diameters ranged from 7.9 cm to 15.7 cm corresponding to tube current values ranging from 0.6 mA to 4.9 mA in order to deliver the same dose as two-view mammography in a 50% glandular breast with a 80 kV x-ray beam and 16.6 second scan time. Conclusion: The present work provides a robust framework for AEC in dedicated bCT using only the width measurements derived from the two orthogonal pre-scan scout views. Future work will investigate how these automatically chosen exposure levels affect the quality of the reconstructed image.

60 APPLIED LIFE SCIENCES↗