Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “radiation detection algorithm”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Free spectator nucleons in ultracentral relativistic heavy-ion collisions as a probe of neutron skin

Besides the yield ratio of free spectator neutrons produced in ultracentral 96 Zr + 96 Zr to 96 Ru + 96 Ru collisions, we propose that the yield ratio N n /N p of free spectator neutrons to protons in a single collision system at the BNL Relativistic Heavy Ion Collider or the CERN Large Hadron Collider can be a more sensitive probe of the neutron-skin thickness Δrnp and the slope parameter L of the symmetry energy. Here, the idea is demonstrated based on the proton and neutron density distributions of colliding nuclei obtained from Skyrme-Hartree-Fock-Bogolyubov calculations, and on a Glauber model that provides information of spectator matter. The final spectator particles are produced from direct emission, clusterization by a minimum spanning tree algorithm or a Wigner function approach, and deexcitation of heavy clusters by gemini. A larger Δr np associated with a larger L value increases the isospin asymmetry of spectator matter and thus leads to a larger N n /N p , especially in ultracentral collisions where the multiplicity of free nucleons are free from the uncertainties of cluster formation and deexcitation. We have further shown that the double ratio of N n /N p in isobaric collision systems or in collisions by isotopes helps to cancel the detecting efficiency for protons. Effects from nuclear deformation and electromagnetic excitation are studied, and they are found to be subdominant compared to the expected sensitivity to Δr np .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Can the GPM IMERG Final Product Accurately Represent MCSs’ Precipitation Characteristics over the Central and Eastern United States?

Mesoscale convective systems (MCSs) play an important role in water and energy cycles as they produce heavy rainfall and modify the radiative profile in the tropics and midlatitudes. An accurate representation of MCSs’ rainfall is therefore crucial in understanding their impact on the climate system. The V06B Integrated Multisatellite Retrievals from Global Precipitation Measurement (IMERG) half-hourly precipitation final product is a useful tool to study the precipitation characteristics of MCSs because of its global coverage and fine spatiotemporal resolutions. However, errors and uncertainties in IMERG should be quantified before applying it to hydrology and climate applications. This study evaluates IMERG performance on capturing and detecting MCSs’ precipitation in the central and eastern United States during a 3-yr study period against the radar-based Stage IV product. The tracked MCSs are divided into four seasons and are analyzed separately for both datasets. IMERG shows a wet bias in total precipitation but a dry bias in hourly mean precipitation during all seasons due to the false classification of nonprecipitating pixels as precipitating. These false alarm events are possibly caused by evaporation under the cloud base or the misrepresentation of MCS cold anvil regions as precipitating clouds by the algorithm. IMERG agrees reasonably well with Stage IV in terms of the seasonal spatial distribution and diurnal cycle of MCSs precipitation. A relative humidity (RH)-based correction has been applied to the IMERG precipitation product, which helps reduce the number of false alarm pixels and improves the overall performance of IMERG with respect to Stage IV.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Marine Boundary Layer Cloud Boundaries and Phase Estimation Using Airborne Radar and In Situ Measurements During the SOCRATES Campaign over Southern Ocean

The Southern Ocean Clouds, Radiation, Aerosol Transport Experimental Study (SOCRATES) was an aircraft-based campaign (15 January–26 February 2018) that deployed in situ probes and remote sensors to investigate low-level clouds over the Southern Ocean (SO). A novel methodology was developed to identify cloud boundaries and classify cloud phases in single-layer, low-level marine boundary layer (MBL) clouds below 3 km using the HIAPER Cloud Radar (HCR) and in situ measurements. The cloud base and top heights derived from HCR reflectivity, Doppler velocity, and spectrum width measurements agreed well with corresponding lidar-based and in situ estimates of cloud boundaries, with mean differences below 100 m. A liquid water content–reflectivity (LWC-Z) relationship, LWC = 0.70Z0.29, was derived to retrieve the LWC and liquid water path (LWP) from HCR profiles. The cloud phase was classified using HCR measurements, temperature, and LWP, yielding 40.6% liquid, 18.3% mixed-phase, and 5.1% ice samples, along with drizzle (29.1%), rain (3.2%), and snow (3.7%) for drizzling cloud cases. The classification algorithm demonstrates good consistency with established methods. This study provides a framework for the boundary and phase detection of MBL clouds, offering insights into SO cloud microphysics and supporting future efforts in satellite retrievals and climate model evaluation.

