Search NASA⌕ Search

SEARCH · Search NASA

Results for “predictive equivalence”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Auto-BLAST (AutoBLAST) [SWR-20-93]

Battery life modeling often involves a lot of manual parameter fitting and is not easy for users to adopt the model and use it. To reduce the difficulties for users to adopt battery lifetime models, an automatic battery lifetime modeling analysis and simulation tool suite, Auto-BLAST, has been developed. Auto-BLAST includes a lithium-loss-base life model, a battery electric model, and an algorithm which automatically fits all key parameters in the electric and life models using user provided data. The models and algorithm are packaged into two user-friendly GUIs, Auto-LifeMod, for easy battery life prognostic model fitting and Auto-LifeSim, for easy battery lifetime simulation. The lithium-loss-based life model adopts a similar model framework that models degradations using aging rate models, using battery cycling data to predict battery degradation and expected lifetime. The electric model is an equivalent circuit model which simulates battery voltage responses based on current/power demand profiles. The auto-fitting algorithm uses the user-input data to generate custom battery life model(s). Two GUIs, wrapping around the life model and the auto-fitting algorithm, provide friendly interfaces for users to generate a life model and use it for case study. One GUI requires summary data from accelerated battery life degradation tests as an input and produces battery life models predicting (a) capacity degradation, and (b) resistance growth of the battery. The GUI displays the electrical model response and the input experimental data against the life model predictions with fitted model parameters and fitting error. The GUI also generates an output file to save the fitted model parameters. The second GUI uses the life model from the first GUI and a user-defined battery cycling profile to predict battery lifetime degradation and expected lifetime. Predicted capacity and resistance vs. time are plotted in the GUI and saved to output file. There is a user guide for both GUIs

Mishra, Partha↗
DOE OSTI2024

Machine Learning and Data Science to Advance Laboratory Earthquake Prediction and Illuminate the Mechanics of Precursors to Failure

Earthquakes represent one of our greatest natural hazards and in recent years human induced seismicity is adding to the threat. Even a modest improvement in the ability to forecast devastating large earthquakes or smaller shallow events associated with fluid injection could save thousands of lives and billions of dollars. Current efforts to forecast earthquakes are limited by knowledge of earthquake physics and hampered by a lack of reliable lab or field observations. However, recent work has provided a critical opportunity for advancement. We have found: 1) clear and consistent precursors prior to earthquake-like failure in the laboratory and 2) that lab earthquakes can be predicted using machine learning (ML). These works show that stick-slip failure events –the lab equivalent of earthquakes– are preceded by a cascade of micro-failure events that radiate elastic energy in a manner that foretells catastrophic failure. Remarkably, ML predicts the fault zone stress state, the failure time and in some cases the magnitude of lab earthquakes. In addition, the observations include clear precursors to failure in the form of changes in fault zone properties prior to lab earthquakes. Precursors have been observed in previous laboratory studies but their origin is poorly understood and their possible connection to ML based earthquake prediction is unknown. The work conducted under our project has dramatically expanded these efforts. We have developed an integrated data science approach to illuminate the physics of earthquake precursors and lab earthquake prediction. Our work has accelerated the development of ML, artificial intelligence (AI), and related data science approaches by providing massive data sets that are tightly connected to critical scientific problems and by bringing together leading subject matter experts and data scientists. Earthquake physics involves phenomena that are far from equilibrium. Our work has leveraged data science methods to illuminate these phenomena and investigate how they relate to earthquake prediction. In addition to a large database with many types of labeled events that is available to everyone, our work has advanced the fundamental understanding of seismic forecasting, earthquake physics, and fault rheology

58 GEOSCIENCES↗
DOE OSTI2020

Automated fault detection of residential air-conditioning systems using thermostat drive cycles

