SEARCH · Engineering Papers

Results for “Iterative methods”

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 595 records · Page 33

All sources NASA NTRS DOE OSTI

DOE OSTI2012

Development and preliminary verification of the 3D core neutronic code: COCO

As the recent blooming economic growth and following environmental concerns (China)) is proactively pushing forward nuclear power development and encouraging the tapping of clean energy. Under this situation, CGNPC, as one of the largest energy enterprises in China, is planning to develop its own nuclear related technology in order to support more and more nuclear plants either under construction or being operation. This paper introduces the recent progress in software development for CGNPC. The focus is placed on the physical models and preliminary verification results during the recent development of the 3D Core Neutronic Code: COCO. In the COCO code, the non-linear Green's function method is employed to calculate the neutron flux. In order to use the discontinuity factor, the Neumann (second kind) boundary condition is utilized in the Green's function nodal method. Additionally, the COCO code also includes the necessary physical models, e.g. single-channel thermal-hydraulic module, burnup module, pin power reconstruction module and cross-section interpolation module. The preliminary verification result shows that the COCO code is sufficient for reactor core design and analysis for pressurized water reactor (PWR). (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2016

WE-FG-207B-06: Plaque Composition Measurement with Dual Energy Computed Tomography

Purpose: To investigate the feasibility of characterizing arterial plaque composition in terms of water, lipid and protein or calcium using dual energy computed tomography. Characterization of plaque composition can potentially help distinguish vulnerable from stable plaques. Methods: Simulations studies were performed by the CT simulator based on ASTRA tomography toolbox. The beam energy for dual energy images was selected to be 80 kVp and 135 kVp. The radiation dose and energy spectrum for the CT simulator were carefully calibrated with respect to a 320-slice CT scanner. A digital chest phantom was constructed using Matlab for calibration and plaque measurement. Pure water, lipid, protein or calcium was used for calibration and a mixture of different volume percentages of these materials were used for validation purposes. Non-calcified plaque was simulated using water, lipid and protein with volumetric percentage range of 35%∼65%, 5%∼60% and 5%∼40%, respectively. Calcified plaque was simulated using water, lipid and calcium with volumetric percentage range of 50%∼80%, 8%∼45% and 3%∼13%, respectively. We employed iterative sinogram processing (ISP) to reduce the beam hardening effect in the simulation to improve the decomposition results. Results: The simulated known composition and dual energy decomposition results were in good agreement. Water, lipid and protein (calcium) mixtures were decomposed into water, lipid and protein (calcium) contents. The RMS errors of volumetric percentage for the water, lipid and protein (non-calcified plaque) decomposition, as compared to known values, were estimated to be approximately 5.74%, 2.54%, and 0.95% respectively. The RMS errors of volumetric percentage for the water, lipid and Calcium (calcified plaque) decomposition, as compared to known values, were estimated to be approximately 7.4%, 8.64%, and 0.08% respectively. Conclusion: The results of this study suggest that the dual energy decomposition can potentially be used to quantify the water, lipid, and protein or calcium composition of a plaque with relatively good accuracy. Grant funding from Toshiba Medical Systems and Philips Medical Systems.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2015

SU-E-T-792: Validation of a Secondary TPS for IROC-H Recalculation of Anthropomorphic Phantoms

Purpose: To validate a secondary treatment planning system (sTPS) for use by the Imaging & Radiation Oncology Core-Houston (IROC-H). The TPS will recalculate phantom irradiations submitted by institutions to IROC-H and compare plan results of the institution to the sTPS. Methods: In-field dosimetric data was collected by IROC-H for numerous linacs at 6, 10, 15, and 18 MV. The data was aggregated and used to define reference linac classes; each class was then modeled in the sTPS (Mobius3D) by matching the in-field characteristics. Fields used to collect IROC-H data were recreated and recalculated using Mobius3D. The same dosimetric points were measured in the recalculation and compared to the initial collection data. Additionally, a 6MV Monte Carlo beam configuration was used to compare penumbrae in the Mobius3D models. Finally, a handful of IROC-H head and neck phantoms were recalculated using Mobius3D. Results: Recalculation and quantification of differences between reference data and Mobius3D values resulted in a relative matching score of 12.45 (0 is a perfect match) for the default 6MV Mobius3D beam configuration. By adjusting beam configuration options, iterations resulted in scores of 8.45, 6.32, and 3.52, showing that customization could have a dramatic effect on beam configuration. After in-field optimization, penumbra was compared between Monte Carlo and Mobius3D for the reference fields. For open jaw fields, FWHM field widths and penumbra widths were different by <0.6 and <1mm respectively; for MLC open fields the penumbra widths were up to 1.5mm different. Phantom recalculations showed good agreement, having an average of 0.6% error per beam. Conclusion: A secondary TPS has been validated for simple irradiation geometries using reference data collected by IROC-H. The beam was customized to the reference data iteratively and resulted in a good match. This system can provide independent recalculation of phantom plans based on independent reference data.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

SU-C-209-03: Anti-Scatter Grid-Line Artifact Minimization for Removing the Grid Lines for Three Different Grids Used with a High Resolution CMOS Detector

Purpose: Demonstrate the effectiveness of an anti-scatter grid artifact minimization method by removing the grid-line artifacts for three different grids when used with a high resolution CMOS detector. Method: Three different stationary x-ray grids were used with a high resolution CMOS x-ray detector (Dexela 1207, 75 µm pixels, sensitivity area 11.5cm × 6.5cm) to image a simulated artery block phantom (Nuclear Associates, Stenosis/Aneurysm Artery Block 76–705) combined with a frontal head phantom used as the scattering source. The x-ray parameters were 98kVp, 200mA, and 16ms for all grids. With all the three grids, two images were acquired: the first for a scatter-less flat field including the grid and the second of the object with the grid which may still have some scatter transmission. Because scatter has a low spatial frequency distribution, it was represented by an estimated constant value as an initial approximation and subtracted from the image of the object with grid before dividing by an average frame of the grid flat-field with no scatter. The constant value was iteratively changed to minimize residual grid-line artifact. This artifact minimization process was used for all the three grids. Results: Anti-scatter grid lines artifacts were successfully eliminated in all the three final images taken with the three different grids. The image contrast and CNR were also compared before and after the correction, and also compared with those from the image of the object when no grid was used. The corrected images showed an increase in CNR of approximately 28%, 33% and 25% for the three grids, as compared to the images when no grid at all was used. Conclusion: Anti-scatter grid-artifact minimization works effectively irrespective of the specifications of the grid when it is used with a high spatial resolution detector. Partial support from NIH Grant R01-EB002873 and Toshiba Medical Systems Corp.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

MO-AB-BRA-09: Development and Evaluation of a Biomechanical Modeling-Assisted CBCT Reconstruction Technique (Bio-Recon)

Purpose: To develop a Bio-recon technique by incorporating the biomechanical properties of anatomical structures into the deformation-based CBCT reconstruction process. Methods: Bio-recon reconstructs the CBCT by deforming a prior high-quality CT/CBCT using a deformation-vector-field (DVF). The DVF is solved through two alternating steps: 2D–3D deformation and finite-element-analysis based biomechanical modeling. 2D–3D deformation optimizes the DVF through an ‘intensity-driven’ approach, which updates the DVF to minimize intensity mismatches between the acquired projections and the simulated projections from the deformed CBCT. In contrast, biomechanical modeling optimizes the DVF through a ‘biomechanical-feature-driven’ approach, which updates the DVF based on the biophysical properties of anatomical structures. In general, Biorecon extracts the 2D–3D deformation-optimized DVF at high-contrast structure boundaries, and uses it as the boundary condition to drive biomechanical modeling to optimize the overall DVF, especially at low-contrast regions. The optimized DVF is fed back into the 2D–3D deformation for further optimization, which forms an iterative loop. The efficacy of Bio-recon was evaluated on 11 lung patient cases, each with a prior CT and a new CT. Cone-beam projections were generated from the new CTs to reconstruct CBCTs, which were compared with the original new CTs for evaluation. 872 anatomical landmarks were also manually identified by a clinician on both the prior and new CTs to track the lung motion, which was used to evaluate the DVF accuracy. Results: Using 10 projections for reconstruction, the average (± s.d.) relative errors of reconstructed CBCTs by the clinical FDK technique, the 2D–3D deformation-only technique and Bio-recon were 46.5±5.9%, 12.0±2.3% and 10.4±1.3%, respectively. The average residual errors of DVF-tracked landmark motion by the 2D–3D deformation-only technique and Bio-recon were 5.6±4.3mm and 3.1±2.4mm, respectively. Conclusion: Bio-recon improved accuracy for both the reconstructed CBCT and the DVF. The accurate DVF can benefit multiple clinical practices, such as image-guided adaptive radiotherapy. We acknowledge funding support from the American Cancer Society (RSG-13-326-01-CCE), from the US National Institutes of Health (R01 EB020366), and from the Cancer Prevention and Research Institute of Texas (RP130109).

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

SU-F-T-61: Treatment Planning Observations for the CivaSheet Directional Brachytherapy Device Using VariSeed 9.0

Purpose: The VariSeed 9.0 brachytherapy TPS is recently available and has new features such as ability to rotate a brachytherapy source away from normal to the imaging plane. Consequently, a dosimetric analysis was performed for a directional brachytherapy source (CivaSheet) with tests of this functionality and experiences from clinical treatment planning were documented. These observations contribute to safe, practical, and accurate use of such new software features. Methods: Several tests were established to evaluate the new rotational feature, specific to the CivaSheet for the first patients treated using this new brachytherapy device. These included suitability of imaging slice-thickness and in-plane resolution, window/level adjustments for brachytherapy source visualization, commissioning the source physical length for performing rotations, and using different planar and 3D window views to identify source orientation. Additional CivaSheet-specific tests were performed to determine the dosimetric influence on target coverage: changing the source tilt angle, source positioning in the treatment plan based on the CivaSheet rectangular array of CivaDots, and influence of prescription depth on the necessary treatment margin for adequate target coverage. Results: Higher imaging-resolution produced better accuracy for source orientation and positioning, with sub-millimeter CT slice-thickness and in-plane resolution preferred. Source rotation was possible only in sagittal or coronal views. The process for validating source orientation required iteratively altering rotations then checking them in the 3D view, which was cumbersome given the absence of quantitative plan documentation to indicate orientation. Given the small Pd-103 source size, influence of source tilt within 30° was negligible for <1.0 cm. Influence of source position was important when the source was positioned in/out of the adjacent source plane, causing changes of 15%, 7%, and 3% at depths of 0.5, 0.7, and 1.0 cm. Conclusion: The new TPS rotational feature worked well, but several issues were identified to improve the treatment planning process. Research supported in part by CivaTech Oncology, Inc. for Dr. Rivard.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

A Prototype Gamma Tomosynthesis System for Molecular Breast Imaging

We are developing an improved gamma camera for molecular breast imaging (MBI) that allows tomographic imaging and improved signal-to-noise over conventional planar methods. We have built a novel gamma camera for MBI that utilizes limited angle tomography, similar to x-ray tomosynthesis, to generate images with depth information and improved signal-to-noise. A unique feature of the gamma tomosynthesis system is a variable-angle slant hole (VASH) collimator. This collimator allows the camera to remain flush against the compression paddle during the tomographic acquisition, which achieves high spatial resolution and simplified detector motion. We constructed a VASH collimator from a stack of tungsten sheets, each sheet containing a matrix of holes created by photo-etching. With the holes of the sheets aligned, a parallel-hole collimator is created; shearing the sheets creates a variable angle slant hole collimator. The shearing is controlled by a motorized mechanism. The collimator was mounted on a commercial MBI camera and projection data were acquired of a set of capillary tube line sources over a 50° range of angles. The projection data were reconstructed using the iterative MLEM method. Results demonstrated the ability of the system to resolve in the depth dimension, and spatial resolution matched expected values. We conclude that the proposed gamma tomosynthesis system has the potential to provide improved signal-to-noise over conventional MBI methods and allow a reduction in the administered radioactivity and patient dose.

Gilland, David↗

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 OSTI2021

Plateau Phenomenon in Gradient Descent Training of RELU Networks: Explanation, Quantification, and Avoidance

The ability of neural networks to provide ‘best in class’ approximation across a wide range of applications is well-documented. Nevertheless, the powerful expressivity of neural networks comes to naught if one is unable to effectively train (choose) the parameters defining the network. In general, neural networks are trained by gradient descent type optimization methods,a stochastic variant thereof. In practice, such methods result in the loss function decreases rapidly at the beginning of training but then, after a relatively small number of steps, significantly slow down. The loss may even appear to stagnate over the period of a large number of epochs, only to then suddenly start to decrease fast again for no apparent reason. This so-called plateau phenomenon manifests itself in many learning tasks. The present work aims to identify and quantify the root causes of plateau phenomenon.analysis is carried out in the setting of univariate ReLU networks. No assumptions are made on the number of neurons relative to the number of training data, and our results hold for both the lazy and adaptive regimes. Here, the main findings are: plateaux correspond to periods during which activation patterns remain constant, where activation pattern refers to the number of data points that activate a given neuron; quantification of convergence of the gradient flow dynamics; and, characterization stationary points in terms solutions of local least squares regression lines over subsets of the training data. Based on these conclusions, we propose a new iterative training method, the Active Neuron Least Squares (ANLS), characterised by the explicit adjustment of the activation pattern at each step, which is designed to enable a quick exit from a plateau. Illustrative numerical examples are included throughout.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

LoTToR: An Algorithm for Missing-Wedge Correction of the Low-Tilt Tomographic 3D Reconstruction of a Single-Molecule Structure

A single-molecule three-dimensional (3D) structure is essential for understanding the thermal vibrations and dynamics as well as the conformational changes during the chemical reaction of macromolecules. Individual-particle electron tomography (IPET) is an approach for obtaining a snap-shot 3D structure of an individual macromolecule particle by aligning the tilt series of electron tomographic (ET) images of a targeted particle through a focused iterative 3D reconstruction method. The method can reduce the influence on the 3D reconstruction from large-scale image distortion and deformation. Due to the mechanical tilt limitation, 3D reconstruction often contains missing-wedge artifacts, presented as elongation and an anisotropic resolution. Here, we report a post-processing method to correct the missing-wedge artifact. This low-tilt tomographic reconstruction (LoTToR) method contains a model-free iteration process under a set of constraints in real and reciprocal spaces. A proof of concept is conducted by using the LoTToR on a phantom, i.e., a simulated 3D reconstruction from a low-tilt series of images, including that within a tilt range of ±15°. The method is validated by using both negative-staining (NS) and cryo-electron tomography (cryo-ET) experimental data. A significantly reduced missing-wedge artifact verifies the capability of LoTToR, suggesting a new tool to support the future study of macromolecular dynamics, fluctuation and chemical activity from the viewpoint of single-molecule 3D structure determination.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Porting hypre to heterogeneous computer architectures: Strategies and experiences

We report that linear systems are occurring in many applications, and solving them can take a large amount of the total simulation time. The high performance library hypre provides a variety of interfaces and linear solvers, including various multigrid methods, that have achieved good scalability on a variety of homogeneous parallel computer architectures. Heterogeneous architectures with nodes that have both CPUs and accelerators provide new challenges, since they require more fine-grained parallelism and reduced data movement between different memories on a single node as well as across nodes. We will discuss our experiences and strategies to port hypre to heterogeneous computers with accelerators, including the design of a new memory model, the use of abstractions, the BoxLoop macros in the structured and semi-structured interfaces, and the restructuring of algebraic multigrid (AMG) into modular components. We present numerical experiments comparing CPU and GPU performance for several test problems.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2016

Thermal analysis of the in-vessel components of the ITER plasma-position reflectometry

The ITER plasma position reflectometry system measures the edge electron density profile of the plasma, providing real-time supplementary contribution to the magnetic measurements of the plasma-wall distance. Some of the system components will be in direct sight of the plasma and therefore subject to plasma and stray radiation, which may cause excessive temperatures and stresses. In this work, thermal finite element analysis of the antenna and adjacent waveguides is conducted with ANSYS V17 (ANSYS® Academic Research, Release 17.0, 2016). Results allow the identification of critical temperature points, and solutions are proposed to improve the thermal behavior of the system.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2016

MO-DE-BRA-06: 3D Image Acquisition and Reconstruction Explained with Online Animations

Purpose: Understanding the principles of 3D imaging and image reconstruction is fundamental to the field of medical imaging. Clinicians, technologists, physicists, patients, students, and inquisitive minds all stand to benefit from greater comprehension of the supporting technologies. To help explain the basic principles of 3D imaging, we developed multi-frame animations that convey the concepts of tomographic imaging. The series of free (gif) animations are accessible online, and provide a multimedia introduction to the main concepts of image reconstruction. Methods: Text and animations were created to convey the principles of analytic tomography in CT, PET, and SPECT. Specific topics covered included: principles of sinograms/image data storage, forward projection, principles of PET acquisitions, and filtered backprojection. A total of 8 animations were created and presented for CT, PET, and digital phantom formats. In addition, a free executable is also provided to allow users to create their own tomographic animations – providing an opportunity for interaction and personalization to help foster user interest. Results: Tutorial text and animations have been posted online, freely available to view or download. The animations are in first position in a google search of “image reconstruction animations”. The website currently receives approximately 200 hits/month, from all over the world, and the usage is growing. Positive feedback has been collected from users. Conclusion: We identified a need for improved teaching tools to help visualize the (temporally variant) concepts of image reconstruction, and have shown that animations can be a useful tool for this aspect of education. Furthermore, posting animations freely on the web has shown to be a good way to maximize their impact in the community. In future endeavors, we hope to expand this animated content, to cover principles of iterative reconstruction, as well as other phenomena relating to imaging.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

SU-G-TeP2-11: Initial Evaluation of a Novel Split-Filter Dual-Energy CT for Use in Radiation Oncology

Purpose: To perform an initial evaluation of a novel split-filter dual-energy CT (DECT) system with the goal of understanding the clinical utility and limitations of the system for radiation therapy. Methods: Several phantoms were imaged using the split-filter DECT technique on the Siemens Edge CT scanner using a range of clinically-relevant doses. The optimum-contrast reconstruction, the mixed reconstruction, and the monoenergetic reconstructions (ranging from 40 keV to 190 keV) were evaluated. Each image was analyzed for CT number accuracy, uniformity, noise, low-contrast visibility (LCV), spatial resolution and geometric distortion. For comparison purposes, all parameters were evaluated on 120 kVp single-energy CT (SECT) scans used for treatment planning, as well as, a sequential-scan DECT technique for corresponding doses. Results: For all DECT reconstructions no observable geometric distortion was found. Both the optimal-contrast and mixed images demonstrated slight improvements in LCV and noise when compared to the SECT, and slight reductions in CT number accuracy and spatial resolution. The CT numbers trended as expected for the monoenergetic reconstructions, with CT number accuracy within 50 HU for materials of density <2 g/cm3. Spatial resolution increased with energy, and for monoenergetic reconstructions >70 keV the spatial resolution exceeded that of the SECT. The noise in the monoenergetic reconstructions increased with decreasing energy. Thus, the image uniformity, signal-to-noise ratio and LCV were diminished at lower energies (70 keV). Applying iterative reconstruction techniques to the low-energy images reduced noise and improved LCV. The signal-to-noise ratio was stable for energies >100 keV. Conclusion: The initial commissioning of the novel split-filter DECT technology demonstrated favorable results for clinical implementation. The mixed reconstruction showed potential as a replacement for the treatment planning SECT. The image parameters for the monoenergetic reconstructions varied appropriately with energy. This work provides an initial understanding of the limitations and potential applications for monoenergetic imaging.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2014

Application of the TITAN Hybrid Deterministic Transport Code to Medical Physics - Paper 118

The TITAN code is a 3-dimensional parallel deterministic transport code that solves the linear Boltzmann equation. The TITAN code has hybrid algorithms that are ideal for medical physics applications: i) SN and characteristics method and ii) SN with a fictitious quadrature set and ray-tracing. These hybrid algorithms have been used in the simulation of computed tomography (CT) and single photon emission computed tomography (SPECT) In this paper, the hybrid algorithms developed for medical physics applications in the TITAN code are discussed. First, the TITAN code's block-oriented approach for utilizing SN and characteristic solvers for CT is presented and results for a simple test case are given. Second, the paper describes a hybrid formulation for SPECT with an SN solver in the phantom, a fictitious quadrature set for directions from the phantom surface passing through the collimator, and a simplified ray-tracing algorithm for particle transport to the detector. Finally, ongoing work is described and concluding remarks are given. The TITAN code's hybrid algorithms allow it to efficiently solve problems with low-scatter regions, which are frequently found in medical physics applications. The TITAN code has been used to solve a simplified CT problem and generate projection images for SPECT simulations. The WCOS technique was developed to improve the TITAN code's collimator representation. For these applications, the TITAN code was shown to require significantly less computation time than MCNP5. Work is currently underway to implement an iterative image reconstruction algorithm in the TITAN code. A script that calls the TITAN code to complete the forward projection step in SPECT has already been used to reconstruction a 2-dimensional phantom. The ability of TITAN to fully model particle transport in the patient in short computation times is expected to be a great advantage to iterative reconstruction. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2014

