Search NASA⌕ Search

SEARCH · Search NASA

Results for “Level set method”

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 1,315 records · Page 73

Application of dynamic metabolomics to examine in vivo skeletal muscle glucose metabolism in the chronically high-fat fed mouse

Rationale: Defects in muscle glucose metabolism are linked to type 2 diabetes. Mechanistic studies examining these defects rely on the use of high fat-fed rodent models and typically involve the determination of muscle glucose uptake under insulin-stimulated conditions. While insightful, they do not necessarily reflect the physiology of the postprandial state. In addition, most studies do not examine aspects of glucose metabolism beyond the uptake process. Here we present an approach to study rodent muscle glucose and intermediary metabolism under the dynamic and physiologically relevant setting of the oral glucose tolerance test (OGTT). Methods and results: In vivo muscle glucose and intermediary metabolism was investigated following oral administration of [U-{sup 13}C] glucose. Quadriceps muscles were collected 15 and 60 min after glucose administration and metabolite flux profiling was determined by measuring {sup 13}C mass isotopomers in glycolytic and tricarboxylic acid (TCA) cycle intermediates via gas chromatography–mass spectrometry. While no dietary effects were noted in the glycolytic pathway, muscle from mice fed a high fat diet (HFD) exhibited a reduction in labelling in TCA intermediates. Interestingly, this appeared to be independent of alterations in flux through pyruvate dehydrogenase. In addition, our findings suggest that TCA cycle anaplerosis is negligible in muscle during an OGTT. Conclusions: Under the dynamic physiologically relevant conditions of the OGTT, skeletal muscle from HFD fed mice exhibits alterations in glucose metabolism at the level of the TCA cycle. - Highlights: • Dynamic metabolomics was used to investigate muscle glucose metabolism in vivo. • Mitochondrial TCA cycle metabolism is altered in muscle of HFD mice. • This defect was not pyruvate dehydrogenase mediated, as has been previously thought. • Mitochondrial TCA cycle anaplerosis in muscle is virtually absent during the OGTT.

60 APPLIED LIFE SCIENCES↗

Carbon Dioxide Conversion to Nanomaterials: Methods, Applications, and Challenges

Anthropogenic carbon dioxide (CO 2 ) emission is one of the biggest global threats that we are confronted with in the 21st century. A lack of prompt action to mitigate the atmospheric CO 2 level could lead to global catastrophic events. As a result of the significance of the aforementioned threat, multiple nations have agreed to enact carbon capture, storage, utilization, and sequestration (CCUS) technologies at academic and industrial scales with the aim of lowering the atmospheric CO 2 level. Although there are many reviews on CCUS technologies, there is a lack of comprehensive review on CO 2 conversion to more economically valuable products, such as nanomaterials. Carbon dioxide capture and conversion into nanomaterials is an interesting research topic in many ways because it contributes to our fundamental knowledge of material conversion, provides alternative ways to reduce atmospheric pollution, and yields lucrative products, which has not been explored broadly. Aside from the positive outcomes, we have a moral obligation to rectify a global threat caused by anthropogenic emission. Herein, we present current methods of converting CO 2 into organic and inorganic nanomaterials with a focus on operating conditions and challenges as well as the potential of nanomaterials in specifically biomedical applications as opposed to extensively reviewed energy-related uses. Furthermore, our review emphasizes the toxicity issues related to using such nanomaterials in biological settings based on literature and other comparable studies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

BioCARS: Synchrotron facility for probing structural dynamics of biological macromolecules

A major goal in biomedical science is to move beyond static images of proteins and other biological macromolecules to the internal dynamics underlying their function. This level of study is necessary to understand how these molecules work and to engineer new functions and modulators of function. Stemming from a visionary commitment to this problem by Keith Moffat decades ago, a community of structural biologists has now enabled a set of x-ray scattering technologies for observing intramolecular dynamics in biological macromolecules at atomic resolution and over the broad range of timescales over which motions are functionally relevant. Many of these techniques are provided by BioCARS, a cutting-edge synchrotron radiation facility built under Moffat leadership and located at the Advanced Photon Source at Argonne National Laboratory. BioCARS enables experimental studies of molecular dynamics with time resolutions spanning from 100 ps to seconds and provides both time-resolved x-ray crystallography and small- and wide-angle x-ray scattering. Structural changes can be initiated by several methods—UV/Vis pumping with tunable picosecond and nanosecond laser pulses, substrate diffusion, and global perturbations, such as electric field and temperature jumps. Studies of dynamics typically involve subtle perturbations to molecular structures, requiring specialized computational techniques for data processing and interpretation. In this review, we present the challenges in experimental macromolecular dynamics and describe the current state of experimental capabilities at this facility. As Moffat imagined years ago, BioCARS is now positioned to catalyze the scientific community to make fundamental advances in understanding proteins and other complex biological macromolecules.

