Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “core model proposal”

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

Uncertainty Quantification of Spatial Correction Factor for Sjoestrand Method due to Cross-section Data

Accelerator-driven System (ADS), which couples a subcritical core with a spallation neutron source, is proposed as a nuclear transmutation technology to reduce the radioactive half-life of high-level waste. To maintain a safe operation in ADS, the core should be kept in subcritical state. Thus, the subcritical state of the core must be quantified through measurement of subcriticality, which is a quantitative margin from critical state. The area ratio method, or the Sjoestrand method, is one of the subcriticality measurement techniques applicable to the ADS. In this method, the absolute value of subcriticality can be obtained by analyzing the temporal variation of neutron density in the core due to the periodic injection of neutron pulses. In practice, the value of subcriticality is measured by the temporal variation of neutron count rate detected by a detector in the core. Time integrals of count rates due to prompt and delayed neutrons A{sub p} and A{sub d}, are measured. Then, the ratio of A{sub p} to A{sub d}, which is called as the area ratio, yields negative reactivity -ρ, i.e. subcriticality, in dollar units: A{sub p}/A{sub d} ≅ -ρ/β{sub eff'} (1) where β{sub eff} is the effective delayed neutron fraction. Equation (1) represents the basic idea of the Sjoestrand method. Note, however, that Eq. (1) is approximately derived using the point kinetics model. In an actual subcritical system, the measured area ratio is not necessarily equal to the exact value of subcriticality in dollar units, due distortion by spatial higher-order modes. Hence, in the actual measurement, this spatial effect is corrected by a spatial correction factor such as the Bell and Glasstone correction factor (the Bell factor), or by extraction of the fundamental mode (so-called the extrapolated area ratio method). In addition, there is an effect due to the neutron detector materials (e.g. {sup 3}He(n,p) and {sup 235}U(n,f)). These effects must be considered in the correction of the measured area ratio. In this study, we focused on the uncertainty of Bell factor. The Bell factor f is defined by the ratio of the area ratio to the subcriticality as follows: f ≡ -[ρ/β{sub eff}]/[A{sub p}/A{sub d}] (2) The Bell factor f can be evaluated by a core analysis. By multiplying the measured area ratio by the Bell factor?, the spatial higher-order mode effect in the area ratio can be corrected. However, the uncertainty of the Bell factor due to uncertainties of cross-section data (covariance data) has not been quantitatively evaluated so far. In this study, we propose a method to quantify the uncertainty of the Bell factor in an actual core analysis. As a feasibility study, the uncertainty quantification is carried out in a simple one-dimensional subcritical core system. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Line formation in collapsing cloud cores with rotation and applications to B335 and IRAS 16293-2422

We present radiative transfer models of collapsing clouds with rotation and apply them to B335 and IRAS 16293-2422, the two best candidates of protostellar collapse known to date. By including rotation in the manner of Tereby, Shu, & Cassen (1984), we can reproduce the profiles of several molecular lines not only toward the center position, but also toward a grid of positions near the center. We find that the model of B335 by Zhou et al. (1993) is not affected significantly by the presence of rotation. By including rotation, we can reproduce the observations of IRAS 16293-2422 by Menden et al. (1987); hence, we support the infall interpretations of spectral lines from IRAS 16293-2422 proposed by Walker et al. (1986). We also observed the large-scale rotation of IRAS 16293-2422 in the C(18)O J = 2 approaches 1 line. The observed rotation rate is a factor of 6 smaller than that required to explain the small-scale CS data of Menden et al. (1987). This probably means the precollapse cloud has differential rotation, possibly due to more efficient magnetic braking on large scales.

Zhou, Shudong↗

THE NUCLEAR SPECTROSCOPIC TELESCOPE ARRAY (NuSTAR) HIGH-ENERGY X-RAY MISSION

The Nuclear Spectroscopic Telescope Array (NuSTAR) mission, launched on 2012 June 13, is the first focusing high-energy X-ray telescope in orbit. NuSTAR operates in the band from 3 to 79 keV, extending the sensitivity of focusing far beyond the {approx}10 keV high-energy cutoff achieved by all previous X-ray satellites. The inherently low background associated with concentrating the X-ray light enables NuSTAR to probe the hard X-ray sky with a more than 100-fold improvement in sensitivity over the collimated or coded mask instruments that have operated in this bandpass. Using its unprecedented combination of sensitivity and spatial and spectral resolution, NuSTAR will pursue five primary scientific objectives: (1) probe obscured active galactic nucleus (AGN) activity out to the peak epoch of galaxy assembly in the universe (at z {approx}< 2) by surveying selected regions of the sky; (2) study the population of hard X-ray-emitting compact objects in the Galaxy by mapping the central regions of the Milky Way; (3) study the non-thermal radiation in young supernova remnants, both the hard X-ray continuum and the emission from the radioactive element {sup 44}Ti; (4) observe blazars contemporaneously with ground-based radio, optical, and TeV telescopes, as well as with Fermi and Swift, to constrain the structure of AGN jets; and (5) observe line and continuum emission from core-collapse supernovae in the Local Group, and from nearby Type Ia events, to constrain explosion models. During its baseline two-year mission, NuSTAR will also undertake a broad program of targeted observations. The observatory consists of two co-aligned grazing-incidence X-ray telescopes pointed at celestial targets by a three-axis stabilized spacecraft. Deployed into a 600 km, near-circular, 6 Degree-Sign inclination orbit, the observatory has now completed commissioning, and is performing consistent with pre-launch expectations. NuSTAR is now executing its primary science mission, and with an expected orbit lifetime of 10 yr, we anticipate proposing a guest investigator program, to begin in late 2014.

79 ASTRONOMY AND ASTROPHYSICS↗

Methodology and Tool for the Physical Security Analysis of Micro and Advanced Reactors

This work proposes a dynamic evaluation methodology to relax the conservatism in physical security evaluation, by leveraging an ongoing work in the Light Water Reactor Sustainability pathway. This methodology is implemented in a dynamic risk assessment tool named Event Modeling Risk Assessment using Linked Diagrams (EMRALD). The work extends EMRALD’s capability to support a sandbox feature where analysts can easily create attack scenarios and modify advanced/small modular reactor (A/SMR) security and safety features using templates. This approach saves time and cost since the analysis does not require creating detailed computer-aided design models, as is commonly required in commercial force-on-force software tools. EMRALD is completely free to use at https://emraldapp.inl.gov. We have developed basic templates including physical barriers, intrusion sensors, physical areas, and safety actions, that can be downloaded from EMRALD’s GitHub site: https://github.com/idaholab/EMRALD. These templates use generic data commonly used for training purposes, which do not reflect any actual operating nuclear reactor. Users may adjust the data in the templates with their own dataset and/or create new templates in EMRALD. The proposed methodology combines security and safety by assessing sabotage effects up to the radiological consequence to the public instead of merely the core damage state. This practice follows the industry standard for advanced non-light-water reactors currently proposed for endorsement by the Nuclear Regulatory Commission. The combination of security and safety is expressed in an achievability-consequence chart. EMRALD can be used to generate data for this chart. A hypothetical case study using a representative sodium-cooled fast reactor (SFR) facility is presented in this report to demonstrate this methodology. This case study does not contain any actual nuclear plant information. This work will benefit A/SMR vendors and utilities to implement security by design during the reactor design iteration phase, such that they do not have to perform upgrades and retrofits to the reactor after it is installed to improve its physical protection system. The tool may also be used to analyze domestic or foreign reactor designs to support the International Nuclear Security Techniques for Advanced Reactors (INSTAR) bilateral missions. Future works are planned to implement the methodology on a reference SFR reactor and a reference high-temperature gas-cooled reactor to obtain insights and lessons-learned for the A/SMR community.

97 MATHEMATICS AND COMPUTING↗

Transformational challenge reactor analysis to inform preconceptual core design decisions: Sensitivity study of transient analysis in a hydride-moderated microreactor

The Transformational Challenge Reactor (TCR) program aims to demonstrate a revolutionary design approach enabled by advanced manufacturing and data analytics in the design of nuclear reactors. This article discusses scoping analyses of preconceptual designs to inform TCR design decisions and the evaluation of sensitivities and uncertainties on postulated transient scenarios. The applicability of the systems codes TRACE and RELAP5-3D to TCR transient analysis are examined, and RELAP5-3D models are used to examine the transient response of two candidate core designs at multiple power levels. Then, the uncertainty quantification code RAVEN is used to quantify the effect of several design parameters on reactivity-initiated accident (RIA) progression at hot zero power (HZP) and hot full power (HFP) as well as to assess the impact of uncertainties in transient parameters for the pressurized loss of forced cooling (PLOFC).When results were compared, TRACE and RELAP5-3D showed good agreement in their ability to predict system behavior, but RELAP5-3D calculations were closer to analytical predictions for the RIA. Furthermore, models for a PLOFC accident in two designs (a UO2 and TRISO core) at multiple power levels showed greater temperature margins for the TRISO core at all power levels. Using this information, along with other scoping analyses and constraints, the TCR design team selected a power level of 3 MWth and a TRISO-based core design. For this design RAVEN was applied to vary RIA parameters in RELAP5-3D models at HZP and HFP to understand the effect on figures of merit such as peak power, fuel and coolant temperature, and energy deposition. This sensitivity study found that the inserted reactivity worth was the most important parameter controlling all figures of merit, but for insertion up to 1.5$ no failure of TCR fuel is anticipated. For constant reactivity insertions, the magnitude of the fuel temperature coefficient was found to have the greatest effect on all figures of merit under most circumstances. These results not only demonstrate the anticipated robust safety of the proposed TCR fuel form but also provide a reference for future metal-hydride moderated systems to understand RIA behavior. In the PLOFC, the impact of heat transfer enhancement due to wavy flow channel effects dominated the variance in peak temperatures, and variations in heat exchanger elevation provided the greatest control on natural circulation flow rate. Additionally, no fuel particle failure is anticipated in the PLOFC.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Migrational Instabilities in Particle Suspensions

This work deals with an instability arising from the shear-induced migration of particles in dense suspensions coupled with a dependence of viscosity on particle concentration. The analysis summarized here treats the inertialess (Re = O) linear stability of homogeneous simple shear flows for a Stokesian suspension model of the type proposed by Leighton and Acrivos (1987). Depending on the importance of shear-induced migration relative to concentration-driven diffusion, this model admits short-wave instability arising from wave-vector stretching by the base flow and evolving into particle-depleted shear bands. Moreover, this instability in the time-dependent problem corresponds to loss of ellipticity in the associated static problem (Re = O, Pe = O). While the isotropic version of the Leighton-Acrivos model is found to be stable with their experimentally determined parameters for simple shear, it is known that the stable model does not give a good quantitative description of particle clustering in the core of pipe flow (Nott and Brady 1994). This leads to the conjecture that an appropriate variant on the above model could explain such clustering as a two-phase bifurcation in the base flow.

Goddard, Joe D.↗

Core-Plasma Refilling and Erosion: Science Justification

Global magnetospheric models do not currently couple dynamic plasmaspheric plasma due to computational limitations that will be overcome by 2050. A true GGCM will need to include a dynamic plasmasphere for both physics and space weather prediction. Such inclusion requires that we answer longstanding, fundamental questions about the life cycle of core (<10 eV, plasmaspheric) plasma: how it gets into the inner magnetosphere, and how it is eroded. Erosion and refilling are processes controlling a fundamental space plasma population. During every geomagnetic disturbance, tens of metric tons of plasma are rapidly eroded away, then slowly and unevenly replenished. After several decades we still do not understand the cross-scale mechanisms proposed to be responsible for these 2 processes. Removal and replacement of this enormous plasma mass is as important to the dynamics of the magnetosphere as solar-wind driving. Our community must dedicate the resources and effort needed to solve this enduring puzzle. Other white papers submitted in response to the HP2050 call are discussing cold plasma and ionospheric outflow in more general contexts. This white paper deals specifically with the science justification for studying refilling and erosion.

Plasmasphere, core-plasma, erosion, refilling↗

When will low-contrast features be visible in a STEM X-ray spectrum image?

When will a small or low-contrast feature, such as an embedded second-phase particle, be visible in a scanning transmission electron microscopy (STEM) X-ray map? This work illustrates a computationally inexpensive method to simulate X-ray maps and spectrum images (SIs), based upon the equations of X-ray generation and detection. To particularize the general procedure, an example of nanostructured ferritic alloy (NFA) containing nm-sized Y 2 Ti 2 O 7 embedded precipitates in ferritic stainless steel matrix is chosen. The proposed model produces physically appearing simulated SI data sets, which can either be reduced to X-ray dot maps or analyzed via multivariate statistical analysis. Comparison to NFA X-ray maps acquired using three different STEM instruments match the generated simulations quite well, despite the large number of simplifying assumptions used. A figure of merit of electron dose multiplied by X-ray collection solid angle is proposed to compare feature detectability from one data set (simulated or experimental) to another. The proposed method can scope experiments that are feasible under specific analysis conditions on a given microscope. As a result, future applications, such as spallation proton–neutron irradiations, core-shell nanoparticles, or dopants in polycrystalline photovoltaic solar cells, are proposed.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Role of Amorphous Chains in Nanoplastic Formation from Semicrystalline Polymers

Semicrystalline polymers, e.g., polyethylene terephthalate (PET), release micro- (100 nm to 1 mm) and nanoplastics (up to 100 nm) [MNPL] when they are degraded under quiescent conditions. However, the exact molecular mechanisms leading to material fragmentation into MNPL are unknown. Here, we monitor the evolution of chain molecular weight and the MNPL production kinetics during hydrolysis of PET pellets. We find that only ~0.6% of the amorphous phase ester bonds are hydrolyzed at the onset of MNPL release. Then, by combining a random scission model with measured amorphous spacings, we estimate that only ~15% of the stress transmitters in the amorphous phase, namely, bridges and entangled loops, have failed by this point. Thus, spontaneous fragmentation of the semicrystalline nanostructure occurs despite significant intercrystalline connectivity. We resolve this apparent contradiction by proposing that fragmentation can only occur when progressive tie-chain scission causes the material to undergo the ductile/brittle transition. Mechanistically, we speculate that the internal stresses responsible for material fragmentation are caused by processing (i.e., residual stresses) and/or sample densification induced by chemi-crystallization. We discuss additional factors that may affect our analysis and thus require further investigation, such as skin-core effects, preferential degradation of stress transmitters and/or recrystallization processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

THREE-DIMENSIONAL BOLTZMANN HYDRO CODE FOR CORE COLLAPSE IN MASSIVE STARS. I. SPECIAL RELATIVISTIC TREATMENTS

We propose a novel numerical method for solving multi-dimensional, special relativistic Boltzmann equations for neutrinos coupled with hydrodynamics equations. This method is meant to be applied to simulations of core-collapse supernovae. We handle special relativity in a non-conventional way, taking account of all orders of v/c. Consistent treatment of the advection and collision terms in the Boltzmann equations has been a challenge, which we overcome by employing two different energy grids: Lagrangian remapped and laboratory fixed grids. We conduct a series of basic tests and perform a one-dimensional simulation of core-collapse, bounce, and shock-stall for a 15 M {sub ☉} progenitor model with a minimum but essential set of microphysics. We demonstrate in the latter simulation that our new code is capable of handling all phases in core-collapse supernova. For comparison, a non-relativistic simulation is also conducted with the same code, and we show that they produce qualitatively wrong results in neutrino transfer. Finally, we discuss a possible incorporation of general relativistic effects into our method.

79 ASTRONOMY AND ASTROPHYSICS↗

Attribute-Aware RBFs: Interactive Visualization of Time Series Particle Volumes Using RT Core Range Queries

Smoothed-particle hydrodynamics (SPH) is a mesh-free method used to simulate volumetric media in fluids, astrophysics, and solid mechanics. Visualizing these simulations is problematic because these datasets often contain millions, if not billions of particles carrying physical attributes and moving over time. Radial basis functions (RBFs) are used to model particles, and overlapping particles are interpolated to reconstruct a high-quality volumetric field; however, this interpolation process is expensive and makes interactive visualization difficult. Existing RBF interpolation schemes do not account for color-mapped attributes and are instead constrained to visualizing just the density field. To address these challenges, we exploit ray tracing cores in modern GPU architectures to accelerate scalar field reconstruction. We use a novel RBF interpolation scheme to integrate per-particle colors and densities, and leverage GPU-parallel tree construction and refitting to quickly update the tree as the simulation animates over time or when the user manipulates particle radii. We also propose a Hilbert reordering scheme to cluster particles together at the leaves of the tree to reduce tree memory consumption. Finally, we reduce the noise of volumetric shadows by adopting a spatially temporal blue noise sampling scheme. Our method can provide a more detailed and interactive view of these large, volumetric, time-series particle datasets than traditional methods, leading to new insights into these physics simulations.

Particle Volumes↗

Hydrogen-poor Superluminous Supernovae with Late-time H α Emission: Three Events From the Intermediate Palomar Transient Factory

We present observations of two new hydrogen-poor superluminous supernovae (SLSN-I), iPTF15esb and iPTF16bad, showing late-time Hα emission with line luminosities of $(1\mbox{-}3)\times {10}^{41}$ erg s -1 and velocity widths of (4000-6000) km s - . Including the previously published iPTF13ehe, this makes up a total of three such events to date. iPTF13ehe is one of the most luminous and the slowest evolving SLSNe-I, whereas the other two are less luminous and fast decliners. We interpret this as a result of the ejecta running into a neutral H-shell located at a radius of ~10 16 cm. This implies that violent mass loss must have occurred several decades before the supernova explosion. Such a short time interval suggests that eruptive mass loss could be common shortly before core collapse, and more importantly helium is unlikely to be completely stripped off the progenitor and could be present in the ejecta. It is a mystery why helium features are not detected, even though nonthermal energy sources, capable of ionizing He, may exist as suggested by the O ii absorption series in the early-time spectra. Our late-time spectra (+240 days) appear to have intrinsically lower [O i] 6300 Å luminosities than that of SN2015bn and SN2007bi, which is possibly an indication of less oxygen (<10 M ⊙ ). The blueshifted Hα emission relative to the hosts for all three events may be in tension with the binary model proposed for iPTF13ehe. Lastly, iPTF15esb has a peculiar light curve (LC) with three peaks separated from one another by ~22 days. The LC undulation is stronger in bluer bands. One possible explanation is ejecta-circumstellar medium interaction.

79 ASTRONOMY AND ASTROPHYSICS↗

The internal structure of the convective envelope of the sun

An axisymmetric hydrodynamic model of the sun consisting of a core at radiant equilibrium and a convective envelope is studied with consideration of gravitational forces, rotation, turbulence effects, and electromagnetic forces. The dimensionless parameters that control the equations of plasma motion within the sun and the boundary conditions are investigated; estimates are given for these parameters on the basis of experimental data and the existing theoretical models. Stability conditions are derived for convective equilibrium, developed cellular convection, and developed turbulent convection with consideration of rotation and turbulent Reynolds stresses. These conditions characterize the possible physical state of the matter in the convective envelope. Using estimates of the dimensionless parameters and models of some of them, a method of successive approximations is proposed for solution of the problem as a whole.

Yavorskaya, I. M.↗

Achieving High Efficiency in Reduced Order Modeling for Large Scale Polycrystal Plasticity Simulations

Reduced order models for the nonlinear response of heterogeneous microstructures typically require a construction (or training) stage to build the reduced order basis. In this manuscript, an efficient model construction strategy for the eigenstrain homogenization method (EHM) is presented. The proposed strategy relies on a parallel, element-by-element, conjugate gradient solver. Near linear scaling has been achieved with respect to the number of degrees of freedom used to resolve the microstructure. Linear scaling with respect to the number of pre-analyses required to construct the reduced order model (ROM) follows from the EHM formulation. Furthermore, a parallel implementation for fast evaluation of the constructed ROM has been developed using shared memory parallelization. It has been shown that for large microstructures with ≈ 10,000 grains, the total computational cost of evaluating the nonlinear response of a polycrystal could be reduced by approximately an order of magnitude using 32 cores with respect to serial ROM simulation. The present methodology has been verified using an additively manufactured polycrystalline microstructure of a nickel-based superalloy, Inconel 625. The capability of the developed framework to construct a ROM for such large microstructures, as well as the ability of the ROM to predict average and local quantities of interest has been demonstrated.

microscale↗

Multiphysics Simulation of Moisture-Graphite Oxidation in MHTGR

A small amount of moisture on the ppm level could present in the primary helium coolant under normal operations of High-Temperature Gas-cooled Reactors (HTGRs). One potential safety issue related to HTGRs is the chronic moisture-graphite oxidation, which could affect the integrity of the fuel blocks and reflectors in the reactor core, and the support columns in the hot plenum. However, it is infeasible to perform chronic moisture-graphite oxidation tests under the prototypic high-pressure, high-temperature conditions over a long period of time that is comparable to the full service time, i.e., 36 months for Modular High Temperature Gas-cooled Reactor (MHTGR). As an alternate, accelerated moisture-graphite oxidation tests at the atmospheric pressure and high moisture concentrations have been performed in the literature. Through these accelerated tests, global reaction rate equations have been proposed, which makes it possible to study the moisture-graphite oxidation numerically. This paper is aimed to establish a multiphysics model that can evaluate the moisture-graphite oxidation under MHTGR normal operation conditions. COMSOL Multiphysics was applied to couple the fluid flow, heat and mass transfer, chemical reaction and material structural changes. To reduce the calculation time, the prototypic three-dimensional (3-D) structures were simplified into a two-dimensional (2-D) simulation domain. The performance of four nuclear grades of graphite (IG-110, 2114, PCEA, and NBG-17) over a 36-month service period was then investigated under the most likely prototypic conditions of MHTGR. The simulation results indicate that most oxidation occurs in the three or four bottom fuel blocks due to their higher temperatures. In general, however, only a thin layer of 1 mm into the graphite will be considerably oxidized even for those three or four bottom fuel blocks.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Inertia-driven particle migration and mixing in a wall-bounded laminar suspension flow

Laminar pressure-driven suspension flows are studied in the situation of neutrally buoyant particles at finite Reynolds number. The numerical method is validated for homogeneous particle distribution (no lateral migration across the channel): the increase of particle slip velocities and particle stress with inertia and concentration is in agreement with former works in the literature. In the case of a two-phase channel flow with freely moving particles, migration towards the channel walls due to the Segré-Silberberg effect is observed, leading to the development of a non-uniform concentration profile in the wall-normal direction (the concentration peaks in the wall region and tends towards zero in the channel core). The particle accumulation in the region of highest shear favors the shear-induced particle interactions and agitation, the profile of which appears to be correlated to the concentration profile. A 1D model predicting particle agitation, based on the kinetic theory of granular flows in the quenched state regime when Stokes number St = O(1) and from numerical simulations when St < 1, fails to reproduce the agitation profile in the wall normal direction. Instead, the existence of secondary flows is clearly evidenced by long time simulations. These are composed of a succession of contra-rotating structures, correlated with the development of concentration waves in the transverse direction. The mechanism proposed to explain the onset of this transverse instability is based on the development of a lift force induced by spanwise gradient of the axial velocity fluctuations. The establishment of the concentration profile in the wall-normal direction therefore results from the combination of the mean flow Segré-Silberberg induced migration, which tends to stratify the suspension and secondary flows which tend to mix the particles over the channel cross section.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Power System Resilience Evaluation Framework and Metric Review

Power system resilience has been an emerging hot topic in recent years to investigate the increasing threats of extreme events, such as natural disasters, severe weather, and cyberattacks. Although much research has been done to define, model, and quantify resilience from different aspects, the lack of universally accepted evaluation methods and resilience metrics makes it difficult to assess and compare resilience across different power systems, such as what is typically done in power system reliability studies. In this paper, first, we review the definitions of resilience, and we summarize two core concepts shared by most of the literature. Then, we develop a new framework to assess power system resilience from two perspectives - i.e., pre-event estimation and post-event evaluation - to capture system resilience performance in both general and specific fashions. We conduct a thorough review of existing resilience metrics and categorize them using the proposed framework, where recommendations are also proposed to capture core concepts of resilience.

power system resilience↗

Model-based reconstruction for enhanced x-ray CT of dense tri-structural isotropic particles

Tri-structural isotropic (TRISO) fuel particles are a key component of next generation nuclear fuels. Using x-ray computed tomography (CT) to characterize TRISO particles is challenging because of the strong attenuation of the x-ray beam by the uranium core, leading to severe photon starvation in a substantial fraction of the measurements. Furthermore, the overall acquisition time for a high-resolution CT scan can be very long when using conventional laboratory-based x-ray systems and reconstruction algorithms. Specifically, when analytic methods such as the Feldkamp–Davis–Kress (FDK) algorithm are used for reconstruction, it results in severe streak artifacts and noise in the corresponding 3D volume, which makes subsequent analysis of the particles challenging. In this paper, we develop and apply model-based image reconstruction (MBIR) algorithms to improve the quality of CT reconstructions for TRISO particles to facilitate better characterization. We demonstrate that the proposed MBIR algorithms can significantly suppress artifacts with minimal pre-processing compared to conventional approaches. We also demonstrate that the proposed MBIR approach can obtain high-quality reconstruction compared to the FDK approach even when using a fraction of the typically acquired measurements, thereby enabling dramatically faster measurement times for TRISO particles.

Venkatakrishnan, Singanallur↗