Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “tomography problems”

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 235 records · Page 13

Fully Arbitrary Control of Frequency-Bin Qubits

Accurate control of two-level systems is a longstanding problem in quantum mechanics. One such quantum system is the frequency-bin qubit: a single photon existing in superposition of two discrete frequency modes. In this Letter, we demonstrate fully arbitrary control of frequency-bin qubits in a quantum frequency processor for the first time. Furthermore, we numerically establish optimal settings for multiple configurations of electro-optic phase modulators and pulse shapers, experimentally confirming near-unity mode-transformation fidelity for all fundamental rotations. Performance at the single-photon level is validated through the rotation of a single frequency-bin qubit to 41 points spread over the entire Bloch sphere, as well as tracking of the state path followed by the output of a tunable frequency beam splitter, with Bayesian tomography confirming state fidelities Fρ>0.98 for all cases. Overall, such high-fidelity transformations expand the practical potential of frequency encoding in quantum communications, offering exceptional precision and low noise in general qubit manipulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Nondestructive Evaluation of Advanced Materials with X-ray Phase Mapping

X-ray radiation has been widely used for imaging applications since Rontgen first discovered X-rays over a century ago. Its large penetration depth makes it ideal for the nondestructive visualization of the internal structure and/or defects of materials unobtainable otherwise. Currently used nondestructive evaluation (NDE) tools, X-ray radiography and tomography, are absorption-based, and work well in heavy-element materials where density or composition variations due to internal structure or defects are high enough to produce appreciable absorption contrast. However, in many cases where materials are light-weight and/or composites that have similar mass absorption coefficients, the conventional absorption-based X-ray methods for NDE become less useful. Indeed, the light-weight and ultra-high-strength requirements for the most advanced materials used or developed for current flight mission and future space exploration pose a great challenge to the standard NDE tools in that the absorption contrast arising from the internal structure of these materials is often too weak to be resolved. In this presentation, a solution to the problem, the use of phase information of X-rays for phase contrast X-ray imaging, will be discussed, along with a comparison between the absorption-based and phase-contrast imaging methods. Latest results on phase contrast X-ray imaging of lightweight Space Shuttle foam in 2D and 3D will be presented, demonstrating new opportunities to solve the challenging issues encountered in advanced materials development and processing.

Hu, Zhengwei↗

An improved specimen preparation method for marine shallow gas-bearing sand sediments and its validations

Instabilities of shallow gas-charged seabed are potential geological hazards in ocean engineering. In practice, the conventional field sampling techniques failed to obtain undisturbed gas-bearing sediments from the seabed for laboratory mechanical testing because of sensitive gas exsolution and escape from sediments. However, preparation of representative remoulded gas-charged specimens is a challenging issue, because it is rather difficult to quantitatively control the gas content and obtain uniform distribution of gas bubbles within the specimen. Given the above problems, this work proposes a reliable approach to reconstitute the high-saturation specimen of gas-charged sediments in the laboratory by an improved multifunction integrated triaxial apparatus (MITA). This apparatus is developed based on an advanced stress path triaxial system by introducing a temperature-controlled system and a wave-monitoring system. The temperature-controlled system is used to accurately mimic the in situ environments of sediments in the seabed. The wave-monitoring system is utilized to identify exsolution point of free gas and examine the disturbance of gas to specimens during gas exsolution. The detailed procedure of gassy specimen preparation is introduced. Then, the quality of prepared specimens using our improved apparatus is validated by the high-resolution micro-X-ray computed tomography (μCT) scanning test, from which bubble occurrence and size distribution within the gassy sand specimen can be obtained; and preliminary mechanical tests on gassy sand specimens with various initial saturation degrees are performed. The proposed specimen preparation procedure succeeds in proving the postulated occurrence state of gas bubbles in coarse-grained sediments and accurately controlling the gas content.

58 GEOSCIENCES↗

Developing Open-Source Training Materials for AI/ML and Space Biological Sciences Using NASA Cloud-Based Data