59 BASIC BIOLOGICAL SCIENCES↗

Spin-orbit coupling induced demagnetization in Ni: Ab initio nonadiabatic molecular dynamics perspective

Spin-orbit coupling (SOC), which can induce spin flip during the relaxation of photoexcited charge carrier, plays a crucial role in spin dynamics. Here in this work, we have used time-domain ab initio nonadiabatic molecular dynamics (NAMD) method to study the SOC induced ultrafast demagnetization in Ni at 300K. The spin-diabatic representation using spin-polarized Kohn-Sham (KS) basis sets and spin-adiabatic representation using spinor basis sets have been applied, and both of them achieve demagnetization in Ni with a timescale around 100fs. The spin-diabatic representation suggests a picture that the electron-phonon coupling (EPC) provides direct energy relaxation channel among the same-spin states, while the SOC can induce spin flip. After photoexcitation, it is found the spin-minority electrons relax to the same-spin states rather than the opposite-spin states, since EPC is larger than SOC by one order of magnitude. By contrast, for the spin-majority electrons, spin flip occurs since there are no empty same-spin states as electron acceptor above the Fermi level. The different relaxation pathways for spin-majority and spin-minority electrons induce the demagnetization. The spin-adiabatic representation provides an Elliott-Yafet spin-phonon scattering picture. The SOC induced reduction of magnetic moment in Ni may induce magnon to drive further demagnetization. The ab initio NAMD simulation provides a critical angle to understand how the SOC and EPC affect demagnetization process in Ni.

36 MATERIALS SCIENCE↗

Estimating local noise power spectrum from a few FBP-reconstructed CT scans

Purpose: Traditional ways to estimate 2D CT noise power spectrum (NPS) involve an ensemble average of the power spectrums of many noisy scans. When only a few scans are available, regions of interest are often extracted from different locations to obtain sufficient samples to estimate the NPS. Using image samples from different locations ignores the nonstationarity of CT noise and thus cannot accurately characterize its local properties. The purpose of this work is to develop a method to estimate local NPS using only a few fan-beam CT scans. Methods: As a result of FBP reconstruction, the CT NPS has the same radial profile shape for all projection angles, with the magnitude varying with the noise level in the raw data measurement. This allows a 2D CT NPS to be factored into products of a 1D angular and a 1D radial function in polar coordinates. The polar separability of CT NPS greatly reduces the data requirement for estimating the NPS. The authors use this property and derive a radial NPS estimation method: in brief, the radial profile shape is estimated from a traditional NPS based on image samples extracted at multiple locations. The amplitudes are estimated by fitting the traditional local NPS to the estimated radial profile shape. The estimated radial profile shape and amplitudes are then combined to form a final estimate of the local NPS. We evaluate the accuracy of the radial NPS method and compared it to traditional NPS methods in terms of normalized mean squared error (NMSE) and signal detectability index. Results: For both simulated and real CT data sets, the local NPS estimated with no more than six scans using the radial NPS method was very close to the reference NPS, according to the metrics of NMSE and detectability index. Even with only two scans, the radial NPS method was able to achieve a fairly good accuracy. Compared to those estimated using traditional NPS methods, the accuracy improvement was substantial when a few scans were available. Conclusions: The radial NPS method was shown to be accurate and efficient in estimating the local NPS of FBP-reconstructed 2D CT images. It presents strong advantages over traditional NPS methods when the number of scans is limited and can be extended to estimate the in-plane NPS of cone-beam CT and multislice helical CT scans.

61 RADIATION PROTECTION AND DOSIMETRY↗

Dynamic Tuning of Seismic Signal Detector Trigger Levels for Local Networks