WE-D-BRF-05: Quantitative Dual-Energy CT Imaging for Proton Stopping Power Computation

Purpose: To extend the two-parameter separable basis-vector model (BVM) to estimation of proton stopping power from dual-energy CT (DECT) imaging. Methods: BVM assumes that the photon cross sections of any unknown material can be represented as a linear combination of the corresponding quantities for two bracketing basis materials. We show that both the electron density (ρe) and mean excitation energy (Iex) can be modeled by BVM, enabling stopping power to be estimated from the Bethe-Bloch equation. We have implemented an idealized post-processing dual energy imaging (pDECT) simulation consisting of monogenetic 45 keV and 80 keV scanning beams with polystyrene-water and water-CaCl2 solution basis pairs for soft tissues and bony tissues, respectively. The coefficients of 24 standard ICRU tissue compositions were estimated by pDECT. The corresponding ρe, Iex, and stopping power tables were evaluated via BVM and compared to tabulated ICRU 44 reference values. Results: BVM-based pDECT was found to estimate ρe and Iex with average and maximum errors of 0.5% and 2%, respectively, for the 24 tissues. Proton stopping power values at 175 MeV, show average/maximum errors of 0.8%/1.4%. For adipose, muscle and bone, these errors result range prediction accuracies less than 1%. Conclusion: A new two-parameter separable DECT model (BVM) for estimating proton stopping power was developed. Compared to competing parametric fit DECT models, BVM has the comparable prediction accuracy without necessitating iterative solution of nonlinear equations or a sample-dependent empirical relationship between effective atomic number and Iex. Based on the proton BVM, an efficient iterative statistical DECT reconstruction model is under development.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2014

