Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Iterative”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

SU-F-18C-02: Evaluations of the Noise Power Spectrum of a CT Iterative Reconstruction Technique for Radiation Therapy

Purpose: To quantitatively assess the noise power spectrum (NPS) of the new, commercially released CT iterative reconstruction technique, iDose{sup 4} from Philips, to compare it with filtered back-projection techniques (FBP), and to provide clinical practice suggestions for radiation therapy. Methods: A uniform phantom was CT imaged with 120kVp tube potential over a range of mAs (250-3333). The image sets were reconstructed using two reconstruction algorithms (FBP and iDose{sup 4} with noise reduction levels 1, 3, and 6) and three reconstruction filters (standard B, smooth A, and sharp C), after which NPS variations were analyzed and compared on region of interest (ROI) sizes (16×16 to 128×128 pixels), ROI radii (0–65 mm), reconstruction algorithms, reconstruction filters, and mAs. Results: The NPS magnitude and shape depended considerably on ROI size and location for both reconstruction algorithms. Regional noise variance became more stationary as ROI size decreased, minimizing NPS artifacts. The optimal 32×32-pixel ROI size balanced the trade-off between stationary noise and adequate sampling. NPS artifacts were greatest at the center of reconstruction space and decreased with increasing ROI distance from the center. The optimal ROI position was located near the phantom's radial midpoint (∼40mm). For sharper filters, the NPS magnitude and the maximum magnitude frequency increased. Higher dose scans yielded lower NPS magnitudes for both reconstruction algorithms and all filters. Compared to FBP, the iDose{sup 4} algorithm reduced the NPS magnitude while preferentially reducing noise at mid-range spatial frequencies, altering noise texture. This reduction was more significant with increasing iDose{sup 4} noise reduction level. Conclusion: Compared to pixel standard deviation, NPS has greater clinical potential for task-based image quality assessment, describing both the magnitude and spatial frequency characteristics of image noise. While iDose{sup 4} reduces the NPS magnitude more dramatically than FBP, further studies are needed to assess the impacts of noise texture alteration on its clinical usage.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

SU-F-18C-06: Prospective Patient Evaluation of Iterative Reconstruction in Radiation Oncology

Purpose: This work incorporates iterative reconstruction (IR) into a dose reduction study to characterize image quality metrics, delineation, and dosimetric assessment, with the goal of reducing imaging dose in Radiation Oncology. Methods: Three-dimensional noise power spectrum (NPS) analysis characterized noise magnitude/texture (120 kVp, 50–200 mAs, IR levels 1–6 yielding noise reduction of 0.89–0.55 compared to filtered backprojection (FBP)). Task-specific Modulation Transfer Functions (MTFtask) were characterized across varied subject contrasts. A prospective dose reduction study (500 to 150 mAs) was conducted for 12 patients (43 inter-fraction CTs) for high-dose rate brachytherapy. Three physicians performed qualitative image assessment between full-dose FBP (FD-FBP, 500 mAs), low-dose FBP (LD-FBP, 150–250 mAs), and low-dose IRL5-6 (LD-IR) scans for image noise, cuff/bladder interface detectability, spatial resolution, texture, and segmentation confidence. Comparisons between LD-FBP and LD-IR were conducted for the following metrics: delineation (bladder and rectum evaluated via overlap indices (OI) and Dice similarity coefficients (DSC)), noise, boundary changes, dose calculation, and digitally reconstructed radiographs (DRRs). Results: NPS showed ∼50% reduction in noise magnitude and ∼0.1 1/mm spatial frequency shift with IRL6. The largest MTFtask decrease between FBP and IR was 0.08 A.U. Qualitative patient image evaluation revealed LD-IR was equivalent or slightly worse than FD-FBP, and superior to LD-FBP for all metrics except low contrast interface and texture. The largest CT number discrepancy from FBP occurred at a bone/tissue interface using IRL6 (−1.2 ± 4.9 HU (range: −17.6 – 12.5 HU)). No significant contour differences (OIs and DSCs = 0.85 – 0.95) and dose calculation discrepancy (<0.02%) were observed. DRRs preserved anatomical detail and demonstrated <2% difference in intensity between LD-FBP and LD-IRL6. Conclusion: While phantom analysis showed slight noise texture differences with IR, patient results revealed that image quality, contouring ability, and dosimetric parameters were not adversely affected, thus support integrating IR into treatment planning. Research supported in part by a grant from Philips HealthCare.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