The quality of automatic signal detections from sensor networks depends on individual detector trigger levels (TLs) from each sensor. The largely manual process of identifying effective TLs is painstaking and does not guarantee optimal configuration settings, yet achieving superior automatic detection of signals and ultimately, events, is closely related to these parameters. In this paper, we present a Dynamic Detector Tuning (DDT) system that automatically adjusts effective TL settings for signal detectors to the current state of the environment by leveraging cooperation within a local neighborhood of network sensors. After a stabilization period, the DDT algorithm can adapt in near-real time to changing conditions and automatically tune a signal detector to identify (detect) signals from only events of interest. Our current work focuses on reducing false signal detections early in the seismic signal processing pipeline, which leads to fewer false events and has a significant impact on reducing analyst time and effort. This system provides an important new method to automatically tune detector TLs for a network of sensors and is applicable to both existing sensor performance boosting and new sensor deployment. Finally, with ground truth on detections from a local neighborhood of seismic sensors within a network monitoring the Mount Erebus volcano in Antarctica, we show that DDT reduces the number of false detections by 18% and the number of missed detections by 11% when compared with optimal fixed TLs for all sensors.

58 GEOSCIENCES↗

Bridging Equipment Reliability Data and Robust Decisions in a Plant Operation Context

In order to reduce operation and maintenance (O&M) costs, nuclear power plants (NPPs) are moving from corrective and periodic maintenance to predictive maintenance strategies. Such transition requires changes on the data that needs to be retrieved and on the type of decision processes to be employed. Advanced monitoring and data analysis technologies are essential to support predictive strategies. They can in fact provide precise information about health of a component, track its degradation trends, and provide information of its expected failure time. With such information, maintenance operations for a component can be performed right before its expected failure time. This dynamic context of O&M operations requires new methods to analyze data, propagate component health information from the component to the system level, and optimize plant resources. In this respect, the risk informed asset management (RIAM) project has been tasked to develop and test this new class of methods into a risk analytics toolset. This toolset consists of data analytics tools coupled with reliability methods designed to manage plant assets and performances in a predictive maintenance context. This report shows the latest improvements on such development and the initial testing of our methods on the three main research areas that the RIAM project is focusing on. These areas are the following: equipment reliability data analytics, system reliability modeling, and plant resources optimization methods. We show how the methods developed in these areas can support predictive maintenance strategies by: 1) analyzing equipment reliability data (either in numeric and textual form), 2) assessing component and system health through an innovative margin-based reliability approach, and 3) identifying the most critical components and set optimal maintenance schedule based on plant economic and operational constraints.

97 MATHEMATICS AND COMPUTING↗

Identification of synthetic consortia from a set of plant-beneficial bacteria

The use of microbial inoculants in agriculture as biofertilisers and/or biopesticides is an appealing alternative to replace or reduce the practice of agrochemicals. Plant microbiota studies are revealing the different bacterial groups which are populating plant microbiomes re-energising the plant probiotic bacteria (PPB) translational research sector. Some single-microbial strain bioinoculants have proven valid in agriculture (e.g., based on Trichoderma, mycorrhiza or rhizobia); however, it is now recommended to consider multistrain consortia since plant-beneficial effects are often a result of community-level interactions in plant microbiomes. A limiting step is the selection of a fitting combination of microbial strains in order to accomplish the best beneficial effect upon plant inoculation. In this study, we have used a subset of 23 previously identified and characterised rice-beneficial bacterial colonisers to design and test a series of associated experiments aimed to identify potential PPB consortia which are able to co-colonise and induce plant growth promotion. Bacterial strains were co-inoculated in vitro and in planta using several different methods and their co-colonisation and co-persistence monitored. Results include the identification of two 5-strain and one 2-strain consortia which displayed plant growth-promoting features. Future practical applications of microbiome research must include experiments aimed at identifying consortia of bacteria which can be most effective as crop amendments.

59 BASIC BIOLOGICAL SCIENCES↗

Path optimization with limited sensing ability

We propose a computational strategy to find the optimal path for a mobile sensor with limited coverage to traverse a cluttered region. The goal is to find one of the shortest feasible paths to achieve the complete scan of the environment. We pose the problem in the level set framework, and first consider a related question of placing multiple stationary sensors to obtain the full surveillance of the environment. By connecting the stationary locations using the nearest neighbor strategy, we form the initial guess for the path planning problem of the mobile sensor. Then the path is optimized by reducing its length, via solving a system of ordinary differential equations (ODEs), while maintaining the complete scan of the environment. Furthermore, we use intermittent diffusion, which converts the ODEs into stochastic differential equations (SDEs), to find an optimal path whose length is globally minimal. To improve the computation efficiency, we introduce two techniques, one to remove redundant connecting points to reduce the dimension of the system, and the other to deal with the entangled path so the solution can escape the local traps. Numerical examples are shown to illustrate the effectiveness of the proposed method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A method for estimating the amplitude response of smartphone built-in microphone sensors below 4 kHz

