Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computerized simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

SU-F-T-122: 4Dand 5D Proton Dose Evaluation with Monte Carlo

Purpose: We evaluated uncertainties in therapeutic proton doses of a lung treatment, taking into account intra-fractional geometry changes, such as breathing, and inter-fractional changes, such as tumor shrinkage and weight loss. Methods: A Monte Carlo study was performed using four dimensional CT image sets (4DCTs) and weekly repeat imaging (5DCTs) to compute fixed RBE (1.1) and variable RBE weighted dose in an actual lung treatment geometry. The MC2 Monte Carlo system was employed to simulate proton energy deposition and LET distributions according to a thoracic cancer treatment plan developed with a 3D-CT in a commercial treatment planning system, as well as in each of the phases of 4DCT sets which were recorded weekly throughout the course of the treatment. A cumulative dose distribution in relevant structures was computed and compared to the predictions of the treatment planning system. Results: Using the Monte Carlo method, dose deposition estimates with the lowest possible uncertainties were produced. Comparison with treatment planning predictions indicates that significant uncertainties may be associated with therapeutic lung dose prediction from treatment planning systems, depending on the magnitude of inter- and intra-fractional geometry changes. Conclusion: As this is just a case study, a more systematic investigation accounting for a cohort of patients is warranted; however, this is less practical because Monte Carlo simulations of such cases require enormous computational resources. Hence our study and any future case studies may serve as validation/benchmarking data for faster dose prediction engines, such as the track repeating algorithm, FDC.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Planning Irreversible Electroporation in the Porcine Kidney: Are Numerical Simulations Reliable for Predicting Empiric Ablation Outcomes?

PurposeNumerical simulations are used for treatment planning in clinical applications of irreversible electroporation (IRE) to determine ablation size and shape. To assess the reliability of simulations for treatment planning, we compared simulation results with empiric outcomes of renal IRE using computed tomography (CT) and histology in an animal model.MethodsThe ablation size and shape for six different IRE parameter sets (70–90 pulses, 2,000–2,700 V, 70–100 µs) for monopolar and bipolar electrodes was simulated using a numerical model. Employing these treatment parameters, 35 CT-guided IRE ablations were created in both kidneys of six pigs and followed up with CT immediately and after 24 h. Histopathology was analyzed from postablation day 1.ResultsAblation zones on CT measured 81 ± 18 % (day 0, p ≤ 0.05) and 115 ± 18 % (day 1, p ≤ 0.09) of the simulated size for monopolar electrodes, and 190 ± 33 % (day 0, p ≤ 0.001) and 234 ± 12 % (day 1, p ≤ 0.0001) for bipolar electrodes. Histopathology indicated smaller ablation zones than simulated (71 ± 41 %, p ≤ 0.047) and measured on CT (47 ± 16 %, p ≤ 0.005) with complete ablation of kidney parenchyma within the central zone and incomplete ablation in the periphery.ConclusionBoth numerical simulations for planning renal IRE and CT measurements may overestimate the size of ablation compared to histology, and ablation effects may be incomplete in the periphery.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2016

MO-FG-BRA-01: 4D Monte Carlo Simulations for Verification of Dose Delivered to a Moving Anatomy

