Search NASA⌕ Search

SEARCH · Search NASA

Results for “Multiscale analysis”

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

Distinguishing Mesoscale Polar Order (Unidirectional vs Bidirectional) of Cellulose Microfibrils in Plant Cell Walls Using Sum Frequency Generation Spectroscopy

Cellulose in plant cell walls are synthesized as crystalline microfibrils with diameters of 3–4 nm and lengths of around 1–10 μm. These microfibrils are known to be the backbone of cell walls, and their multiscale three-dimensional organization plays a critical role in cell wall functions including plant growth and recalcitrance to degradation. The mesoscale organization of microfibrils over a 1–100 nm range in cell walls is challenging to resolve because most characterization techniques investigating this length scale suffer from low spatial resolution, sample preparation artifacts, or inaccessibility of specific cell types. Here, we report a sum frequency generation (SFG) study determining the mesoscale polarity of cellulose microfibrils in intact plant cell walls. SFG is a nonlinear optical spectroscopy technique sensitive to the molecular-to-mesoscale order of noncentrosymmetric domains in amorphous matrices. However, the quantitative theoretical model to unravel the effect of polarity in packing of noncentrosymmetric domains on SFG spectral features has remained unresolved. In this work, we show how the phase synchronization principle of the SFG process is used to predict the relative intensities of vibrational modes with different polar angles from the noncentrosymmetric domain axis. Applying this model calculation for the first time and employing SFG microscopy, we found that cellulose microfibrils in certain xylem cell walls are deposited unidirectionally (or biased in one direction) instead of the bidirectional polarity which was believed to be dominant in plant cell walls from volume-averaged characterizations of macroscopic samples. Here, with this advancement in SFG analysis, one can now determine the relative polarity of noncentrosymmetric domains such as crystalline biopolymers interspersed in amorphous polymer matrices, which will open opportunities to study new questions that have not been conceived in the past.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Computer-aided detection of clustered microcalcifications in multiscale bilateral filtering regularized reconstructed digital breast tomosynthesis volume

Purpose: Develop a computer-aided detection (CADe) system for clustered microcalcifications in digital breast tomosynthesis (DBT) volume enhanced with multiscale bilateral filtering (MSBF) regularization. Methods: With Institutional Review Board approval and written informed consent, two-view DBT of 154 breasts, of which 116 had biopsy-proven microcalcification (MC) clusters and 38 were free of MCs, was imaged with a General Electric GEN2 prototype DBT system. The DBT volumes were reconstructed with MSBF-regularized simultaneous algebraic reconstruction technique (SART) that was designed to enhance MCs and reduce background noise while preserving the quality of other tissue structures. The contrast-to-noise ratio (CNR) of MCs was further improved with enhancement-modulated calcification response (EMCR) preprocessing, which combined multiscale Hessian response to enhance MCs by shape and bandpass filtering to remove the low-frequency structured background. MC candidates were then located in the EMCR volume using iterative thresholding and segmented by adaptive region growing. Two sets of potential MC objects, cluster centroid objects and MC seed objects, were generated and the CNR of each object was calculated. The number of candidates in each set was controlled based on the breast volume. Dynamic clustering around the centroid objects grouped the MC candidates to form clusters. Adaptive criteria were designed to reduce false positive (FP) clusters based on the size, CNR values and the number of MCs in the cluster, cluster shape, and cluster based maximum intensity projection. Free-response receiver operating characteristic (FROC) and jackknife alternative FROC (JAFROC) analyses were used to assess the performance and compare with that of a previous study. Results: Unpaired two-tailedt-test showed a significant increase (p < 0.0001) in the ratio of CNRs for MCs with and without MSBF regularization compared to similar ratios for FPs. For view-based detection, a sensitivity of 85% was achieved at an FP rate of 2.16 per DBT volume. For case-based detection, a sensitivity of 85% was achieved at an FP rate of 0.85 per DBT volume. JAFROC analysis showed a significant improvement in the performance of the current CADe system compared to that of our previous system (p = 0.003). Conclusions: MBSF regularized SART reconstruction enhances MCs. The enhancement in the signals, in combination with properly designed adaptive threshold criteria, effective MC feature analysis, and false positive reduction techniques, leads to a significant improvement in the detection of clustered MCs in DBT.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2021

