Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model uncertainty”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,171 records · Page 65

PVDeg: Development of a Streamlined Tool for PV Degradation Modeling

The photovoltaic (PV) industry constantly aims for lower costs, higher-efficiency cells, and improved module designs. These trends lead to using new materials, designs, and manufacturing processes, resulting in a continually changing technological landscape. These changes can potentially introduce new, unknown degradation mechanisms and failure modes that are difficult to diagnose, analyze, test, and model. This introduces uncertainty into the expected lifetime of PV modules of 25 to 50 years. research efforts aim to achieve this while keeping performance degradation at a minimum for decades. This puts considerable pressure on improving the accuracy of long-term durability and reliability assessments. There is a need to organize the existing degradation data into an accessible format and to provide industry relevant tools for extrapolation from laboratory to field conditions. Because the core of this type of analysis involves calculations that are complicated but ubiquitous for many degradation processes, an enhanced predictive modeling framework will facilitate the analysis to help researchers keep up with the rapid pace of technological changes. In this work, we present an online tool that can be used to search for and analyze degradation information and extrapolate PV module performance and durability to field exposure. The tool will simplify many of the routine computational operations that are common to many degradation studies. The prediction tool will be built modular and published as open source, enabling users to expand on the existing framework. This repository will contain various degradation models and material parameters suitable for the reliability and durability assessment of materials and components deployed outdoors.

degradation↗

U.S. Efforts in Support of Examinations at Fukushima Daiichi (November 2021 Meeting Notes with Updated Information Requests)

Information obtained from Fukushima Daiichi Nuclear Power Station (Daiichi) is required to inform future Decontamination and Decommissioning (D&D) activities, improving the ability of the Tokyo Electric Power Company Holdings, Incorporated (TEPCO Holdings) to characterize potential hazards and to ensure the safety of workers involved with cleanup activities. This information also has important implications for the safety and operation of U.S. commercial nuclear power plants. This document summarizes results from the Fiscal Year 2022 (FY2022) U.S. effort to review Daiichi information and extract insights to enhance the safety of existing and future nuclear power plant designs. This U.S. effort, which was initiated in 2014 by the Department of Energy Office of Nuclear Energy, is completed by a group of experts in reactor safety and plant operations that identify examination needs and evaluate recent Daiichi examination data to address these needs. Fukushima-related information and associated discussions during forensics meetings benefit operating, new, and advanced reactors. Significant safety insights are being obtained in several areas: system and component performance, radionuclide surveys and sampling, debris end-state location, combustible gas effects, and plant operations and maintenance. In addition to reducing uncertainties related to severe accident modeling progression, these insights have and continue to be used to update guidance for severe accident prevention, mitigation, and emergency planning. As discussed in this document, revised operator guidance was successfully used to improve operator response during a loss of off-site power event at the Duane Arnold Energy Center plant. Reduced uncertainties in modeling the events at Daiichi improve the realism of reactor safety evaluations that inform future D&D activities. U.S. evaluations of information from Fukushima and input regarding future examinations are of interest to several organizations within Japan. Meeting presentations by Japan describe how comments and recommendations documented in prior U.S. forensics effort reports, including consensus information requests developed by forensics effort participants, are considered in future Fukushima D&D activities. As discussed in this report, TEPCO Holdings considered these information requests in their D&D planning activities. An updated list of consensus information requests is included in this FY2022 report.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Southern Ocean Aerosols Field Campaign Report

