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 811 records · Page 45

Development of Fusion Neutronics Methodology, Software, and Testing Facility in FDS Team

Fusion neutronics including methodology, software, experiments, etc. is studied concerning the generation, utilizing and protection of neutrons in fusion facilities. Achievements have been made in the past decades to solve the new problems fusion neutronics. However there is still some way to go such as following. Automatic modeling method based on CAD has been developed to avoid the heavy burden of manual geometry description. The pre-processing costs much time to simplify the CAD model for converting to Constructive Solid Geometry (CSG), 2D or 3D mesh based geometry. This problem has become a new focus. Deep penetration of radiation shielding is notable problem for fusion. Adaptive variance reduction methods such as Consistent Adjoint Driven Importance Sampling (CADIS) were developed to accelerate convergence. However the gap still exists between engineering requirements and current methodology. The geometry simplification, homogenizing based tally, error in nuclear data and error propagation in coupled calculation, etc., causes considerable uncertainty, which leads to 'safety factor' applied in simulation results to achieve required safety margin. It can significantly raise the cost and even makes it impractical to initiate the fusion facilities. Uncertainty qualification, error analysis and error propagation problem should be studied. New methodology and software should be validated. For example the (n,p), (n,a), (n,t) cross section of Cu, Fe, W, Ni which is important for evaluation of material radiation damage is absent. The validation of neutronics properties of special components such as first wall, blanket and coil are needed. The D-T fusion neutron generator can directly represent the neutron environment in fusion reactors. However the neutron yield is several orders of magnitudes lower than requirements. To resolve obstacles in fusion neutronics, broad improvements of methodology, software and testing facility have been achieved by FDS Team in recent years and the main progress was introduced in this paper.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Spectroscopy of states in 136 Ba using the 138 Ba(p,t) reaction

Background: The 136 Ba isotope is the daughter nucleus in 136 Xe $ββ$ decay. It also lies in a shape transitional region of the nuclear chart, making it a suitable candidate to test a variety of nuclear models. Purpose: To obtain spectroscopic information on states in 136 Ba, which will allow a better understanding of its low-lying structure. These data may prove useful to constrain future 136 Xe → 136 Ba neutrinoless $ββ$ decay matrix element calculations. Methods: Here, a 138 Ba(p, t) reaction was used to populate states in 136 Ba up to approximately 4.6 MeV in excitation energy. The tritons were detected using a high-resolution Q3D magnetic spectrograph. A distorted wave Born approximation analysis was performed for the measured triton angular distributions. Results: 102 excited states in 136 Ba were observed, out of which 52 are reported for the first time. Definite spin-parity assignments are made for 26 newly observed states, while previously ambiguous assignments for ten other states are resolved. Together with other available data, the results are used to determine level densities in 136 Ba. These were compared with theory predictions, obtained using shell model calculations with Hamiltonians previously used for 136 Xe neutrinoless $ββ$ decay matrix element evaluations. Conclusions: The shell model predicted level densities agree reasonably well for the two Hamiltonians. However the results for theory and experiment are found to agree only at lower energies, diverging from one another for the higher lying states, with the discrepancy increasing with energy. This is presumably because of lower production cross sections for a majority of the higher-lying predicted states and the experimental limitations in resolving a large number of nearly degenerate states predicted by the theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Analytical capabilities for iodine detection: Review of possibilities for different applications

This Review summarizes a range of analytical techniques that can be used to detect, quantify, and/or distinguish between isotopes of iodine (e.g., long-lived 129 I, short-lived 131 I, stable 127 I). One reason this is of interest is that understanding potential radioiodine release from nuclear processes is crucial to prevent environmental contamination and to protect human health as it can incorporate into the thyroid leading to cancer. It is also of interest for evaluating iodine retention performances of next-generation iodine off-gas capture materials and long-term waste forms for immobilizing radioiodine for disposal in geologic repositories. Depending upon the form of iodine (e.g., molecules, elemental, and ionic) and the matter state (i.e., solid, liquid, and gaseous), the available options can vary. In addition, several other key parameters vary between the methods discussed herein, including the destructive vs nondestructive nature of the measurement process (including in situ vs ex situ measurement options), the analytical data collection times, and the amount of sample required for analysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The role of material modeling on strain range estimation for elevated temperature cyclic life evaluation