Residential air conditioning equipment comprises a significant portion of the total energy consumption of a home. Unfortunately, air-conditioning systems can be susceptible to faulty operation either from installation errors or faults that accrue over the equipment’s lifetime. This paper presents a novel automated fault detection algorithm for residential air-conditioning systems that can alert the homeowner of the presence of these faults. The proposed algorithm utilizes only the home’s thermostat and outside air temperature to perform automated fault detection over the course of the equipment’s lifetime, including immediately after installation. The algorithm uses an extended Kalman filter approach to identify a three-resistor, two-capacitor (3R2C) electrical equivalent thermodynamic model. The identified 3R2C model is used to predict cooling times during a testing period comprising of a series of thermostat drive-cycle experiments. We tested the algorithm on an EnergyPlus™ model of a typical residential building in Orlando, Florida. Duct faults, indoor airflow faults, and refrigerant undercharge faults were introduced into the building model one at a time. The algorithm was able to accurately determine duct-leak faults, 40% airflow faults, 40% undercharge faults, and no-fault cases with an accuracy of 70%, 77%, 82%, and 87%, respectively.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Effects of elastic softening and helium accumulation kinetics on surface morphological evolution of plasma-facing tungsten

Based on a continuous-domain model, capable of accessing the spatiotemporal scales relevant to fuzz formation on the surface of plasma-facing component (PFC) tungsten, we report self-consistent simulation results that elucidate the effects of elastic softening and helium (He) accumulation kinetics on the surface morphological response of PFC tungsten. The model accounts for the softening of the elastic moduli in the near-surface region of PFC tungsten, including both thermal softening at high temperature and softening due to He accumulation upon He implantation. The dependence of the elastic moduli on the He content follows an exponential scaling relation predicted by molecular-dynamics simulations, while the He content in the near-surface region of PFC tungsten evolves according to a first-order saturation kinetics, consistent with experimental and simulation results reported in the literature. We establish that this elastic softening accelerates both nanotendril growth on the PFC surface and the onset of fuzz formation. We also explore the role of the rate of He accumulation to a saturation level in the near-surface region of irradiated tungsten in the onset of fuzz formation. For PFC tungsten surfaces such as W(110) where, under typical irradiation conditions, the characteristic time scale for stress-driven surface diffusion is comparable to the characteristic time scale for He accumulation, we find that accelerating the rate of He accumulation accelerates the growth rate of nanotendrils emanating from the surface. Additionally, we present a systematic parametric study of the PFC surface morphological response to explore its dependence on the He accumulation kinetics that is controlled by the irradiation conditions for low-energy implantation. Finally, we introduce an incubation time for nanotendril growth on the PFC surface, a concept equivalent to that of incubation fluence discussed in the literature, to predict and explain the minimum exposure time required to observe fuzz formation on PFC tungsten surfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Application of Cell Painting for chemical hazard evaluation in support of screening-level chemical assessments

‘Cell Painting’ is an imaging-based high-throughput phenotypic profiling (HTPP) method in which cultured cells are fluorescently labeled to visualize subcellular structures (i.e., nucleus, nucleoli, endoplasmic reticulum, cytoskeleton, Golgi apparatus / plasma membrane and mitochondria) and to quantify morphological changes in response to chemicals or other perturbagens. HTPP is a high-throughput and cost-effective bioactivity screening method that detects effects associated with many different molecular mechanisms in an untargeted manner, enabling rapid in vitro hazard assessment for thousands of chemicals. Here, 1201 chemicals from the ToxCast library were screened in concentration-response up to ~100 μM in human U-2 OS cells using HTPP. A phenotype altering concentration (PAC) was estimated for chemicals active in the tested range. PACs tended to be higher than lower bound potency values estimated from a broad collection of targeted high-throughput assays, but lower than the threshold for cytotoxicity. In vitro to in vivo extrapolation (IVIVE) was used to estimate administered equivalent doses (AEDs) based on PACs for comparison to human exposure predictions. AEDs for 18/412 chemicals overlapped with predicted human exposures. Phenotypic profile information was also leveraged to identify putative mechanisms of action and group chemicals. Of 58 known nuclear receptor modulators, only glucocorticoids and retinoids produced characteristic profiles; and both receptor types are expressed in U-2 OS cells. Thirteen chemicals with profile similarity to glucocorticoids were tested in a secondary screen and one chemical, pyrene, was confirmed by an orthogonal gene expression assay as a novel putative GR modulating chemical. Most active chemicals demonstrated profiles not associated with a known mechanism-of-action. However, many structurally related chemicals produced similar profiles, with exceptions such as diniconazole, whose profile differed from other active conazoles. Overall, the present study demonstrates how HTPP can be applied in screening-level chemical assessments through a series of examples and brief case studies.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Simulated X-Ray Radiographic Performance of a Bismuth-Loaded PVT Array