Atmospheric processes over the Southern Ocean have a profound influence on regional and global climate. This is a part of the world where global climate models perform particularly poorly, with models persistently overpredicting the amount of sunlight reaching the Earth's surface (Regayre et al. 2020). A major challenge when trying to improve the representation of atmospheric processes in this part of the globe is the scarcity of observations capable of constraining our understanding. During 2024/25 a U.S. Department of Energy Atmospheric Radiation Measurement (ARM) User Facility field campaign (CAPE-k) deployed a suite of instrumentation at the kennaook/Cape Grim atmospheric monitoring station in Tasmania, Australia to provide detailed cloud-aerosol observations in this region in order to address the greatest source of uncertainty in climate models – atmospheric aerosols and how they influence cloud formation (Carslaw et al. 2013). In support of the ARM observations, an Australian Research Council (ARC)-funded project “Southern Ocean aerosols: sources, sinks and impact on cloud properties,” led by the Queensland University of Technology, complemented the ARM measurements with a suite of chemical measurements at the kennaook/Cape Grim site in northwestern Tasmania to provide more information on the processes controlling aerosol formation and growth. The deployment of the University of York LIF-SO2 instrument (Temple et al. 2025) in the ARM mobile facility container was part of this chemistry-focused deployment. Increased observational efforts are critical for understanding sources and sinks of Southern Ocean aerosols, in particular the role of marine micro-organisms (phytoplankton, algae) on aerosols formation, their properties and growth, cloud droplet formation, and the conversion of cloud droplets to ice crystals and precipitation. The observational focus of the ARC-funded project is aerosol chemical composition and their gaseous precursors. The remote nature of the Southern Ocean poses a significant analytical challenge when studying gaseous precursors, as low concentrations are often below the limits of detection of commercial instrumentation. Historically this has meant that observations have had to be time-averaged over days to weeks to achieve the limits of detection required. Although useful when considering overall average levels of aerosol precursors, this time-averaging obscures temporal variability that can provide insight into the controlling processes. Over recent years advances in analytical techniques have made high-time-resolution measurements of key gas-phase aerosol precursors achievable with sufficiently low limits of detection.

54 ENVIRONMENTAL SCIENCES↗

Estimating the color lifetime of energetic quarks

Using a simple geometric framework with a realistic nuclear density distribution, we fit published HERMES data to determine fundamental properties of hadronization using the nuclear medium as a spatial analyzer. Our approach uses a fit to the transverse momentum broadening observable and the hadronic multiplicity ratio; the simultaneous fit to two different observables strongly constrains the outcome. Using the known sizes of the target nuclei, we extract the color lifetime, finding a z h -dependent range of values from 2 to 8 fm/c for these data. We also extract estimates for the qˆ transport coefficient characterizing the strength of the interaction between the quark and the cold nuclear medium, finding an average value of 0.035±0.011 GeV 2 /fm. With a three-parameter model we obtain satisfactory fits to the data with a goodness-of-fit parameter χ 2 /dof of 1.1 or less. In a secondary fit of the results from that model we independently find a value for the Lund String Model string tension of 1.00±0.05 GeV/fm. We evaluated the sensitivity for extracting quark energy loss and effective in-medium hadronic cross sections using four-parameter variants of the model, finding large uncertainties in both cases. Our results suggest that hadronic interaction of forming hadrons in the nuclear medium is the primary dynamical cause of meson attenuation in the HERMES data, with quark energy loss playing a more minor role.

Electron-ion collider↗

Industry Facing PV Degradation Prediction Tool and Database to Enable a 50-Year Life Module

The Photovoltaic (PV) industry constantly aims for lower costs, higher-efficiency cells, and improved module designs. These trends lead to using new materials, designs, and manufacturing processes, resulting in a continually changing technological landscape. These changes can potentially introduce new, unknown degradation mechanisms and failure modes that are difficult to diagnose, analyze, test, and model. This introduces uncertainty into the expected lifetime of PV modules of 25 years. Furthermore, research efforts aim for up to 50 years of service life while keeping performance degradation at a minimum for decades of outdoor weathering - putting additional pressure on improving the accuracy of long-term durability and reliability assessments. Here there is a need to organize the existing degradation data into an accessible format and to provide industry relevant tools for extrapolation from laboratory to field conditions. Because the core of this type of analysis involves calculations that are complicated but ubiquitous for many degradation processes, an enhanced predictive modeling framework will facilitate the analysis to help researchers keep up with the rapid pace of technological changes. In this work, we present an online tool that can be used to search for and analyze degradation information and extrapolate PV module performance and durability to field exposure. A graphical user interface will aid in the understanding of the results. The prediction tool will be built modular and published as open source, enabling users to expand on the existing framework. We use an integration pipeline approach that allows us to leverage weather data from the National Solar Radiation Database to perform geospatial degradation analysis in the US and worldwide. Our repository will contain various degradation models and material parameters suitable for the reliability and durability assessment of materials and components deployed outdoors. We hope to become a repository that can be used for weathering and degradation analysis for various applications beyond the PV industry.