Artificial Intelligence (AI) and Machine Learning (ML) has gained significant traction in the biological and biomedical research fields in the last two decades, in part thanks to an increasing culture of open data sharing and reuse. Due to its capability for identifying complex relationships and patterns, AI/ML methodology is particularly well suited to recognize and predict biological patterns from high-dimensional next-generation sequencing data (e.g. whole genome sequencing, transcriptomic sequencing), as well as from biological or medical imaging data (e.g. microscopy, computed tomography, ultrasound, magnetic resonance imaging, radiography). These methodologies hold particular promise for space biosciences research and automated space health monitoring systems. However, there are many key considerations for properly training, validating, and testing a machine learning model in biological research or clinical application. Even with the positive culture of Open Science and data sharing, inexperienced researchers working quickly without proper checks can produce models that perform poorly outside of the immediate training dataset. Lessons learned from biological AI/ML research indicate that Open Science principles such as data sharing and open-source code must go hand-in-hand with publicly available, high-quality training curricula in best practices, with modules centered on real-life scientific use cases and data so future AI/ML practitioners gain experience on real problems. Here we present the development of open-source training materials for AI/ML and space biosciences, as part of the NASA Transform to Open Science Training (TOPST) initiative. We develop 4 independent training programs, focused on the following topics: 1) Fundamentals of Machine Learning and Space Biosciences Domain, 2) Open Science, Artificial Intelligence, and Ethical Best Practices for Data Sharing and Analysis, 3) Using AI/ML Classification to Identify Gene Networks Affected By Space Exposure in Mouse Liver, and 4) Using Neural Networks to Find DNA Damage Patterns in Immune Cells after Radiation. All programs leverage cloud-based NASA biological datasets. The curriculum we present will enable worldwide access to training in AI/ML and scientific analysis.

James Andrew Casaletto↗

Inflation and Dark Energy from spectroscopy at $\mathcal{z} > 2$

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $\mathcal{z} < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. Most of the processes involved during Inflation impact observations on the very largest spatial scales. Traditionally, these have been accessed through observations of the Cosmic Microwave Background (CMB). While very powerful, the CMB originates from a 2D surface and the finite number of modes that it contains will largely be measured by experiments over the next decade. Observations of large 3D volumes with large-scale structure (LSS) access similar scales and will dramatically increase the number of available modes. For example, LSS observations in the range $2 \lesssim z \lesssim 5$ can more than triple the volume surveyed at $z \lesssim 2$, and, together with the sufficiently high galaxy number in this interval, strongly motivates a future spectroscopic survey that exploits this opportunity. In addition, tomography allows mapping the growth of structure with redshift, which provides robust constraints on Dark Energy and neutrino masses while relaxing restrictive assumptions such as a power-law primordial power spectrum. Finally, cross-correlation with external tracers, such as CMB lensing, Intensity Mapping or the Lyman-$\alpha$ forest, immunises the constraints to the systematics that make measurement challenging and further improves the precision through ‘sample variance cancellation’ and degeneracy breaking.

79 ASTRONOMY AND ASTROPHYSICS↗

Evaluation of Instability Phenomena in Sands: Plane Strain Versus Triaxial Conditions

Extensive research was carried out in the 1950s on theories of plasticity to extend the concepts developed for metals to materials that failed according to the Mohr-Coulomb criterion. The new ideas made it possible to merge the two distinct concepts (strength and deformation techniques) into one that relies on better understanding of plasticity and resulted in a rapid growth in the field of constitutive modeling of soil behavior. At the same time advanced experimental apparatuses and laboratory procedures were developed to calibrate the models. However, most laboratory experiments on granular materials are performed under Conventional Triaxial Conditions (CTC) for the purposes of evaluating constitutive behavior and stability properties, whereas most geotechnical field problems are closer to the Plane Strain (PS) condition. The triaxial tests performed in most laboratories comprise a simplification over in situ states and allow easier and robust experimentation. Most landslide problems, failure of soils beneath shallow and deep foundations, and failure of retaining structures, are cases that can generally be considered as plane strain. Strength and deformation characteristics of granular materials loaded in plane strain may be considerably different from those observed in CTC. Most studies on sands were limited to evaluating the constitutive behavior and in some cases extended to briefly describing the associated instability phenomena. This paper presents the results of a series of PS and CTC experiments performed on fine uniform silica sand known as F-75 Ottawa sand. Advanced analysis techniques were used to study the instability phenomena, which yielded very accurate measurements of shear bands occurrences and patterns. Destructive thin-sectioning technique along with monitoring the specimen surface deformation was used in the PS experiments and Computed Tomography (CT) was used to investigate the progress of primary and secondary shear bands in specimens subjected to CTC. Comparison between the two cases will be presented and discussed.