Tissue scale agent-based simulation of premalignant progressions in Barrett’s esophagus

Barrett’s esophagus (BE) is a benign condition of the distal esophagus that initiates a multistage pathway to esophageal adenocarcinoma (EAC). Short of frequent intrusive (and costly) surveillance, effective screening for neoplasia in BE populations is yet to be established since progressors are rare and virtually undetectable without routine biopsies, which often sample only a small portion of the BE tissue. As a result, reliable estimation of the true prevalence of dysplasia in a BE population and evidence-based optimization of screening for at-risk individuals is challenging. Data-driven microsimulations, i.e., model-generated instances of disease history in a predefined virtual population, have found utility in the EAC screening literature as low-overhead alternatives to real-world hypothesis testing of optimal interventions for dysplasia. Despite the successes, computational limitations, paucity of knowledge and data on Barrett’s dysplasia, and the complexities of disease progression as a multiscale multiphysics process have hindered the treatment of disease progression in BE as a spatial process. Agent-based modeling of nucleation and proliferation processes in dysplasia warrants exploration in this context as an approximation that operates at a trade-off between computational tractability and precise representation of the composition and physics of the substrate (tissue). In this study, we describe spatially resolved simulations of premalignant progression toward EAC in a coarse-grained model of Barrett’s tissue that resolves the metaplastic tissue at a length scale of 0.42 mm (~3300 crypts/mm 2 ). Finally, the model is calibrated to reproduce historical high-grade dysplasia prevalence when model-generated patients are screened using the Seattle protocol.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Fractal Characterization of Multimodal, Multiscale Images of Shale Rock Fracture Networks

An array of multimodal and multiscale images of fractured shale rock samples and analogs was collected with the aim of improving the numerical representation of fracture networks. 2D and 3D reconstructions of fractures and matrices span from 10 –6 to 100 m. The origin of the fracture networks ranged from natural to thermal maturation to hydraulically induced maturation. Images were segmented to improve fracture identification. Then, the fractal dimension and length distribution of the fracture networks were calculated for each image dataset. The resulting network connectivity and scaling associations are discussed at length on the basis of scale, sample and order of magnitude. Fracture network origin plays an important role in the resulting fracture systems and their scaling. Rock analogs are also evaluated using these descriptive tools and are found to be faithful depictions of maturation-induced fracture networks.

58 GEOSCIENCES↗
DOE OSTI2024

United States Multi-Sector Dynamics land use and land cover base maps to support Human-Earth System Modeling

Datasets are land use and land cover (LULC) rasterized base maps at 30-m resolution for the conterminous United States (CONUS) for the years 2008, 2011, 2016, and 2019. Separate base maps are provided where LULC classifications are thematically congruent with Community Land Model (CLM), Land Use Harmonization (LUH2), and Global Change Analysis Model (GCAM), and a detailed decomposition of all combined land classes into a Multisector Dynamics (MSD) LULC product. Base maps were developed using empirically derived satellite (National Land Cover Dataset, MODIS) and combined observation datasets (Crop Data Layer, Protected Areas Database) and represent the most up-to-date accurate information on LULC in the CONUS. The four datasets encompass four different landcover classification systems: MSD Layers - The raw landcover classes obtained from reclassifying NLCD and USDA Crop data layers into a respective landcover class GCAM Layers - The MSD classes mosaiced, reclassified, and combined into the respective GCAM landcover classes CLM Layers - Similar process to GCAM layers, but mosaiced, reclassified, and combined MSD layers to their respective PFT classes LUH2 Layers - Similar process to both GCAM and CLM Layers, but mosaiced, reclassified and combined the MSD layers to align with the respective states

Food↗
DOE OSTI2024

Polysulfide-incompatible additive suppresses spatial reaction heterogeneity of Li-S batteries

