Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “suitability”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Assessment the leaching characteristics and long-term leaching behavior of some radionuclides from synthesized zeolite cement matrix

Highlights: • Synthesis and characterization of zeolite cement synthesized from fly ash • Leaching tests were used to assess the release rate of {sup 85}Sr and {sup 60}Co radionuclides. • Simple mathematical models were used to detect leaching mechanisms. • Long-term leaching behavior of {sup 85}Sr and {sup 60}Co from cementitious matrix was assessed. The development of cementitious materials remains a vital goal to produce valuable products with good mechanical, physical and chemical properties suitable for the safe disposal of concentrated radionuclides resulted from the treatment of contaminated solutions. In this paper, zeolite cement with good properties was laboratory synthesized from industrial by-products fly ash and characterized using various characterization techniques. The leaching characteristics of strontium and cobalt radionuclides from the synthesized zeolite cement were scrutinized conforming with the standard leaching methodology of International Atomic Energy Authority (IAEA). A mechanical strength assessment was executed to characterize the extent of immobilization process of the solidified matrices. The cumulative leaching fraction (φ) of the two studied radionuclides was found to be less than 5% in all examined conditions, which implies the applicability of the IAEA-recommended methodology for estimating the diffusion coefficient. The experimental leaching data were regressed nonlinearly to various mathematical kinetic models to assess the controlling leaching mechanism and to determine the leaching parameters. The regression results indicated that strontium ({sup 85}Sr) and cobalt ({sup 60}Co) leaching resulted from two succeeding mechanisms: as first order kinetic reaction, and then diffusion. The calculated values of leachability indices signify that the performance of each of the studied matrices is within an acceptable range. A simplified mathematical model, rooted in the first order reaction and diffusion mechanisms, was simulated to predict the radionuclides leaching rates from zeolite cement matrices. By comparing the synthesized zeolite cement with other cementitious materials, it can be concluded that the synthesized material can be classified as an efficient material suitable to immobilize {sup 85}Sr and {sup 60}Co from radioactive wastes. The acquired findings demonstrated that the studied immobilized waste matrices have acceptable mechanical effectiveness.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Analysis of persistent contaminants and personal care products by dispersive liquid-liquid microextraction using hydrophobic magnetic deep eutectic solvents

Here, in this work, hydrophobic magnetic deep eutectic solvents (HMDESs) were used in the development of a simple and rapid dispersive liquid-liquid microextraction (DLLME) approach coupled to high performance liquid chromatography with UV detection (HPLC-UV) for the determination of ten organic contaminants including five polycyclic aromatic hydrocarbons, four UV filters, and a pesticide from water at trace levels. The HMDESs were prepared by mixing a hydrogen bond acceptor, metal halide salt, and hydrogen bond donor in suitable molar ratios. Two HMDESs, 2 tetraoctylammonium bromide ([N 8888 + ][Br - ]): cobalt chloride (CoCl 2 ): 4 octanoic acid (OA) and 3 trioctylphosphine oxide (TOPO): neodymium chloride (NdCl 3 ): 3 OA, offered the highest analyte extraction efficiency overall and were chosen as suitable solvents for validation of the microextraction method. Under optimized extraction conditions, the method required 30 µL of HMDES as extraction solvent, acetone (87.5 µL) as disperser solvent, a NaCl concentration of 30% (w/v), and an extraction time of 120 s at 20°C. Enrichment factors of the analytes ranged from 44.6 for 3-(4-methylbenzylindene) camphor to 66.0 for 2-ethylhexyl-4-(dimethyl)aminobenzoate. The method provided low limits of detection (LODs) ranging from 0.5 to 4.5 µg L -1 , and acceptable precision, with RSD values lower than 9.6%. Furthermore, the validated method was successfully applied for tap and lake water analysis, resulting in relative recoveries of spiked samples ranging between 94.7 and 119.2%.

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

Physics-Informed Machine Learning Models for Predicting the Progress of Reactive-Mixing

