Search NASA⌕ Search

SEARCH · Search NASA

Results for “Nuclear Data Evaluation 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.

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

At least 793 records · Page 44

Numerical Study of Vortex Core in Turbulent Swirling Flow

Swirling flow has attracted a great deal of research interest. An application of swirling flow can be found in many engineering areas; such as gas cyclone and heat transfer enhancement. In nuclear engineering, swirling flow is used as heat transfer enhancement in rod bundle by installing grid spacer and mixing vane. Because of its extensive applications, there has been a great number of experimental and numerical investigations of swirling flow. McClusky et al., experimentally studied mapping of the lateral velocity in rod bundle using Laser Doppler Velocimetry (LDV) and Conner et al. predicted the flow field in rod bundle by CFD. To understand the characteristics of swirling flow such as tangential velocity and vortex core swirling flow, some researchers investigated swirling flow in a pipe instead of in rod bundle. This is due to the similarity between swirling flow in a pipe and rod bundle. In recent study, flow mapping of 2D velocity field of single-phase swirling flow had been conducted experimentally by Particle Image Velocimetry (PIV). However, this method needs calibration and complicated. The new measurement technique for flow mapping in swirling flow was done by phased array Ultrasonic Velocity Profiler (PAUVP). Hamdani et al., successfully plotted velocity field and determined the swirling core region (vortex core). This technique is nonintrusive and no need calibration. However, there was a limitation of the measurement because phased array UVP cannot capture 360 deg. velocity field inside the pipe. In addition, tangential velocity in CFD results was not agree well with experiment data due asymmetric problem. To evaluate this limitation of experimental method, the numerical simulations using Computational Fluid Dynamics (CFD) needs to be performed. The objective of this study is to use the commercial CFD code Fluent to model the developing turbulent swirling flow induced by twisted tape inside a straight pipe, and to gain an understanding of such a flow and the turbulence models involved in the simulation. Straight pipe is used instead of rod bundle for a simplification. In this work, the Renormalized Group (RNG) k-{epsilon} turbulence model and Reynolds Stress Model (RSM) turbulence model, are performed to study the phenomena of vortex core in a pipe with twisted tape inserts. Effects of the Reynolds number and asymmetric of twisted tape (inclination, {alpha} = 5 deg. and {alpha} = 20 deg.) on swirling core position are examined using water as the testing fluid. Furthermore, the Navier-Stokes equations will be solved to obtain results that can be validated with the previous experimental results. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Nuclear structure studies in the seaborgium region at SHIP

New decay data for the isotopes {sup 259}Sg and {sup 255}Rf were obtained at the velocity filter SHIP using an α-decay spectroscopy measurement. Both isotopes were produced and studied via a one neutron evaporation channel in the compound fusion reaction {sup 54}Cr+{sup 208}Pb. New isomeric states were observed and the single-particle level systematics for isotones with 151 and 153 neutrons were extended. A change of the ground-state configuration for the heaviest N = 151 isotones was observed. Detailed Monte-Carlo simulation for the α decay of {sup 259}Sg applying the GEANT4 toolkit was performed and compared with experimental data.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Impact of new prompt fission neutron spectrum evaluations and the associated cross-correlations on critical assemblies using MCNP6 - 14385

Two recent re-evaluations of the neutron-induced prompt fission neutron spectrum (PFNS) of uranium and plutonium actinides have included some cross-correlated uncertainties that have not been released in a new nuclear data library. The first re-evaluation included cross-isotope correlations of the PFNS within each of the uranium and plutonium suites of actinides below the second-chance fission threshold. The second, more recent re-evaluation includes incident energy-dependent correlations of the {sup 239}Pu(n,f) PFNS. Both of these new evaluations and their cross-correlations, given the importance of these actinide quantities and because the impact of the incident energy-dependent correlations has never been studied previously, could have an impact on a variety of criticality safety applications. For example, the ke uncertainty due to the uranium PFNS uncertainties in the Big Ten critical assembly is underpredicted by ∼10% when the cross-isotope correlations in the uranium PFNS are ignored. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Lunar and Planetary Science XXXV: Meteorites