Purpose: To validate 4D Monte Carlo (MC) simulations of dose delivery by an Elekta Agility linear accelerator to a moving phantom. Methods: Monte Carlo simulations were performed using the 4DdefDOSXYZnrc/EGSnrc user code which samples a new geometry for each incident particle and calculates the dose in a continuously moving anatomy. A Quasar respiratory motion phantom with a lung insert containing a 3 cm diameter tumor was used for dose measurements on an Elekta Agility linac with the phantom in stationary and moving states. Dose to the center of tumor was measured using calibrated EBT3 film and the RADPOS 4D dosimetry system. A VMAT plan covering the tumor was created on the static CT scan of the phantom using Monaco V.5.10.02. A validated BEAMnrc model of our Elekta Agility linac was used for Monte Carlo simulations on stationary and moving anatomies. To compare the planned and delivered doses, linac log files recorded during measurements were used for the simulations. For 4D simulations, deformation vectors that modeled the rigid translation of the lung insert were generated as input to the 4DdefDOSXYZnrc code as well as the phantom motion trace recorded with RADPOS during the measurements. Results: Monte Carlo simulations and film measurements were found to agree within 2mm/2% for 97.7% of points in the film in the static phantom and 95.5% in the moving phantom. Dose values based on film and RADPOS measurements are within 2% of each other and within 2σ of experimental uncertainties with respect to simulations. Conclusion: Our 4D Monte Carlo simulation using the defDOSXYZnrc code accurately calculates dose delivered to a moving anatomy. Future work will focus on more investigation of VMAT delivery on a moving phantom to improve the agreement between simulation and measurements, as well as establishing the accuracy of our method in a deforming anatomy. This work was supported by the Ontario Consortium of Adaptive Interventions in Radiation Oncology (OCAIRO), funded by the Ontario Research Fund Research Excellence program.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Monte Carlo simulations of the dose from imaging with GE eXplore 120 micro-CT using GATE

Purpose: Small animals are increasingly used as translational models in preclinical imaging studies involving microCT, during which the subjects can be exposed to large amounts of radiation. While the radiation levels are generally sublethal, studies have shown that low-level radiation can change physiological parameters in mice. In order to rule out any influence of radiation on the outcome of such experiments, or resulting deterministic effects in the subjects, the levels of radiation involved need to be addressed. The aim of this study was to investigate the radiation dose delivered by the GE eXplore 120 microCT non-invasively using Monte Carlo simulations in GATE and to compare results to previously obtained experimental values. Methods: Tungsten X-ray spectra were simulated at 70, 80, and 97 kVp using an analytical tool and their half-value layers were simulated for spectra validation against experimentally measured values of the physical X-ray tube. A Monte Carlo model of the microCT system was set up and four protocols that are regularly applied to live animal scanning were implemented. The computed tomography dose index (CTDI) inside a PMMA phantom was derived and multiple field of view acquisitions were simulated using the PMMA phantom, a representative mouse and rat. Results: Simulated half-value layers agreed with experimentally obtained results within a 7% error window. The CTDI ranged from 20 to 56 mGy and closely matched experimental values. Derived organ doses in mice reached 459 mGy in bones and up to 200 mGy in soft tissue organs using the highest energy protocol. Dose levels in rats were lower due to the increased mass of the animal compared to mice. The uncertainty of all dose simulations was below 14%. Conclusions: Monte Carlo simulations proved a valuable tool to investigate the 3D dose distribution in animals from microCT. Small animals, especially mice (due to their small volume), receive large amounts of radiation from the GE eXplore 120 microCT, which might alter physiological parameters in a longitudinal study setup.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

A 3.55 keV line from DM →a→γ: predictions for cool-core and non-cool-core clusters

We further study a scenario in which a 3.55 keV X-ray line arises from decay of dark matter to an axion-like particle (ALP), that subsequently converts to a photon in astrophysical magnetic fields. We perform numerical simulations of Gaussian random magnetic fields with radial scaling of the magnetic field magnitude with the electron density, for both cool-core 'Perseus' and non-cool-core 'Coma' electron density profiles. Using these, we quantitatively study the resulting signal strength and morphology for cool-core and non-cool-core clusters. Our study includes the effects of fields of view that cover only the central part of the cluster, the effects of offset pointings on the radial decline of signal strength and the effects of dividing clusters into annuli. We find good agreement with current data and make predictions for future analyses and observations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

SU-E-T-644: QuAArC: A 3D VMAT QA System Based On Radiochromic Film and Monte Carlo Simulation of Log Files