This paper presents a physics-informed machine learning (ML) framework to construct reduced-order models (ROMs) for reactive-transport quantities of interest (QoIs) based on high-fidelity numerical simu-lations. QoIs include species decay, product yield, and degree of mixing. The ROMs for QoIs are applied to quantify and understand how the chemical species evolve over time. First, high-resolution datasets for constructing ROMs are generated by solving anisotropic reaction-di?usion equations using a non-negative finite element formulation for di?erent input parameters. The reactive-mixing model input parameters are: time-scale associated with flipping of velocity, spatial-scale controlling small/large vortex structures of velocity, perturbation parameter of the vortex-based velocity, anisotropic dispersion strength/contrast, and molecular diffusion. Second, random forests, F-test, and mutual information criterion are used to evaluate the importance of model inputs/features with respect to QoIs. We observed that anisotropic dispersion strength/contrast is the most important feature and time-scale associated with flipping of velocity is the least important feature. Third, Support Vector Machines (SVM) and Support Vector Regression (SVR) are used to construct ROMs based on the model inputs. The constructed SVR-ROMs are then used to predict scaling of QoIs. We also present estimates and inequalities on the QoIs, which inform that the species decay, mix, and produce in an exponential fashion. These inequalities also inform that a radial basis function is the most suitable kernel for the SVM/SVR models for QoIs. It is observed that R2-score for SVR-ROMs on unseen data is greater than 0.9, implying that the SVR-ROMs are able to predict the reaction-diffusion system state reasonably well. Finally, in terms of the computational cost, the proposed SVM-ROMs are O(107) times faster than running a high-fidelity finite element simulation for evaluating QoIs. This makes the proposed ML-based ROMs attractive for reactive-transport sensing and real-time monitoring applications as they are significantly faster yet reasonably accurate.

Mudunuru, Maruti K.↗
DOE OSTI2021

Kink-antikink collisions and multi-bounce resonance windows in higher-order field theories

In this work, we study collisions of coherent structures in higher-order field-theoretic models, such as the Φ 8 , Φ 10 and Φ 12 ones. The main distinguishing feature, of the example models considered herein, is that the collision arises due to the long-range interacting algebraic tails of these solitary waves. We extend the approach to suitably initialize the relevant kinks, in the additional presence of finite initial velocity, in order to minimize the dispersive wave radiation potentially created by their slow spatial decay. We find that, when suitably initialized, these models still feature the multi-bounce resonance windows earlier found in models in which the kinks bear exponential tails, such as the Φ 4 and Φ 6 field theories among others. Also present is the self-similar structure of the associated windows with three- and more-bounce windows at the edges of two- and lower-bounce ones. Moreover, phenomenological but highly accurate (and predictive), scaling relations are derived for the dependence of the time between consecutive collisions and, e.g., the difference in kinetic energy between the incoming one and the critical one for one-bounces. Such scalings are traced extensively over two-bounce collision windows throughout the three models, hinting at the possibility of an analytical theory in this direction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Evaluation of a nonlinear variational multiscale method for fluid transport problems

Diverse transport problems, especially those based on fluid flow models, are intrinsically multiscale and nonlinear, characteristics that often lead to intricate dynamics such as the development of instabilities and turbulence. Computational simulations that resolve all scales in these problems are often unfeasible, prompting to coarse-grained simulation strategies in which small-scale features are modeled instead of resolved. Variational Multiscale (VMS) methods, and particularly residual-based Large-Eddy Simulation (LES) approaches, have proven effective and robust for the coarse-grained simulation of complex transport problems. VMS methods avoid the assumption of separable nonlinearity and the reliance on empirical small-scale models by using a variational decomposition of scales together with a residual-based approximation of the small-scales. Evaluation of a nonlinear VMS approach, denoted as VMS n , is presented for the coarse-grained simulation of transient-advective-diffusive-reactive (TADR) transport problems arising from fluid flow models. In contrast to classical VMS approaches that neglect the effect of the small scales on the transport operator, VMS n treats the inter-dependence between large- and small-scales upfront. The treatment of inter-scale coupling involves the solution of a local algebraic nonlinear system describing the evolution of the small-scales. The VMS n approach is complemented with two algebraic approximations of the small-scales: one based on the main diagonal of the transport matrices and another that preserves transport fluxes and is suitable for generic TADR systems. The suitability of the VMS n approach for handling general TADR problems and regimes is evaluated with benchmark incompressible, compressible, and magnetohydrodynamic laminar flow problems, the incompressible Taylor-Green vortex flow, the turbulent free jet, and the two-temperature arc in crossflow. Furthermore, simulation results show that VMS n leads to minor improvements in accuracy with respect to the classical VMS for the laminar flow problems, but to significantly greater accuracy for the turbulent flows and the unsteady plasma flow problems, while using the same cohesive numerical formulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Simulating low-energy neutrino interactions with MARLEY