Recent material advancements in organic plastic scintillators enable marked increases in material detection efficiency, light yield, and pulse shape discrimination properties. These advances may resolve significant capability gaps for lowcost, portable, and durable dual-particle imaging (DPI) systems for nuclear safety, security, and safeguard purposes. One such material, a 21% bismuth-loaded polyvinyl toluene (BiPVT), is computationally evaluated as a small, pixelated radiographic array using Monte Carlo N-Particle (MCNP) and Zemax OpticStudio, and it is compared to identical evaluations of EJ-200 and EJ-256 arrays. MCNP software enables estimates of particle interaction and energy deposition, while OpticStudio computes optical light transport within each material. Here, computational estimates of spatial resolution and relative light collection at 370 kVp are found to agree with experimental results for both EJ-200 and EJ-256 arrays, thereby validating predictions of the same for the BiPVT array. As such, for equivalent exposures at 370 kVp, a BiPVT array may provide ~20× the light collection expected from EJ-200 and ~10× that expected from EJ-256. Similar comparisons of estimated light collection are also computed at 150 and 270 kVp, and these results suggest that BiPVT will provide significantly improved performance over EJ-200 and EJ-256 across all energies practical for portable X-ray radiography.

42 ENGINEERING↗
DOE OSTI2023

Hybrid Tandem Photovoltaics

Tandem solar cell structures are the only strategy demonstrated to surpass the detailed balance efficiency limit of high-quality single-junction solar cells. To continue to improve the efficiencies of cost-effective terrestrial solar power, hybrid tandems of dissimilar subcells are being considered by many around the world, especially designs that incorporate silicon solar cells as a bottom subcell. In this project, we studied a wide variety of tandem design possibilities including those with three-terminal (3T) and four-terminal (4T) configurations. The use of 3T and 4T designs could be useful for efficient and economical hybrid tandem designs that utilize the best available subcell materials such as emerging perovskite materials. Three-terminal configurations, in particular, have not been sufficiently studied previously. We have laid the foundational groundwork in this project for understanding the operation of 3T tandems: developing a taxonomy for naming, a methodology for measuring and interconnecting, and models for simply characterizing 3T tandems. Electrical and optical subcell coupling between the subcells was also measured and modeled. An important part of this work was the fabrication of novel example tandem structures, including 4T GaAs/Si, 3T GaInP/Si, 3T GaAs/Si, and 3T GaInP/GaAs devices. Using these high-quality tandem cells, we have been able to clearly demonstrate the achievability of high-efficiencies, and subtle physical effects such as photon recycling and luminescent coupling. We have developed and demonstrated essential building-block tools such as transparent conductive adhesives (TCA) and 3T silicon bottom cells with interdigitated back contacts (IBC) that can also be used in many other tandem designs. We have tested the reliability of these tools and devices under standardized testing and outdoor measurements. We have found 4T GaAs/Si tandems to be relatively straightforward to fabricate and robust in real-world outdoor conditions. While we have demonstrated working hybrid 3T III-V/TCA/Si IBC tandems, we experienced low yields even with our best process flows yet. Further work is still needed to improve the processing yield of these devices. We therefore also created tandem cells using an all-III-V 3T tandem process which was very robust with high yields, allowing for the creation of voltage-matched strings in many different configurations using 8 nearly identical 3T tandems. Using these robust 3T tandem examples, we were able measure and precisely characterize 3T tandem behaviors to predict their operation under changing spectrum and temperature. The optoelectronic equivalent-circuit model was shown to be very general and applicable to hybrid tandems, and encompassed the operation 3T Si IBC cells. This general model has been distributed to the public in as open-source Python-based software called PVcircuit. We have calculated the implications of these new tandem device designs on the real-world energy production and shown how the relative performance of different tandem configurations is situational and can be engineered using the tools developed here.