WE-D-18A-04: How Iterative Reconstruction Algorithms Affect the MTFs of Variable-Contrast Targets in CT Images

Purpose: To determine how filtered back-projection (FBP), adaptive statistical (ASiR), and model based (MBIR) iterative reconstruction algorithms affect the measured modulation transfer functions (MTFs) of variable-contrast targets over a wide range of clinically applicable dose levels. Methods: The Catphan 600 CTP401 module, surrounded by an oval, fat-equivalent ring to mimic patient size/shape, was scanned on a GE HD750 CT scanner at 1, 2, 3, 6, 12 and 24 mGy CTDIvol levels with typical patient scan parameters: 120kVp, 0.8s, 40mm beam width, large SFOV, 2.5mm thickness, 0.984 pitch. The images were reconstructed using GE's Standard kernel with FBP; 20%, 40% and 70% ASiR; and MBIR. A task-based MTF (MTFtask) was computed for six cylindrical targets: 2 low-contrast (Polystyrene, LDPE), 2 medium-contrast (Delrin, PMP), and 2 high-contrast (Teflon, air). MTFtask was used to compare the performance of reconstruction algorithms with decreasing CTDIvol from 24mGy, which is currently used in the clinic. Results: For the air target and 75% dose savings (6 mGy), MBIR MTFtask at 5 lp/cm measured 0.24, compared to 0.20 for 70% ASiR and 0.11 for FBP. Overall, for both high-contrast targets, MBIR MTFtask improved with increasing CTDIvol and consistently outperformed ASiR and FBP near the system's Nyquist frequency. Conversely, for Polystyrene at 6 mGy, MBIR (0.10) and 70% ASiR (0.07) MTFtask was lower than for FBP (0.18). For medium and low-contrast targets, FBP remains the best overall algorithm for improved resolution at low CTDIvol (1–6 mGy) levels, whereas MBIR is comparable at higher dose levels (12–24 mGy). Conclusion: MBIR improved the MTF of small, high-contrast targets compared to FBP and ASiR at doses of 50%–12.5% of those currently used in the clinic. However, for imaging low- and mediumcontrast targets, FBP performed the best across all dose levels. For assessing MTF from different reconstruction algorithms, task-based MTF measurements are necessary.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

WE-G-18A-03: Cone Artifacts Correction in Iterative Cone Beam CT Reconstruction

Purpose: For iterative reconstruction (IR) in cone-beam CT (CBCT) imaging, data truncation along the superior-inferior (SI) direction causes severe cone artifacts in the reconstructed CBCT volume images. Not only does it reduce the effective SI coverage of the reconstructed volume, it also hinders the IR algorithm convergence. This is particular a problem for regularization based IR, where smoothing type regularization operations tend to propagate the artifacts to a large area. It is our purpose to develop a practical cone artifacts correction solution. Methods: We found it is the missing data residing in the truncated cone area that leads to inconsistency between the calculated forward projections and measured projections. We overcome this problem by using FDK type reconstruction to estimate the missing data and design weighting factors to compensate the inconsistency caused by the missing data. We validate the proposed methods in our multi-GPU low-dose CBCT reconstruction system on multiple patients' datasets. Results: Compared to the FDK reconstruction with full datasets, while IR is able to reconstruct CBCT images using a subset of projection data, the severe cone artifacts degrade overall image quality. For head-neck case under a full-fan mode, 13 out of 80 slices are contaminated. It is even more severe in pelvis case under half-fan mode, where 36 out of 80 slices are affected, leading to inferior soft-tissue delineation. By applying the proposed method, the cone artifacts are effectively corrected, with a mean intensity difference decreased from ∼497 HU to ∼39HU for those contaminated slices. Conclusion: A practical and effective solution for cone artifacts correction is proposed and validated in CBCT IR algorithm. This study is supported in part by NIH (1R01CA154747-01)

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