Purpose: VMAT involves two main sources of uncertainty: one related to the dose calculation accuracy, and the other linked to the continuous delivery of a discrete calculation. The purpose of this work is to present QuAArC, an alternative VMAT QA system to control and potentially reduce these uncertainties. Methods: An automated MC simulation of log files, recorded during VMAT treatment plans delivery, was implemented in order to simulate the actual treatment parameters. The linac head models and the phase-space data of each Control Point (CP) were simulated using the EGSnrc/BEAMnrc MC code, and the corresponding dose calculation was carried out by means of BEAMDOSE, a DOSXYZnrc code modification. A cylindrical phantom was specifically designed to host films rolled up at different radial distances from the isocenter, for a 3D and continuous dosimetric verification. It also allows axial and/or coronal films and point measurements with several types of ion chambers at different locations. Specific software was developed in MATLAB in order to process and evaluate the dosimetric measurements, which incorporates the analysis of dose distributions, profiles, dose difference maps, and 2D/3D gamma index. It is also possible to obtain the experimental DVH reconstructed on the patient CT, by an optimization method to find the individual contribution corresponding to each CP on the film, taking into account the total measured dose, and the corresponding CP dose calculated by MC. Results: The QuAArC system showed high reproducibility of measurements, and consistency with the results obtained with the commercial system implemented in the verification of the evaluated treatment plans. Conclusion: A VMAT QA system based on MC simulation and high resolution dosimetry with film has been developed for treatment verification. It shows to be useful for the study of the real VMAT capabilities, and also for linac commissioning and evaluation of other verification devices.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

Robust primary modulation-based scatter estimation for cone-beam CT

Purpose: Scattered radiation is one of the major problems facing image quality in flat detector cone-beam computed tomography (CBCT). Previously, a new scatter estimation and correction method using primary beam modulation has been proposed. The original image processing technique used a frequency-domain-based analysis, which proved to be sensitive to the accuracy of the modulator pattern both spatially and in amplitude as well as to the frequency of the modulation pattern. In addition, it cannot account for penumbra effects that occur, for example, due to the finite focal spot size and the scatter estimate can be degraded by high-frequency components of the primary image. Methods: In this paper, the authors present a new way to estimate the scatter using primary modulation. It is less sensitive to modulator nonidealities and most importantly can handle arbitrary modulator shapes and changes in modulator attenuation. The main idea is that the scatter estimation can be expressed as an optimization problem, which yields a separation of the scatter and the primary image. The method is evaluated using simulated and experimental CBCT data. The scattering properties of the modulator itself are analyzed using a Monte Carlo simulation. Results: All reconstructions show strong improvements of image quality. To quantify the results, all images are compared to reference images (ideal simulations and collimated scans). Conclusions: The proposed modulator-based scatter reduction algorithm may open the field of flat detector-based imaging to become a quantitative modality. This may have significant impact on C-arm imaging and on image-guided radiation therapy.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

Detectability limitations with 3-D point reconstruction algorithms using digital radiography

The estimated impact of pores in clusters on component fatigue will be highly conservative when based on 2-D rather than 3-D pore positions. To 3-D position and size defects using digital radiography and 3-D point reconstruction algorithms in general require a lower inspection time and in some cases work better with planar geometries than X-ray computed tomography. However, the increase in prior assumptions about the object and the defects will increase the intrinsic uncertainty in the resulting nondestructive evaluation output. In this paper this uncertainty arising when detecting pore defect clusters with point reconstruction algorithms is quantified using simulations. The simulation model is compared to and mapped to experimental data. The main issue with the uncertainty is the possible masking (detectability zero) of smaller defects around some other slightly larger defect. In addition, the uncertainty is explored in connection to the expected effects on the component fatigue life and for different amount of prior object-defect assumptions made.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Technical Note: Phantom study to evaluate the dose and image quality effects of a computed tomography organ-based tube current modulation technique