Alshibli, Khalid A.↗

Earthquake Locations

The location of an earthquake is a classic inverse problem. Arrival times of various seismic phases are combined with the location of the instrument to form the data vector. The inverse operator is applied in an iterative manor until the solution converges and the resultant solution is the hypocentral location of the earthquake in space and time. In this presentation, I will present how seismologists obtain local earthquake hypocenters. I will present how the generalized inverse operator is built using the simple distance = rate x time equation and how the inverse operator is used in an iterative manor to converge on the hypocenter of the earthquake. Once an earthquake location is determined, seismologists can then use the location to infer other information about the medium through which the waves traveled. One possible application of the earthquake locations is a branch of seismology known as tomography. In tomography, source and station locations are used to compute the travel times of seismic phases and these phases are then inverted to obtain a velocity model of the medium. Seismic tomography can be applied in both 2-D and 3-D settings and in this presentation, I will present examples from both. The first example comes from a technique known as ambient noise tomography. Rayleigh waves are obtained by cross correlating the background noise between two stations. Because we know the interstation distance between the two stations and the arrival times of the Rayleigh waves, we can use tomography to build a group velocity map of the study region. I will present results from a study of the Eastern Tennessee Seismic Zone using the EarthScope’s Transportable Array together with different seismic networks in the region, figure 1. The second example uses P and S-waves from 412 earthquakes which occurred during the 1994 Draney Peak swarm to calculate a 3-D velocity model of the area. Preliminary results show a strong velocity contrast at a depth of about 4 km where the majority of the earthquakes occur, figure 2.

58 GEOSCIENCES↗

Dynamic in-situ imaging of methane hydrate formation and self-preservation in porous media

Here we present the results of dynamic in-situ 3D X-ray imaging of methane hydrates microstructure during methane hydrate formation and dissociation in sand samples. Short scanning times and high resolution provided by synchrotron X-rays allowed for better understanding of water movement and different types of gas-hydrate formation. Complementing previous observations, we conclude that the process of gas-hydrate formation is accompanied by the water movements caused by cryogenic water suction that happens in sequences of short fast movements with longer equilibrium states in between (when the water is immobile). Based on the 3D microstructure we identified two distinct types of gas-hydrate formation: (i) into the gas pockets and (ii) inside water volumes. For both mechanisms we do not see problems in gas or water supply to support the gas-hydrate formation. The rate of dissociation in the self-preservation mode (pressure drop at negative temperatures) appears to be different for these two types of gas hydrates. This means that the history of the gas-hydrate formation may influence its behaviour at the dissociation stage (e.g. gas-hydrate production).

58 GEOSCIENCES↗

A Voxel-Based Approach to Explore Local Dose Differences Associated With Radiation-Induced Lung Damage

Purpose: To apply a voxel-based (VB) approach aimed at exploring local dose differences associated with late radiation-induced lung damage (RILD). Methods and Materials: An interinstitutional database of 98 patients who were Hodgkin lymphoma (HL) survivors treated with postchemotherapy supradiaphragmatic radiation therapy was analyzed in the study. Eighteen patients experienced late RILD, classified according to the Radiation Therapy Oncology Group scoring system. Each patient's computed tomographic (CT) scan was normalized to a single reference case anatomy (common coordinate system, CCS) through a log-diffeomorphic approach. The obtained deformation fields were used to map the dose of each patient into the CCS. The coregistration robustness and the dose mapping accuracy were evaluated by geometric and dose scores. Two different statistical mapping schemes for nonparametric multiple permutation inference on dose maps were applied, and the corresponding P<.05 significance lung subregions were generated. A receiver operating characteristic (ROC)-based test was performed on the mean dose extracted from each subregion. Results: The coregistration process resulted in a geometrically robust and accurate dose warping. A significantly higher dose was consistently delivered to RILD patients in voxel clusters near the peripheral medial-basal portion of the lungs. The area under the ROC curves (AUC) from the mean dose of the voxel clusters was higher than the corresponding AUC derived from the total lung mean dose. Conclusions: We implemented a framework including a robust registration process and a VB approach accounting for the multiple comparison problem in dose-response modeling, and applied it to a cohort of HL survivors to explore a local dose–RILD relationship in the lungs. Patients with RILD received a significantly greater dose in parenchymal regions where low doses (∼6 Gy) were delivered. Interestingly, the relation between differences in the high-dose range and RILD seems to lack a clear spatial signature.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Multiscale evolution of charmed particles in a nuclear medium