degradation↗

Neutrino Flavor Model Building and the Origins of Flavor and CP Violation

The neutrino sector offers one of the most sensitive probes of new physics beyond the Standard Model of Particle Physics (SM). The mechanism of neutrino mass generation is still unknown. The observed suppression of neutrino masses hints at a large scale, conceivably of the order of the scale of a rand unified theory (GUT), which is a unique feature of neutrinos that is not shared by the charged fermions. The origin of neutrino masses and mixing is part of the outstanding puzzle of fermion masses and mixings, which is not explained ab initio in the SM. Flavor model building for both quark and lepton sectors is important in order to gain a better understanding of the origin of the structure of mass hierarchy and flavor mixing, which constitute the dominant fraction of the SM parameters. Recent activities in neutrino flavor model building based on non-Abelian discrete flavor symmetries and modular flavor symmetries have been shown to be a promising direction to explore. The emerging models provide a framework that has a significantly reduced number of undetermined parameters in the flavor sector. In addition, such a framework affords a novel origin of CP violation from group theory due to the intimate connection between physical CP transformation and group theoretical properties of non-Abelian discrete groups. Model building based on non-Abelian discrete flavor symmetries and their modular variants enables the particle physics community to interpret the current and anticipated upcoming data from neutrino experiments. Non-Abelian discrete flavor symmetries and their modular variants can result from compactification of a higher-dimensional theory. Pursuit of flavor model building based on such frameworks thus also provides the connection to possible UV completions: in particular, to string theory. We emphasize the importance of constructing models in which the uncertainties of theoretical predictions are smaller than, or at most compatible with, the error bars of measurements in neutrino experiments. While there exist proof-of-principle versions of bottom-up models in which the theoretical uncertainties are under control, it is remarkable that the key ingredients of such constructions were discovered first in top-down model building. We outline how a successful unification of bottom-up and top-down ideas and techniques may guide us towards a new era of precision flavor model building in which future experimental results can give us crucial insights into the UV completion of the SM.

Almumin, Yahya (ORCID:0000000196396772)↗

Moisture availability mediates the relationship between terrestrial gross primary production and solar-induced chlorophyll fluorescence: Insights from global-scale variations

Effective use of solar-induced chlorophyll fluorescence (SIF) to estimate and monitor gross primary production (GPP) in terrestrial ecosystems requires a comprehensive understanding and quantification of the relationship between SIF and GPP. To date, this understanding is incomplete and somewhat controversial in the literature. Here we derived the GPP/SIF ratio from multiple data sources as a diagnostic metric to explore its global-scale patterns of spatial variation and potential climatic dependence. We found that the growing season GPP/SIF ratio varied substantially across global land surfaces, with the highest ratios consistently found in boreal regions. Spatial variation in GPP/SIF was strongly modulated by climate variables. The most striking pattern was a consistent decrease in GPP/SIF from cold-and-wet climates to hot-and-dry climates. We propose that the reduction in GPP/SIF with decreasing moisture availability may be related to stomatal responses to aridity. Furthermore, we show that GPP/SIF can be empirically modeled from climate variables using a machine learning (random forest) framework, which can improve the modeling of ecosystem production and quantify its uncertainty in global terrestrial biosphere models. Finally, our results point to the need for targeted field and experimental studies to better understand the patterns observed and to improve the modeling of the relationship between SIF and GPP over broad scales.

59 BASIC BIOLOGICAL SCIENCES↗

Evaluation and projection of long period return values of extreme daily precipitation in the CMIP5 and CMIP6 models

Using a non-stationary Generalized Extreme Value statistical method, we calculate selected extreme daily precipitation indices and their 20 year return values from the CMIP5 and CMIP6 climate models over the historical and future periods. We evaluate model performance of these indices and their return values in replicating similar quantities calculated from multiple gridded observational products. Difficulties in interpreting model quality in the context of observational uncertainties are discussed. Projections are framed in terms of specified global warming target temperatures rather than at specific times and under specific emissions scenarios. The change in framing shifts projection uncertainty due to differences in model climate sensitivity from the values of the projections to the timing of the global warming target. At their standard resolutions, we find there are no meaningful differences between the two generations of models in their quality or projections of simulated extreme daily precipitation.