Purpose: This technical note quantifies the dose and image quality performance of a clinically available organ-dose-based tube current modulation (ODM) technique, using experimental and simulation phantom studies. The investigated ODM implementation reduces the tube current for the anterior source positions, without increasing current for posterior positions, although such an approach was also evaluated for comparison. Methods: Axial CT scans at 120 kV were performed on head and chest phantoms on an ODM-equipped scanner (Optima CT660, GE Healthcare, Chalfont St. Giles, England). Dosimeters quantified dose to breast, lung, heart, spine, eye lens, and brain regions for ODM and 3D-modulation (SmartmA) settings. Monte Carlo simulations, validated with experimental data, were performed on 28 voxelized head phantoms and 10 chest phantoms to quantify organ dose and noise standard deviation. The dose and noise effects of increasing the posterior tube current were also investigated. Results: ODM reduced the dose for all experimental dosimeters with respect to SmartmA, with average dose reductions across dosimeters of 31% (breast), 21% (lung), 24% (heart), 6% (spine), 19% (eye lens), and 11% (brain), with similar results for the simulation validation study. In the phantom library study, the average dose reduction across all phantoms was 34% (breast), 20% (lung), 8% (spine), 20% (eye lens), and 8% (brain). ODM increased the noise standard deviation in reconstructed images by 6%–20%, with generally greater noise increases in anterior regions. Increasing the posterior tube current provided similar dose reduction as ODM for breast and eye lens, increased dose to the spine, with noise effects ranging from 2% noise reduction to 16% noise increase. At noise equal to SmartmA, ODM increased the estimated effective dose by 4% and 8% for chest and head scans, respectively. Increasing the posterior tube current further increased the effective dose by 15% (chest) and 18% (head) relative to SmartmA. Conclusions: ODM reduced dose in all experimental and simulation studies over a range of phantoms, while increasing noise. The results suggest a net dose/noise benefit for breast and eye lens for all studied phantoms, negligible lung dose effects for two phantoms, increased lung dose and/or noise for eight phantoms, and increased dose and/or noise for brain and spine for all studied phantoms compared to the reference protocol.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

WE-AB-303-03: A Scheme for Real-Time, Marker-Less and Low-Dose Tumor Tracking Using Scattering Imaging: Monte Carlo Simulation Study

Purpose: In-treatment tumor localization is critical for the management of tumor motion in lung cancer radiotherapy. Conventional tumor-tracking methods using a kV or MV x-ray projection has limited contrast. To facilitate real-time, marker-less and low-dose in-treatment image tumor tracking, we propose a novel scheme using Compton scatter imaging. This study reports Monte Carlo (MC) simulations on this scheme for the purpose of proof-of-principle. Methods: A slit x-ray beam along the patient superior-inferior (SI) direction is directed to the patient, intersecting the patient lung at a 2D plane containing majority part of the tumor motion trajectory. X-ray photons are scattered due to Compton effect from this plane, which are spatially collimated by, e.g., a pinhole, on one side of the plane and then captured by a detector behind it. The captured image, after correcting for x-ray attenuation and scatter angle variation, reflects the electron density, which allows visualization of the instantaneous anatomy on this plane. We performed MC studies on a phantom and a patient case for the initial test of this proposed method. Results: In the phantom case, the contrast-resolution calculated using tumor/lung as foreground/background for kV fluoroscopy, cone-beam CT, and scattering image were 0.0625, 0.6993, and 0.5290, respectively. In the patient case, tumor motion can be clearly observed in the scatter images. Compared to fluoroscopy, scattering imaging also significantly reduced imaging dose because of its narrower beam design. Conclusion: MC simulation studies demonstrated the potential of the proposed scheme in terms of capturing the instantaneous anatomy of a patient on a 2D plane. Clear visualization of the tumor will probably facilitate ‘marker-less’ and ‘real-time’ tumor tracking with low imaging dose. NIH (1R01CA154747-01, 1R21CA178787-01A1 and 1R21EB017978-01A1)

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Application of Geant4 Monte Carlo simulation in dose calculations for small radiosurgical fields