The session "Meteorites" included the following reports:Description of a New Stony Meteorite Find from Bulloch County, Georgia; Meteorite Ablation Derived from Cosmic Ray Track Data Dhofar 732: A Mg-rich Orthopyroxenitic Achondrite Halogens, Carbon and Sulfur in the Tagish Lake Meteorite: Implications for Classification and Terrestrial Alteration; Electromagnetic Scrape of Meteorites and Probably Columbia Tiles; Pre-Atmospheric Sizes and Orbits of Several Chondrites; Research of Shock-Thermal History of the Enstatite Chondrites by Track, Thermoluminescence and Neutron-Activation (NAA) Methods; Radiation and Shock-thermal History of the Kaidun CR2 Chondrite Glass Inclusions; On the Problem of Search for Super-Heavy Element Traces in the Meteorites: Probability of Their Discovery by Three-Prong Tracks due to Nuclear Spontaneous Fission Trace Element Abundances in Separated Phases of Pesyanoe, Enstatite Achondrite; Evaluation of Cooling Rate Calculated by Diffusional Modification of Chemical Zoning: Different Initial Profiles for Diffusion Calculation; Mineralogical Features and REE Distribution in Ortho- and Clinopyroxenes of the HaH 317 Enstatite Chondrite Dhofar 311, 730 and 731: New Lunar Meteorites from Oman; The Deuterium Content of Individual Murchison Amino Acids; Clues to the Formation of PV1, an Enigmatic Carbon-rich Chondritic Clast from the Plainview H-Chondrite Regolith Breccia ;Numerical Simulations of the Production of Extinct Radionuclides and ProtoCAIs by Magnetic Flaring.

Source record↗

Experimental study of single-electron detachment for Cu{sup −}, Ag{sup −}, Au{sup −}, Co{sup −}, and Ir{sup −} in collisions with Kr and Xe

The single electron detachment (SED) cross sections for Cu{sup −}, Ag{sup −}, Au{sup −}, Co{sup −}, and Ir{sup −} in collisions with Kr and Xe in the energy range of 10keV≤E≤30keV are measured by using the growth rate method in our experiments. The SED cross sections increase as the incident velocity increases. These data are compared with other available experimental data for the corresponding negative ions colliding with Ar. The affinities of the negative ions and the mass relation between the projectile and the target seem important in the electron detachment process.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A Bayesian Learning Approach to Wireless Outdoor Heatmap Construction using Deep Gaussian Process

We present a novel Bayesian learning approach to outdoor radio heatmap construction utilizing deep Gaussian process (GP). The proposed approach employs a two-layer hierarchy which consists of two cascaded Gaussian processes that are capable of modeling more complex input-output relations than standard single-layer Gaussian processes. Since deriving the exact model likelihood is challenging, a lower bound is optimized instead so that gradient descent-based methods can be performed to find out the optimal model parameters. Typically, inducing points are used in GPs to facilitate low-rank approximation of covariance (kernel) matrices for computation speedup. However, the inaccuracy induced by inducing points can accumulate when stacking multiple layers of GP which may hinder the performance of deep GP. Moreover, since inducing points need to be learned, having them at all layers of deep GP also incurs computational burden. To overcome the above challenges, in contrast to the canonical deep GP model, we use a modified architecture where a full standard GP resides in the first layer and inducing points are only introduced for the second layer. This modified architecture strikes a balance between model accuracy and training complexity. In the proposed model, the noise parameter of the first GP layer is also eliminated to improve the training efficiency as the noise parameter at the output of the second layer suffices to model the uncertainty in the output. The proposed approach is evaluated on real-world datasets, in the form of location-Received Signal Strength (RSS) pairs, collected from the Platform for Open Wireless Data-driven Experimental Research (POWDER) located at the campus of the University of Utah. Experiment results show that the proposed approach can achieve smaller prediction errors on various training and testing data configurations than DNN-based and GP-based methods.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Scale effects on core design, fuel costs, and spent fuel volume of pressurized water reactors