14 SOLAR ENERGY↗
DOE OSTI2025

Nondestructive Evaluation of Concrete: Elastic Property Imaging Through Full Waveform Inversion

Concrete is a major construction material worldwide and plays a crucial role in the nuclear industry. The elastic properties of concrete are prone to change and degrade while in service, as it is often subjected to extreme operational and environmental conditions. An accurate evaluation of concrete's elastic properties is thus essential to ensure structural integrity and safety. This is especially true for concrete in nuclear power plants, where irradiation effects significantly impact concrete mechanical properties. There are various methods to assess these properties, with ultrasound-based techniques showing high potential due to their nondestructive nature, cost-effectiveness, and safety. While several nondestructive evaluation methods exist, most rely on idealizations such as assuming homogeneous material and plane wavefronts. In this work, we address these issues by introducing an ultrasound-based nondestructive method aimed at reconstructing spatially varying images of concrete mechanical properties. By accurately modeling wave physics, including scattering and reflection, we overcome several of the aforementioned idealizations and aim to utilize the full waveform for imaging material properties through depth, resulting in more reliable images. Full waveform inversion (FWI) was first introduced by geophysicists to reconstruct subsurface elastic property images. The goal is to minimize the difference between simulated and recorded wavefield signals, often through gradient-based optimization algorithms. While FWI is primarily conducted using the acoustic approximation of the wave equation, few works focus on elastic FWI, where the goal is to reconstruct images of not only the pressure wave speed but also the shear wave speed and density (or their equivalents). This work explores the potential of using elastic FWI to predict concrete mechanical properties as an initial effort for a more accurate monitoring of concrete conditions in service. Reconstructing images of different elastic parameters enables more specificity and accurate condition diagnosis. This paper will detail this approach and provide examples demonstrating the effectiveness of elastic FWI in reconstructing comprehensive maps of concrete mechanical properties.

42 - ENGINEERING↗
DOE OSTI2021

Predicting transcriptional responses to cold stress across plant species

Although genome-sequence assemblies are available for a growing number of plant species, gene-expression responses to stimuli have been cataloged for only a subset of these species. Many genes show altered transcription patterns in response to abiotic stresses. However, orthologous genes in related species often exhibit different responses to a given stress. Accordingly, data on the regulation of gene expression in one species are not reliable predictors of orthologous gene responses in a related species. Here, we trained a supervised classification model to identify genes that transcriptionally respond to cold stress. A model trained with only features calculated directly from genome assemblies exhibited only modest decreases in performance relative to models trained by using genomic, chromatin, and evolution/diversity features. Models trained with data from one species successfully predicted which genes would respond to cold stress in other related species. Cross-species predictions remained accurate when training was performed in cold-sensitive species and predictions were performed in cold-tolerant species and vice versa. Models trained with data on gene expression in multiple species provided at least equivalent performance to models trained and tested in a single species and outperformed single-species models in cross-species prediction. These results suggest that classifiers trained on stress data from well-studied species may suffice for predicting gene-expression patterns in related, less-studied species with sequenced genomes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Influence of tower shadow on downwind flexible rotors: Field tests and simulations