Monte Carlo event generators are a critical tool for the interpretation of data obtained by neutrino experiments. Several modern event generators are available which are well-suited to the GeV energy scale used in studies of accelerator neutrinos. However, theoretical modeling differences make their immediate application to lower energies difficult. In this paper, I present a new event generator, MARLEY, which is designed to better address the simulation needs of the low-energy (tens of MeV and below) neutrino community. The code is written in C++14 with an optional interface to the popular ROOT data analysis framework. The current release of MARLEY (version 1.2.0) emphasizes simulations of the reaction 40 Ar ( ν e , e − ) 40 K ⁎ but is extensible to other channels with suitable user input. This paper provides detailed documentation of MARLEY's implementation and usage, including guidance on how generated events may be analyzed and how MARLEY may be interfaced with external codes such as Geant4. Further information about MARLEY is available on the official website at http://www.marleygen.org. Program title:MARLEY 1.2.0 CPC Library link to program files:https://doi.org/10.17632/4v7zxnc8j3.1 Developer's respository link:http://github.com/MARLEY-MC/marley Code Ocean capsule:https://codeocean.com/capsule/9868179 Licensing provisions: GNU General Public License 3.0 Programming language: C++14 External routines/libraries used: GNU Scientific Library [1,2] (required), ROOT [3,4] (optional) Nature of problem: Simulation of neutrino-nucleus scattering events at energies of tens-of-MeV and below Solution method: Initial two-to-two scattering kinematics are sampled using the allowed approximation differential cross section and tables of precomputed nuclear matrix elements. Subsequent de-excitations of the remnant nucleus are simulated using a Monte Carlo implementation of the Hauser-Feshbach statistical model and tabulated γ-ray decay schemes for discrete nuclear levels. Additional comments including restrictions and unusual features: Input data are provided with the code that are suitable for producing simulations of the charged-current reaction 40 Ar ( ν e , e − ) 40 K ⁎ , coherent elastic neutrino-nucleus scattering on spin-zero target nuclei, and neutrino-electron elastic scattering on any atomic target. Preparation of new reaction input files (whose format is documented in Appendix B) would enable other reaction channels and nuclear targets to be handled by the existing code framework. Although there is no maximum neutrino energy enforced by the code itself, realistic neutrino-nucleus scattering events may be generated up to roughly 50 MeV. Above this energy, the effects of forbidden nuclear transitions, which are neglected in the current treatment of the cross sections (see section 2.1), become increasingly important.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

A bi-layer barrier design for 122-type iron-based superconducting wires and tapes

Iron-based superconducting wires and tapes hold great promise for high-field magnet applications. A promising design for 122-type wires and tapes based on the powder-in-tube method is using silver and copper double-layer sheaths. For this design a heat treatment temperature below ~ 779 °C is required to prevent Ag-Cu liquid formation. However, this may be below the optimal heat treatment temperature for the critical current density, and still cannot prevent Ag-Cu interdiffusion occurring in the solid state. In this work we propose adding a niobium or tantalum or vanadium (or their alloys) barrier layer between the Ag and Cu to solve the Ag-Cu interdiffusion issue, given that the group-VB metals (vanadium, niobium, tantalum) are relatively inert to both Ag and Cu. To investigate the effectiveness of this design, BaFe 1.84 Co 0.16 As 2 wires and tapes with Ag/Cu and Ag/Ta/Cu sheaths, as well as Ba 0.6 K 0.4 Fe 2 As 2 wires and tapes with Ag/Cu and Ag/Nb/Cu sheaths, were fabricated. It was found that both the Ta and Nb layers kept integral after wire drawing, but after a large flat-rolling reduction the Ta layer broke while the Nb layer kept integral. In the tapes with Ag/Cu sheaths (without the Ta or Nb layer) Cu diffused through the Ag layer and into the powder cores during 740 °C heat treatment, while in the tapes with Ag/Nb/Cu sheaths the Nb layer effectively blocked Ag-Cu interdiffusion even at 900 °C. Further, this work demonstrates that Ta is a suitable barrier material for 122-type wires, while Nb is suitable for both wires and tapes. In this design using Ag/Nb (or Ta)/Cu sheaths, we can regard the outer Cu as the conductor matrix while the Ag and Nb (or Ta) serve as two layers of barriers that suppress reactions between the components. Thus, we call this design a “bi-layer barrier” design for 122-type wires and tapes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Evaluating performance of different generative adversarial networks for large-scale building power demand prediction