The Geant4 toolkit was used to develop a Monte Carlo (MC)-based engine for accurate dose calculations in small radiation field sizes. The Geant4 toolkit (version 10.1.p02) was used to simulate 6-MV photon beam of a Varian2100C linear accelerator that is being used for stereotactic radiosurgery (SRS) treatment with small radiation fields. Geometric models of 3 in-house designed radiosurgical divergent cones, with the diameters of their projections at the isocenter being 10, 20, and 30 mm, were simulated. The accuracy of the MC simulation technique was examined by reproducing several different simulated dosimetric parameters of the primary beams with the experimental data. The dose distributions are first checked for single beams for each cone, then standard multiple field (SMF) techniques are applied. A sample set of DICOM files from computed tomography (CT) scan imaging of a patient's head was converted to the Geant4 geometry format to implement MC-based engine for a clinical test. To validate the accuracy of the MC-based calculations for SMF arrangements, the isodose lines from MC simulation in water phantom were compared with the measured isodose lines using EBT3 Gafchromic film in Solid Water phantoms. Agreements between measured and simulated depth dose values and beam profiles for SRS cones were generally within 2%/2 mm. For output factors, the largest discrepancy was observed for 10 mm SRS cone, which was 1.7%. For SMF techniques, in SRS cones, the MC simulation and EBT3 Gafchromic film dosimetry were in acceptable agreement (5%/5 mm). Excellent agreement between the results of the MC-based and measured dose values for both single and SMF techniques in SRS cones indicates the ability of the Geant4 toolkit to be applied as the platform for treatment planning of advanced radiotherapy techniques.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

Analytic reconstruction algorithms for triple-source CT with horizontal data truncation

Purpose: This paper explores a triple-source imaging method with horizontal data truncation to enlarge the field of view (FOV) for big objects. Methods: The study is conducted by using theoretical analysis, mathematical deduction, and numerical simulations. The proposed algorithms are implemented in c + + and MATLAB. While the basic platform is constructed in MATLAB, the computationally intensive segments are coded in c + +, which are linked via a MEX interface. Results: A triple-source circular scanning configuration with horizontal data truncation is developed, where three pairs of x-ray sources and detectors are unevenly distributed on the same circle to cover the whole imaging object. For this triple-source configuration, a fan-beam filtered backprojection-type algorithm is derived for truncated full-scan projections without data rebinning. The algorithm is also extended for horizontally truncated half-scan projections and cone-beam projections in a Feldkamp-type framework. Using their method, the FOV is enlarged twofold to threefold to scan bigger objects with high speed and quality. The numerical simulation results confirm the correctness and effectiveness of the developed algorithms. Conclusions: The triple-source scanning configuration with horizontal data truncation cannot only keep most of the advantages of a traditional multisource system but also cover a larger FOV for big imaging objects. In addition, because the filtering is shift-invariant, the proposed algorithms are very fast and easily parallelized on graphic processing units.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Detecting electron neutrinos from solar dark matter annihilation by JUNO

We explore the electron neutrino signals from light dark matter (DM) annihilation in the Sun for the large liquid scintillator detector JUNO. In terms of the spectrum features of three typical DM annihilation channels χχ→νν-bar,τ{sup +}τ{sup −},bb-bar, we take two sets of selection conditions to calculate the expected signals and atmospheric neutrino backgrounds based on the Monte Carlo simulation data. Then the JUNO sensitivities to the spin independent DM-nucleon and spin dependent DM-proton cross sections are presented. It is found that the JUNO projected sensitivities are much better than the current spin dependent direct detection experimental limits for the νν-bar and τ{sup +}τ{sup −} channels. In the spin independent case, the JUNO will give the better sensitivity to the DM-nucleon cross section than the LUX and CDMSlite limits for the νν-bar channel with the DM mass lighter than 6.5 GeV. If the νν-bar or τ{sup +}τ{sup −} channel is dominant, the future JUNO results are very helpful for us to understand the tension between the DAMA annual modulation signal and other direct detection exclusions.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

Intrinsic Efficiency Calibration Considering Geometric Factors in Gamma-ray Computed Tomography for Radioactive Waste Assay