Abstract As wind turbine rotors become larger, the blades become more flexible, requiring extra stiffness and cost to avoid the risk of tower strike. Wind turbines in a downwind configuration have a reduced risk of tower strike because the rotor thrust acts away from the tower. However, downwind blades pass through the wake of the tower, and the resulting load variation may contribute to blade fatigue. To date, there have been no field tests to quantify this tower shadow effect on unsteady blade moments. The present study reports on the first field testing of a flexible, downwind, coned rotor and compares the experimental data against simulations run in OpenFAST. The tower shadow effect is simulated using the conventional Powles model and a new Eames model (developed herein), which includes the influence of upstream turbulence. Both models reasonably predict the blade root out‐of‐plane bending moment data and the tower shadow dip magnitude when compared to field test data in Region 3. Tower shadow was found to increase the short‐term Damage Equivalent Loads (DELs) by less than 10% compared to other effects (gravity, shear, and turbulence), and the predictions were consistent with experiments. These results indicate that the tower shadow effect can be reasonably modeled with the simpler Powles model and that the tower shadow effect can be small compared to the effect of turbulence. However, long‐term fatigue due to tower shadow should be included in detailed structural analysis and design of the rotor and the tower.

Simpson, Juliet G.↗
DOE OSTI2024

Inverse prediction of PuO2 processing conditions using Bayesian seemingly unrelated regression with functional data

Over the past decade, a variety of innovative methodologies have been developed to better characterize the relationships between processing conditions and the physical, morphological, and chemical features of special nuclear material (SNM). Different processing conditions generate SNM products with different features, which are known as “signatures” because they are indicative of the processing conditions used to produce the material. These signatures can potentially allow a forensic analyst to determine which processes were used to produce the SNM and make inferences about where the material originated. This article investigates a statistical technique for relating processing conditions to the morphological features of PuO 2 particles. We develop a Bayesian implementation of seemingly unrelated regression (SUR) to inverse-predict unknown PuO 2 processing conditions from known PuO 2 features. Model results from simulated data demonstrate the usefulness of the technique. Applied to empirical data from a bench-scale experiment specifically designed with inverse prediction in mind, our model successfully predicts nitric acid concentration, while results for Pu concentration and precipitation temperature were equivalent to a simple mean model. Our technique compliments other recent methodologies developed for forensic analysis of nuclear material and can be generalized across the field of chemometrics for application to other materials.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Improbability of Nuclear Criticality in Compacted Criticality Control Overpacks after Room Closure by Salt Creep at Waste Isolation Pilot Plant

Based on the rationale presented, nuclear criticality is improbable after salt creep causes compaction of criticality control overpacks (CCOs) disposed at the Waste Isolation Pilot Plant, an operating repository in bedded salt for the disposal of transuranic (TRU) waste from atomic energy defense activities. For most TRU waste, the possibility of post-closure criticality is exceedingly small either because the salt neutronically isolates TRU waste canisters or because closure of a disposal room from salt creep does not sufficiently compact the low mass of fissile material. The criticality potential has been updated here because of the introduction of CCOs, which may dispose up to 380 fissile gram equivalent plutonium-239 in each container. The criticality potential is evaluated through high-fidelity geomechanical modeling of a disposal room filled with CCOs during two representative conditions: (1) large salt block fall, and (2) gradual salt compaction (without brine seepage and subsequent gas generation to permit maximum room closure). Geomechanical models of rock fall demonstrate three tiers of CCOs are not greatly disrupted. Geomechanical models of gradual room closure from salt creep predict irregular arrays of closely packed CCOs after 1000 years, when room closure has asymptotically approached maximum compaction. Criticality models of spheres and cylinders of 380 fissile gram equivalent of plutonium (as oxide) at the predicted irregular spacing demonstrate that an array of CCOs is not critical when surrounded by salt and magnesium oxide, provided the amount of hydrogenous material shipped in the CCO (usually water and plastics) is controlled or boron carbide (a neutron poison) is mixed with the fissile contents.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

Improbability of Post-Closure Criticality in Compacted Criticality Control Overpacks after Room Closure at Waste Isolation Pilot Plant

