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.

SourceResearch connectionFrom yearThrough year

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

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

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 ↗

DOE OSTI2015

Speckle contrast diffuse correlation tomography of complex turbid medium flow

Purpose: Developed herein is a three-dimensional (3D) flow contrast imaging system leveraging advancements in the extension of laser speckle contrast imaging theories to deep tissues along with our recently developed finite-element diffuse correlation tomography (DCT) reconstruction scheme. This technique, termed speckle contrast diffuse correlation tomography (scDCT), enables incorporation of complex optical property heterogeneities and sample boundaries. When combined with a reflectance-based design, this system facilitates a rapid segue into flow contrast imaging of larger, in vivo applications such as humans. Methods: A highly sensitive CCD camera was integrated into a reflectance-based optical system. Four long-coherence laser source positions were coupled to an optical switch for sequencing of tomographic data acquisition providing multiple projections through the sample. This system was investigated through incorporation of liquid and solid tissue-like phantoms exhibiting optical properties and flow characteristics typical of human tissues. Computer simulations were also performed for comparisons. A uniquely encountered smear correction algorithm was employed to correct point-source illumination contributions during image capture with the frame-transfer CCD and reflectance setup. Results: Measurements with scDCT on a homogeneous liquid phantom showed that speckle contrast-based deep flow indices were within 12% of those from standard DCT. Inclusion of a solid phantom submerged below the liquid phantom surface allowed for heterogeneity detection and validation. The heterogeneity was identified successfully by reconstructed 3D flow contrast tomography with scDCT. The heterogeneity center and dimensions and averaged relative flow (within 3%) and localization were in agreement with actuality and computer simulations, respectively. Conclusions: A custom cost-effective CCD-based reflectance 3D flow imaging system demonstrated rapid acquisition of dense boundary data and, with further studies, a high potential for translatability to real tissues with arbitrary boundaries. A requisite correction was also found for measurements in the fashion of scDCT to recover accurate speckle contrast of deep tissues.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2016

SU-F-I-06: Evaluation of Imaging Dose for Modulation Layer Based Dual Energy Cone-Beam CT

Purpose: Dual energy cone beam CT system is finding a variety of promising applications in diagnostic CT, both in imaging of endogenous materials and exogenous materials across a range of body sites. Dual energy cone beam CT system to suggest in this study acquire image by rotating 360 degree with half of the X-ray window covered using copper modulation layer. In the region that covered by modulation layer absorb the low energy X-ray by modulation layer. Relative high energy X-ray passes through the layer and contributes to image reconstruction. Dose evaluation should be carried out in order to utilize such an imaging acquirement technology for clinical use. Methods: For evaluating imaging dose of modulation layer based dual energy cone beam CT system, Prototype cone beam CT that configured X-ray tube (D054SB, Toshiba, Japan) and detector (PaxScan 2520V, Varian Medical Systems, Palo Alto, CA) is used. A range of 0.5–2.0 mm thickness of modulation layer is implemented in Monte Carlo simulation (MCNPX, ver. 2.6.0, Los Alamos National Laboratory, USA) with half of X-ray window covered. In-house phantom using in this study that has 3 cylindrical phantoms configured water, Teflon air with PMMA covered for verifying the comparability the various material in human body and is implemented in Monte Carlo simulation. The actual dose with 2.0 mm copper covered half of X-ray window is measured using Gafchromic EBT3 film with 5.0 mm bolus for compared with simulative dose. Results: Dose in phantom reduced 33% by copper modulation layer of 2.0 mm. Scattering dose occurred in modulation layer by Compton scattering effect is 0.04% of overall dose. Conclusion: Modulation layer of that based dual energy cone beam CT has not influence on unnecessary scatter dose. This study was supported by the Radiation Safety Research Programs (1305033) through the Nuclear Safety and Security Commission.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

Tomography of atomic number and density of materials using dual-energy imaging and the Alvarez and Macovski attenuation model

Dual-energy computed tomography and the Alvarez and Macovski [Phys. Med. Biol. 21, 733 (1976)] transmitted intensity (AMTI) model were used in this study to estimate the maps of density (ρ) and atomic number (Z) of mineralogical samples. In this method, the attenuation coefficients are represented [Alvarez and Macovski, Phys. Med. Biol. 21, 733 (1976)] in the form of the two most important interactions of X-rays with atoms that is, photoelectric absorption (PE) and Compton scattering (CS). This enables material discrimination as PE and CS are, respectively, dependent on the atomic number (Z) and density (ρ) of materials [Alvarez and Macovski, Phys. Med. Biol. 21, 733 (1976)]. Dual-energy imaging is able to identify sample materials even if the materials have similar attenuation coefficients at single-energy spectrum. We use the full model rather than applying one of several applied simplified forms [Alvarez and Macovski, Phys. Med. Biol. 21, 733 (1976); Siddiqui et al., SPE Annual Technical Conference and Exhibition (Society of Petroleum Engineers, 2004); Derzhi, U.S. patent application 13/527,660 (2012); Heismann et al., J. Appl. Phys. 94, 2073–2079 (2003); Park and Kim, J. Korean Phys. Soc. 59, 2709 (2011); Abudurexiti et al., Radiol. Phys. Technol. 3, 127–135 (2010); and Kaewkhao et al., J. Quant. Spectrosc. Radiat. Transfer 109, 1260–1265 (2008)]. This paper describes the tomographic reconstruction of ρ and Z maps of mineralogical samples using the AMTI model. The full model requires precise knowledge of the X-ray energy spectra and calibration of PE and CS constants and exponents of atomic number and energy that were estimated based on fits to simulations and calibration measurements. The estimated ρ and Z images of the samples used in this paper yield average relative errors of 2.62% and 1.19% and maximum relative errors of 2.64% and 7.85%, respectively. Furthermore, we demonstrate that the method accounts for the beam hardening effect in density (ρ) and atomic number (Z) reconstructions to a significant extent.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2015