We report as an unsupervised-learning data-driven model, Generative Adversarial Networks (GANs) have recently attracted a lot of attention for various applications. There is potential to apply GANs for large-scale building power demand prediction, which is needed for power grid operation. However, there are many GAN variations and it is unclear which GAN is suitable for this application. To answer this question, this paper identifies five promising GANs (Original GAN, cGAN, SGAN, InfoGAN, and ACGAN) and evaluates their performance for predicting building power demand at a large scale. Physics-based building energy models are developed to generate training and reference data. A new evaluation indicator that combines accuracy and reproducibility is proposed to evaluate the performance of different GANs in predicting building power demand. The results show that SGAN and InfoGAN are not suitable because they cannot control the number of generated building samples for different building types. The prediction performance among the Original GAN, cGAN, and ACGAN can vary depending on training sample sizes and number of building types. If the training sample size is sufficiently large, Original GAN and cGAN can predict building power demand more accurately than ACGAN with the same number of samples. If training samples are limited, Original GAN provides better accuracy than cGAN and ACGAN. When the number of building types increase, the prediction accuracy increases for cGAN, decreases for ACGAN, and remains the same for Original GAN. As a result, cGAN and Original GAN are recommended for large-scale building power demand prediction.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Spatial and temporal characterization of municipal solid waste based on resource recovery pathways

This study presents a two-year, quarterly assessment of MSW across four source sectors (residential, schools, restaurants, and grocery stores) from sixteen sites across five U.S. states. MSW was manually sorted into 27 categories and aggregated into pathway fractions: high-moisture (HM) organics, low-moisture (LM) organics, recyclable (RC) materials, and residuals for disposal. Organics represented 89 % of the MSW stream. The largest fraction was HM organics consisting of food waste (31 %) and yard waste (3 %), with large coefficient of variations (CV), 79 and 278 %, respectively, reflecting high seasonal and site variability that varied significantly (p < 0.01) across sampling periods. The HM fraction showed properties favorable for anaerobic digestion, with moisture content ranging from 56 to 95 % and volatile solids ranges of 86-95 %. In contrast, the LM and RC fractions remained more stable (plastics CV = 41 %; paper CV = 53 %) with heating values up to 26.9 MJ/kg across sources, reflecting suitability for gasification. Microstructural analysis revealed less porosity in residential waste sampled at the landfill, which can influence preprocessing efficiency and microbial accessibility. Pathway informed allocations showed that 35 % of MSW is suitable for anaerobic digestion, 36 % for gasification, and 18 % for recycling, leaving 11 % requiring landfill disposal. These results provide quantitative evidence to determine feedstock allocation, waste-to-energy system design, and the development of data-driven sustainability and resource recovery strategies within a circular bioeconomy.

09 BIOMASS FUELS↗
DOE OSTI2025

Enhancing building resilience in cold climates: Integrating heat pump technologies with renewable energy

As electrification advances and Cold Climate Heat Pump technology progresses, ensuring grid stability becomes increasingly critical for effective heating in cold climates. However, natural disasters, especially during winter, pose significant threats to grid stability, impacting the reliability of air-source heat pumps. Despite these challenges, the integration of renewable energy sources and storage solutions in heating systems has not been extensively studied within the context of resilience. Here, this paper delves into the literature on renewable-powered heat pumps to assess their potential in enhancing building resilience in U.S. cold climate zones, which are particularly susceptible to extreme weather and grid disruptions. By leveraging renewable sources—solar, geothermal, and water—in conjunction with heat pump technology and supported by thermal or battery storage, this approach aims to provide a dependable solution for maintaining indoor heating during grid failures. Our analysis begins with a review of various renewable energy sources suitable for heat pumps, followed by an exploration of their application in cold climate regions across the U.S., and discussions on potential integration strategies with heat pump systems. This study highlights the advantages and suitability of solar irradiance and geothermal resources, emphasizing the importance of tailored, site-specific assessments to maximize energy efficiency and resilience. Additionally, it outlines the economic and environmental considerations necessary for implementing such systems and identifies potential challenges and areas for future research to facilitate the broader integration of renewable energy in heating solutions for enhanced resilience.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Whole-tree dormant season nitrogen pools for different species receiving combinations of fertilization and irrigation after one short rotation