High temperature nuclear reactors operating in the creep regime are designed to withstand numerous cyclic events. Current ASME code rules provide two basic paths for evaluating creep fatigue and ratcheting under these conditions; one based on full inelastic analysis intended to provide a representative stress and strain history and the other based on elastic material models with adjustments of varying complexity to account for inelastic stress and strain redistribution. More recent developments have used elastic-perfectly plastic analysis to bound the effects of cyclic service. However, these methods still rely on the separate evaluation of fatigue and creep damage utilizing a damage interaction diagram. There is a procedure under current development that uses creep-fatigue data from key feature test articles directly without the use of the damage interaction diagram. However, it requires a reasonable representation of the strain range in a structure as an input. This work develops a simplified procedure based on elastic perfectly-plasticity analysis that can be used to represent the strain range in a structure in the steady state under cyclic loading conditions.

Messner, M. C.↗

SU-F-T-266: Dynalogs Based Evaluation of Different Dose Rate IMRT Using DVH and Gamma Index

Purpose: This work investigates the impact of low and high dose rate on IMRT through Dynalogs by evaluating Gamma Index and Dose Volume Histogram. Methods: The Eclipse™ treatment planning software was used to generate plans on prostate and head and neck sites. A range of dose rates 300 MU/min and 600 MU/min were applied to each plan in order to investigate their effect on the beam ON time, efficiency and accuracy. Each plan had distinct monitor units per fraction, delivery time, mean dose rate and leaf speed. The DVH data was used in the assessment of the conformity and plan quality.The treatments were delivered on Varian™ Clinac 2100C accelerator equipped with 120 leaf millennium MLC. Dynalogs of each plan were analyzed by MATLAB™ program. Fluence measurements were performed using the Sun Nuclear™ 2D diode array and results were assessed, based on Gamma analysis of dose fluence maps, beam delivery statistics and Dynalogs data. Results: Minor differences found by adjusted R-squared analysis of DVH’s for all the plans with different dose rates. It has been also found that more and larger fields have greater time reduction at high dose rate and there was a sharp decrease in number of control points observed in dynalog files by switching dose rate from 300 MU/min to 600 MU/min. Gamma Analysis of all plans passes the confidence limit of ≥95% with greater number of passing points in 300 MU/min dose rate plans. Conclusion: The dynalog files are compatible tool for software based IMRT QA. It can work perfectly parallel to measurement based QA setup and stand-by procedure for pre and post delivery of treatment plan.

60 APPLIED LIFE SCIENCES↗

THE ROLE OF MATERIAL MODELING ON STRAIN RANGE ESTIMATION FOR ELEVATED TEMPERATURE CYCLIC LIFE EVALUATION

High temperature nuclear reactors operating in the creep regime are designed to withstand numerous cyclic events. Current ASME code rules provide two basic paths for evaluating creep fatigue and ratcheting under these conditions; one based on full inelastic analysis intended to provide a representative stress and strain history and the other based on elastic material models with adjustments of varying complexity to account for inelastic stress and strain redistribution. More recent developments have used elastic-perfectly plastic analysis to bound the effects of cyclic service. However, these methods still rely on the separate evaluation of fatigue and creep damage utilizing a damage interaction diagram. There is a procedure under current development that uses creep-fatigue data from key feature test articles directly without the use of the damage interaction diagram. However, it requires a reasonable representation of the strain range in a structure as an input. This work develops a simplified procedure based on elastic perfectly-plasticity analysis that can be used to represent the strain range in a structure in the steady state under cyclic loading conditions.

Messner, Mark C.↗

Automated Waterbox Inspection for Nuclear Power Plants Using Computer Vision - Based Change Detection

Nuclear power plant waterboxes require regular inspection for leaks, missing components, and structural damage during maintenance outages. Traditional manual inspection is time-consuming and poses safety risks from confined space entry. We developed an automated computer vision system for drone-based waterbox inspection in partnership with Florida Light and Power. Our approach uses feature detection and matching to identify critical changes between baseline and current inspection images, automatically flagging additions (leaks/debris), removals (missing plugs), and translations (displaced components) while compensating for drone movement and environmental variations. We systematically evaluated six feature matching methods, from classical approaches (SIFT+BF) to state-of-the-art neural networks (SuperPoint+SuperGlue), using both standard benchmarks (HPatches) and waterbox-specific validation with real-world augmentations. SuperPoint+SuperGlue achieved superior performance with 7.82 pixels RMSE and 100% success rate—2.8x better accuracy than our baseline. While the pre-trained model has commercial licensing restrictions for nuclear deployment, our findings validate this architecture for custom training. We implemented a real-time GUI demonstrating the SIFT+BF approach for immediate deployment, processing drone feeds at 30 FPS with color-coded change visualization. Future work includes training a custom SuperPoint+SuperGlue model on waterbox data and integrating Vision-Language Models for automated reporting and maintenance guidance.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗

Update on Covariance Data Testing Strategy at LANL [Slides]

LANL is working towards an ENDF/B-VIII.0-based Covariance Library, with several key goals and work processes outlined. This includes processing through NJOY’s ERRORR module, identifying and correcting mathematical and physical deficiencies, communicating across pipeline from evaluator to end user, understanding use cases and interpreting results, and releasing to customers. Their testing approach includes interaction, processing, checks (mathematical properties, constraints, and physical bounds), and error propagation.

97 MATHEMATICS AND COMPUTING↗

Off-Gas Adsorption Model Capabilities and Recommendations

Off-gas treatment is required to reduce emissions from aqueous fuel reprocessing. Evaluating the products of innovative gas adsorption research requires increased computational simulation capability to more effectively transition from fundamental research to operational design. Early modeling efforts produced the Off-Gas SeParation and REcoverY (OSPREY) model that, while efficient in terms of computation time, was of limited value for complex systems. However, the computational and programming lessons learned in development of the initial model were used to develop Discontinuous Galerkin OSPREY (DGOSPREY), a more effective model. Initial comparisons between OSPREY and DGOSPREY show that, while OSPREY does reasonably well to capture the initial breakthrough time, it displays far too much numerical dispersion to accurately capture the real shape of the breakthrough curves. DGOSPREY is a much better tool as it utilizes a more stable set of numerical methods. In addition, DGOSPREY has shown the capability to capture complex, multispecies adsorption behavior, while OSPREY currently only works for a single adsorbing species. This capability makes DGOSPREY ultimately a more practical tool for real world simulations involving many different gas species. While DGOSPREY has initially performed very well, there is still need for improvement. The current state of DGOSPREY does not include any micro-scale adsorption kinetics and therefore assumes instantaneous adsorption. This is a major source of error in predicting water vapor breakthrough because the kinetics of that adsorption mechanism is particularly slow. However, this deficiency can be remedied by building kinetic kernels into DGOSPREY. Another source of error in DGOSPREY stems from data gaps in single species, such as Kr and Xe, isotherms. Since isotherm data for each gas is currently available at a single temperature, the model is unable to predict adsorption at temperatures outside of the set of data currently available. Thus, in order to improve the predictive capabilities of the model, there is a need for more single-species adsorption isotherms at different temperatures, in addition to extending the model to include adsorption kinetics. This report provides background information about the modeling process and a path forward for further model improvement in terms of accuracy and user interface.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Total Mercury Analysis Comparison: Deployment of Analytical Method for the Savannah River Site Liquid Waste System

A direct mercury analyzer (DMA; EPA 7473) and Cold Vapor Atomic Absorption Spectroscopy (CVAAS; EPA 245.2) have been used to analyze samples high in methylmercury and other forms of mercury with oxidation states of Hg o (elemental), Hg 1+ (mercurous), and/or Hg 2+ (mercuric) to determine analytical data equivalency for total mercury. The DMA instrument has been determined to be the preferred method based on both simulated and radioactive sample testing due to analytical results and ease of use in a radiological hood. SRNL made a number of modifications to the DMA instrument and analysis method to limit the amount of time personnel are required to work in a radiological environment. Equivalency is based on comparative statistical evaluation of the results from: spike solutions/simulants, Tank 50 samples, Mixed Analyte Performance Evaluation Program (MAPEP) samples, and Modular Caustic-Side Solvent Extraction Unit (MCU) solvent samples. SRNL performed a cross-comparison of results using CVAAS, DMA, and Inductively Coupled Plasma Mass Spectrometer (ICP-MS), as well as independent analysis from an offsite laboratory (using cold vapor sample processing with atomic fluorescence detection, EPA 1631). All data were within the combined uncertainty (2σ ± 20%) of the CVAAS and DMA methods

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

KilonovaNet : Surrogate models of kilonova spectra with conditional variational autoencoders