A method for estimating the amplitude response of built-in smartphone microphone sensors is pre- sented. The method is intended to be accessible to the general public, and comparison calibration measurements are performed in a non-isolated environment using conventional consumer products. A double reference sensor scheme is set up, consisting of a MB3 digital microbarometer (reference over the 0.5–2 Hz octave bands) and the iPrecisionMic smartphone microphone (reference above the 2 Hz octave band). The amplitude response of the iPrecisionMic sensor is first evaluated over the 1–2000 Hz octave bands. The amplitude response of three Samsung Galaxy S8 built-in smart- phone microphone sensors is then measured over the 0.5–2000 Hz octave bands. The Redvox Infrasound Recorder application (app) for Android is used to measure acoustic signals with the built-in smartphone microphone sensors. Amplitude response models in terms of digital gain are estimated for the test sensors based on the results. Last, self-noise levels for the iPrecisionMic and Samsung Galaxy S8 microphones are estimated and compared to infrasound ambient noise models. Results show an experimental capability for estimating the amplitude response of built-in smartphone microphone sensors in a non-isolated environment with conventional consumer products.

47 OTHER INSTRUMENTATION↗

Solvent hold tank sample results for MCU-15-914-915-916. December 2015 Monthly sample

Savannah River National Laboratory (SRNL) received one set of Solvent Hold Tank (SHT) samples (MCU-15-914-915-916), pulled on 12/22/2015 for analysis. The samples were combined and analyzed for composition. Analysis of the composite sample MCU-15-914-915-916 indicated the TiDG, Isopar™L, and MaxCalix are at nominal levels. The modifier concentration is 3% below its nominal concentration. This analysis confirms the addition of TiDG, MaxCalix, and modifier to the solvent in November 2015. Based on the current monthly sample, the levels of TiDG, Isopar™L, MaxCalix, and modifier are sufficient for continuing operation but are expected to decrease with time. Periodic characterization and trimming additions to the solvent are recommended. No impurities above the 1000 ppm level were found in this solvent by the Semi-Volatile Organic Analysis (SVOA). No impurities were observed in the Hydrogen Nuclear Magnetic Resonance (HNMR). However, the Fourier transform infra-red spectroscopy (FTIR) method detected trace levels (a few ppm) of amides (more indicative of bacteria than a possible degradation product of TiDG). In addition, up to 18 ± 4 micrograms of mercury per gram of solvent (or 14.8 μg/mL) was detected in this sample. The current gamma concentration level (8.48E4 dpm/mL) confirmed that the gamma concentration has returned to the previous level where the process operated normally as expected. The laboratory will continue to monitor the quality of the solvent in particular for any new impurities or degradation of the solvent components.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

SU-F-R-18: Updates to the Computational Environment for Radiological Research for Image Analysis

Purpose: To present new tools in CERR for Texture Analysis and Visualization. Method: (1) Quantitative Image Analysis: We added the ability to compute Haralick texture features based on local neighbourhood. The Texture features depend on many parameters used in their derivation. For example: (a) directionality, (b) quantization of image, (c) patch-size for the neighborhood, (d) handling of the edge voxels within the region of interest, (e) Averaging co-occurance matrix vs texture features for different directions etc. A graphical user interface was built to set these parameters and then visualize their impact on the resulting texture maps. The entire functionality was written in Matlab. Array indexing was used to speed up the texture calculation. The computation speed is very competitive with the ITK library. Moreover, our implementation works with multiple CPUs and the computation time can be further reduced by using multiple processor threads. In order to reduce the Haralick texture maps into scalar features, we propose the use of Texture Volume Histograms. This lets users make use of the entire distribution of texture values within the region of interest rather than using just the mean and the standard deviations. (2) Qualitative/Visualization tools: The derived texture maps are stored as a new scan (derived) within CERR’s planC data structure. A display that compares various scans was built to show the raw image and the derived texture maps side-by-side. These images are positionally linked and can be navigated together. CERR’s graphics handling was updated and sped-up to be compatible with the newer Matlab versions. As a result, the users can use (a) different window levels and colormaps for different viewports, (b) click-and-drag or use mouse scroll-wheel to navigate slices. Results: The new features and updates are available via https://www.github.com/adityaapte/cerr . Conclusion: Features added to CERR increase its utility in Radiomics and Outcomes modeling.