Our understanding of the accretion and distribution of tree tissue nitrogen (N) pools across a variety of species and genotypes suitable for short-rotation woody crop (SRWC) production in response to water and N availability remains limited. We measured dormant-season, rotation-length, whole-tree N pools for five tree genotypes from four species (two eastern cottonwood, Populus deltoides Bartr., genotypes; American sycamore, Platanus occidentalis L.; American sweetgum, Liquidambar styraciflua L.; and loblolly pine, Pinus taeda L.) receiving irrigation (I), fertilization (F), their combination (IF), or no resources manipulation (C). Our results demonstrate that foliar nitrogen concentration [N] responded to fertilization but was constrained within genotype-specific ranges and varied temporally. Tree genotypes differed in their composite and component tissue N content (N C ), and these differences mostly reflected tissue mass (i.e., larger components and trees resulted in higher N C ). Resource amendments (I, F, IF) resulted in up to 3.8-fold increases in N C compared with C, which were most pronounced for sycamore and sweetgum with F and IF treatments, respectively. By the end of the rotation, forest stands accumulated 73 to 452 kg N ha –1 in tree tissues with 40 – 78 % distributed aboveground and 22 – 60 % distributed belowground. A critical difference between genotypes was that all hardwoods exhibited larger belowground N distributions than the evergreen conifer. Our results stress the importance of belowground N pools and highlight differences among genotypes. Furthermore, our study underscores valuable information about N pools across genotypes suitable for SRWC production, which can be leveraged to inform fertilization plans and devise sustainable nutrient management as production expands across marginal lands.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Assessment of oil and gas fields in California as potential CO 2 storage sites

California's total annual greenhouse gas (GHG) emissions (425.3 MtCO 2 e) in 2018 were about 6.4% of the US total (6,677 MtCO 2 e) and around 1% of global emissions. About 39% of 2018 GHG emissions in California were from the industrial and electrical sectors. Many of these emissions were from large stationary point sources and were suitable for carbon capture retrofit with subsequent storage of the captured carbon dioxide (CO 2 ) in geological formations. Previous studies of California found suitable geology and CO 2 storage resource. This study refines and furthers prior work using a three-stage screening process of oil fields, gas fields, and underground natural gas storage (UGS) sites by combining criteria from previous studies while excluding sites that pose technical risk or are located in regions with surface restrictions including sensitive habitats and dense populations. In the first stage, 129 CO 2 storage sites in California were identified using qualification criteria based upon formation properties including geological conditions and pore pressure. The second stage identified sensitive sites by applying conservative screens including seismic activity, faulting, population density, restricted lands, and sensitive habitats. During the third stage, 61 CO 2 potential storage sites were identified by subtraction of stage 2 areas from stage 1. The potential storage volume in the third stage ranged from 1.0 to 2.0 GtCO 2 . Finally, we applied a scoring system with seven parameters to rank the 61 potential sites based on subsurface technical criteria. The scored sites are classified as high priority, medium priority, and sites for future study. Prospective CO 2 storage sites with high and moderate priority were selected and linked to CO 2 sources. There are 14 prospective sites (above 20 MtCO 2 storage resource per site) with a total storage resource of 1024 MtCO 2 distributed in Northern and Southern California. Of these sites, there are 9 potential CO 2 -EOR sites and 1 depleted oil field with a total estimated CO 2 storage volume of ~800 MtCO 2 in the Southern San Joaquin and Ventura Basin. These 10 prospective sites with a storage resource greater than 20 MtCO 2 could potentially deliver more than 20 years of storage with an average injection rate of 40 MtCO 2 /year. The remaining 4 highly prospective sites are in Northern California. Additionally, study results suggest that saline formations should be re-evaluated in concert with storage in oil, gas, and natural gas storage reservoirs.

58 GEOSCIENCES↗
DOE OSTI2022

GCS site selection in saline Miocene formations in South Louisiana