MBL clouds over Southern Ocean↗
DOE OSTI2024

Minijet clustering algorithm using transverse-momentum seeds in high-energy nuclear collisions

We propose an algorithm to detect mini-jet clusters in high-energy nuclear collisions, by selecting a high-transverse-momentum (pT) particle as a seed and assigning a clustering radius (R) in the pseudorapidity and azimuthal-angle space. Our PYTHIA simulations for p+p collisions show that a scheme with a seeding p T of around 0.5 GeV/c and R of approximately 0.6 satisfactorily identifies mini-jet clusters. The correlation between clusters obtained in PYTHIA calculations using the algorithm exhibits the proper behavior of hard-scattering-like processes, suggesting its usefulness in isolating mini-jet-like clusters from non-hard-scattering soft processes when applied to actual nuclear-collision data, thereby allowing a closer examination of both the mini-jet and the soft mechanisms.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Neutrino event selection in the MicroBooNE liquid argon time projection chamber using Wire-Cell 3D imaging, clustering, and charge-light matching

An accurate and efficient event reconstruction is required to realize the full scientific capability of liquid argon time projection chambers (LArTPCs). The current and future neutrino experiments that rely on massive LArTPCs create a need for new ideas and reconstruction approaches. Wire-Cell, proposed in recent years, is a novel tomographic event reconstruction method for LArTPCs. The Wire-Cell 3D imaging approach capitalizes on charge, sparsity, time, and geometry information to reconstruct a topology-agnostic 3D image of the ionization electrons prior to pattern recognition. A second novel method, the many-to-many charge-light matching, then pairs the TPC charge activity to the detected scintillation light signal, thus enabling a powerful rejection of cosmic-ray muons in the MicroBooNE detector. A robust processing of the scintillation light signal and an appropriate clustering of the reconstructed 3D image are fundamental to this technique. In this paper, we describe the principles and algorithms of these techniques and their successful application in the MicroBooNE experiment. A quantitative evaluation of the performance of these techniques is presented. Using these techniques, a 95% efficient pre-selection of neutrino charged-current events is achieved with a 30-fold reduction of non-beam-coincident cosmic-ray muons, and about 80% of the selected neutrino charged-current events are reconstructed with at least 70% completeness and 80% purity.

3D imaging↗
DOE OSTI2022

Synergistic Retrieval of Temperature and Humidity Profiles from Space-Based and Ground-Based Infrared Sounders Using an Optimal Estimation Method

The atmospheric temperature and humidity profiles of the troposphere are generally measured by radiosondes and satellites, which are essential for analyzing and predicting weather. Nevertheless, the insufficient observation frequencies and low detection accuracy of the boundary layer restricts the description of atmospheric state changes by the temperature and humidity profiles. Therefore, this work focus on retrieving the temperature and humidity profiles using observations of the FengYun-4 (FY-4) Geostationary Interferometric Infrared Sounder (GIIRS) combined with ground-based infrared spectral observations from the Atmospheric Emitted Radiance Interferometer (AERI), which are more accurate than space-based individual retrieval results and have a wider effective retrieval range than ground-based individual retrieval results. Based on the synergistic observations, which are made by matching the space-based and ground-based data with those of different spatial and temporal resolutions, a synergistic retrieval process is proposed to obtain the temperature and humidity profiles at a high frequency under clear-sky conditions based on the optimal estimation method. In this research, using the line-by-line radiative transfer model (LBLRTM) as the forward model for observing simulations, a retrieval experiment was carried out in Qingdao, China, where an AERI is situated. Taking radiosonde data as a reference for comparing the retrieval results of the temperature and humidity profiles of the troposphere, the root-mean-square error (RMSE) of the synergistic retrieval algorithm below 400 hPa is within 2 K for temperature and within 12% for relative humidity. Compared with the GIIRS individual retrieval, the RMSE of temperature and relative humidity for the synergistic method is reduced by 0.13–1.5 K and 2.7–4.4% at 500 hPa, and 0.13–2.1 K and 2.5–7.2% at 900 hPa. Moreover, the forecast index (FI) calculated from the retrieval results shows reasonable consistency with the FIs calculated from the ERA5 reanalysis and from radiosonde data. The synergistic retrieval results have higher temporal resolution than space-based retrieval results and can reflect the changes in the atmospheric state more accurately. Overall, the results demonstrated the promising potential of the synergistic retrieval of temperature and humidity profiles at high accuracy and high temporal resolution under clear-sky conditions from FY-4/GIIRS and AERI.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Detection of Isotopes in Urban Source Search Low-Count Gamma Spectra Using Hopfield Neural Networks