The desire to improve the economic competitiveness and deployment pace of nuclear energy through modularization, manufacturing, and series production had led to the development of smaller size reactors. As the standard 17x17 fuel technology is mainly maintained in the pressurized water reactors (PWRs) category, this translates into a lower number of fuel assemblies in the core and sometimes a reduced fuel height. To assess the impact of such scale change in core design on fuel cycle cost and spent fuel volume, a scoping analysis tool is developed based on infinite lattice calculations, leakage, fuel management reduced models, and levelized unit cost of electricity (LCOE) estimate. As such, cost dynamics driven by fuel specific power, burnup, core leakage, feed, cycle length, fuel assembly height as well as uranium market data are captured with consistent set of assumptions and analysis methods. A selection of 5 reactor designs representative of leading PWR developers is assessed and compared. Pursuing higher specific powers and optimal burnups are highlighted as the main fuel cost reduction drivers, nevertheless, practical limitations and opportunities must be evaluated to establish the feasibility of such enhanced fuel operation. In consequence, a detailed core design is performed using SIMULATE3 code for 5 PWR variations including natural and forced coolant circulation modes, two reactor scales, power densities of 73, 112, and 123 kW/l and higher discharge burnups. Design and optimization are performed at the lattice level, for the reflector, and at the core loading level. Satisfactory steady-state operation including power distribution, coolant operating limits, and reactivity requirements are analyzed and reported in this paper. The fuel economics of the detailed core designs confirm the scoping analysis findings. Despite the unlocked power uprates in small PWRs, the achievable burnup for a given fuel specific power requires more enrichment and shorter fuel height results in higher fabrication costs per mass of fuel, which makes scaling down core size a more expensive endeavor on the fuel cycle front. Spent fuel volumes are reported for the PWRs designed in this paper. Furthermore, these volumes are driven by the core average discharge burnup regardless of the scale in consideration. Additional cost and core performance aspects related to heavy reflector gains, fuel-reflector substitution, and disposal cost policy in the U.S. are examined.

42 ENGINEERING↗

In-situ Calibration of Detectors using Muon-induced Neutrons

In this work we investigate a method that confirms the operability of neutron detectors requiring neither radiological sources nor radiation-generating devices. This is desirable when radiological sources are not available, but confidence in the functionality of the instrument is required. The “source”, based on the production of neutrons in high-Z materials by muons, provides a tagged, low-background and consistent rate of neutrons that can be used to check the functionality of or calibrate a detector. Using a Monte Carlo guided optimization, an experimental apparatus was designed and built to evaluate the feasibility of this technique. Through a series of trial measurements in a variety of locations we show that gated muon-induced neutrons appear to provide a consistent source of neutrons (35.9 ± 2.3 measured neutrons/10,000 muons in the instrument) under normal environmental variability (less than one statistical standard deviation for 10,000 muons) with a combined environmental + statistical uncertainty of ~18% for 10,000 muons. This is achieved in a single 21-22 minute measurement at sea level.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

The upgrade of the Thomson scattering system for measurement on the C-2/C-2U devices

The C-2/C-2U Thomson scattering system has been substantially upgraded during the latter phase of C-2/C-2U program. A Rayleigh channel has been added to each of the three polychromators of the C-2/C-2U Thomson scattering system. Onsite spectral calibration has been applied to avoid the issue of different channel responses at different spots on the photomultiplier tube surface. With the added Rayleigh channel, the absolute intensity response of the system is calibrated with Rayleigh scattering in argon gas from 0.1 to 4 Torr, where the Rayleigh scattering signal is comparable to the Thomson scattering signal at electron densities from 1 × 10{sup 13} to 4 × 10{sup 14} cm{sup −3}. A new signal processing algorithm, using a maximum likelihood method and including detailed analysis of different noise contributions within the system, has been developed to obtain electron temperature and density profiles. The system setup, spectral and intensity calibration procedure and its outcome, data analysis, and the results of electron temperature/density profile measurements will be presented.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Using choice architecture to exploit a university Distinct Urban Mine