WE-G-18A-04: 3D Dictionary Learning Based Statistical Iterative Reconstruction for Low-Dose Cone Beam CT Imaging

Purpose: To develop a 3D dictionary learning based statistical reconstruction algorithm on graphic processing units (GPU), to improve the quality of low-dose cone beam CT (CBCT) imaging with high efficiency. Methods: A 3D dictionary containing 256 small volumes (atoms) of 3x3x3 voxels was trained from a high quality volume image. During reconstruction, we utilized a Cholesky decomposition based orthogonal matching pursuit algorithm to find a sparse representation on this dictionary basis of each patch in the reconstructed image, in order to regularize the image quality. To accelerate the time-consuming sparse coding in the 3D case, we implemented our algorithm in a parallel fashion by taking advantage of the tremendous computational power of GPU. Evaluations are performed based on a head-neck patient case. FDK reconstruction with full dataset of 364 projections is used as the reference. We compared the proposed 3D dictionary learning based method with a tight frame (TF) based one using a subset data of 121 projections. The image qualities under different resolutions in z-direction, with or without statistical weighting are also studied. Results: Compared to the TF-based CBCT reconstruction, our experiments indicated that 3D dictionary learning based CBCT reconstruction is able to recover finer structures, to remove more streaking artifacts, and is less susceptible to blocky artifacts. It is also observed that statistical reconstruction approach is sensitive to inconsistency between the forward and backward projection operations in parallel computing. Using high a spatial resolution along z direction helps improving the algorithm robustness. Conclusion: 3D dictionary learning based CBCT reconstruction algorithm is able to sense the structural information while suppressing noise, and hence to achieve high quality reconstruction. The GPU realization of the whole algorithm offers a significant efficiency enhancement, making this algorithm more feasible for potential clinical application. A high zresolution is preferred to stabilize statistical iterative reconstruction. This work was supported in part by NIH(1R01CA154747-01), NSFC((No. 61172163), Research Fund for the Doctoral Program of Higher Education of China (No. 20110201110011), China Scholarship Council.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

TH-C-18A-01: Is Automatic Tube Current Modulation Still Necessary with Statistical Iterative Reconstruction?

Purpose: Automatic tube current modulation (TCM) has been widely used in modern multi-detector CT to reduce noise spatial nonuniformity and streaks to improve dose efficiency. With the advent of statistical iterative reconstruction (SIR), it is expected that the importance of TCM may diminish, since SIR incorporates statistical weighting factors to reduce the negative influence of photon-starved rays. The purpose of this work is to address the following questions: Does SIR offer the same benefits as TCM? If yes, are there still any clinical benefits to using TCM? Methods: An anthropomorphic CIRS chest phantom was scanned using a state-of-the-art clinical CT system equipped with an SIR engine (Veo™, GE Healthcare). The phantom was first scanned with TCM using a routine protocol and a low-dose (LD) protocol. It was then scanned without TCM using the same protocols. For each acquisition, both FBP and Veo reconstructions were performed. All scans were repeated 50 times to generate an image ensemble from which noise spatial nonuniformity (NSN) and streak artifact levels were quantified. Monte-Carlo experiments were performed to estimate skin dose. Results: For FBP, noise streaks were reduced by 4% using TCM for both routine and LD scans. NSN values were actually slightly higher with TCM (0.25) than without TCM (0.24) for both routine and LD scans. In contrast, for Veo, noise streaks became negligible (<1%) with or without TCM for both routine and LD scans, and the NSN was reduced to 0.10 (low dose) or 0.08 (routine). The overall skin dose was 2% lower at the shoulders and more uniformly distributed across the skin without TCM. Conclusion: SIR without TCM offers superior reduction in noise nonuniformity and streaks relative to FBP with TCM. For some clinical applications in which skin dose may be a concern, SIR without TCM may be a better option. K. Li, W. Zhao, D. Gomez-Cardona: Nothing to disclose; G.-H. Chen: Research funded, General Electric Company Research funded, Siemens AG Research funded, Varian Medical Systems, Research funded, Hologic, Inc.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Towards the clinical implementation of iterative low-dose cone-beam CT reconstruction in image-guided radiation therapy: Cone/ring artifact correction and multiple GPU implementation

Purpose: Compressed sensing (CS)-based iterative reconstruction (IR) techniques are able to reconstruct cone-beam CT (CBCT) images from undersampled noisy data, allowing for imaging dose reduction. However, there are a few practical concerns preventing the clinical implementation of these techniques. On the image quality side, data truncation along the superior–inferior direction under the cone-beam geometry produces severe cone artifacts in the reconstructed images. Ring artifacts are also seen in the half-fan scan mode. On the reconstruction efficiency side, the long computation time hinders clinical use in image-guided radiation therapy (IGRT). Methods: Image quality improvement methods are proposed to mitigate the cone and ring image artifacts in IR. The basic idea is to use weighting factors in the IR data fidelity term to improve projection data consistency with the reconstructed volume. In order to improve the computational efficiency, a multiple graphics processing units (GPUs)-based CS-IR system was developed. The parallelization scheme, detailed analyses of computation time at each step, their relationship with image resolution, and the acceleration factors were studied. The whole system was evaluated in various phantom and patient cases. Results: Ring artifacts can be mitigated by properly designing a weighting factor as a function of the spatial location on the detector. As for the cone artifact, without applying a correction method, it contaminated 13 out of 80 slices in a head-neck case (full-fan). Contamination was even more severe in a pelvis case under half-fan mode, where 36 out of 80 slices were affected, leading to poorer soft tissue delineation and reduced superior–inferior coverage. The proposed method effectively corrects those contaminated slices with mean intensity differences compared to FDK results decreasing from ∼497 and ∼293 HU to ∼39 and ∼27 HU for the full-fan and half-fan cases, respectively. In terms of efficiency boost, an overall 3.1 × speedup factor has been achieved with four GPU cards compared to a single GPU-based reconstruction. The total computation time is ∼30 s for typical clinical cases. Conclusions: The authors have developed a low-dose CBCT IR system for IGRT. By incorporating data consistency-based weighting factors in the IR model, cone/ring artifacts can be mitigated. A boost in computational efficiency is achieved by multi-GPU implementation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

Clinical evaluation of the iterative metal artifact reduction algorithm for CT simulation in radiotherapy

Purpose: To clinically evaluate an iterative metal artifact reduction (IMAR) algorithm prototype in the radiation oncology clinic setting by testing for accuracy in CT number retrieval, relative dosimetric changes in regions affected by artifacts, and improvements in anatomical and shape conspicuity of corrected images. Methods: A phantom with known material inserts was scanned in the presence/absence of metal with different configurations of placement and sizes. The relative change in CT numbers from the reference data (CT with no metal) was analyzed. The CT studies were also used for dosimetric tests where dose distributions from both photon and proton beams were calculated. Dose differences and gamma analysis were calculated to quantify the relative changes between doses calculated on the different CT studies. Data from eight patients (all different treatment sites) were also used to quantify the differences between dose distributions before and after correction with IMAR, with no reference standard. A ranking experiment was also conducted to analyze the relative confidence of physicians delineating anatomy in the near vicinity of the metal implants. Results: IMAR corrected images proved to accurately retrieve CT numbers in the phantom study, independent of metal insert configuration, size of the metal, and acquisition energy. For plastic water, the mean difference between corrected images and reference images was −1.3 HU across all scenarios (N = 37) with a 90% confidence interval of [−2.4, −0.2] HU. While deviations were relatively higher in images with more metal content, IMAR was able to effectively correct the CT numbers independent of the quantity of metal. Residual errors in the CT numbers as well as some induced by the correction algorithm were found in the IMAR corrected images. However, the dose distributions calculated on IMAR corrected images were closer to the reference data in phantom studies. Relative spatial difference in the dose distributions in the regions affected by the metal artifacts was also observed in patient data. However, in absence of a reference ground truth (CT set without metal inserts), these differences should not be interpreted as improvement/deterioration of the accuracy of calculated dose. With limited data presented, it was observed that proton dosimetry was affected more than photons as expected. Physicians were significantly more confident contouring anatomy in the regions affected by artifacts. While site specific preferences were detected, all indicated that they would consistently use IMAR corrected images. Conclusions: IMAR correction algorithm could be readily implemented in an existing clinical workflow upon commercial release. While residual errors still exist in IMAR corrected images, these images present with better overall conspicuity of the patient/phantom geometry and offer more accurate CT numbers for improved local dosimetry. The variety of different scenarios included herein attest to the utility of the evaluated IMAR for a wide range of radiotherapy clinical scenarios.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-D-206-01: Employing a Novel Consensus Optimization Strategy to Achieve Iterative Cone Beam CT Reconstruction On a Multi-GPU Platform

Purpose: While compressed sensing-based cone-beam CT (CBCT) iterative reconstruction techniques have demonstrated tremendous capability of reconstructing high-quality images from undersampled noisy data, its long computation time still hinders wide application in routine clinic. The purpose of this study is to develop a reconstruction framework that employs modern consensus optimization techniques to achieve CBCT reconstruction on a multi-GPU platform for improved computational efficiency. Methods: Total projection data were evenly distributed to multiple GPUs. Each GPU performed reconstruction using its own projection data with a conventional total variation regularization approach to ensure image quality. In addition, the solutions from GPUs were subject to a consistency constraint that they should be identical. We solved the optimization problem with all the constraints considered rigorously using an alternating direction method of multipliers (ADMM) algorithm. The reconstruction framework was implemented using OpenCL on a platform with two Nvidia GTX590 GPU cards, each with two GPUs. We studied the performance of our method and demonstrated its advantages through a simulation case with a NCAT phantom and an experimental case with a Catphan phantom. Result: Compared with the CBCT images reconstructed using conventional FDK method with full projection datasets, our proposed method achieved comparable image quality with about one third projection numbers. The computation time on the multi-GPU platform was ∼55 s and ∼ 35 s in the two cases respectively, achieving a speedup factor of ∼ 3.0 compared with single GPU reconstruction. Conclusion: We have developed a consensus ADMM-based CBCT reconstruction method which enabled performing reconstruction on a multi-GPU platform. The achieved efficiency made this method clinically attractive.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-P-45: Clinical Experience with Radiation Dose Reduction of CT Examinations Using Iterative Reconstruction Algorithms

Purpose: Iterative reconstruction (IR) algorithms have been adopted by medical centers in the past several years. IR has a potential to substantially reduce patient dose while maintaining or improving image quality. This study characterizes dose reductions in clinical settings for CT examinations using IR. Methods: We retrospectively analyzed dose information from patients who underwent abdomen/pelvis CT examinations with and without contrast media in multiple locations of our Healthcare system. A total of 743 patients scanned with ASIR on 64 slice GE lightspeed VCTs at three sites, and 30 patients scanned with SAFIRE on a Siemens 128 slice Definition Flash in one site was retrieved. For comparison, patient data (n=291) from a GE scanner and patient data (n=61) from two Siemens scanners where filtered back-projection (FBP) was used was collected retrospectively. 30% and 10% ASIR, and SAFIRE Level 2 was used. CTDIvol, Dose-length-product (DLP), weight and height from all patients was recorded. Body mass index (BMI) was calculated accordingly. To convert CTDIvol to SSDE, AP and lateral dimensions at the mid-liver level was measured for each patient. Results: Compared with FBP, 30% ASIR reduces dose by 44.1% (SSDE: 12.19mGy vs. 21.83mGy), while 10% ASIR reduced dose by 20.6% (SSDE 17.32mGy vs. 21.83). Use of SAFIRE reduced dose by 61.4% (SSDE: 8.77mGy vs. 22.7mGy). The geometric mean for patients scanned with ASIR was larger than for patients scanned with FBP (geometric mean is 297.48 mmm vs. 284.76 mm). The same trend was observed for the Siemens scanner where SAFIRE was used (geometric mean: 316 mm with SAFIRE vs. 239 mm with FBP). Patient size differences suggest that further dose reduction is possible. Conclusion: Our data confirmed that in clinical practice IR can significantly reduce dose to patients who undergo CT examinations, while meeting diagnostic requirements for image quality.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

TH-CD-207B-09: Task-Driven Fluence Field Modulation Design for Model-Based Iterative Reconstruction in CT

Purpose: There has been increasing interest in integrating fluence field modulation (FFM) devices with diagnostic CT scanners for dose reduction purposes. Conventional FFM strategies, however, are often either based on heuristics or the analysis of filtered-backprojection (FBP) performance. This work investigates a prospective task-driven optimization of FFM for model-based iterative reconstruction (MBIR) in order to improve imaging performance at the same total dose as conventional strategies. Methods: The task-driven optimization framework utilizes an ultra-low dose 3D scout as a patient-specific anatomical model and a mathematical formation of the imaging task. The MBIR method investigated is quadratically penalized-likelihood reconstruction. The FFM objective function uses detectability index, d’, computed as a function of the predicted spatial resolution and noise in the image. To optimize performance throughout the object, a maxi-min objective was adopted where the minimum d’ over multiple locations is maximized. To reduce the dimensionality of the problem, FFM is parameterized as a linear combination of 2D Gaussian basis functions over horizontal detector pixels and projection angles. The coefficients of these bases are found using the covariance matrix adaptation evolution strategy (CMA-ES) algorithm. The task-driven design was compared with three other strategies proposed for FBP reconstruction for a calcification cluster discrimination task in an abdomen phantom. Results: The task-driven optimization yielded FFM that was significantly different from those designed for FBP. Comparing all four strategies, the task-based design achieved the highest minimum d’ with an 8–48% improvement, consistent with the maxi-min objective. In addition, d’ was improved to a greater extent over a larger area within the entire phantom. Conclusion: Results from this investigation suggests the need to re-evaluate conventional FFM strategies for MBIR. The task-based optimization framework provides a promising approach that maximizes imaging performance under the same total dose constraint.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Movement of the Melt Metal Layer under Conditions Typical of Transient Events in ITER

During the operation of ITER, protective coatings of the divertor and the first wall will be exposed to significant plasma heat loads which may cause a huge erosion. One of the major failure mechanisms of metallic armor is diminution of their thickness due to the melt layer displacement. New experimental data are required in order to develop and validate physical models of the melt layer movement. The paper presents the experiments where metal targets were irradiated by a plasma stream at the quasi-stationary high-current plasma accelerator QSPA-T. The obtained data allow one to determine the velocity and acceleration of the melt layer at various distances from the plasma stream axis. The force causing the radial movement of the melt layer is shown to create an acceleration whose order of magnitude is 1000g. The pressure gradient is not responsible for creating this large acceleration. To investigate the melt layer movement under a known force, the experiment with a rotating target was carried out. The influence of centrifugal and Coriolis forces led to appearance of curved elongated waves on the surface. The surface profile changed: there is no hill in the central part of the erosion crater in contrast to the stationary target. The experimental data clarify the trends in the melt motion that are required for development of theoretical models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Specific features of measuring the isotopic composition of hydrogen ions in ITER plasma by using neutral particle diagnostics under neutral beam injection conditions

Results of numerical simulation of signals from neutral particle analyzers under injection of the heating and diagnostic neutral beams in different operating modes of the ITER tokamak are presented. The distribution functions of fast ions in plasma are simulated, and the corresponding neutral particle fluxes escaping from the plasma along the line of sight of the analyzers are calculated. It is shown that the injection of heating deuterium (D{sup 0}) beams results in the appearance of an intense background signal hampering measurements of the ratio between the densities of deuterium and tritium fuel ions in plasma in the thermal energy range. The injection of a diagnostic hydrogen (H{sup 0}) beam does not affect measurements owing to the high mass resolution of the analyzers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

PolyCheck: Dynamic Verification of Iteration Space Transformations on Affine Programs

High-level compiler transformations, especially loop transformations, are widely recognized as critical optimizations to restructure programs to improve data locality and expose parallelism. Guaranteeing the correctness of program transformations is essential, and to date three main approaches have been developed: proof of equivalence of affine programs, matching the execution traces of programs, and checking bit-by-bit equivalence of the outputs of the programs. Each technique suffers from limitations in either the kind of transformations supported, space complexity, or the sensitivity to the testing dataset. In this paper, we take a novel approach addressing all three limitations to provide an automatic bug checker to verify any iteration reordering transformations on affine programs, including non-affine transformations, with space consumption proportional to the original program data, and robust to arbitrary datasets of a given size. We achieve this by exploiting the structure of affine program control- and data-flow to generate at compile-time lightweight checker code to be executed within the transformed program. Experimental results assess the correctness and effectiveness of our method, and its increased coverage over previous approaches.

affine programs↗
DOE OSTI2020

Vertical forces during VDEs in an ITER plasma and the role of halo currents

Vertical displacement events (VDEs) can occur in elongated tokamaks causing large currents to flow in the vessel and other adjacent metallic structures. To better understand the potential magnitude of the associated forces and the role of the so called halo currents'' on them, we have used the M3D-C1 code to simulate potential VDEs in ITER. We used actual values for the vessel resistivity and pre-quench temperatures and, unlike most of the previous studies, the halo region is naturally formed by triggering the thermal quench with an increase in the plasma thermal conductivity. We used the 2D non-linear version of the code and vary the post-thermal quench thermal conductivity profile as well as the boundary temperature in order to generate a wide range of possible cases that could occur in the experiment. We also show that, for a similar condition, increasing the halo current does not increase the total force on the wall since it is offset by a decrease in the toroidal contribution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Source term evaluation for accident transients in the experimental fusion facility ITER

We have studied the transport and chemical speciation of radio-toxic and toxic species for an event of water ingress in the vacuum vessel of experimental fusion facility ITER with the ASTEC code. In particular our evaluation takes into account an assessed thermodynamic data for the beryllium gaseous species. This study shows that deposited beryllium dusts of atomic Be and Be(OH){sub 2} are formed. It also shows that Be(OT){sub 2} could exist in some conditions in the drain tank. (authors)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Block-Iterative Methods for 3D Constant-Coefficient Stencils on GPUs and Multicore CPUs

Block iterative methods are extremely important as smoothers for multigrid methods, as preconditioners for Krylov methods, and as solvers for diagonally dominant linear systems. Developing robust and efficient smoother algorithms suitable for current and evolving GPU and multicore CPU systems is a significant challenge. We address this issue in the case of constant-coefficient stencils arising in the solution of elliptic partial differential equations on structured 3D uniform and adaptively refined block structured grids. Robust, highly parallel implementations of block Jacobi and chaotic block Gauss-Seidel algorithms with exact inversion of the blocks are developed using different parallelization techniques. Experimental results for NVIDIA Fermi/Kepler GPUs and AMD multicore systems are presented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Progress Update on Iterative Reconstruction of Neutron Tomographic Images

This report satisfies the fiscal year 2016 technical deliverable to report on progress in development of fast iterative reconstruction algorithms for project OR16-3DTomography-PD2Jb, "3D Tomography and Image Processing Using Fast Neutrons." This project has two overall goals. The first of these goals is to extend associated-particle fast neutron transmission and, particularly, induced-reaction tomographic imaging algorithms to three dimensions. The second of these goals is to automatically segment the resultant tomographic images into constituent parts, and then extract information about the parts, such as the class of shape and potentially shape parameters. This report addresses of the component of the project concerned with three-dimensional (3D) image reconstruction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