Parton energy-momentum exchange with the quark gluon plasma (QGP) is a multiscale problem. In this work, we calculate the interaction of charm quarks with the QGP within the higher twist formalism at high virtuality and high energy using the Modular All Twist Transverse-scattering Elastic-drag and Radiation (MATTER) model, while the low-virtuality and high-energy portion is treated via a linearized Boltzmann transport formalism. Coherence effect that reduces the medium-induced emission rate in the MATTER model is also taken into account through a virtuality-dependent qˆ, leaving the simultaneous dependence of qˆ on heavy quark mass and virtuality for future studies. The interplay between these two formalisms is studied phenomenologically and used to produce a first description of the D-meson and charged hadron nuclear modification factor R AA across multiple centralities. As a result, all calculations were carried out utilizing the Jet Energy-loss Tomography with a Statistically and Computationally Advanced Program Envelope framework.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Machine Learning for Predictive Performance Analysis in Charged Particle Beam Tools

Imaging methods driven by probes, electrons, and ions have played a dominant role in modern science and engineering. Opportunities for machine vision and AI that focus on consumer problems like driving and feature recognition, are now presenting themselves for automating aspects of the scientific processes. This proposal aims to enable and drive discovery in ultra-low energy implantation by taking advantage of faster processing, flexible control and detection methods, and architecture-agnostic workflows that will result in higher efficiency and shorter scientific development cycles. Custom microscope control, collection and analysis hardware will provide a framework for conducting novel in situ experiments revealing unprecedented insight into surface dynamics at the nanoscale. Ion implantation is a key capability for the semiconductor industry. As devices shrink, novel materials enter the manufacturing line, and quantum technologies transition to being more mainstream. Traditional implantation methods fall short in terms of energy, ion species, and positional precision. Here we demonstrate 1 keV focused ion beam Au implantation into Si and validate the results via atom probe tomography. We show the Au implant depth at 1 keV is 0.8 nm and that identical results for low energy ion implants can be achieved by either lowering the column voltage, or decelerating ions using bias – while maintaining a sub-micron beam focus. We compare our experimental results to static calculations using SRIM and dynamic calculations using binary collision approximation codes TRIDYN and IMSIL. A large discrepancy between the static and dynamic simulation is found that is due to lattice enrichment with high stopping power Au and surface sputtering. Additionally, we demonstrate how model details are particularly important to the simulation of these low-energy heavy-ion implantations. Finally, we discuss how our results pave a way to much lower implantation energies, while maintaining high spatial resolution.

47 OTHER INSTRUMENTATION↗

Baylor University Campus-Wide Deep Dive

In January 2020, staff members from the Engagement and Performance Operations Center (EPOC) and the Lonestar Education And Research Network (LEARN) met with researchers and staff at Baylor University for the purpose of a Campus-Wide Deep Dive into research drivers. The goal of this meeting was to help characterize the requirements for five campus research use cases and to enable cyberinfrastructure support staff to better understand the needs of the researchers they support. Profiled scientific use cases included: - Experimental High Energy Physics (HEP) - Proton Computed Tomography (pCT) - Nutrition and Relation to Digestive Microbiome - Baylor University Core Research Facilities - Molecular Quantum-dot Cellular Automata (QCA), and Material Science of Quantum Computing - Modeling and Simulation of Low-Dimensional and Nano-Structured Materials - Computational Fluid Dynamics Material for this event included the written documentation from each of the research areas at Baylor University, documentation about the current state of technology support, and a write-up of the discussion that took place in person. The Case Studies highlighted the ongoing challenges that Baylor University has in supporting a cross-section of established and emerging research use cases. Each Case Study mentioned unique challenges which were summarized into common needs. These included: - Tradeoffs for network/software security, and usability of the resulting infrastructure. Better communication to set expectations and understand realities is required. - Computation use on campus is widespread and healthy. While no major problems were uncovered, upgrades to maintain current usage patterns and encourage growth will be required. - Storage is a critical need for enterprise use cases and research. In particular, a campus wide ‘storage architecture’ to support research use cases (e.g. instruments, data sharing) is required in the 2-5 year time window. - Instrumentation on campus is healthy and expanding. Technology must scale with this in the form of computation and storage. - Working with LEARN to upgrade network capacity (in multiples of 10G, or upgrades to 100G) will be required in the 1-3 year time frame. - Network monitoring and visibility will help to establish external science use cases. - Data sharing via portal systems is not currently a critical need, but growing in scope. EPOC can assist Baylor with options.