Rational electrolyte engineering for practical pouch cells remains elusive because the correlation between the cathode/solid-electrolyte interphase layer and cell-level reaction behavior is poorly understood. Here, by combining multiscale characterization and computational modeling, we show that—counter to the conventional perception of polysulfide-incompatible additives—the spontaneous reaction of sparingly solvated polysulfides with Lewis acid additives (LAAs) can induce in situ formation of a homogeneous interphase on thick and tortuous S cathode. Multiscale synchrotron X-ray characterization consistently affirms that such interface design could effectively eliminate the notorious problems of polysulfide shuttle and lithium corrosion and, more importantly, provide an interconnected “ion transport highway” to alleviate the uneven ion transport within the tortuous S cathode. Hence, this design dramatically reduces the reaction heterogeneity of lithium-sulfur (Li-S) pouch cells under lean electrolyte conditions. Further, this work resolves controversy around the role of polysulfide-incompatible additives in high-energy Li-S pouch cells and highlights the importance of suppressing reaction heterogeneity for practical batteries.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Understanding the dynamic evolution of atomically dispersed Cu catalyst for CO 2 electrochemical conversion using integrated XANES analysis and mechanistic studies

Direct electrochemical conversion of CO 2 to ethanol (CH 3 CH 2 OH) offers a promising strategy to lower CO 2 emission while storing energy from renewable electricity. Our recent study reported a carbon-supported atomically dispersed Cu catalyst that achieved the highest reported selectivity for CH 3 CH 2 OH formation (91%) at a relatively low potential (-0.6 V), however, the active site structure that is responsible for such high activity and selectivity has yet to be understood. Here, in this paper, we demonstrate a computational investigation combining X-ray absorption near edge structure (XANES) simulations and a mechanistic study via density functional theory (DFT) to understand the catalyst structures of this Cu catalyst during electrocatalysis and the corresponding reaction mechanisms of the key products. An integrated computational and experimental XANES analysis depicted the dynamic evolution of the catalytic site during electrocatalysis. The as-prepared, atomically dispersed Cu catalyst aggregates and forms metallic clusters/nanoparticles under electrochemical condition, which then break down to smaller oxidized clusters after electrocatalysis. The formed Cu clusters/nanoparticles showed distinct catalytic activity and selectivity as a function of particle size based on the mechanistic investigation using DFT, which is consistent with experimental observations for catalyst samples with different Cu loadings. This comprehensive study which combines experimental and computational XANES investigation, mechanistic study via DFT calculations, and experimental performance of the catalysts, provides unprecedented dynamic and mechanistic insights into the supported atomically dispersed metal catalysts for CO 2 reduction. Such strategy and details gained can further guide discovery of novel catalyst materials for CO 2 electrochemical reduction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Multiscale Electricity Modeling for Evaluating Carbon Capture and Sequestration Technologies (Final Report)

Carbon capture and sequestration (CCS) technologies that can operate with a high degree of operating flexibility could provide necessary electric grid flexibility in a system with high shares of variable renewables. This project examines the deployment and dispatch potential of twelve unique flexible CCS (FLECCS) technologies that encompass post-combustion carbon dioxide (CO 2 ) capture designs, concepts using a storage media to enable energy arbitrage, and hybrid processes that integrate CCS with direct air capture (DAC) for flexibility with net zero or negative CO 2 emissions. FLECCS technology potential is explored with a multi-model, multi-scale framework including the Regional Energy Deployment System (ReEDS) electric sector capacity expansion model (CEM) and the PLEXOS production cost model (PCM). Innovative methods were developed to represent FLECCS technology operating modes, performance, and cost in the two models. ReEDS was then used to simulate nine scenarios for each FLECCS technology, three CO 2 emissions price futures reaching $\$$150, $\$$225, and $\$$300/tCO 2 in 2050; and three scenarios for FLECCS technology deployment favorability relative to competing technologies. For each CO 2 price and reference FLECCS favorability, the 2050 infrastructures from ReEDS model are downscaled and implemented in PLEXOS to examine hourly dispatch under detailed operational constraints that are not included in ReEDS. FLECCS technologies exhibited a wide range of deployment potential ranging from none to several hundred gigawatts of capacity, with outcomes highly sensitive to input cost and performance parameters that are inherently highly uncertain. When deployed, FLECCS tended to displace a combination of wind, solar, and natural gas-based technologies rather than supporting increased renewable deployment. As a result, CO 2 emissions reductions facilitated by FLECCS deployment tended to come with higher overall system costs and electricity prices. When economically competitive, FLECCS technologies can contribute significant flexible generation and firm capacity to the grid, but continued technology development and an expanded analytical scope are necessary to fully understand FLECCS deployment potential its impact on the electric power sector. Follow-on analysis incorporating captured CO 2 tax credit value from the Inflation Reduction Act (IRA) and other potential policy scenarios could be particularly valuable, as this policy can substantially change the relative competitiveness of FLECCS technologies.