54 ENVIRONMENTAL SCIENCES↗

Incorporating the Impacts of Climate Change on Hydrology in a Performance Assessment Model - 20403

The evidence of climate change is increasingly well-documented and impacts should be incorporated in performance assessment studies. The current climate literature provides both observational evidence and climate model projections of climate trends and/or climate change in the late 20. and early 21. centuries for North America and the northeast United States. Probabilistic modeling is a core requirement for quantifying uncertainty and evaluating its impacts. Not evaluating future climate states in a performance assessment because of the existence of uncertainty is contradictory to good modeling practices - the most uncertain issues and parameters require the most attention in effective probabilistic modeling. Excluding climate change limits development of modeling information that could aid in effective decision making. In this work we develop methods to use the output from hydrologic models and analysis of historical aerial imagery to quantify and implement the impacts of climate change on hydrologic processes at a nuclear waste site in West Valley, New York. Specifically, we used the HELP (Hydraulic Performance of Landfill Performance) model to characterize key hydrologic processes under both current and future climate conditions to assess the impacts of changing climate on hydrology. A suite of previous climatic models were reviewed and synthesized to produce a cohesive representation of the current state of knowledge of the impact of climate change on important model inputs such as precipitation. Output from the HELP simulations was coupled to the GoldSim model that was used to develop the Probabilistic Performance Assessment (PPA) approach through the application of a novel 'nearest neighbor' technique. First, several thousand realizations were generated from the HELP model using a Latin Hypercube experimental design to ensure adequate coverage of the parameter space of explanatory variables used to drive HELP. For example, porosity is a physical parameter that is used as an input to both HELP and the GoldSim PA model. We then ran sensitivity analysis (SA) algorithms on the output of HELP for each of the responses of interest. For each predictor, each time we build an SA model we get a different value for the sensitivity index (SI). From the collection of all the SA models, the average was computed among all of the SIs to represent the predictor within the context of the nearest neighbor approach. That is, we conduct SA on each HELP outcome for each scenario. This gives us parameter sensitivity indices for the outcomes. We average the parameter sensitivity indices across the outcomes to get the average SI for a scenario. For each realization that is generated from the Goldsim PA model, Goldsim generates random values for physical/empirical parameters that HELP uses as well. For each vector of physical/empirical parameters that Goldsim generates, the vector from the 5,000 HELP runs that is most 'similar' to the Goldsim vector is computed using the nearest neighbor approach. In this context 'similar' means minimization of the SA-weighted sum of the absolute differences among the 5,000 values computed for this statistic, where each value corresponds to a different HELP realization. In order to account for the impacts of climate change, this process was repeated using the spatially downscaled future climate projections. For each of the key parameters of interest, it was assumed that a linear change depicted the relationship between the values for the present day and those for 2100. In this way, the climatically-driven changes in key parameters used to inform the GoldSim model are quantified and incorporated into the PA model output for the future. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

GPLaSDI: Gaussian Process Latent Space Dynamics Identification

GPLaSDI is a software that introduces GPLaSDI, a novel LaSDI-based framework that relies on Gaussian process (GP) for latent space ODE interpolations. Using GPs offers two significant advantages. First, it enables the quantification of uncertainty over the reduced order model (ROM) predictions. Second, leveraging this prediction uncertainty allows for efficient adaptive training through a greedy selection of additional training data points. This approach does not require prior knowledge of the underlying PDEs. Consequently, GPLaSDI is inherently non-intrusive and can be applied to problems without a known PDE or its residual. We demonstrate the effectiveness of our approach on the Burgers equation, Vlasov equation for plasma physics, and a rising thermal bubble problem. Our proposed method achieves between 200 and 100,000 times speed-up, with up to 7% relative error.

Choi, Youngsoo↗

MACHINE LEARNING-ENABLED PREDICTION OF TRANSIENT INJECTION MAP IN AUTOMOTIVE INJECTORS WITH UNCERTAINTY QUANTIFICATION