99 GENERAL AND MISCELLANEOUS↗

SQMS Quantum R&D in Machine Learning, Optimization and Sensing beyond Fundamental Physics Applications

This newly formed team at SQMS under the Ecosystem Thrust is looking to develop capabilities impacting societal advances outside the core domain of HEP and condensed matter physics. We explicitly leverage the experimental and algorithmic innovations developed across all groups as well as connect to broad-scope external projects of the diverse team of PIs. As the inaugural set of projects, we are studying numerically quantum machine learning models inspired by efficiently trainable echo-state and orthogonal neural networks and developing designs for related experiments to be performed on quantum processors based on SQMS SRF cQED technology and Rigetti s transmon arrays. Investigated models exploit ideas and lessons learned from multiple prior work by SQMS team members in a variety of internal and external activities [R1]. Target initial applications include noisy signal processing, potentially captured by quantum sensors or noisy QPUs, as well as simulation and classification of healthcare data. For instance, image reconstruction of the brain s electrical properties by solving the inverse Maxwell equation problem with uncertainty [R2] through a hybrid quantum-classical physics-informed architecture for time-dependent processes [R3]. The group is also investigating the application and development of novel quantum sensors based on magnetic levitation of a superconducting sphere coupled to a superconducting qubit. This coupling enables high-precision measurements of the position of the sphere, which can be used for sensitive detection of forces, enabling practical applications such as gravimetry for geophysics analysis, or accelerometry for GPS-denied navigation [R4] [R1] Rieffel, Eleanor G., Ata Akbari Asanjan, M. Sohaib Alam, Namit Anand, David E. Bernal Neira, Sophie Block, Lucas T. Brady et al. "Assessing and advancing the potential of quantum computing: A NASA case study." Future Generation Computer Systems (2024). [R2] Yu, X., Serrall s, J.E., Giannakopoulos, I.I., Liu, Z., Daniel, L., Lattanzi, R. and Zhang, Z., 2023. Pifon-ept: Mr-based electrical property tomography using physics-informed fourier networks. IEEE Journal on Multiscale and Multiphysics Computational Techniques. [R3] Wudarski, Filip, Daniel OConnor, Shaun Geaney, Ata Akbari Asanjan, Max Wilson, Elena Strbac, P. Aaron Lott, and Davide Venturelli. "Hybrid quantum-classical reservoir computing for simulating chaotic systems." arXiv preprint arXiv:2311.14105 (2023). [R4] Higgins, Gerard, Saarik Kalia, and Zhen Liu. "Maglev for dark matter: Dark-photon and axion dark matter sensing with levitated superconductors." Physical Review D 109.5 (2024): 055024.

Venturelli, Davide↗

Elucidating the Onset of Plasticity in Sliding Contacts Using Differential Computational Orientation Tomography

Abstract Depending on the mechanical and thermal energy introduced to a dry sliding interface, the near-surface regions of the mated bodies may undergo plastic deformation. In this work, we use large-scale molecular dynamics simulations to generate “differential computational orientation tomographs” (dCOT) and thus highlight changes to the microstructure near tribological FCC alloy surfaces, allowing us to detect subtle differences in lattice orientation and small distances in grain boundary migration. The analysis approach compares computationally generated orientation tomographs with their undeformed counterparts via a simple image analysis filter. We use our visualization method to discuss the acting microstructural mechanisms in a load- and time-resolved fashion, focusing on sliding conditions that lead to twinning, partial lattice rotation, and grain boundary-dominated processes. Extracting and laterally averaging the color saturation value of the generated tomographs allows us to produce quantitative time- and depth-resolved maps that give a good overview of the progress and severity of near-surface deformation. Corresponding maps of the lateral standard deviation in the color saturation show evidence of homogenization processes occurring in the tribologically loaded microstructure, frequently leading to the formation of a well-defined separation between deformed and undeformed regions. When integrated into a computational materials engineering framework, our approach could help optimize material design for tribological and other deformation problems. Graphic Abstract .