03 NATURAL GAS↗
DOE OSTI2024

Determining the Crosslinking and Degradation Reaction Kinetics in Photovoltaic Encapsulants Using Accelerated Aging: Preprint

Degradation of photovoltaic (PV) module encapsulant mechanical characteristics that lead to embrittlement and delamination remains a cause of failure in solar installations. A multiscale reliability model based on detailed molecular degradation reaction kinetics was previously published, connecting the encapsulant mechanical properties (elastic modulus, yield strength, and adhesion energy) to the degraded molecular structure and interfacial bond density to adjacent solar cell and glass substrates. The model, developed primarily for poly(ethylene-co-vinyl acetate) (EVA) encapsulants, remains to be experimentally validated. Determining the degradation and crosslinking kinetics of alternative encapsulants, such as polyolefin elastomer (POE) and EVA/POE/EVA composites (EPE), can further generalize the model. Importantly, the activation energy for crosslinking of fully cured PV encapsulant products from temperature or UV is presently unknown or unavailable. In this work, we subject EVA, POE, and EPE encapsulants to a set of accelerated aging conditions, varying the temperature, UV intensity, and relative humidity. We use DSC, FTIR, and Soxhlet extraction (gel content) to characterize the encapsulants' changing molecular structure. This allows for determination of the photochemical degradation and crosslinking kinetics of the encapsulants. Preliminary results show an increase in gel content (crosslinking) and a decrease in crystallinity of EVA, POE, and EPE encapsulants under hot-aerobic (90 deg C, 22% RH) and hot-anaerobic (90oC, sealed in N2 air) accelerated aging, even in the absence of UV and crosslinking initiating agents. For a least squares regression with an assumed first-order crosslinking kinetics, the crosslinking rate parameters were computed for the three encapsulants under hot-aerobic and hot-anaerobic aging conditions. FTIR analysis showed insignificant encapsulant degradation for hot-aerobic and hot-anaerobic aging, demonstrating the critical role of UV and moisture in accelerating degradation. Aging conditions with UV exposure and elevated humidity are presently in progress.

adhesion↗
DOE OSTI2021

Cross-scale dynamics and the evolutionary emergence of infectious diseases

When emerging pathogens encounter new host species for which they are poorly adapted, they must evolve to escape extinction. Pathogens experience selection on traits at multiple scales, including replication rates within host individuals and transmissibility between hosts. We analyze a stochastic model linking pathogen growth and competition within individuals to transmission between individuals. Our analysis reveals a new factor, the cross-scale reproductive number of a mutant virion, that quantifies how quickly mutant strains increase in frequency when they initially appear in the infected host population. This cross-scale reproductive number combines with viral mutation rates, single-strain reproductive numbers, and transmission bottleneck width to determine the likelihood of evolutionary emergence, and whether evolution occurs swiftly or gradually within chains of transmission. We find that wider transmission bottlenecks facilitate emergence of pathogens with short-term infections, but hinder emergence of pathogens exhibiting cross-scale selective conflict and long-term infections. Furthermore, our results provide a framework to advance the integration of laboratory, clinical, and field data in the context of evolutionary theory, laying the foundation for a new generation of evidence-based risk assessment of emergence threats.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

A multiscale model of the action of a capsid assembly modulator for the treatment of chronic hepatitis B

Chronic hepatitis B virus (HBV) infection is strongly associated with increased risk of liver cancer and cirrhosis. While existing treatments effectively inhibit the HBV life cycle, viral rebound frequently occurs following treatment interruption. Consequently, functional cure rates of chronic HBV infection remain low and there is increased interest in a novel treatment modality, capsid assembly modulators (CAMs). Here, we develop a multiscale mathematical model of CAM treatment in chronic HBV infection. By fitting the model to participant data from a phase I trial of the first-generation CAM vebicorvir, we estimate the drug’s dose-dependent effectiveness and identify the physiological mechanisms that drive the observed biphasic decline in HBV DNA and RNA, and mechanistic differences between HBeAg-positive and negative infection. Finally, we demonstrate analytically and numerically that the relative change of HBV RNA more accurately reflects the antiviral effectiveness of a CAM than the relative change in HBV DNA.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Simulations of grain growth in tungsten armor materials under ARC plasma edge operation conditions using an integrated plasma-edge/materials model