Accurate prediction of injection profiles is a critical aspect of linking injector operation with engine performance and emissions. However, highly resolved injector simulations can take one to two weeks of wall-clock time, which is incompatible with engine design cycles with desired turnaround times of less than a day. Hence, it is important to reduce the time-to-solution of the internal flow simulations by several orders of magnitude to make it compatible with engine simulations. This work demonstrates a data-driven approach for tackling the computational overhead of injector simulations, whereby the transient injection profiles are emulated for a side-oriented, single-hole diesel injector using a Bayesian machine-learning framework. First, an interpretable Bayesian learning strategy was employed to understand the effect of design parameters on the total void fraction field. Then, autoencoders are utilized for efficient dimensionality reduction of the flowfields. Gaussian process models are finally used to predict the spatiotemporal void fraction field at the injector exit for unknown operating conditions. The Gaussian process models produce principled uncertainty estimates associated with the emulated flowfields, which provide the engine designer with valuable information of where the data-driven predictions can be trusted in the design space. The Bayesian flowfield predictions are compared with the corresponding predictions from a deep neural network, which has been transfer-learned from static needle simulations from a previous work by the authors. The emulation framework can predict the void fraction field at the exit of the orifice within a few seconds, thus achieving a speed-up factor of up to 38 x 10(6) over the traditional simulation-based approach of generating transient injection maps.

machine learning↗

Quantifying mean, variability, and uncertainty in indoor radon exposure in Pennsylvania using random forest and quantile regression forest models

Radon is a naturally occurring radioactive gas that poses a serious health risk as the primary cause of lung cancer in non-smokers. Despite the well-known adverse association with health outcomes, current radon exposure assessments are limited to county-level or average-level estimates, which fail to capture regional variability. This study uses Machine Learning models, including Random Forest (RF) and Quantile Regression Forest (QRF), to estimate the indoor radon concentrations at the ZCTA (Zip code tabulation area)-level and characterize uncertainties in model estimates. Incorporating geological, meteorological, and building-specific data, the models aim to improve radon risk assessment by capturing mean exposure, variability, and extreme concentration levels. Processed radon test data (n = 718,111) were analyzed using average, variability, and quantile prediction methods. Models that estimate the average radon exposure at the ZCTA-level can yield promising model-fit results, but they do not capture the underlying variability of indoor radon exposure within a ZCTA. We utilize volatility analyses to identify characteristics indicative of high variability of indoor radon exposure. We also show that a QRF model can be used to estimate upper quantiles of residential radon exposure, thereby uncovering localized areas of elevated exposure that were not apparent in mean estimates. The results highlighted the need for a deep characterization of exposure risk and show that regions with moderate average exposure levels could still harbor extreme outliers with implications for evaluating health risks. Utilizing multiple radon exposure models allows for a deeper characterization of radon risk within a geographic area and can better identify high-risk areas. The results from this study provide a foundation for developing mitigation strategies and examining associations between radon exposure and health outcomes at fine scales. Future research should extend the geographic scope and incorporate additional environmental risk factors to establish a comprehensive framework for risk assessment.

Lee, Heechan [ORNL]↗

Structural uncertainty assessment for fire-engulfed objects in crosswind: Establishing credibility for a multiphysics wall-modeled large-eddy simulation paradigm