Source search campaigns involve measurements of background gamma-ray spectra with a mobile detector-spectrometer traveling along arbitrarily chosen trajectories over a wide screening area. Radiation counts are typically measured with a tellurium-doped sodium iodide [NaI(Tl)] scintillator detector-spectrometer in short acquisition intervals, usually 1 s. The objective is to detect orphan isotopes with half-lives shorter than those of the isotopes in the natural background. In principle, radioisotopes can be identified by their unique gamma emission spectrum. However, detecting orphan isotopes in search data is challenging because low counts measured in short acquisition intervals result in incomplete spectral lines. In this study, we investigate the performance of a Hopfield neural network (HNN) that implements an auto-associative memory for the detection of isotopes of interest in an urban search campaign. The HNN is trained on one example of gamma spectra with well-resolved spectral lines of each isotope of interest. During testing, the auto-associative memory implementation of the HNN processes low-count gamma spectra with partially complete isotopic lines by matching incoming measurements to the closest one of its memory-stored patterns. The testing database consisted of almost 10 000 1-s gamma spectra, including measurements of orphan isotopes 137 Cs, 241 Am, and 131 I, obtained during two urban search surveys with a NaI(Tl) detector. The performance of the HNN detection algorithm was evaluated using precision, recall, and F1 scores, and benchmarked with a multiple linear regression (MLR) identification algorithm. In conclusion, the test results demonstrate that HNN outperforms MLR in the detection of all the isotopes of interest.

Auto associative memory↗
DOE OSTI2021

Holes in Optical Lightning Flashes: Identifying Poorly Transmissive Clouds in Lightning Imager Data

Space-based optical lightning sensors including the lightning imaging sensor (LIS) and geostationary lightning mapper (GLM) are pixelated imagers that detect lightning as transient increases in cloud top illumination. Detection requires optical emissions to escape the cloud top to space with sufficient energy to trigger a pixel on the imaging array. Through scattering and absorption, certain clouds are able to block most light from reaching the instrument, causing a reduction in detection efficiency (DE) and possibly location accuracy (LA). Radiant lightning emissions that illuminate large cloud top areas are used to examine scenarios where clouds block light from reaching orbit. In some cases, these anomalies in the spatial radiance distribution from the lightning pulse lead to “holes” in the optical lightning flash where certain pixels fail to trigger. Such holes are identified algorithmically in the Tropical Rainfall Measuring Mission satellite LIS record and the microphysical properties of the coincident storm region are queried. We find that holes primarily occur in tall (IR T b < 235 K) convection (87%) and overhanging anvil clouds (10%). The remaining 3% of holes occur in moderate-to-weak convection or in clear air breaks between stormclouds. We further demonstrate how an algorithm that assesses the spatial radiance patterns from energetic lightning pulses might be used to construct an optical transmission gridded stoplight product for GLM that could help operators identify clouds with a potentially reduced DE and LA.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Real-Time, Adaptive Radiological Anomaly Detection and Isotope Identification Using Non-Negative Matrix Factorization

Spectroscopic anomaly detection and isotope identification algorithms are integral components in nuclear nonproliferation applications such as search operations. The task is especially challenging in the case of mobile detector systems because the observed gamma-ray background changes more than for a static detector system, and a pretrained background model can easily find itself out of domain. The result is that algorithms may exceed their intended false alarm rate or sacrifice detection sensitivity to maintain the desired false alarm rate. Non-negative matrix factorization (NMF) is a powerful tool for spectral anomaly detection and identification, but, like many similar algorithms that rely on data-driven background models, in its conventional implementation, it is unable to update in real time to account for environmental changes that affect the background spectroscopic signature. Here, we have developed a novel NMF-based algorithm that periodically updates its background model to accommodate changing environmental conditions. The adaptive NMF algorithm involves fewer assumptions about its environment, making it more generalizable than existing NMF-based methods while maintaining or exceeding detection performance on simulated and real-world datasets.