An integrated model of grain growth deuterium-exposed tungsten polycrystals, consisting of a two-dimensional vertex dynamics model fitted to atomistic data, has been developed to assess the grain growth kinetics of deuterium-exposed polycrystalline tungsten (W). The model tracks the motion of grain boundaries under the effect of driving forces stemming from grain boundary curvature and differential deuterium concentration accumulation. Here, we apply the model to experimentally synthesized tungsten polycrystals under deuterium-saturated conditions relevant to the ARC concept design. The results indicate rapid grain growth kinetics in the near-surface region adjacent to the plasma, where the temperature reaches 1400 K, whereas the microstructure remains stable deeper in the material with the lower temperature of 1000 K. The combined modeling and analysis further reveal that monolithic tungsten produced via conventional fabrication routes is highly susceptible to grain coarsening at temperatures exceeding 1000 K, largely independent of the magnitude of the applied driving force. Moreover, the accumulation of deuterium near grain boundaries has a pronounced inhibitory effect on grain boundary migration. High-angle grain boundaries ( > 50°) contribute more significantly to the overall grain growth process.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Non-autoregressive time-series methods for stable parametric reduced-order models

Advection-dominated dynamical systems, characterized by partial differential equations, are found in applications ranging from weather forecasting to engineering design where accuracy and robustness are crucial. There has been significant interest in the use of techniques borrowed from machine learning to reduce the computational expense and/or improve the accuracy of predictions for these systems. These rely on the identification of a basis that reduces the dimensionality of the problem and the subsequent use of time series and sequential learning methods to forecast the evolution of the reduced state. Often, however, machine-learned predictions after reduced-basis projection are plagued by issues of stability stemming from incomplete capture of multiscale processes as well as due to error growth for long forecast durations. To address these issues, we have developed a non-autoregressive time series approach for predicting linear reduced-basis time histories of forward models. In particular, we demonstrate that non-autoregressive counterparts of sequential learning methods such as long short-term memory (LSTM) considerably improve the stability of machine-learned reduced-order models. Further, we evaluate our approach on the inviscid shallow water equations and show that a non-autoregressive variant of the standard LSTM approach that is bidirectional in the principal component directions obtains the best accuracy for recreating the nonlinear dynamics of partial observations. Moreover-and critical for many applications of these surrogates-inference times are reduced by three orders of magnitude using our approach, compared with both the equation-based Galerkin projection method and the standard LSTM approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A multiresolution adaptive wavelet method for nonlinear partial differential equations

We report the multiscale complexity of modern problems in computational science and engineering can prohibit the use of traditional numerical methods in multi-dimensional simulations. Therefore, novel algorithms are required in these situations to solve partial differential equations (PDEs) with features evolving on a wide range of spatial and temporal scales. To meet these challenges, we present a multiresolution wavelet algorithm to solve PDEs with significant data compression and explicit error control. We discretize in space by projecting fields and spatial derivative operators onto wavelet basis functions. We provide error estimates for the wavelet representation of fields and their derivatives. Then, our estimates are used to construct a sparse multiresolution discretization which guarantees the prescribed accuracy. Additionally, we embed a predictor-corrector procedure within the temporal integration to dynamically adapt the computational grid and maintain the accuracy of the solution of the PDE as it evolves. We present examples to highlight the accuracy and adaptivity of our approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Final Technical Report -- Bridging the PSI Knowledge Gap: A Multiscale Approach