As part of its periodic re-certification of the Waste Isolation Pilot Plant (WIPP), an operating repository in bedded salt for the disposal of transuranic (TRU) waste from atomic energy defense activities, the United States Environmental Protection Agency expects a re-evaluation of features, events, and processes, such as post-closure nuclear criticality. Although salt creep beneficially encapsulates the TRU waste in the closed WIPP repository, the spacing between an array of waste packages is disrupted as the salt creep closes disposal rooms and containers lose structural integrity. For most TRU waste, the possibility of post-closure criticality is exceedingly small either because the salt neutronically isolates TRU waste canisters or because closure of a disposal room from salt creep does not sufficiently compact the low mass of fissile material. The criticality evaluation was updated, however, because of the introduction of criticality control overpack (CCO) containers, which may dispose up to 380 fissile gram equivalent plutonium-239 in each container. The criticality potential is evaluated through high-fidelity geomechanical modeling of a disposal room filled with CCO containers during two representative conditions: (1) large salt block fall, and (2) gradual disposal room closure from salt creep. Geomechanical models of roof fall demonstrate three tiers of CCO containers are not greatly disrupted. Geomechanical models of gradual room closure from salt creep (without brine seepage and subsequent gas generation to permit maximum room closure) were used to predict irregular arrays of closely packed CCOs after 1000 years, when room closure has asymptotically approached maximum compaction. Models of spheres or cylinders with 380 fissile gram equivalent of plutonium (as oxide) at the predicted irregular compacted spacing demonstrate that an array of CCO containers is not critical when surrounded by salt and magnesium oxide, provided the mass of hydrogenous material shipped in CCO containers (usually plastics) is controlled or boron carbide (a neutron poison) is mixed with the fissile contents.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Using a physics-informed neural network and fault zone acoustic monitoring to predict lab earthquakes

Abstract Predicting failure in solids has broad applications including earthquake prediction which remains an unattainable goal. However, recent machine learning work shows that laboratory earthquakes can be predicted using micro-failure events and temporal evolution of fault zone elastic properties. Remarkably, these results come from purely data-driven models trained with large datasets. Such data are equivalent to centuries of fault motion rendering application to tectonic faulting unclear. In addition, the underlying physics of such predictions is poorly understood. Here, we address scalability using a novel Physics-Informed Neural Network (PINN). Our model encodes fault physics in the deep learning loss function using time-lapse ultrasonic data. PINN models outperform data-driven models and significantly improve transfer learning for small training datasets and conditions outside those used in training. Our work suggests that PINN offers a promising path for machine learning-based failure prediction and, ultimately for improving our understanding of earthquake physics and prediction.

42 ENGINEERING↗
DOE OSTI2025

Global Evaluation of Process Conditions and Wave Modes in a Rotating Detonation Engine

Rotating Detonation Engines (RDEs) show significant promise for enhancing the efficiency of gas turbine engines while maintaining low 𝑁𝑂𝑥 emissions. This work investigates the predictability of wave modes in a water-cooled RDE under varying operational conditions. Experimental data comprising over 6,700 samples was collected, including parameters such as flow rates, temperatures, pressures, and equivalence ratios. A machine learning approach using the XGBoost library was used to build a multi-class classifier, predicting wave modes based on these inputs. The model achieved a high accuracy of 97%, demonstrating that wave modes are not random but deterministic based on the process conditions. SHAP analysis was used to identify the most influential parameters affecting wave mode prediction. The results show that for the water-cooled NETL RDE, wave mode is determinant and predictable based on the process parameters.

Weber, Justin [NETL] (ORCID:0000000218487035)↗
DOE OSTI2021

Hybrid Large Eddy Simulation/Reynolds-Averaged Navier–Stokes Analysis of a Premixed Ethylene-Fueled Dual-Mode Scramjet Combustor