Anomaly detection↗
DOE OSTI2020

Radiation Source Localization Using Surrogate Models Constructed from 3-D Monte Carlo Transport Physics Simulations

Recent research has focused on the development of surrogate models for radiation source localization in a simulated urban domain. We employ the Monte Carlo N-Particle (MCNP) code to provide high- delity simulations of radiation transport within an urban domain. The model is constructed to employ a source location (x, y, z) as input and return the estimated count rate for a set of speci ed detector locations. Because MCNP simulations are computationally expensive, we develop e cient and accurate surrogate models of the detector responses. We construct surrogate models using Gaussian processes (GP) and neural networks (NN) that we train and verify using the MCNP simulations. The trained surrogate models provide an e cient framework for Bayesian inference and experimental design. We employ Delayed Rejection Adaptive Metropolis (DRAM), a Markov Chain Monte Carlo (MCMC) algorithm, to infer the location and intensity of an unknown source. The DRAM results yield a posterior probability distribution for the source's location conditioned on the observed detector count rates. The posterior distribution exhibits regions of high and low probability within the simulated environment identifying potential source locations. In this manner, we can quantify the source location to within at least one of these regions of high probability in the considered cases. Employing these methods, we are able to reduce the space of potential source locations by at least 60%.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Synthesis of ARM User Facility Surface Rainfall Datasets to Construct a Best Estimate Value Added Product (PrecipBE)