Highlights: • More than USD 11 million lie within this UK DUM for small household appliances. • Different collection systems should be established for working or broken appliances. • Consumer behaviour should inform choice architecture. •DUM evaluations should inform future collection system investments. - Abstract: There are widespread concerns regarding the potential future scarcity of ferrous and non-ferrous materials. However, there are already potentially rich reserves of secondary materials via high ownership of Electrical and Electronic Equipment (EEE) in economically-developed nations. Young people are particularly high consumers of EEE, thus university students and campuses may present an opportunity to harness this potential. University Distinct Urban Mines (DUM) may be used to exemplify how potential reserves of secondary metals may be exploited, and could contribute to the transition from a linear to a circular economy. This study aimed to evaluate small household appliances (SHA) DUM from a UK university, with the objectives to identify and quantify student households’ SHA ownership, WEEE recycling, stockpiling and discarding habits amongst student households, assess and evaluate the monetary potential of SHA DUM at UK level, and propose methods to exploit DUM for universities in the UK. To this purpose, a quantitative survey was undertaken to measure students’ ownership and discarding behaviour with respect to SHA. The amounts of ferrous and non-ferrous materials were then estimated and converted to monetary values from secondary materials market data to appraise the SHA DUM overall value. Thirty-five per cent of SHA are discarded in the general refuse. Broken personal care appliances (PCA) tend to be discarded due to hygiene and small size factors. When in working order, SHA tend to be equally reused, recycled or stockpiled. We conclude that a total of 189 tonnes of ferrous and non-ferrous materials were available via discarding or being stockpiled at the University of Southampton. Extrapolated to UK higher education level, discarded and stockpiled SHA represent a potential worth ∼USD 11 million. To initiate DUM exploitation within Higher Education campuses, we suggest improving users’ choice architecture by providing collection methods specific to broken SHA.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

L3:RTM.MCH.P19.06: Neutronics UQ

This report summarizes some of the work performed in the last year on the nuclear data uncertainty tasks using the multipole representation. This approach was proposed to provide a simple and integrated way of evaluating the nuclear data uncertainty in a Monte Carlo simulation and evaluating their impact on multigroup cross sections. As previously described, the windowed multipole format allows for a resonance based representation of the cross sections with direct Doppler broadening. This representation also lends itself nicely to the evaluation of sensitivities of resonance parameters and the direct perturbation of resonance parameters in the simulation environment. Both of these approaches rely on the need of conversion of the covariance data (File 32 of the ENDF library) into a pole and residues form: i.e. into a multipole covariance. Section 5 details the current pathways that have been used so far to achieve this. In order to verify the conversion process of the covariance data, the sensitivies to the resonance parameters and the sampling from these parameters, an analytical benchmark of the slowing process was developed. The benchmark provides a closed-form solution for the flux in the presence of Hydrogen and multiple resonant nuclides. Section 2 details the analytical benchmark and provides the closed-form solutions for the flux and its sensitivity to resonance parameters. A verification of the proposed methodology is performed and compared with the analytic solution to demonstrate the accuracy of the benchmark and the applicability of our embedded resonance perturbation approach in a small Monte Carlo code. Papers are currently being drafted documenting the benchmark and the generation of sensitivities. As paths to best generate multipole covariance matrices are being compared, we also explored new methods to propagate nuclear data uncertainty across the Monte Carlo neutron transport calculation. Section 4 studies the use of multipole covariance matrices to perform sensitivity analysis. Section 3.2 documents the results of our investigations into a new way of propagating nuclear data uncertainty across Monte Carlo transport simulations by sampling new epistemic cross sections for each new neutron.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Risk Model for EM Decision Support Toolsets