The Plasma Surface Interactions (PSI) Science Center formed by the grant undertook a multidisciplinary set of studies on the complex interface between the plasma and solid states of matter. The strategy of the center was to combine and integrate the experimental, diagnostic and modeling toolkits from multiple institutions towards specific PSI problems. In this way the Center could tackle integrated science issues which were not addressable by single institutions, as well as evolve the underlying science of the PSI in a more general way than just for fusion applications. The overall strategy proved very successful. The research result and highlights of the MIT portion of the Center are primarily described. A particular highlight is the study of tungsten nano-tendril growth in the presence of helium plasmas. The Center research provided valuable new insights to the mechanisms controlling the nano-tendrils by developing coupled modeling and in situ diagnostic methods which could be directly compared. For example, the role of helium accumulation in tungsten distortion in the surface was followed with unique in situ helium concentration diagnostics developed. These depth-profiled, time-resolved helium concentration measurements continue to challenge the numerical models of nano-tendrils. The Center team also combined its expertise on tungsten nano-tendrils to demonstrate for the first time the growth of the tendrils in a fusion environment on the Alcator C-Mod fusion experiment, thus having significant impact on the broader fusion research effort. A new form of isolated nano-tendril “columns” were identified which are now being used to understand the underlying mechanisms controlling the tendril growth. The Center also advanced PSI science on a broader front with a particular emphasis on developing a wide range of in situ PSI diagnostic tools at the DIONISOS facility at MIT. For example the strong suppression of sputtering by the certain combination of light-species plasmas and metals was experimentally studied with independent measurement methods across the Center. This surprising result challenges the universal use of the binary-collision approximation in sputtering predictions and continues to be the subject of study. In order to address this issue MIT developed a new in situ erosion measurement technique based on ion beam analysis which can be used at elevated material temperatures. This exciting new technique is now being used to study material erosion in high performance plasma thrusters for space exploration and is being adopted to fusion experimental devices. This is an indicator of the positive synergies that arise from such a Center, with the research having impact beyond the initial area of study. The Center also served successfully as an organizing force for communication to the science community. The MIT members of the Center provided many high-profile overview presentations at prestigious international conferences and national workshops. The research resulted in three student theses and 24 peer-reviewed publications. PSI research continues to be identified as a critical area for fusion energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Micro-cantilever beam experiments and modeling in porous polycrystalline UO 2