ABSTRACT Detailed radiative transfer simulations of kilonova spectra play an essential role in multimessenger astrophysics. Using the simulation results in parameter inference studies requires building a surrogate model from the simulation outputs to use in algorithms requiring sampling. In this work, we present kilonovanet, an implementation of conditional variational autoencoders (cVAEs) for the construction of surrogate models of kilonova spectra. This method can be trained on spectra directly, removing overhead time of pre-processing spectra, and greatly speeds up parameter inference time. We build surrogate models of three state-of-the-art kilonova simulation data sets and present in-depth surrogate error evaluation methods, which can in general be applied to any surrogate construction method. By creating synthetic photometric observations from the spectral surrogate, we perform parameter inference for the observed light-curve data of GW170817 and compare the results with previous analyses. Given the speed with which kilonovanet performs during parameter inference, it will serve as a useful tool in future gravitational wave observing runs to quickly analyse potential kilonova candidates.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Probing Interactions between Adsorbed Polymers and γ-Al 2 O 3 Using Solid-State NMR

The interactions between polymers and metal oxide surfaces affect the physical properties of polymer composites and the selectivity of catalysts used for plastic conversion. Gauging the strengths of these interactions and how they depend on adsorption conformations and polymer structure is nevertheless quite challenging. Herein, we evaluate whether a recently introduced dynamic nuclear polarization-based surface-to-atom distance measurement method can be used to gauge the specific interactions between five vinyl polymers and an oxide support (alumina). The adsorption conformations of polymer monolayers were studied using site-specific 13 C{ 27 Al} rotational-echo saturation-pulse double-resonance (RESPDOR) experiments. The method is effective in revealing the average distance between the adsorbed polymer and the surface, as well as the fraction of the monomers that are adsorbed. Polymers that possess larger side chains are more weakly adsorbed, as evidenced by their formation of more numerous loop and tail structures. In conclusion, the measured RESPDOR dephasing data further reveals that polymer adsorption conformations are largely disordered and involve few hydrogen bonding interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simulations of a PSD Plastic Neutron Collar for Assaying Fresh Fuel

The potential performance of a notional active coincidence collar for assaying uranium fuel based on segmented detectors constructed from the new PSD plastic fast organic scintillator with pulse shape discrimination capability was investigated in simulation. Like the International Atomic Energy Agency's present Uranium Neutron Collar for LEU (UNCL), the PSD plastic collar would also function by stimulating fission in the 235 U content of the fuel with a moderated 241 Am/Li neutron source and detecting instances of induced fission via neutron coincidence counting. In contrast to the moderated detectors of the UNCL, the fast time scale of detection in the scintillator eliminates statistical errors due to accidental coincidences that limit the performance of the UNCL. However, the potential to detect a single neutron multiple times historically has been one of the properties of organic scintillator detectors that has prevented their adoption for international safeguards applications. Consequently, as part of the analysis of simulated data, a method was developed by which true neutron-neutron coincidences can be distinguished from inter-detector scatter that takes advantage of the position and timing resolution of segmented detectors. Then, the performance of the notional simulated coincidence collar was evaluated for assaying a variety of fresh fuels, including some containing burnable poisons and partial defects. In these simulations, particular attention was paid to the analysis of fast mode measurements. In fast mode, a Cd liner is placed inside the collar to shield the fuel from the interrogating source and detector moderators, thereby eliminating the thermal neutron flux that is most sensitive to the presence of burnable poisons that are ubiquitous in modern nuclear fuels. The simulations indicate that the predicted precision of fast mode measurements is similar to what can be achieved by the present UNCL in thermal mode. For example, the statistical accuracy of a ten-minute measurement of fission coincidences collected in fast mode will be approximately 1% for most fuels of interest, yielding a ~1.4% error after subtraction of a five minute measurement of the spontaneous fissions from 238 U in the fuel, a ~2% error in analyzed linear density after accounting for the slope of the calibration curve, and a ~2.9% total error after addition of an assumed systematic error of 2%.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NSRD-16: Computational Capability to Substantiate DOE-HDBk-3010 Data