SU-E-T-401: Evaluation of TG-43 Dose Calculation Accuracy for SAVI-Based Accelerated Partial Breast Irradiation (APBI) Via Monte Carlo Simulations

Purpose: The current standard TG-43 dose calculation method for SAVI-based Accelerated Partial Breast Irradiation (APBI) assumes an ideal geometry of infinite homogeneous water. However, in SAVI treatments, the air cavity inside the device and the short source-to-skin distance raise concerns about the dose accuracy of the TG-43 method. This study is to evaluate TG-43 dose calculation accuracy in SAVI treatments using Monte Carlo (MC) simulations. Methods: We recalculated the dose distributions of 15 APBI patients treated with SAVI devices, including five cases with a size of 6–1, five with 8−1 and five with 10−1, using our in-house developed fast MC dose package for HDR brachytherapy (gBMC). A phase-space file was used to model the Ir-192 HDR source. For each case, the patient CT was converted into a voxelized phantom and the dwell positions and times were extracted from treatment plans for MC dose calculations. Clinically relevant dosimetric parameters of the recalculated dose were compared to those computed via the TG-43 approach. Results: A systematic overestimation of doses was found for the 15 cases in TG-43 results, with D90, V150, and V200 for PTV-eval 2.8±1.8%, 2.0±2.2%, and 1.8±3.5% higher than MC results. TG-43 also overestimated the dose to skin with the maximum dose 4.4±8.4% higher on average. The relatively large standard deviation seen in the difference of maximum skin dose is partially ascribed to the statistical uncertainty of MC simulations when computing the maximum dose. It took gBMC ∼1 minute to compute dose for a SAVI plan. Conclusion: The high efficiency of our gBMC package facilitated the studies with a relatively large number of cases. An overestimation of TG-43 doses was found when using this MC package to recompute doses in SAVI cases. Clinical utilization of TG-43 dose calculation method in this scenario should be aware of this fact.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

Registration of human skull computed tomography data to an ultrasound treatment space using a sparse high frequency ultrasound hemispherical array

Purpose: Transcranial focused ultrasound (FUS) shows great promise for a range of therapeutic applications in the brain. Current clinical investigations rely on the use of magnetic resonance imaging (MRI) to monitor treatments and for the registration of preoperative computed tomography (CT)-data to the MR images at the time of treatment to correct the sound aberrations caused by the skull. For some applications, MRI is not an appropriate choice for therapy monitoring and its cost may limit the accessibility of these treatments. An alternative approach, using high frequency ultrasound measurements to localize the skull surface and register CT data to the ultrasound treatment space, for the purposes of skull-related phase aberration correction and treatment targeting, has been developed. Methods: A prototype high frequency, hemispherical sparse array was fabricated. Pulse-echo measurements of the surface of five ex vivo human skulls were made, and the CT datasets of each skull were obtained. The acoustic data were used to rigidly register the CT-derived skull surface to the treatment space. The ultrasound-based registrations of the CT datasets were compared to the gold-standard landmark-based registrations. Results: The results show on an average sub-millimeter (0.9 ± 0.2 mm) displacement and subdegree (0.8° ± 0.4°) rotation registration errors. Numerical simulations predict that registration errors on this scale will result in a mean targeting error of 1.0 ± 0.2 mm and reduction in focal pressure of 1.0% ± 0.6% when targeting a midbrain structure (e.g., hippocampus) using a commercially available low-frequency brain prototype device (InSightec, 230 kHz brain system). Conclusions: If combined with ultrasound-based treatment monitoring techniques, this registration method could allow for the development of a low-cost transcranial FUS treatment platform to make this technology more widely available.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2019

Computational Investigation of Melt Pool Process Dynamics and Pore Formation in Laser Powder Bed Fusion