Understanding the impact of microstructure on the thermo-mechanical behavior of oxide nuclear fuels is vital to predicting their performance through multiscale models. Evaluating the mechanical properties at the sub-grain length scale is key to developing these multiscale models. In this work, 3D finite element (FE) models were constructed to simulate the micrometer-scale bending of micro-cantilever beams fabricated using porous polycrystalline uranium dioxide (UO 2 ) and tested at room temperature. Here, the results showed that the porosity and elastic anisotropy of individual grains can play a significant role in determining the effective mechanical properties of the material deduced from the tests. Specifically, the porosity had a non-negligible effect, given that the pore size was of the same order of magnitude as the dimensions of the micro-beams. Correlations between load-deflection data, pore location, and elastic properties (effective Young's modulus) were investigated using UO 2 micro-beam FE models, where pore clusters were included and placed at different locations along the length of the beam. Results indicated that the presence of pore clusters near the substrate, i.e., the clamp of the micro-cantilever beam, has the strongest effect on the load-deflection behavior, with the porosity leading to a reduction of stiffness that is the largest for any location of the pore clusters. Furthermore, it was also found that pore clusters located towards the middle of the span and close to the end of the beam have a comparatively small effect on the load-deflection behavior. Therefore, it is concluded that accurate estimates of Young's modulus can be obtained from micro-cantilever experiments after accounting for porosity on the one third of the beam length close to the clamp. This, in turn, provides an avenue to improve microscale experiments and their analysis in porous, anisotropic elastic materials.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

MOOSE ProbML: Parallelizable Probabilistic Machine Learning and Uncertainty Quantification Capabilities

The Multiphysics Object Oriented Simulation Environment (MOOSE) is a widely used open- source finite element software for performing multiphysics multiscale simulations in a massively parallel fashion. Recently, the computational team at Idaho National Laboratory (INL) has implemented Probabilistic Machine Learning (ProbML) capabilities in MOOSE—in a parallelized fashion—and enable active learning with large-scale computational models for tasks such as surrogate model development, scale bridging, forward/inverse uncertainty quantification (UQ), Bayesian optimization, etc. This presentation summarizes these developments in MOOSE along with demonstrations on several real applications relevant to nuclear energy. At the fundamental level, samplers like Monte Carlo/Latin Hypercube, variance reduction, parallelized Markov Chain Monte Carlo (MCMC) support uncertainty propagation in both forward and inverse settings. These samplers can be integrated with the Gaussian processes (GP) suite in MOOSE, which offer several variants like scalar GPs, multi-output GPs, and deep GPs, to enable active learning. These GPs can be tuned using gradient-based optimization methods like Adam and its variants or gradient-free methods like the elliptical slice sampler (a variant of MCMC adept under Gaussian settings) for more complex covariance kernels or likelihoods whose gradient computations can be cumbersome. A variety of batch acquisition functions permit parallelized evaluation of the computational model and support different learning objectives with high efficiency like Bayesian inference, global surrogate development, optimization, etc. Furthermore, libtorch integration supports training, evaluation, and re-training of neural networks and other complex machine learning models in active learning settings. The impacts of these developments are shown on several real applications: (1) nuclear fuel inverse UQ and model inadequacy assessment using the Kennedy O’Hagan framework; (2) uncertainty aware surrogate modeling for additive manufacturing to predict field quantities; (3) nuclear reactor rare events analysis; and (4) complex fluid flow prediction using a global surrogate with quantified prediction uncertainty. Finally, the outlook of MOOSE ProbML is discussed for both outer-loop and inner-loop computations in the broad view to accelerate fuels and materials qualification, address gaps in knowledge and data, and assess new reactor/fuel systems.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Multiphysics and Multiscale Modeling of Coupled Transport of Chloride Ions in Concrete

Chloride ions (Cl−)-induced corrosion is one of the main degradation mechanisms in reinforced concrete (RC) structures. In most situations, the degradation initiates with the transport of Cl− from the surface of the concrete towards the reinforcing steel. The accumulation of Cl− at the steel-concrete interface could initiate reinforcement corrosion once a threshold Cl− concentration is achieved. An accurate numerical model of the Cl− transport in concrete is required to predict the corrosion initiation in RC structures. However, existing numerical models lack a representation of the heterogenous concrete microstructure resulting from the varying environmental conditions and the indirect effect of time dependent temperature and relative humidity (RH) on the water adsorption and Cl− binding isotherms. In this study, a numerical model is developed to study the coupled transport of Cl− with heat, RH and oxygen (O2) into the concrete. The modeling of the concrete microstructure is performed using the Virtual Cement and Concrete Testing Laboratory (VCCTL) code developed by the U.S. National Institute of Standards and Technology (NIST). The concept of equivalent maturation time is utilized to eliminate the limitation of simulating concrete microstructure using VCCTL in specific environmental conditions such as adiabatic. Thus, a time-dependent concrete microstructure, which depends on the hydration reactions coupled with the temperature and RH of the environment, is achieved to study the Cl− transport. Additionally, Cl− binding isotherms, which are a function of the pH of the concrete pore solution, Cl− concentration, and weight fraction of mono-sulfate aluminate (AFm) and calcium-silicate-hydrate (C-S-H), obtained from an experimental study by the same authors are utilized to account for the Cl− binding of cement hydration products. The temperature dependent RH diffusion was considered to account for the transport of Cl− with moisture transport. The temperature and RH diffusion in the concrete domain, composite theory, and Cl− binding and water adsorption isotherms are used in combination, to estimate the ensuing Cl− diffusion field within the concrete. The coupled transport process of heat, RH, Cl−, and O2 is implemented in the Multiphysics Object-Oriented Simulation Environment (MOOSE) developed by the U.S. Idaho National Laboratory (INL). The model was verified and validated using data from multiple experimental studies with different concrete mixture proportions, curing durations, and environmental conditions. Additionally, a sensitivity analysis was performed to identify that the water-to-cement (w/c) ratio, the exposure duration, the boundary conditions: temperature, RH, surface Cl− concentration, Cl− diffusion coefficient in the capillary water, and the critical RH are the important parameters that govern the Cl− transport in RC structures. In a case study, the capabilities of the developed numerical model are demonstrated by studying the complex 2D diffusion of Cl− in a RC beam located in two different climatic regions: warm and humid weather in Galveston, Texas, and cold and dry weather in North Minnesota, Minnesota, subjected to time varying temperature, RH, and surface Cl− concentrations.

composite theory↗