36 MATERIALS SCIENCE↗

Observing crystal nucleation in four dimensions using atomic electron tomography

Nucleation plays a critical role in many physical and biological phenomena that range from crystallization, melting and evaporation to the formation of clouds and the initiation of neurodegenerative diseases. However, nucleation is a challenging process to study in experiments, especially in the early stage when several atoms or molecules start to form a new phase from a parent phase. A number of experimental and computational methods have been used to investigate nucleation processes, but experimental determination of the three-dimensional atomic structure and the dynamics of early stage nuclei has been unachievable. Here we use atomic electron tomography to study early stage nucleation in four dimensions (4D: that is, including time) at atomic resolution. Using FePt nanoparticles as a model system, we find that early stage nuclei are irregularly shaped, each has a core of one to a few atoms with the maximum order parameter, and the order parameter gradient points from the core to the boundary of the nucleus. We capture the structure and dynamics of the same nuclei undergoing growth, fluctuation, dissolution, merging and/or division, which are regulated by the order parameter distribution and its gradient. These experimental observations are corroborated by molecular dynamics simulations of heterogeneous and homogeneous nucleation in liquid–solid phase transitions of Pt. Our experimental and molecular dynamics results differ from classical nucleation theory, indicating that a theory beyond this is needed to describe early stage nucleation at the atomic scale. Looking forward, we anticipate that the reported approach will open the door to the study of many fundamental problems in materials science, nanoscience, condensed matter physics and chemistry, such as phase transition, atomic diffusion, grain boundary dynamics, interface motion, defect dynamics and surface reconstruction with 4D atomic resolution.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Synthetic modeling of soft x-ray emissivity for magnetic island analysis in LTX-β

We present a synthetic soft x-ray (SXR) forward-modeling framework that characterizes the emissivity structure of rotating magnetic islands in the Lithium Tokamak eXperiment-β across space and time. Magnetic islands associated with tearing modes produce modulations in line-integrated SXR brightness. Traditional tomographic methods struggle to resolve these structures in devices with limited sightlines, resulting in an under-determined inversion problem, or otherwise necessitate equilibrium reconstruction data, preventing their use in active control. Our approach removes this challenge by combining a fully three-dimensional ray-tracing model with a time-dependent emissivity prescription of a $m/n$ = 2/1 magnetic island. The model is validated through helical island geometry, rotation measured by magnetic diagnostics, and equilibrium constraints from PSI-Tri reconstructions. Synthetic brightness signals are generated for each photodiode sightline from a modeled emissivity profile and directly compared with experimental data from a tangential midplane SXR array. By fitting the synthetic diagnostic output to the observed brightness across time, we simultaneously infer key geometric parameters-such as island radial location, island width, and rotation frequency-without requiring full tomographic reconstruction, all with an average absolute deviation of less than 3%. This work demonstrates that forward modeling with a single tangential array can extract key island parameters in a spherical tokamak, providing a computationally efficient alternative to conventional SXR tomography and providing a pathway toward real-time magnetic-island characterization in future devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quantitative imaging and automated fuel pin identification for passive gamma emission tomography

Compliance of member States to the Treaty on the Non-Proliferation of Nuclear Weapons is monitored through nuclear safeguards. The Passive Gamma Emission Tomography (PGET) system is a novel instrument developed within the framework of the International Atomic Energy Agency (IAEA) project JNT 1510, which included the European Commission, Finland, Hungary and Sweden. The PGET is used for the verification of spent nuclear fuel stored in water pools. Advanced image reconstruction techniques are crucial for obtaining high-quality cross-sectional images of the spent-fuel bundle to allow inspectors of the IAEA to monitor nuclear material and promptly identify its diversion. In this work, we have developed a software suite to accurately reconstruct the spent-fuel cross sectional image, automatically identify present fuel rods, and estimate their activity. Unique image reconstruction challenges are posed by the measurement of spent fuel, due to its high activity and the self-attenuation. While the former is mitigated by detector physical collimation, we implemented a linear forward model to model the detector responses to the fuel rods inside the PGET, to account for the latter. The image reconstruction is performed by solving a regularized linear inverse problem using the fast-iterative shrinkage-thresholding algorithm. We have also implemented the traditional filtered back projection (FBP) method based on the inverse Radon transform for comparison and applied both methods to reconstruct images of simulated mockup fuel assemblies. Higher image resolution and fewer reconstruction artifacts were obtained with the inverse-problem approach, with the mean-square-error reduced by 50%, and the structural-similarity improved by 200%. We then used a convolutional neural network (CNN) to automatically identify the bundle type and extract the pin locations from the images; the estimated activity levels finally being compared with the ground truth. The proposed computational methods accurately estimated the activity levels of the present pins, with an associated uncertainty of approximately 5%.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Anatomical pulmonary magnetic resonance imaging segmentation for regional structure-function measurements of asthma