In the laser powder bed fusion additive manufacturing process, the presence of porosity may result in cracks and significantly affects the part performance. A comprehensive understanding of the melt pool process dynamics and porosity evolution can help to improve build quality. In this study, a novel multi-physics computational fluid dynamics (CFD) model has been applied to investigate the fluid dynamics in melt pools and resultant pore defects. To accurately capture the melting and solidification process, major process physics, such as the surface tension, evaporation as well as laser multi-reflection, have been considered in the model. A discrete element method is utilized to model the generation of powder spreading upon build plate by additional numerical simulations. Multiple single track experiments have been performed to obtain the melt pool shape and cross-sectional dimension information. The predicted melt pool dimensions were found to have a reasonable agreement with experimental measurements, e.g., the errors are in the range of 1.3 to 10.6% for melt pool width, while they are between 1.4 and 15.9% for melt depth. Pores are captured by both CFD simulation and x-ray computed tomography measurement for the case with a laser power of 350 W and laser speed of 100 mm/s. The formation of keyholes maybe related to the melt pool front wall angle, and it is found that the front wall angle increases with the increase in laser line energy density. In addition, a larger laser power or smaller scanning speed can help to generate keyhole-induced pores; they also contribute to produce larger sized pores.

36 MATERIALS SCIENCE↗

DOE OSTI2015

TU-EF-204-01: Accurate Prediction of CT Tube Current Modulation: Estimating Tube Current Modulation Schemes for Voxelized Patient Models Used in Monte Carlo Simulations

Purpose: Most patient models used in Monte Carlo-based estimates of CT dose, including computational phantoms, do not have tube current modulation (TCM) data associated with them. While not a problem for fixed tube current simulations, this is a limitation when modeling the effects of TCM. Therefore, the purpose of this work was to develop and validate methods to estimate TCM schemes for any voxelized patient model. Methods: For 10 patients who received clinically-indicated chest (n=5) and abdomen/pelvis (n=5) scans on a Siemens CT scanner, both CT localizer radiograph (“topogram”) and image data were collected. Methods were devised to estimate the complete x-y-z TCM scheme using patient attenuation data: (a) available in the Siemens CT localizer radiograph/topogram itself (“actual-topo”) and (b) from a simulated topogram (“sim-topo”) derived from a projection of the image data. For comparison, the actual TCM scheme was extracted from the projection data of each patient. For validation, Monte Carlo simulations were performed using each TCM scheme to estimate dose to the lungs (chest scans) and liver (abdomen/pelvis scans). Organ doses from simulations using the actual TCM were compared to those using each of the estimated TCM methods (“actual-topo” and “sim-topo”). Results: For chest scans, the average differences between doses estimated using actual TCM schemes and estimated TCM schemes (“actual-topo” and “sim-topo”) were 3.70% and 4.98%, respectively. For abdomen/pelvis scans, the average differences were 5.55% and 6.97%, respectively. Conclusion: Strong agreement between doses estimated using actual and estimated TCM schemes validates the methods for simulating Siemens topograms and converting attenuation data into TCM schemes. This indicates that the methods developed in this work can be used to accurately estimate TCM schemes for any patient model or computational phantom, whether a CT localizer radiograph is available or not. Funding Support: NIH Grant R01-EB017095; Disclosures - Michael McNitt-Gray: Institutional Research Agreement, Siemens AG; Research Support, Siemens AG; Consultant, Flaherty Sensabaugh Bonasso PLLC; Consultant, Fulbright and Jaworski; Disclosures - Cynthia McCollough: Research Grant, Siemens Healthcare.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

SU-C-BRB-06: Dosimetric Impact of Breast Contour Reconstruction Errors in GammaPod Stereotactic Radiotherapy

Purpose: The first GammaPod™ unit, a dedicated prone stereotactic treatment device for early stage breast cancer, has been installed and commissioned at University of Maryland School of Medicine. The objective of this study was to investigate potential dosimetric impact of inaccurate breast contour. Methods: In GammaPod treatments, patient’s beast is immobilized by a breast cup device (BCID) throughout the entire same-day imaging and treatment procedure. 28 different BICD sizes are available to accommodate patients with varying breast sizes. A mild suction helps breast tissue to conform to the shape of the cup with selected size. In treatment planning, dose calculation utilizes previously calculated dose distributions for available cup geometry rather than the breast shape from CT image. Patient CT images with breast cups indicate minor geometric discrepancy between the matched shape of the cup and the breast contour, i.e., the contour size is larger or smaller. In order to investigate the dosimetric impact of these discrepancies, we simulated such discrepancies and reassessed the dose to target as well as skin. Results: In vicinity of skin, hot/cold spots were found when matched cup size was smaller/larger than patient’s breast after comparing the corrected dose profiles from Monte Carlo simulation with the planned dose from TPS. The overdosing/underdosing of target could yield point dose differences as large as 5% due to these setup errors (D95 changes within 2.5%). Maximal skin dose was overestimated/underestimated up to 25%/45% when matched cup size was larger/smaller than real breast contour. Conclusion: The dosimetric evaluation suggests substantial underdosing/overdosing with inaccurate cup geometry during planning, which is acceptable for current clinical trial. Further studies are needed to evaluate such impact to treating small volume close to skin.

60 APPLIED LIFE SCIENCES↗

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 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 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 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↗