The Government Office of Accountability (GAO) has published several reports identifying the need for the Department of Energy (DOE) Office of Environmental Management (EM) to address mounting costs for DoE's cleanup program. DOEEM could greatly benefit from independent decision tool-sets/models that would allow them to evaluate options and inform business decision at the enterprise level considering site-specific life cycle costs and system plans. Program goals: Develop a tool-set that enables EM to independently evaluate alternatives, assess outcomes from different contracting strategies, and inform critical decisions for the enterprise. Project goals: Adaption of a Risk Model for integration with complimentary tool-sets for project level decision making. Operational Events: Discrete event model that evaluates operational variables (e.g. capacity, throughput, maintenance constraints) and identify bottlenecks. Identifying and Bounding Project Risk: Identify risks that impact confidence in meeting goals (e.g. production). Life cycle Cost: Evaluate impacts of staffing levels, inventory, and capital investments on life cycle costs. Methods and approach: Monte Carlo Analysis is being executed to generate results: Input derives from Risk Register Data; Tied to Projected and Target Schedules; Incorporates float duration within the model. Assumes associated risk mitigation actions being completed within a timeline of five fiscal years. Metrics include: Confidence in Meeting Production Goal; Confidence in Safety Standards; Confidence in Continuous Operation; Other Metrics can be added as appropriate regarding specific site needs. The adapted risk model can be used as a stand alone decision tool or can be integrated complementary tool-sets (i.e. process and cost models) for project-specific decisions. These support tool-sets can then be integrated with others for site- and complex- level evaluations. Future work includes designing an adaptable and modular framework that would allow integration of multiple tool-sets for holistic and/or targeted evaluation of alternative strategies for decision making that could lead to risk and cost reduction across the enterprise.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Evaluating 239 Pu(n,f) cross sections via machine learning using experimental data, covariances, and measurement features

In this paper, the neutron-induced 239 Pu fission cross section, 239 Pu(n,f), is evaluated from 1–20 MeV using experimental data and associated covariances while also considering information on the measurement, termed features here. For instance, methods to determine the background, sample backing material, or impurities in the sample, are explicitly taken into account in the evaluation process. To this end, outliers in the experimental data are identified with a modified version of the Hybrid Robust Support Vector Machine. In a second step, two machine learning methods (logistic regression with elastic net regularization and random forest regression with SHAP feature importance metric) are used to highlight measurement features that are common among many of the outlying data points. Based on this analysis, penalty uncertainties are added to the experimental covariances of outlying data points that have outlier measurement features and are put through the generalized-least-squares evaluation. The resulting evaluated mean values and covariances differ distinctly from those data evaluated without the penalty uncertainties. These results highlight that certain measurement features should be more closely examined.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

nCTEQ15 - Global analysis of nuclear parton distributions with uncertainties

We present the first official release of the nCTEQ nuclear parton distribution functions with errors. The main addition to the previous nCTEQ PDFs is the introduction of PDF uncertainties based on the Hessian method. Another important addition is the inclusion of pion production data from RHIC that give us a handle on constraining the gluon PDF. This contribution summarizes our results from arXiv:1509.00792 and concentrates on the comparison with other groups providing nuclear parton distributions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Absolute mass calibration of fission product distributions measured with the E-$v$ method

The mass calibration of fission product distributions measured with the energy–velocity (E-v) method constitutes a technical challenge. The energy loss of the fission fragments in the various dead layers of the spectrometer and other sources of pulse-height defects of the energy detectors are a significant source of systematic uncertainty in the mass calculations used to determine the measured fission product yields. Here, in this study, the absolute calibration of the fission mass distributions is accomplished by measuring prompt γ-rays in coincidence with the fission fragments. This allows for the direct calibration of the mass spectra and eliminates any dependence on potentially complicated energy corrections. The first test measurement was performed at the Los Alamos Neutron Science Center employing the SPectrometer for Ion DEtermination in fission Research (SPIDER). SPIDER is a 2E-2v spectrometer designed for measuring independent fission product yields from neutron-induced fission. In this test, the single-arm SPIDER (E-v) system and an array of 252 Cf sources were used. The single-arm system consisted of two time pick-off detectors for measuring the time-of-flight of the fission fragments and a double-sided silicon strip detector (DSSD) for measuring the kinetic energy. Characteristic γ-rays from fission fragments were detected using three high-purity germanium (HPGe) detectors. For the mass calibration, γ–mass coincidence events from twelve product isotopes were used. The measured FPYs from 252 Cf spontaneous fission were found to be in excellent agreement with the evaluated data after applying the absolute mass calibration. From the γ–mass coincidence events, the mass resolution of the system was also extracted. An average mass resolution of ~1.4 AMU (FWHM) for the light fragments and ~2.4 AMU (FWHM) for the heavy fragments were found. This was the first in-situ calibration and characterization of the SPIDER spectrometer, which paves the way for high-quality FPY measurements with this instrument.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