Purpose: Pulmonary magnetic-resonance-imaging (MRI) and x-ray computed-tomography have provided strong evidence of spatially and temporally persistent lung structure-function abnormalities in asthmatics. This has generated a shift in their understanding of lung disease and supports the use of imaging biomarkers as intermediate endpoints of asthma severity and control. In particular, pulmonary {sup 1}H MRI can be used to provide quantitative lung structure-function measurements longitudinally and in response to treatment. However, to translate such biomarkers of asthma, robust methods are required to segment the lung from pulmonary {sup 1}H MRI. Therefore, their objective was to develop a pulmonary {sup 1}H MRI segmentation algorithm to provide regional measurements with the precision and speed required to support clinical studies. Methods: The authors developed a method to segment the left and right lung from {sup 1}H MRI acquired in 20 asthmatics including five well-controlled and 15 severe poorly controlled participants who provided written informed consent to a study protocol approved by Health Canada. Same-day spirometry and plethysmography measurements of lung function and volume were acquired as well as {sup 1}H MRI using a whole-body radiofrequency coil and fast spoiled gradient-recalled echo sequence at a fixed lung volume (functional residual capacity + 1 l). We incorporated the left-to-right lung volume proportion prior based on the Potts model and derived a volume-proportion preserved Potts model, which was approximated through convex relaxation and further represented by a dual volume-proportion preserved max-flow model. The max-flow model led to a linear problem with convex and linear equality constraints that implicitly encoded the proportion prior. To implement the algorithm, {sup 1}H MRI was resampled into ∼3 × 3 × 3 mm{sup 3} isotropic voxel space. Two observers placed seeds on each lung and on the background of 20 pulmonary {sup 1}H MR images in a randomized dataset, on five occasions, five consecutive days in a row. Segmentation accuracy was evaluated using the Dice-similarity-coefficient (DSC) of the segmented thoracic cavity with comparison to five-rounds of manual segmentation by an expert observer. The authors also evaluated the root-mean-squared-error (RMSE) of the Euclidean distance between lung surfaces, the absolute, and percent volume error. Reproducibility was measured using the coefficient of variation (CoV) and intraclass correlation coefficient (ICC) for two observers who repeated segmentation measurements five-times. Results: For five well-controlled asthmatics, forced expiratory volume in 1 s (FEV{sub 1})/forced vital capacity (FVC) was 83% ± 7% and FEV{sub 1} was 86 ± 9%{sub pred}. For 15 severe, poorly controlled asthmatics, FEV{sub 1}/FV C = 66% ± 17% and FEV{sub 1} = 72 ± 27%{sub pred}. The DSC for algorithm and manual segmentation was 91% ± 3%, 92% ± 2% and 91% ± 2% for the left, right, and whole lung, respectively. RMSE was 4.0 ± 1.0 mm for each of the left, right, and whole lung. The absolute (percent) volume errors were 0.1 l (∼6%) for each of right and left lung and ∼0.2 l (∼6%) for whole lung. Intra- and inter-CoV (ICC) were <0.5% (>0.91%) for DSC and <4.5% (>0.93%) for RMSE. While segmentation required 10 s including ∼6 s for user interaction, the smallest detectable difference was 0.24 l for algorithm measurements which was similar to manual measurements. Conclusions: This lung segmentation approach provided the necessary and sufficient precision and accuracy required for research and clinical studies.

60 APPLIED LIFE SCIENCES↗