One of the key requirements of a commercial scale geologic carbon sequestration (GCS) site selection process is identifying injection sites to safely inject and retain injected CO 2 over extended time scales. This requires thorough evaluation of site-specific surface and subsurface features. A site screening and selection framework was developed and used to identify potential sites in South Louisiana, incorporating surface data including site proximity to large CO 2 stationary sources, vicinity of large population centres and state lands, land usage, legacy oil and gas wells, and local/regional geologic trends. Regionally the Miocene sediments thicken and dip towards the south, reflecting a prograding deltaic environment during the Miocene which resulted in the deposition of thick packages of fine-grained sandstones suitable for CO 2 storage. The results of three example selected sites that are spatially separated with varying subsurface and surface features are presented in this study. The multiple stacked sand zones at these sites provide ample amount of pore volume for CO 2 storage with varying degree of stratigraphic/structural traps for long terms storage containment. Here, the screening framework presented in this study and the detailed evaluation of three specific sites in South Louisiana provides a process to screen and characterize sites for GCS projects and the suitability of these Miocene sediments for commercial scale CO 2 storage projects. It is envisioned that the analysis described in this paper will be beneficial for a large number of audiences involved in GCS work and aid in the site selection process.

58 GEOSCIENCES↗
DOE OSTI2022

Effect of secondary gas-phase reactions (SGR) in pyrolysis of carbon feedstocks for anisotropic carbon materials production – 1: Controlling SGR to modify intermediate coal tar species to improve pitch anisotropy

To meet the increasing demand for graphitizable carbon products, such as needle coke and carbon fiber, more carbon feedstocks capable of forming anisotropy should be utilized. Non-coking coals are widely available but are not typically suitable for producing anisotropic carbons due to lacking proper coal chemistry. This work used secondary gas-phase reactions (SGR) during coal pyrolysis to improve the coal tar chemistry of a non-coking coal for anisotropic carbon production. SGR pyrolysis temperatures and residence times were varied (T = 800–900°C and τ = 0–2.5 s), and analysis of the intermediate coal tar products showed that as these SGR pyrolysis conditions increased, the oxygen and aliphatic concentrations decreased, whereas aromatic contents and molecular weights increased. Without any SGR, microscopy of the thermally-treated coal tar pitch product revealed that the coal tar pitch was isotropic; however, upon using increased SGR pyrolysis conditions, the resulting coal tar pitch samples substantially increased in the percentage and quality of anisotropy formation. Further, the products' analyses show clear trends of modified chemical properties in the intermediate coal tar and improved anisotropy results. Thus, the results presented in this work show that by controlling the SGR during pyrolysis of a non-coking coal, this approach can modify the coal tar chemistry towards a precursor more suitable for quality anisotropic carbon material production.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

Slot-die-coating operability windows for polymer electrolyte membrane fuel cell cathode catalyst layers

Roll-to-roll (R2R) slot-die coating of polymer electrolyte membrane fuel cell (PEMFC) catalyst layers represents a scalable deposition method for producing 10–20 m 2 ·min -1 of catalyst-coated gas diffusion layers (GDLs). This high-throughput production technique will help lower the cost of PEMFC catalyst layers. The uniformity of the wet layer applied by slot die deposition is affected by process parameters such as substrate speed, vacuum pressure applied at the upstream meniscus, gap between the slot die lips and substrate, ink rheology, and other ink and substrate properties. The set of conditions for producing a defect-free coating with a dilute ink typically requires little to no upstream vacuum pressure, so suitable operating conditions can be found easily through trial and error and operator intuition. However, the higher viscosity of more concentrated inks dramatically shifts the range of settings that result in a homogeneous coating to higher vacuum levels, which are harder to find through hit or miss. A predictive model showing the range of operable conditions decreases material wastage inherent in experimentally searching for suitable parameters. In this study, the defect-free coating parameter window is explored experimentally and theoretically for two concentrations of PEFC cathode inks. Both a full capillary hydrodynamic model and a computationally cheaper viscocapillary model successfully predict the experimentally determined coating window within the experimental and model uncertainty limits for inks with 5.3 wt% and 12.0 wt% solids ink while maintaining the 0.1 mg Pt ·cm -2 Pt areal loading target. Here, this paper demonstrates a viable pathway for meeting the $30/kW net ultimate cost target of the United States Department of Energy (U.S. DOE) Hydrogen Fuel Cells Technologies Office (HFTO). The concentrated ink lowers the thermal energy and capital expenditure (CapEx) budget of the coating process by decreasing the amount of time, energy, and floorspace required for drying the coating.