Safety basis analysts throughout the U.S. Department of Energy (DOE) complex rely heavily on the information provided in the DOE Handbook, DOE-HDBK-3010, Airborne Release Fractions/Rates and Respirable Fractions for Nonreactor Nuclear Facilities, to determine radionuclide source terms from postulated accident scenarios. In calculating source terms, analysts tend to use the DOE Handbook's bounding values on airborne release fractions (ARFs) and respirable fractions (RFs) for various categories of insults (representing potential accident release categories). This is typically due to both time constraints and the avoidance of regulatory critique. Unfortunately, these bounding ARFs/RFs represent extremely conservative values. Moreover, they were derived from very limited small-scale bench/laboratory experiments and/or from engineered judgment. Thus, the basis for the data may not be representative of the actual unique accident conditions and configurations being evaluated. The goal of this research is to develop a more accurate and defensible method to determine bounding values for the DOE Handbook using state-of-art multi-physics-based computer codes. This enables us to better understand the fundamental physics and phenomena associated with the types of accidents in the handbook. In this fourth year, we improved existing computational capabilities to better model fragmentation situations to capture small fragments during an impact accident. In addition, we have provided additional new information for various sections of Chapters 4 and 5 of the Handbook on free fall powders and impacts of solids, and have provided the damage ratio simulations for containers (7A drum and standard waste box) for various drops and impact scenarios. Thus, this work provides a low-cost method to establish physics-justified safety bounds by considering specific geometries and conditions that may not have been previously measured and/or are too costly to perform during an experiment.

61 RADIATION PROTECTION AND DOSIMETRY↗

Lost and Found Opportunities Around the Chlorine Worth Study

Los Alamos National Laboratory performed a series of critical experiments in 2021 to examine the worth of chlorine in plutonium-fueled systems. This series of experiments has been dubbed the “Chlorine Worth Study,” and the evaluation of the experiments was presented to the International Criticality Safety Benchmark Evaluation Project (ICSBEP) Technical Review Group in April, 2023. The primary purpose of these experiments was to enable validation of aqueous solutions crediting neutron absorption in 35 Cl. An external, independent view of the events leading up to the design and execution of these experiments indicates a missed opportunity to leverage sensitivity/uncertainty (S/U) analysis to assert validation without the experiments by taking an additional margin for the lack of direct validation of chlorine. On the other hand, the execution of these experiments also presents a rare opportunity to examine the efficacy of the S/U approach and extract useful information about the evaluated chlorine covariance data. TSUNAMI-1D models of representative application solutions were created and used to generate sensitivity data. Varying plutonium and chlorine concentrations were considered to examine the impact of these differences on the chlorine sensitivities and uncertainties. The data-induced uncertainty in k eff resulting from chlorine was calculated directly from uncertainty information calculated in the TSUNAMI-1D sequence. In all cases, this uncertainty was less than 0.1 %Δk. This result could potentially be used to justify a reactivity margin to account for the validation gap related to chlorine in the validation set. On the other hand, given that the experiments were performed, the community should endeavor to extract as much value from them and their results as possible. The results can be used to examine the actual bias associated with chlorine in these systems once the evaluations have been released. These data can be compared with the data-induced uncertainty margin discussed above to test the sufficiency of the validation gap penalty. This result will provide an indication of the performance of the chlorine covariance data specifically and the S/U validation approach generally. More advanced S/U techniques may also be employed to determine reactivity sensitivities associated with the chlorine in the experiments, potentially generating a more robust test of the chlorine covariance data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Infrared spectroscopic method for uranium isotopic analysis

A high performance infrared (HPIR) system was developed and demonstrated for the infrared absorption analysis of {sup 235}U and {sup 238}U isotopes in uranium hexafluoride gas samples. Sweeping the quantum cascade laser light source over the spectral range and sampling via a high-rate analog-to-digital converter provided 0.0005 cm-1 spectral resolution, which allowed for high-precision measurements of the isotopic peak shift. A data analysis method was developed using principal component analysis to predict the isotope weight % content of {sup 235}U. The HPIR precision, accuracy, and error were evaluated for a wide range of isotopic ratio samples (0.287 – 93.7 weight % {sup 235}U), and the results were compared to the International Target Values (ITVS) set forth by the International Atomic Energy Agency (IAEA) for non-destructive and destructive analytical techniques. The method meets or surpasses the IAEA ITVs for non-destructive analysis of samples with isotopic content of depleted to highly enriched. The results also demonstrated the capability of the HPIR system to correctly predict the {sup 235}U weight % content of a mislabeled sample whose isotopic distribution was validated by mass spectroscopic measurements. The HPIR measurement is nondestructive and, thus, allows for confirmatory analyses of the exact sample at a designated IAEA lab if higher-resolution or a certified analysis is needed. (author)

07 ISOTOPE AND RADIATION SOURCES↗

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 ↗