Surface precipitation measurements are essential for Earth system model (ESM) evaluation and understanding cloud processes. An ever-growing need for robust, temporally evolving, and easy-to-use statistical datasets provides motivation for a baseline ground-based precipitation properties data product. The U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility operates an extensive suite of precipitation instruments with various sensitivities and operating mechanisms, which render the decision of which instrument to use based on one or more fixed thresholds challenging and prone to errors and bias. Using a long-term instrument inter-comparison from a unique per-precipitation event perspective, rather than instantaneous sample comparison, we demonstrate that ARM rainfall-measuring instruments are generally consistent with each other at the statistical level. Inter-instrument deviations at the single event level can be large, especially for specific rainfall event properties such as maximum precipitation rates. A machine-learning (ML) analysis using a random forest regressor indicates that in some cases, depending on instrument, local site climatology, and/or specific deployment configuration, certain atmospheric state variables influence the measured quantities in an unpredictable manner. Thus, a-priori weighting of different instruments does not necessarily lead to more accurate and less biased synthesis of instrument data. These results motivate the design of the ARM precipitation best-estimate (PrecipBE) value-added product, which incorporates all valid precipitation data while considering data quality and other instrument limitations. PrecipBE consists of time series and tabular statistics datasets in an easy-to-use and insightful per-precipitation event format. It provides a large set of precipitation event properties supplemented with ancillary data from ARM datasets that correspond to the detected precipitation events. We describe the PrecipBE algorithm and demonstrate its use via the examination of a single-day output as well as a long-term trend analysis of precipitation events at the ARM Southern Great Plains (SGP) site, covering more than 30 years of data. The trend analysis tentatively suggests a long-term temporal tendency for mainly shorter and less intense precipitation events at the SGP site, but a long-term increase in annual rainfall by more than 36 mm (5 %) per decade. This rainfall trend is catalyzed primarily by more extreme event properties of relatively rare, intense precipitation events, with event total and 1 min maximum precipitation rate at a 1 year timeframe increasing up to 5 mm and 9 mm h −1 (several percent) per decade, respectively. While the currently available PrecipBE datasets (at https://adc.arm.gov/discovery/, last access: 8 December 2025) cover rainfall from multiple ARM deployments up to March 2025, PrecipBE is planned to be expanded to include solid-phase precipitation and will soon become an operational product with a several-day lag from real-time. We invite the ARM user community to leverage this new product and welcome user feedback to further enhance the dataset.

Silber, Israel [Pacific Northwest National Laborat↗
DOE OSTI2024

Deployment of Water-based Liquid Scintillator in the Accelerator Neutrino Neutron Interaction Experiment

The Accelerator Neutrino Neutron Interaction Experiment (ANNIE) is a 26-ton water Cherenkov neutrino detector installed on the Booster Neutrino Beam (BNB) at Fermilab. Its main physics goals are to perform a measurement of the neutron yield from neutrino-nucleus interactions, as well as a measurement of the charged-current cross section of muon neutrinos. An equally important focus is the research and development of new detector technologies and target media. Specifically, water-based liquid scintillator (WbLS) is of interest as a novel detector medium, as it allows for the simultaneous detection of Cherenkov light and scintillation. This paper presents the deployment of a 366 L WbLS vessel in ANNIE in March 2023 and the subsequent detection of both Cherenkov light and scintillation from the WbLS. Finally, this proof-of-concept allows for the future development of reconstruction and particle identification algorithms in ANNIE, as well as dedicated analyses within the WbLS volume, such as the search for neutral-current events and the hadronic scintillation component.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

High-Energy Collision of Quarks and Mesons in the Schwinger Model: From Tensor Networks to Circuit QED

With the aim of studying nonperturbative out-of-equilibrium dynamics of high-energy particle collisions on quantum simulators, we investigate the scattering dynamics of lattice quantum electrodynamics in 1 + 1 dimensions. Working in the bosonized formulation of the model and in the thermodynamic limit, we use uniform-matrix-product-state tensor networks to construct multiparticle wave-packet states, evolve them in time, and detect outgoing particles post collision. This facilitates the numerical simulation of scattering experiments in both confined and deconfined regimes of the model at different energies, giving rise to rich phenomenology, including inelastic production of quark and meson states, meson disintegration, and dynamical string formation and breaking. We obtain elastic and inelastic scattering cross sections, together with time-resolved momentum and position distributions of the outgoing particles. Furthermore, we propose an analog circuit-QED implementation of the scattering process that is native to the platform, requires minimal ingredients and approximations, and enables practical schemes for particle wave-packet preparation and evolution. Furthermore, this study highlights the role of classical and quantum simulation in enhancing our understanding of scattering processes in quantum field theories in real time.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

The large-area hybrid-optics CLAS12 RICH: First years of data-taking

The CLAS12 deep-inelastic scattering experiment at the upgraded 12 GeV continuous electron beam accelerator facility of Jefferson Lab conjugates luminosity and wide acceptance to study the 3D nucleon structure in the yet poorly explored valence region, and to perform precision measurements in hadron spectroscopy. A large area ring-imaging Cherenkov detector has been designed to achieve the required hadron identification in the momentum range from 3 GeV/c to 8 GeV/c, with the kaon rate about one order of magnitude lower than the rate of pions and protons. The adopted solution comprises aerogel radiator and composite mirrors in a novel hybrid optics design, where either direct or reflected light could be imaged in a high-packed and high-segmented photon detector. The first RICH module was assembled during the second half of 2017 and installed at the beginning of January 2018, in time for the start of the experiment. The second RICH module, planned with the goal to be ready for the beginning of the operation with polarized targets, has been timely built despite the complications caused by the pandemic crisis and successfully installed in June 2022. Finally, the detector performance is here discussed with emphasis on the operation and stability during the data-taking, calibration and alignment procedures, reconstruction and pattern recognition algorithms, and particle identification.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Supernova electron-neutrino interactions with xenon in the nEXO detector

Electron-neutrino charged-current interactions with xenon nuclei were modeled in the nEXO neutrinoless double-𝛽 decay detector (∼5 metric ton, 90% 136 Xe, 10% 134 Xe) to evaluate its sensitivity to supernova neutrinos. Predictions for event rates and detectable signatures were modeled using the Model of Argon Reaction Low Energy Yields (MARLEY) event generator. We find good agreement between MARLEY’s predictions and existing theoretical calculations of the inclusive cross sections at supernova neutrino energies. The interactions modeled by MARLEY were simulated within the nEXO simulation framework and were run through an example reconstruction algorithm to determine the detector’s efficiency for reconstructing these events. The simulated data, incorporating the detector response, were used to study the ability of nEXO to reconstruct the incident electron-neutrino spectrum and these results were extended to a larger xenon detector of the same isotope enrichment. We estimate that nEXO will be able to observe electron-neutrino interactions with xenon from supernovae as far as 5–8 kpc from Earth, while the ability to reconstruct incident electron-neutrino spectrum parameters from observed interactions in nEXO is limited to closer supernovae.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

The picasso map-making code: application to a simulation of the QUIJOTE northern sky survey

ABSTRACT Map-making is an important step for the data analysis of cosmic microwave background (CMB) experiments. It consists of converting the data, which are typically a long, complex, and noisy collection of measurements, into a map, which is an image of the observed sky. We present in this paper a new map-making code named picasso (Polarization and Intensity CArtographer for Scanned Sky Observations), which was implemented to construct intensity and polarization maps from the Multi Frequency Instrument (MFI) of the QUIJOTE (Q-U-I Joint TEnerife) CMB polarization experiment. picasso is based on the destriping algorithm, and is suited to address specific issues of ground-based microwave observations, with a technique that allows the fit of a template function in the time domain, during the map-making step. This paper describes the picasso code, validating it with simulations and assessing its performance. For this purpose, we produced realistic simulations of the QUIJOTE-MFI survey of the northern sky (approximately ∼20 000 deg2), and analysed the reconstructed maps with picasso, using real and harmonic space statistics. We show that, for this sky area, picasso is able to reconstruct, with high fidelity, the injected signal, recovering all the scales with ℓ > 10 in TT, EE, and BB. The signal error is better than 0.001 per cent at 20 < ℓ < 200. Finally, we validated some of the methods that will be applied to the real wide-survey data, like the detection of the CMB anisotropies via cross-correlation analyses. Despite that the implementation of picasso is specific for QUIJOTE-MFI data, it could be adapted to other experiments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Lipkin model on a quantum computer

Atomic nuclei are important laboratories for exploring and testing new insights into the universe, such as experiments to directly detect dark matter or explore properties of neutrinos. The targets of interest are often heavy, complex nuclei that challenge our ability to reliably model them (as well as quantify the uncertainty of those models) with classical computers. Hence there is great interest in applying quantum computation to nuclear structure for these applications. As an early step in this direction, especially with regards to the uncertainties in the relevant quantum calculations, we develop circuits to implement variational quantum eigensolver (VQE) algorithms for the Lipkin-Meshkov-Glick model, which is often used in the nuclear physics community as a testbed for many-body methods. Here, we present quantum circuits for VQE for two and three particles and discuss the construction of circuits for more particles. Implementing the VQE for a two-particle system on the IBM Quantum Experience, we identify initialization and two-qubit gates as the largest sources of error. We find that error mitigation procedures reduce the errors in the results significantly, but additional quantum hardware improvements are needed for quantum calculations to be sufficiently accurate to be competitive with the best current classical methods.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Inverse Method for Noninvasive Material Discrimination with Multienergy X-Ray Radiography

Active x-ray inspections have great potential for noninvasive detection of illicit materials. Dual-energy x-ray radiography, in particular, has been used to determine material composition by utilizing the energy-dependence of the x-ray attenuation coefficients. However, current implementations of this method are limited in their ability to determine the material composition of composite objects. Here, in this paper, we present an inverse algorithm that uses multienergy x-ray radiography data to noninvasively reconstruct the contents of a container. A critical feature of the current contribution is that material identification can be performed using conventional detectors and x-ray spectra with different endpoint energies produced by varying nothing more than the tube voltage. Adaptive regularization is used to increase the accuracy of material estimations from multienergy data sets. The utility of these methods is demonstrated with experimentally acquired radiographs obtained using a tunable x-ray source that produces spectra with endpoint energies of 100–450 keV. The object inspected is a scale model of a nuclear materials storage container composed of three-dimensional printed plastic and stainless-steel spheres inside a thin-walled steel container. Reconstructions of the steel sphere thicknesses are within a root-mean-square error of 0.37 cm.

47 OTHER INSTRUMENTATION↗