A structural uncertainty validation study for a large-scale, fire-engulfed, elevated object subjected to crosswind is presented to establish the credibility of a high-fidelity, low-Mach, turbulent reacting flow wall-modeled large-eddy simulation (WMLES) approach that includes multiphysics coupling to participating media radiation and conjugate heat transfer. To establish that WMLES can accurately predict surface quantities including drag and pressure coefficient in the low-Mach crosswind regime, a foundational elevated isothermal cylinder validation case is presented at a similar gap-to-diameter ratio of 0.25, spanning the subcritical to supercritical drag regime (Re 𝐷 = 1.1 × 10 5 and 4.3 × 10 5 , respectively). Here, this study exercised both static and dynamic coefficient LES (Smagorinsky and 𝑘 sgs ) with both local and exchange-based velocity sampling. Results showcase that the drag crisis (or the sudden drop in drag coefficient at increased Re 𝐷 ) is well captured when using an exchange-based dynamic coefficient WMLES methodology, while noting lack of mesh convergence and overall drag and pressure coefficient predictively when using a static coefficient, local velocity sampling WMLES. For the 𝒪⁡(10) m JP-8 liquid pool fire crosswind validation study presented, two experimental crosswind configurations (2 m/s and 9.5 m/s) are showcased for a fire-engulfed mock fuselage roughly 4 m in diameter. Using the best model-form practices identified in the isothermal study, dynamic coefficient 𝑘 sgs exchange-based WMLES fire validation findings demonstrate accurate peak irradiation and skin temperature predictions as a function of crosswind magnitude. Excessive yaw in the low-crosswind fuselage configuration, consistent with experimental findings, captured a significant predicted asymmetry in flame attachment and heat flux toward the downwind cylindrical cap—indicative of axial vortex structures transporting the flame along the upper and lower fuselage leeward surface. All fire mesh resolution simulations captured the experimental finding that as crosswind increased, predicted flame shape and peak irradiation magnitude onto the fuselage transitioned from a windward to a leeward cylinder location due to the migration of the upper- to lower-shear fuel/air mixing layer thereby demonstrating the novelty, significance, and credibility of this high-fidelity WMLES reacting flow framework.