SU-E-J-191: Automated Detection of Anatomic Changes in H'N Patients

Purpose: To develop a novel statistics-based method for automated detection of anatomical changes using cone-beam CT data. A method was developed that can provide a reliable and automated early warning system that enables a “just-in-time” adaptation of the treatment plan. Methods: Anatomical changes were evaluated by comparing the original treatment planning CT with daily CBCT images taken prior treatment delivery. The external body contour was computed on a given CT slice and compared against the corresponding contour on the daily CBCT. In contrast to threshold-based techniques, a statistical approach was employed to evaluate the difference between the contours using a given confidence level. The detection tool used the two-sample Kolmogorov-Smirnov test, which is a non-parametric technique that compares two samples drawn from arbitrary probability distributions. 11 H'N patients were retrospectively selected from a clinical imaging database with a total of 186 CBCT images. Six patients in the database were confirmed to have anatomic changes during the course of radiotherapy. Five of the H'N patients did not have significant changes. The KS test was applied to the contour data using a sliding window analysis. The confidence level of 0.99 was used to moderate false detection. Results: The algorithm was able to correctly detect anatomical changes in 6 out of 6 patients with an excellent spatial accuracy as early as at the 14th elapsed day. The algorithm provided a consistent and accurate delineation of the detected changes. The output of the anatomical change tool is easy interpretable, and can be shown overlaid on a 3D rendering of the patient's anatomy. Conclusion: The detection method provides the basis for one of the key components of Adaptive Radiation Therapy. The method uses tools that are readily available in the clinic, including daily CBCT imaging, and image co-registration facilities.

60 APPLIED LIFE SCIENCES↗

Use of Percutaneous Aspiration Thrombectomy vs. Anticoagulation Therapy to Treat Acute Iliofemoral Venous Thrombosis: 1-year Follow-up Results of a Randomised, Clinical Trial

PurposeThe purpose of this study was to compare the efficacy of percutaneous aspiration thrombectomy (PAT) followed by standard anticoagulant therapy, with anticoagulation therapy alone, for the treatment of acute proximal lower extremity deep vein thrombosis.MethodsIn this randomised, prospective study, 42 patients with acute proximal iliofemoral deep vein thrombosis documented via Doppler ultrasound examination, were separated into an interventional treatment group (16 males, 5 females, average age 51 years) and a medical treatment group (13 males, 8 females, average age 59 years). In the interventional group, PAT with large-lumen 9-F diameter catheterisation was applied, after initiation of standard anticoagulant therapy. Balloon angioplasty (n 19) and stent implementation (n: 14) were used to treat patients with residual stenosis (>50 %) after PAT. Prophylactic IVC filters were placed in two patients. The thrombus clearance status of the venous system was evaluated by venography. In both the medical and interventional groups, venous patency rates and clinical symptom scores were evaluated at months 1, 3, and 12 after treatment.ResultsDeep venous systems became totally cleared of thrombi in 12 patients treated with PAT. The venous patency rates in month 12 were 57.1 and 4.76 % in the interventional and medical treatment groups, respectively. A statistically significant improvement was observed in clinical symptom scores of the interventional group (PAT) with or without stenting (4.23 ± 0.51 before treatment; 0.81 ± 0.92 at month 12) compared with the medical treatment group (4.00 ± 0.63 before treatment; 2.43 ± 0.67 at month 12). During follow-up, four patients in the medical treatment and one in the interventional group developed pulmonary embolisms.ConclusionsFor treatment of acute deep vein thrombosis, PAT with or without stenting is superior to anticoagulant therapy alone in terms of both ensuring venous patency and improving clinical symptoms. PAT is a safe, inexpensive, and easily performed method of endovascular treatment with a low rate of major complications. Our present findings and literature data suggest that PAT can be used as first-line treatment in proximal deep vein thrombosis patients, especially when thrombolytic treatment is contraindicated.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

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 ↗