In radioactive waste assay with gamma-ray computed tomography, calibration for intrinsic efficiency of the system is important to the reconstruction of radioactivity distribution. Due to the geometric characteristics of the system, the non-uniformity of intrinsic efficiency for gamma-rays with different incident positions and directions are often un-negligible. Intrinsic efficiency curves versus geometric parameters of incident gamma-ray are obtained by Monte-Carlo simulation, and two intrinsic efficiency models are suggested to characterize the intrinsic efficiency determined by relative source-detector position and system geometry in the system matrix. Monte-Carlo simulation is performed to compare the different intrinsic efficiency models. Better reconstruction results of radioactivity distribution are achieved by both suggested models than by the uniform intrinsic efficiency model. And compared to model based on detector position, model based on point response increases reconstruction accuracy as well as complexity and time of calculation. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2015

SU-E-J-148: Tools for Development of 4D Proton CT

Purpose: To develop tools for performing 4D proton computed tomography (CT). Methods: A suitable patient with a tumor in the right lower lobe was selected from a set of 4D CT scans. The volumetric CT images formed the basis for calculating the parameters of a breathing model that allows reconstruction of a static reference CT and CT images in each breathing phase. The images were imported into the TOPAS Monte Carlo simulation platform for simulating an experimental proton CT scan with 45 projections spaced by 4 degree intervals. Each projection acquired data for 2 seconds followed by a gantry rotation for 2 seconds without acquisition. The scan covered 180 degrees with individual protons passing through a 9-cm slab of the patient’s lung covering the moving tumor. An initial proton energy sufficient for penetrating the patient from all directions was determined. Performing the proton CT simulation, TOPAS provided output of the proton energy and coordinates registered in two planes before and after the patient, respectively. The set of projection data was then used with an iterative reconstruction algorithm to generate a volumetric proton CT image set of the static reference image and the image obtained under breathing motion, respectively. Results: An initial proton energy of 230 MeV was found to be sufficient, while for an initial energy of 200 MeV a substantial number of protons did not penetrate the patient. The reconstruction of the static reference image set provided sufficient detail for treatment planning. Conclusion: We have developed tools to perform studies of proton CT in the presence of lung motion based on the TOPAS simulation toolkit. This will allow to optimize 4D reconstruction algorithms by synchronizing the acquired proton CT data with a breathing signal and utilizing a breathing model obtained prior to the proton CT scan. This research has been supported by the National Institute Of Biomedical Imaging And Bioengineering of the National Institutes of Health under Award Number R01EB013118.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Generalised form factor dark matter in the Sun

We study the effects of energy transport in the Sun by asymmetric dark matter with momentum and velocity-dependent interactions, with an eye to solving the decade-old Solar Abundance Problem. We study effective theories where the dark matter-nucleon scattering cross-section goes as v{sub rel}{sup 2n} and q{sup 2n} with n=−1,0,1 or 2, where v{sub rel} is the dark matter-nucleon relative velocity and q is the momentum exchanged in the collision. Such cross-sections can arise generically as leading terms from the most basic nonstandard DM-quark operators. We employ a high-precision solar simulation code to study the impact on solar neutrino rates, the sound speed profile, convective zone depth, surface helium abundance and small frequency separations. We find that the majority of models that improve agreement with the observed sound speed profile and depth of the convection zone also reduce neutrino fluxes beyond the level that can be reasonably accommodated by measurement and theory errors. However, a few specific points in parameter space yield a significant overall improvement. A 3–5 GeV DM particle with σ{sub SI}∝q{sup 2} is particularly appealing, yielding more than a 6σ improvement with respect to standard solar models, while being allowed by direct detection and collider limits. We provide full analytical capture expressions for q- and v{sub rel}-dependent scattering, as well as complete likelihood tables for all models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

A Monte Carlo study on {sup 223}Ra imaging for unsealed radionuclide therapy