60 APPLIED LIFE SCIENCES↗

Effects of knock intensity measurement technique and fuel chemical composition on the research octane number (RON) of FACE gasolines: Part 1 – Lambda and knock characterization

The Research and Motor Octane Number (RON and MON) rate the knock propensity of gasoline in the Cooperative Fuel Research (CFR) engine by comparing the knock intensity of sample fuels relative to that of primary reference fuels (PRF), a binary blend of iso-octane and n-heptane. Important differences exist between standard octane testing and automotive spark ignition (SI) engine knock testing including speed, load, air-to-fuel equivalence ratio (lambda), and knock characterization, which lead to a discrepancy between a fuel’s RON rating and its knock resistance characterized on an automotive SI engine based on knock-limited spark advance. This publication is the first of a set of three publications which modify operating parameters of the RON test method (ASTM D2699) to investigate the effects of these differences with automotive SI engine knock-limited spark advance testing. A fuel’s standard RON is evaluated at the lambda of the highest knock intensity, whereas automotive SI engines typically operate at stoichiometry. Here, we analyze the effects of a stoichiometric RON rating methodology. Additionally, the knock intensity response from the standard CFR knockmeter system is compared to a cylinder pressure oscillation-based knock intensity at several lambda settings. All experiments were performed with a set of seven Coordinating Research Council (CRC) Fuels for Advanced Combustion Engines (FACE) gasolines with approximately 95 RON. The fuel chemical composition impacted the lambda of the highest knock intensity, which resulted in fuel-specific offsets between the standard and stoichiometric RON ratings. The knock system comparison showed significant offsets between cylinder pressure-based and knockmeter-based knock intensity levels.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Human factors and modeling methods in the development of control room modernization concepts - 315

Control room modernization is one of the most challenging and complex upgrade projects that a nuclear power plant can undertake. It can even have a bigger impact on operations than, for example, turbine replacement. The challenges of migrating an analog control system to a distributed control system are already well known and a number of nuclear utilities have embarked upon various levels of effort to upgrade some of the systems in the control room. When planning for control room upgrades, plants have to deal with a multitude of engineering and operational impacts. This will inevitably include several human factors considerations, including physical ergonomics of workstations, viewing angles, lighting, seating, new communication requirements, and new concepts of operation. In helping nuclear power utilities to deal with these challenges, the Idaho National Laboratory (INL) has developed human factors design and evaluation methods to be used in the development of an end-state control room concept and to manage the various phases of the upgrade life cycle. This included interactive sessions with operators in INL's Human System Simulation Laboratory, and three-dimensional modeling to visualize control board changes. To make the various upgrade phases as concrete and as visible as possible, the end-state concept includes a set of visual representations of the control room before and after various upgrade phases to provide the context and a framework within which to consider the various options in the upgrade. This includes the various control systems, human-system interfaces to be replaced, and possible changes to operator workstations. This paper describes how this framework helps to ensure an integrated and cohesive outcome that is consistent with human factors engineering principles and also provide substantial improvement in operator performance. The paper further describes the application of this integrated approach in the strategic modernization program at a nuclear power plant where legacy systems are upgraded to advanced digital technologies through a systematic process that links human factors principles to the systems engineering process. This approach results in an integrated control room architecture beyond what is possible for individual subsystem upgrades alone. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Extended Finite Element Based Approach in Additive Manufacturing Modeling for Optimizing Highly Complex Manifold in Protonic Ceramic Electrochemical Cells

The objective of this project is to develop a novel numerical approach based on the extended finite element method (XFEM) for modeling moving boundaries of material deposited in additive manufacturing (AM) process with improved accuracy and reduced computational cost. One of the major challenges in simulating AM processes is that the boundaries of the computational domain need to change as material is deposited. Previous approaches have achieved this by activating new finite elements, but this requires a high level of mesh refinement to capture small movement of the boundary. XFEM allows the solution boundary to move independently of the mesh, permitting smooth representation of the evolution of the boundary as material is deposited. As a result, this new method will significantly improve the simulation accuracy while reducing the computational cost compared with existing approaches. Furthermore, this project will bring significant improvements for the design and AM process optimization in an iterative fashion among numerical simuation, parameter optimization, and experimental validation. Specifically, this project will support the protonic ceramic electrochemical cells (PCEC) stack development at INL by providing insights into the density, residual stress, thermal and mechanical properties of the PCEC manifold and interconnect, which is critical in enhancing the system lifetime and reducing the overall cost. Upon the success of proposed development and validation of the proposed approach, simulation will be applied to determine an optimal set of AM process parameters for the PCEC manifold production. This project will consolidate Idaho National Laboratory (INL)’s simulation capabilities provided by the MOOSE framework and the Valhalla AM simulation application, with encouraging expansion to emerging PCEC applications and AM technologies at INL and beyond.