SU-E-T-406: Use of TrueBeam Developer Mode and API to Increase the Efficiency and Accuracy of Commissioning Measurements for the Varian EDGE Stereotactic Linac

Purpose: The Varian EDGE machine is a new stereotactic platform, combining Calypso and VisionRT localization systems with a stereotactic linac. The system includes TrueBeam DeveloperMode, making possible the use of XML-scripting for automation of linac-related tasks. This study details the use of DeveloperMode to automate commissioning tasks for Varian EDGE, thereby improving efficiency and measurement consistency. Methods: XML-scripting was used for various commissioning tasks,including couch model verification,beam-scanning,and isocenter verification. For couch measurements, point measurements were acquired for several field sizes (2×2,4×4,10×10cm{sup 2}) at 42 gantry angles for two couch-models. Measurements were acquired with variations in couch position(rails in/out,couch shifted in each of motion axes) compared to treatment planning system(TPS)-calculated values,which were logged automatically through advanced planning interface(API) scripting functionality. For beam scanning, XML-scripts were used to create custom MLC-apertures. For isocenter verification, XML-scripts were used to automate various Winston-Lutz-type tests. Results: For couch measurements, the time required for each set of angles was approximately 9 minutes. Without scripting, each set required approximately 12 minutes. Automated measurements required only one physicist, while manual measurements required at least two physicists to handle linac positions/beams and data recording. MLC apertures were generated outside of the TPS,and with the .xml file format, double-checking without use of TPS/operator console was possible. Similar time efficiency gains were found for isocenter verification measurements Conclusion: The use of XML scripting in TrueBeam DeveloperMode allows for efficient and accurate data acquisition during commissioning. The efficiency improvement is most pronounced for iterative measurements, exemplified by the time savings for couch modeling measurements(approximately 10 hours). The scripting also allowed for creation of the files in advance without requiring access to TPS. The API scripting functionality enabled efficient creation/mining of TPS data. Finally, automation reduces the potential for human error in entering linac values at the machine console,and the script provides a log of measurements acquired for each session. This research was supported in part by a grant from Varian Medical Systems, Palo Alto, CA.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2015

A parallel multi-domain solution methodology applied to nonlinear thermal transport problems in nuclear fuel pins

This paper describes an efficient and nonlinearly consistent parallel solution methodology for solving coupled nonlinear thermal transport problems that occur in nuclear reactor applications over hundreds of individual 3D physical subdomains. Efficiency is obtained by leveraging knowledge of the physical domains, the physics on individual domains, and the couplings between them for preconditioning within a Jacobian Free Newton Krylov method. Details of the computational infrastructure that enabled this work, namely the open source Advanced Multi-Physics (AMP) package developed by the authors is described. Details of verification and validation experiments, and parallel performance analysis in weak and strong scaling studies demonstrating the achieved efficiency of the algorithm are presented. Furthermore, numerical experiments demonstrate that the preconditioner developed is independent of the number of fuel subdomains in a fuel rod, which is particularly important when simulating different types of fuel rods. Finally, we demonstrate the power of the coupling methodology by considering problems with couplings between surface and volume physics and coupling of nonlinear thermal transport in fuel rods to an external radiation transport code.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