Purpose: Radium-223 ({sup 223}Ra), an α-emitting radionuclide, is used in unsealed radionuclide therapy for metastatic bone tumors. The demand for qualitative {sup 223}Ra imaging is growing to optimize dosimetry. The authors simulated {sup 223}Ra imaging using an in-house Monte Carlo simulation code and investigated the feasibility and utility of {sup 223}Ra imaging. Methods: The Monte Carlo code comprises two modules, HEXAGON and NAI. The HEXAGON code simulates the photon and electron interactions in the tissues and collimator, and the NAI code simulates the response of the NaI detector system. A 3D numeric phantom created using computed tomography images of a chest phantom was installed in the HEXAGON code. {sup 223}Ra accumulated in a part of the spine, and three x-rays and 19 γ rays between 80 and 450 keV were selected as the emitted photons. To evaluate the quality of the {sup 223}Ra imaging, the authors also simulated technetium-99m ({sup 99m}Tc) imaging under the same conditions and compared the results. Results: The sensitivities of the three photopeaks were 147 counts per unit of source activity (cps MBq{sup −1}; photopeak: 84 keV, full width of energy window: 20%), 166 cps MBq{sup −1} (154 keV, 15%), and 158 cps MBq{sup −1} (270 keV, 10%) for a low-energy general-purpose (LEGP) collimator, and those for the medium-energy general-purpose (MEGP) collimator were 33, 13, and 8.0 cps MBq{sup −1}, respectively. In the case of {sup 99m}Tc, the sensitivity was 55 cps MBq{sup −1} (141 keV, 20%) for LEGP and 52 cps MBq{sup −1} for MEGP. The fractions of unscattered photons of the total photons reflecting the image quality were 0.09 (84 keV), 0.03 (154 keV), and 0.02 (270 keV) for the LEGP collimator and 0.41, 0.25, and 0.50 for the MEGP collimator, respectively. Conversely, this fraction was approximately 0.65 for the simulated {sup 99m}Tc imaging. The sensitivity with the LEGP collimator appeared very high. However, almost all of the counts were because of photons that penetrated or were scattered in the collimator; therefore, the proportions of unscattered photons were small. Conclusions: Their simulation study revealed that the most promising scheme for {sup 223}Ra imaging is an 84-keV window using an MEGP collimator. The sensitivity of the photopeaks above 100 keV is too low for {sup 223}Ra imaging. A comparison of the fractions of unscattered photons reveals that the sensitivity and image quality are approximately two-thirds of those for {sup 99m}Tc imaging.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

A promising limited angular computed tomography reconstruction via segmentation based regional enhancement and total variation minimization

X-ray computed tomography (CT) is a powerful and common inspection technique used for the industrial non-destructive testing. However, large-sized and heavily absorbing objects cause the formation of artifacts because of either the lack of specimen penetration in specific directions or the acquisition of data from only a limited angular range of views. Although the sparse optimization-based methods, such as the total variation (TV) minimization method, can suppress artifacts to some extent, reconstructing the images such that they converge to accurate values remains difficult because of the deficiency in continuous angular data and inconsistency in the projections. To address this problem, we use the idea of regional enhancement of the true values and suppression of the illusory artifacts outside the region to develop an efficient iterative algorithm. This algorithm is based on the combination of regional enhancement of the true values and TV minimization for the limited angular reconstruction. In this algorithm, the segmentation approach is introduced to distinguish the regions of different image knowledge and generate the support mask of the image. A new regularization term, which contains the support knowledge to enhance the true values of the image, is incorporated into the objective function. Then, the proposed optimization model is solved by variable splitting and the alternating direction method efficiently. A compensation approach is also designed to extract useful information from the initial projections and thus reduce false segmentation result and correct the segmentation support and the segmented image. The results obtained from comparing both simulation studies and real CT data set reconstructions indicate that the proposed algorithm generates a more accurate image than do the other reconstruction methods. The experimental results show that this algorithm can produce high-quality reconstructed images for the limited angular reconstruction and suppress the illusory artifacts caused by the deficiency in valid data.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