Domino, Stefan Paul [Sandia National Laboratories ↗

Toward $\nu_{\mu}$CC $\eta$ Meson Cross Sections on Argon in SBND

The modeling of resonant neutrino interactions on argon is critical for achieving precision neutrino oscillation measurements and searching for physics beyond the Standard Model with both DUNE and the Short-Baseline Neutrino Program. Higher mass resonances beyond the $\Delta$(1232) baryon are particularly poorly constrained experimentally, and yet contribute a significant rate in multi-GeV interactions. The measurement of $\eta$(548) meson production from the decay of these resonances provides a unique probe, as well as the opportunity to measure observables sensitive to final state interactions with the nucleus. This poster will present progress towards the measurement of meson production with the Short-Baseline Near Detector (SBND) at Fermilab, leveraging the clean $\eta \rightarrow \gamma \gamma$ decay channel. The large statistics of the SBND dataset will provide high-precision measurements of this channel capable of constraining uncertainties in the modeling of higher order resonant interactions. This analysis further provides important handles on electromagnetic shower reconstruction performance through evaluation of the 548 MeV $\eta$ invariant mass peak.

Antonakis, Alexander [UC, Santa Barbara]↗

Searching for beyond the Standard Model physics using the improved description of 100 Mo $2\nu \beta \beta$ decay spectral shape with CUPID-Mo

The current experiments searching for neutrinoless double-β ($0\nu \beta \beta$) decay also collect large statistics of Standard Model allowed two-neutrino double-β ($2\nu \beta \beta$ ) decay events. These can be used to search for Beyond Standard Model (BSM) physics via $2\nu \beta \beta$ decay spectral distortions. 100 Mo has a natural advantage due to its relatively short half-life, allowing higher $2\nu \beta \beta$ decay statistics at equal exposures compared to the other isotopes. We demonstrate the potential of the dual read-out bolometric technique exploiting a 100 Mo exposure of 1.47 kg years, acquired in the CUPID-Mo experiment at the Modane underground laboratory (France). We set limits on $0\nu \beta \beta$ decays with the emission of one or more Majorons, on $2\nu \beta \beta$ decay with Lorentz violation, and $2\nu \beta \beta$ decay with a sterile neutrino emission. In this analysis, we investigate the systematic uncertainty induced by modeling the $2\nu \beta \beta$ decay spectral shape parameterized through an improved model, an effect never considered before. This work motivates searches for BSM processes in the upcoming CUPID experiment, which will collect the largest amount of $2\nu \beta \beta$ decay events among the next-generation experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Uncertainty Quantification for Multiphase Computational Fluid Dynamics Closure Relations with a Physics-Informed Bayesian Approach

Multiphase Computational Fluid Dynamics (MCFD) based on the two-fluid model is considered a promising tool to model complex two-phase flow systems. MCFD simulation can predict local flow features without resolving interfacial information. As a result, the MCFD solver relies on closure relations to describe the interaction between the two phases. Those empirical or semi-mechanistic closure relations constitute a major source of uncertainty for MCFD predictions. In this paper, we leverage a physics-informed uncertainty quantification (UQ) approach to inversely quantify the closure relations’ model form uncertainty in a physically consistent manner. This proposed approach considers the model form uncertainty terms as stochastic fields that are additive to the closure relation outputs. Combining dimensionality reduction and Gaussian processes, the posterior distribution of the stochastic fields can be effectively quantified within the Bayesian framework with the support of experimental measurements. As this UQ approach is fully integrated into the MCFD solving process, the physical constraints of the system can be naturally preserved in the UQ results. Here, in a case study of adiabatic bubbly flow, we demonstrate that this UQ approach can quantify the model form uncertainty of the MCFD interfacial force closure relations, thus effectively improving the simulation results with relatively sparse data support.

42 ENGINEERING↗

Fluidization of Group A Glass Particles: Experiments and Preliminary Validation

This work is aimed at providing reliable high-quality data from fluidization experiments. Geldart Group A glass beads are used as bed material. The fluidizing medium is air at atmospheric pressure which enters the system at 40% relative humidity. The air inlet velocity is varied over a broad range from 2 U mf to 26 U mf , where the minimum fluidization velocity, U mf , is 0.608 cm/s. There is not a significant change in the mean values of differential pressure; however, there are noticeable trends in their standard deviation values. Preliminary validation studies are performed using the open-source software Multiphase Flow with Interphase eXchanges (MFiX) developed by the National Energy Technology Laboratory (NETL). The results from Eulerian-Eulerian (MFiX-TFM, two-fluid model) and Eulerian-Lagrangian (MFiX-PIC, particle-in-cell) formulations are compared. The work presented in this report highlights the applicability of MFiX-PIC for large-scale applications, and reveals that PIC simulation offers a suitable alternative to TFM simulation when effects due to multiple components or polydispersity are significant. PIC modeling provides a readily extendible framework for such systems. However, PIC models do rely on an empirical closure for inter-particle stress. Some solution discrepancies were evident in MFiX-PIC predictions which require further analysis. As such, the study underlines the need for a systematic approach to quantify the sources of numerical uncertainty within the PIC model. However, even with known empiricism, MFiX-PIC has demonstrated considerable potential in analyzing gas-solid systems. The results obtained are encouraging and warrant further investigation to improve the predictive capability of MFiX-PIC.

36 MATERIALS SCIENCE↗

Heavy-Duty Low-Temperature and Diesel Combustion & Heavy-Duty Combustion Modeling (FY 2018 Annual Progress Report)

Regulatory drivers and market demands for lower pollutant emissions, lower carbon dioxide emissions, and lower fuel consumption motivate the development of clean and fuel-efficient engine operating strategies. Most current production engines use a combination of both in-cylinder and exhaust emission control strategies to achieve these goals. The emissions and efficiency performance of in-cylinder strategies depend strongly on flow and mixing processes associated with fuel injection. Both heavy- and light-duty engine/vehicle manufacturers use multiple-injection strategies to reduce noise, emissions, and fuel consumption. For both conventional and low-temperature diesel combustion, the state of knowledge and modeling tools for multiple injections are far less advanced than for single-injection strategies. Engine efficiency is limited to some degree by tradeoffs that must be accepted to meet particulate matter (including soot) emissions limits. Recent work on this project has filled some knowledge gaps on soot oxidation with multiple injections, and the current work for Fiscal Year (FY) 2018 addresses knowledge gaps on soot formation for multiple injections. While total in-cylinder soot is readily measured, discerning formation from oxidation is difficult. The FY 2018 experiments are designed to create in-cylinder conditions at the threshold of soot formation, where processes that affect soot formation can be more readily discerned. Soot formation pathways under such conditions are fraught with uncertainties, and soot models significantly overpredict polyaromatic hydrocarbon (PAH) and soot, so experimental data at these conditions will provide much needed data for improvements to PAH and soot models.

02 PETROLEUM↗