Hydrocarbon fuels offer optimal high energy per volume for scramjet applications for sustained hypersonic flight but require additional residence time due to slower ignition delays (compared with hydrogen fuel). The injection of ethylene at the start of the isolator of a dual-mode scramjet combustor, operating in ramjet mode, allows sufficient mixing to achieve efficient premixed turbulent combustion. A cavity flameholder anchors the flame, supplying sufficient radicals to sustain a stable flame in the high-speed environment. This work investigates flame structure and stabilization limits of a new configuration with a scaled-down cavity embedded in the flow path of the combustor with a strut and insert. The cavity is reduced in size by one-third to enable eventual direct numerical simulations of the flame stabilization process. This work, however, focuses on modeling the full isolator/combustor geometry using a hybrid large eddy simulation/Reynolds-averaged Navier–Stokes simulation strategy. Additionally, particle image velocimetry, planar laser-induced fluorescence, coherent anti-Stokes Raman spectroscopy, and pressure measurements are compared with numerical predictions to analyze and characterize the conditions within the combustor, including flame structure, flow velocities, species composition, and wall pressure. An adjustable air throttle in the extender is able to control the placement of the shock train in the isolator and maintain a stable flame at various equivalence ratios. The simulations show reasonably good agreement with the experimental scalar and velocity data and predict a flame angle consistent with premixed turbulent flame-speed correlations.

42 ENGINEERING↗
DOE OSTI2022

FutureTense

Protective vaccines and reliable diagnostics are essential tools for controlling viral diseases. However, the efficacy of these tools can be diminished by mutations in viral genomes. The delay between the emergence of new viral strains and the redesign of vaccines and diagnostics allows for continued viral transmission. Is it possible to address this challenge by computationally predicting viral genome sequence evolution? Can we “future-proof” vaccines and diagnostics by targeting both current and anticipated future sequence variants? While predicting viral evolution is still an unsolved, “grand challenge” problem in biology, the large, and rapidly growing, number of SARS-CoV-2 genome sequences provide an opportunity to quantify the ability of machine learning to predict viral genome sequence evolution. Towards this end, we have developed a simple computational model for predicting viral evolution at the level of individual nucleotides. The key metric for quantifying the per-base, prediction accuracy for viral evolution is the Mann-Whitney U statistic (or, equivalently, the area under the receiver operator curve). Since the Mann-Whitney U statistic is not a differentiable function, existing deep leaning packages (like Pytorch and Keras/TensorFlow) are not useful, as they require that the accuracy metric/objective function be analytically differentiable with respect to the model parameters. To overcome this challenge, we have implemented custom software, “FutureTense”, that can train a machine learning model by maximizing the non-differentiable Mann-Whitney U statistic. This software trains a machine learning model by exploring along the direction of the discrete gradient of the Mann-Whitney U statistic in the model parameter space. Parallel computing and genome sequence-specific optimizations are used to accelerate model training. The resulting machine learning model learns the observed high C->U mutation rates in the SARS-CoV-2 genome (which are potentially induced by host defenses) and provides prediction accuracies that are significantly better than one would expect from random chance. While predicting viral evolution is still quite far from a solved problem, the surprising performance of this simple model gives hope that the accuracy of predicting viral genome evolution can be further increased by more sophisticated approaches.

Gans, Jason↗
DOE OSTI2022

Accelerated radiation transport modeling techniques for pencil beam computed tomography using gamma rays

Monte Carlo radiation transport modeling studies were performed for a compact, and high-resolution gamma-ray computed tomography system designed for imaging irradiated nuclear fuel. The system comprises a 60 Co source – chosen for its highly penetrating 1173 keV and 1332 keV gamma rays – a pair of high-aspect-ratio pencil beam collimators, and an inorganic scintillator detector. Here two acceleration methods are proposed to rapidly model a transmission type gamma-ray tomography system. The first, a variance reduction technique, is based on performing Monte Carlo simulations with a monodirectionally-biased source, sampled from a characteristic sub-volume of the full source volume. The second acceleration method is based on the deterministic calculations using the Beer–Lambert law and detector response characteristics. Comparison of simulations using acceleration approaches with analog simulations of the fully isotropic, full-volume equivalent, show that the Monte Carlo variance reduction technique gives quantitatively accurate predictions for large collimator aspect ratios while the deterministic calculations are semi-quantitative but converge close to the correct result as the collimator aspect ratio increases. As such, these techniques can be used to reduce the computational cost in generating simulated radiographs and tomographs by several orders of magnitude. Experimental validation efforts are currently underway and will be demonstrated in future work.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