97 MATHEMATICS AND COMPUTING↗

Development and Comparison of Purgeable Mercury Values in SRR Samples Measured by SRNL and Eurofins FGS

Savannah River Remediation (SRR) requested the development of mercury speciation capabilities at the Savannah River National Laboratory (SRNL) to support the Liquid Waste Operations at SRS. As part of that method development, SRR requested that SRNL Analytical Development (AD) compare their results with those obtained from their outside contract laboratory, Eurofins Frontier Global Sciences (FGS). This document reports on this method development work performed at SRNL as well as the comparative analyses conducted between the two laboratories. Development, optimization, and validation were undertaken at SRNL to produce a method for the species-specific analysis of dissolved gaseous mercury (DGM) in samples from SRR. It was determined that DGM present in SRR samples represented a combination of various soluble and volatile mercury species- and is thus referred to in this document as purgeable mercury. A method has been developed and analytically validated for the analysis of purgeable mercury in caustic nuclear waste. Using multivariate optimization, an efficient and sensitive method was produced. Total recovery of spiked purgeable mercury from deionized water was 99.4 ±10.6% and 91.6 ±14.9% from spiked high-activity tank waste samples. Calibration, by external calibration and standard addition calibration, maintained a mean accuracy of 101 ±5% and 99.9 ±1.7%, respectively. An interlaboratory comparison was performed between SRNL and Eurofins FGS in the co-analysis of three quarterly SRR Tank 50 samples. The statistical evaluations presented in this document support the following conclusions: for internal validation data, the two laboratories demonstrate biases of -4.3% for SRNL and 4.4% for Eurofins FGS, though neither are statistically significant at the 5% level. While the percent relative standard deviation for each of the laboratories is rather large (21.5% for SRNL and 18.4% for Eurofins FGS), there is no indication of a statistically significant difference between the two laboratories in the precision of these results. For the three sets of SRR quarterly Tank 50 samples, there is no indication of a statistically significant bias between the two laboratories at the 5% significance level. These results indicate that SRNL may provide data to the customer of statistically similar quality to Eurofins FGS.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A Machine Learning Method for the Forensics Attribution of Separated Plutonium

Plutonium (Pu) source attribution would be a powerful tool to support nuclear nonproliferation efforts. This capability to find the source of a Pu sample would act as a deterrent to smuggling efforts, and also help regulatory agencies verify declared nuclear activities. Work at Texas A&M University yielded a nuclear forensics methodology, which is capable of determining separated Pu’s reactor of origin, fuel burnup, and the time since irradiation (TSI)—three parameters of interest. The methodology used a set of ten intra-element isotopic ratios found in separated Pu, which was compared to a library of isotopic ratio values produced using neutronics simulations for reactors of interest. By calculating the probability that unknown Pu sample’s isotopic ratio set matched a set in the library, the methodology could predict the three parameters of interest of the sample. One shortcoming of this methodology was an inability to correctly attribute spoofed Pu, where Pu sourced from two different reactors or two different fuel burnup levels are mixed. A new methodology to rectify this vulnerability using machine learning (ML) technique is developed, instead of the maximum likelihood calculation previously used and the results are satisfactory. The ML approach leverages the existing simulated data for training the algorithm, but use them efficiently by only using intra-element isotope ratios that contribute to the attribution one of the three parameters at a time. Previously, all isotope ratios were used to attribute all three parameters together. The new methodology attributes the Pu parameters in three steps, one for each parameter, rather than resolving all of the three parameters simultaneously like the previous maximum likelihood approach. First, a support vector machine classifier with a set of seven isotopic ratios finds the reactor of origin and a set of regression models trained using gaussian process predicts the burnup with a different set of seven isotopes. Finally, TSI is calculated analytically using decay equations. Thus far, the new methodology is capable of attributing pure Pu samples and has been validated using experimental data. The next step will to be augment the classifier training data set with spoofed Pu data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