25 ENERGY STORAGE↗
DOE OSTI2021

B-PINNs: Bayesian Physics-informal Neural Networks for Forward and Inverse PDE Problems with Noisy Data

We propose a Bayesian physics-informed neural network (B-PINN) to solve both forward and inverse nonlinear problems described by partial differential equations (PDEs) and noisy data. In this Bayesian framework, the Bayesian neural network (BNN) combined with a PINN for PDEs serves as the prior while the Hamiltonian Monte Carlo (HMC) or the variational inference (VI) could serve as an estimator of the posterior. B-PINNs make use of both physical laws and scattered noisy measurements to provide predictions and quantify the aleatoric uncertainty arising from the noisy data in the Bayesian framework. Compared with PINNs, in addition to uncertainty quantification, B-PINNs obtain more accurate predictions in scenarios with large noise due to their capability of avoiding overfitting. We conduct a systematic comparison between the two different approaches for the B-PINNs posterior estimation (i.e., HMC or VI), along with dropout used for quantifying uncertainty in deep neural networks. Our experiments show that HMC is more suitable than VI with mean field Gaussian approximation for the B-PINNs posterior estimation, while dropout employed in PINNs can hardly provide accurate predictions with reasonable uncertainty. Finally, we replace the BNN in the prior with a truncated Karhunen-Loève (KL) expansion combined with HMC or a deep normalizing flow (DNF) model as posterior estimators. The KL is as accurate as BNN and much faster but this framework cannot be easily extended to high-dimensional problems unlike the BNN based framework.

Non linear PDEs, Noisy data, Bysian physics inform↗
DOE OSTI2021

Robust verification of stochastic simulation codes

We introduce a robust verification tool for computational codes, which we call Stochastic Robust Extrapolation based Error Quantification (StREEQ). Unlike the prevalent Grid Convergence Index (GCI) [1] method, our approach is suitable for both stochastic and deterministic computational codes and is generalizable to any number of discretization variables. Building on ideas introduced in the Robust Verification [2] approach, we estimate the converged solution and orders of convergence with uncertainty using multiple fits of a discretization error model. In contrast to Robust Verification, we perform these fits to many bootstrap samples yielding a larger set of predictions with smoother statistics. Here, bootstrap resampling is performed on the lack-of-fit errors for deterministic code responses, and directly on the noisy data set for stochastic responses. This approach lends a degree of robustness to the overall results, capable of yielding precise verification results for sufficiently resolved data sets, and appropriately expanding the uncertainty when the data set does not support a precise result. For stochastic responses, a credibility assessment is also performed to give the analyst an indication of the trustworthiness of the results. Furthermore, this approach is suitable for both code and solution verification, and is particularly useful for solution verification of high-consequence simulations..

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A survey of unsupervised learning methods for high-dimensional uncertainty quantification in black-box-type problems

Constructing surrogate models for uncertainty quantification (UQ) on complex partial differential equations (PDEs) having inherently high-dimensional O(10 n ), n ≥ 2, stochastic inputs (e.g., forcing terms, boundary conditions, initial conditions) poses tremendous challenges. The “curse of dimensionality” can be addressed with suitable unsupervised learning techniques used as a pre-processing tool to encode inputs onto lower-dimensional subspaces while retaining its structural information and meaningful properties. In this work, we review and investigate thirteen dimension reduction methods including linear and nonlinear, spectral, blind source separation, convex and non-convex methods and utilize the resulting embeddings to construct a mapping to quantities of interest via polynomial chaos expansions (PCE). Here, we refer to the general proposed approach as manifold PCE (m-PCE), where manifold corresponds to the latent space resulting from any of the studied dimension reduction methods. To investigate the capabilities and limitations of these methods we conduct numerical tests for three physics-based systems (treated as black-boxes) having high-dimensional stochastic inputs of varying complexity modeled as both Gaussian and non-Gaussian random fields to investigate the effect of the intrinsic dimensionality of input data. We demonstrate both the advantages and limitations of the unsupervised learning methods and we conclude that a suitable m-PCE model provides a cost-effective approach compared to alternative algorithms proposed in the literature, including recently proposed expensive deep neural network-based surrogates and can be readily applied for high-dimensional UQ in stochastic PDEs.

42 ENGINEERING↗