Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “models”

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 811 records · Page 45

DiffESM: Conditional Emulation of Temperature and Precipitation in Earth System Models With 3D Diffusion Models

Earth system models (ESMs) are essential for understanding the interaction between human activities and the Earth's climate. However, the computational demands of ESMs often limit the number of simulations that can be run, hindering the robust analysis of risks associated with extreme weather events. While low-cost climate emulators have emerged as an alternative to emulate ESMs and enable rapid analysis of future climate, many of these emulators only provide output on at most a monthly frequency. This temporal resolution is insufficient for analyzing events that require daily characterization, such as heat waves or heavy precipitation. We propose using diffusion models, a class of generative deep learning models, to effectively downscale ESM output from a monthly to a daily frequency. Trained on a handful of ESM realizations, reflecting a wide range of radiative forcings, our DiffESM model takes monthly mean precipitation or temperature as input, and is capable of producing daily values with statistical characteristics close to ESM output. Combined with a low-cost emulator providing monthly means, this approach requires only a small fraction of the computational resources needed to run a large ensemble. We evaluate model behavior using a number of extreme metrics, showing that DiffESM closely matches the spatio-temporal behavior of the ESM output it emulates in terms of the frequency and spatial characteristics of phenomena such as heat waves, dry spells, or rainfall intensity.

54 ENVIRONMENTAL SCIENCES↗

The GFDL‐CM4X Climate Model Hierarchy, Part I: Model Description and Thermal Properties

We present the GFDL‐CM4X (Geophysical Fluid Dynamics Laboratory Climate Model version 4X) coupled climate model hierarchy. The primary application for CM4X is to investigate ocean and sea ice physics as part of a realistic coupled Earth climate model. CM4X utilizes an updated MOM6 (Modular Ocean Model version 6) ocean physics package relative to CM4.0, and there are two members of the hierarchy: one that uses a horizontal grid spacing of 0.25° (referred to as CM4X‐p25) and the other that uses a grid 0.125° (CM4X‐p125). CM4X also refines its atmospheric grid from the nominally 100 km (cubed sphere C96) of CM4.0–50 km (C192). Finally, CM4X simplifies the land model to allow for a more focused study of the role of ocean changes to global mean climate. CM4X‐p125 reaches a global ocean area mean heat flux imbalance of -0.02 W m -2 within $\mathcal{O}$ (150) years in a pre‐industrial simulation, and retains that thermally equilibrated state over the subsequent centuries. This 1850 thermal equilibrium is characterized by roughly 400 ZJ less ocean heat than present‐day, which corresponds to estimates for anthropogenic ocean heat uptake between 1870 and present‐day. CM4X‐p25 approaches its thermal equilibrium only after more than 1000 years, at which time its ocean has roughly 1100 ZJ more heat than its early 21st century ocean initial state. Furthermore, the root‐mean‐square sea surface temperature bias for historical simulations is roughly 20% smaller in CM4X‐p125 relative to CM4X‐p25 (and CM4.0). We offer the mesoscale dominance hypothesis for why CM4X‐p125 shows such favorable thermal equilibration properties.

54 ENVIRONMENTAL SCIENCES↗

Calibration of the 2-Phase Bubble Tracking Model for Liquid Mercury Target Simulation with Machine Learning Surrogate Models

The Spallation Neutron Source (SNS) at Oak Ridge National Laboratory is one of the most powerful accelerator-driven neutron sources in the world. The intense protons strike on SNS's mercury target to provide bright neutron beams, which also leads to severe fluid-structure interactions inside the target. Prediction of resultant loading on the target is difficult particularly when helium gas is injected into mercury to reduce the loading and mitigate the pitting damage on vessel walls. A 2-phase material model that incorporates the Rayleigh-Plesset (R-P) model is expected to address this multi-physics problem. However, several uncertain parameters in the R-P model require intensive simulations to determine their optimal values. With the help of machine learning and the measured target strain, we have studied the major uncertain parameters in this R-P model and developed a framework to identify optimal parameters that significantly reduce the discrepancy between simulations and experimental strains. The preliminary results show the possibility of using this mercury/helium mixture and surrogate models to predict a better match of target strain response when the helium gas is injected.

Lin, Lianshan↗

Masked Particle Modeling on Sets: Towards Self-Supervised High Energy Physics Foundation Models

Abstract We propose masked particle modeling (MPM) as a self-supervised method for learning generic, transferable, and reusable representations on unordered sets of inputs for use in high energy physics (HEP) scientific data. This work provides a novel scheme to perform masked modeling based pre-training to learn permutation invariant functions on sets. More generally, this work provides a step towards building large foundation models for HEP that can be generically pre-trained with self-supervised learning and later fine-tuned for a variety of down-stream tasks. In MPM, particles in a set are masked and the training objective is to recover their identity, as defined by a discretized token representation of a pre-trained vector quantized variational autoencoder. We study the efficacy of the method in samples of high energy jets at collider physics experiments, including studies on the impact of discretization, permutation invariance, and ordering. We also study the fine-tuning capability of the model, showing that it can be adapted to tasks such as supervised and weakly supervised jet classification, and that the model can transfer efficiently with small fine-tuning data sets to new classes and new data domains.

Heinrich, Lukas (ORCID:0000000240487584)↗

Correlations between azimuthal anisotropy Fourier harmonics in PbPb collisions at sNN =2.76 TeV in the HYDJET++ model and in the multiphase transport model

Correlations between azimuthal anisotropy Fourier harmonics vn (n=2,3,4) are studied using the events from PbPb collisions at sNN=2.76 TeV generated by the HYDJET++ and multiphase transport (AMPT) models, and compared to the corresponding experimental results obtained by the ATLAS Collaboration. The Fourier harmonics vn are measured over a wide centrality range using the two-particle azimuthal correlation method. The slopes of the v2-v3 correlation from both models are in a good agreement with the ATLAS data. The HYDJET++ model predicts a stronger slope for the v2-v4 and v3-v4 correlations than the ones experimentally measured, while the results from the AMPT model are in a rather good agreement with the experimental results. In contrast to the HYDJET++ predictions, the AMPT model predicts a boomeranglike shape in the structure of the correlations as found in the experimental data. © 2020 American Physical Society.

60 APPLIED LIFE SCIENCES↗

Model-independent versus model-dependent interpretation of the SDSS-III BOSS power spectrum: Bridging the divide

The traditional clustering analyses of galaxy redshift surveys compress the clustering data into a set of late-time physical variables in a model-independent way. This approach has recently been extended by an additional shape variable encoding early-time physics information. Here we apply this new technique, ShapeFit, to SDSS-III BOSS data and show that it matches the constraining power of alternative, model-dependent approaches, which directly constrain the model’s parameters adopting a cosmological model ab initio. ShapeFit is ~30 times faster, model independent, naturally splits early- and late-time variables, and enables a better control of observational systematics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Modeling thin layers of analytes on substrates for spectral analysis: use of solid/liquid n and k values to model reflectance spectra

Since solids are only sometimes seen en masse in a pure bulk form, and for liquids other than water almost never, a capability to model reflectance spectra from analytes deposited on various substrates would be highly advantageous. If available, the real, n ( ν ) , and imaginary, k ( ν ) , components of the complex refractive index, n n + i k , can be used to simulate infrared spectra, accounting for reflection, refraction, and absorption phenomena as a function of wavelength. We focus on using the Pacific Northwest National Laboratory (PNNL) derived n / k vectors for solid and liquid analytes deposited as thin layers on different types of substrates including conductors, such as aluminum, and inorganic dielectrics, such as glass. The model is an adaptation of the Monte Carlo ray trace modeling program, TracePro, extended through the use of its macrolanguage. The model is tested using thin films of organic liquids including silicone oil and no. 2 diesel fuel, as well as organic solids such as caffeine and acetaminophen on aluminum and glass. The predicted spectra for the solid films were compared to experimental hemispherical reflectance data measured using a Fourier transform spectrometer with an integrating sphere. The thickness of the calculated layer is a parameter for predicting the (transflectance) spectra and is obtained using the areal density measured from gravimetric methods to generate the thin-layer samples. Comparison of the calculated spectra with experimental hemispherical reflectance data shows excellent agreement, indicating promise for the use of measured n / k data to synthesize reference spectral data. In particular, accounting for the inhomogeneity of the deposits greatly improved the match with experimental data. Finally, the theoretical modeling shows that for thicker layers (ca. 20 to 100 μm) of typical organics possessing moderately strong k values, the longwave infrared features are often saturated and better spectral contrast is obtained from the overtone/combination bands in the shortwave infrared.

42 ENGINEERING↗

Comparative analysis of model-free and model-based HVAC control for residential demand response

In this paper, we present a comparative analysis of model-free reinforcement learning (RL) and model predictive control (MPC) approaches for intelligent control of heating, ventilation, and air-conditioning (HVAC). Deep-Q-network (DQN) is used as a candidate for model-free RL algorithm. The two control strategies were developed for residential demand-response (DR) HVAC system. We considered MPC as our golden standard to compare DQN's performance. The question we tried to answer through this work was, What % of MPC's performance can be achieved by model-free RL approach for intelligent HVAC control?. Based on our test result, RL achieved an average of ≈ 62% daily cost saving of MPC. Considering the pure optimization and model-based nature of MPC methods, the RL showed very promising performance. We believe that the interpretations derived from this comparative analysis provide useful insights to choose from various DR approaches and further enhance the performance of the RL-based methods for building energy managements.

Kurte, Kuldeep↗

Modeling Plant Metabolism: From Network Reconstruction to Mechanistic Models

Mathematical modeling of plant metabolism allows the plant science community to understand the organization of plant metabolism, obtain quantitative insights into metabolic functions, and derive engineering strategies for manipulation of metabolism. Among the various modeling approaches, metabolic pathway analysis can dissect the basic functional modes of subsections of core metabolism, such as photorespiration, and reveal how classical definitions of metabolic pathways have overlapping functionality. In current research using constraint-based modeling in plants, numerous computational tools are currently available to analyze large-scale and genome-scale metabolic networks. For 13 C-metabolic flux analysis, principles of isotopic steady state have been used to study heterotrophic plant tissues, while nonstationary isotope labeling approaches are amenable to the study of photoautotrophic and secondary metabolism. Enzyme kinetic models explore pathways in mechanistic detail, and we discuss different approaches to determine or estimate kinetic parameters. In this review, we describe recent advances and challenges in modeling plant metabolism.

59 BASIC BIOLOGICAL SCIENCES↗

Model-Data for Joint Estimation of Biogeochemical Model Parameters from Multiple Experiments: A Bayesian Approach Applied to Mercury Methylation

This modeling archive supports the manuscript submitted for publication in the Environmental Modeling and Software. This study is supported by ORNL-SFA and IDEAS-Watershed. This study aims to improve calibration of complex biogeochemical models using datasets from multiple experiments targeting specific subprocesses. The proposed Bayesian joint-fitting scheme calibrates the entire biogeochemical model in one go using all the available datasets and estimate parameter uncertainties using Markov Chain Monte Carlo (MCMC). This allows for complete propagation of uncertainties and utilization of the information shared between different datasets. Mapping joint distribution of parameters guides model improvement by identifying null spaces in the parameter space. This archive contains files used to perform MCMC, post-process outputs and visualize results.

East Fork Poplar Creek↗

Parallel simulated annealing with embedded machine learning and multifidelity models for reactor core design

This paper presents extensions to a penalty-free, parallel simulated annealing (SA) algorithm for multi-constrained combinatorial optimization with the aim of embedding multi-fidelity physics models into the annealing procedure. The method uses a low-fidelity, quickly executing model for rapid design space exploration and a high-fidelity model for detailed constraint resolution and on-the-fly bias correction. Machine learning models updated within the annealing procedure were used to bridge the gap between the multi-fidelity models, which led to accurate rapid exploration and efficient detailed constraint resolution. A software implementation of the new multi-fidelity optimization methods, called ML-PSA, was demonstrated on a continuous multi-fidelity optimization problem and a constrained combinatorial PWR lattice design problem. These problems demonstrate some of the features, parallel performance characteristics, and extensible nature of the multi-fidelity SA methods. This paper shows that the developed software and procedure are a general optimization tool that can be applied to a wide variety of scientific and engineering design optimization applications. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

BISON BWR Fuel Modeling Capability: Material Models (FY2020)

Until this year the BISON fuel performance code has been developed primarily for pressurized water reactor (PWR), accident tolerant fuel (ATF), and advanced reactor (e.g., tri-structural isotropic (TRISO), metallic) analyses. The main focus of this work is to develop and document required material models that would aid predicting cladding failure with a specific hydride distribution to be investigated for boiling water reactor (BWR) applications in the next fiscal year. In terms of materials, the primary differences between a BWR and PWR is the use of gadolinia doping of the UO 2 fuel, the use of Zircaloy-2 for the cladding, and the addition of a liner on the inner surface of the cladding for improved pellet-clad mechanical interaction (PCMI)performance. The liner is typically made of pure zirconium or a lower tin content zirconium alloy (e.g. Zr-0.3%Sn). This work details the development and documentation of material and behavioral models for BWR analysis such as the addition of material and behavior models for gadolinia doped UO 2 , pure zirconium, Zr-0.3%Sn, and Zircaloy-2 as well as the ease-of-use additions to allow internal meshing capabilities to include liners. Another important aspect affecting cladding performance is hydrogen diffusion and hydride precipitation in the cladding. A model was added previously for non-lined, homogeneous claddings. In this work, instructions are provided on how to apply the model to BWR liner cladding. Representative experimental values are provided for the hydrogen diffusion parameters and precipitation-dissolution kinetics for typical liner materials of interest.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

FECM/NETL CO 2 Transport Cost Model (2022): Model Overview [Slides]

The FECM/NETL CO 2 Transport Cost Model (CO 2 _T_COM) is an Excel spreadsheet model that calculates the cost of transporting CO 2 from the beginning to the end of a pipeline. This document provides an overview of the model. This document is intended to provide a user with a briefer description of the model and how to use the model than the user’s manual.

42 ENGINEERING↗

The Schwarz Alternating Method for the Seamless Coupling of Nonlinear Reduced Order Models and Full Order Models

Projection-based model order reduction allows for the parsimonious representation of full order models (FOMs), typically obtained through the discretization of a set of partial differential equations (PDEs) using conventional techniques (e.g., finite element, finite volume, finite difference methods) where the discretization may contain a very large number of degrees of freedom. As a result of this more compact representation, the resulting projection-based reduced order models (ROMs) can achieve considerable computational speedups, which are especially useful in real-time or multi-query analyses. One known deficiency of projection-based ROMs is that they can suffer from a lack of robustness, stability and accuracy, especially in the predictive regime, which ultimately limits their useful application. Another research gap that has prevented the widespread adoption of ROMs within the modeling and simulation community is the lack of theoretical and algorithmic foundations necessary for the “plug-and-play” integration of these models into existing multi-scale and multi-physics frameworks. This paper describes a new methodology that has the potential to address both of the aforementioned deficiencies by coupling projection-based ROMs with each other as well as with conventional FOMs by means of the Schwarz alternating method [41]. Leveraging recent work that adapted the Schwarz alternating method to enable consistent and concurrent multiscale coupling of finite element FOMs in solid mechanics [35, 36], we present a new extension of the Schwarz framework that enables FOM-ROM and ROM-ROM coupling, following a domain decomposition of the physical geometry on which a PDE is posed. In order to maintain efficiency and achieve computation speed-ups, we employ hyper-reduction via the Energy-Conserving Sampling and Weighting (ECSW) approach [13]. We evaluate the proposed coupling approach in the reproductive as well as in the predictive regime on a canonical test case that involves the dynamic propagation of a traveling wave in a nonlinear hyper-elastic material.

97 MATHEMATICS AND COMPUTING↗

Integrated Model of Models for Global Flood Alerting

A dramatic increase in frequency of minor to major flooding since 2000 has caused significant economic losses across the world. To mitigate and recover from these losses, actions have been taken to build resilient communities and infrastructures, specifically, by providing situational awareness in near real-time about flood impacts to enhance response and recovery efforts. Several hydrologic and hydraulic flood models are available at various spatial and temporal resolutions to forecast flood events at regional to global scale. Given the global coverage of two operational flood models – GloFAS (Global Flood Awareness System) and GFMS (Global Flood Monitoring System), the purpose of this project is to implement a Model of Models (MoM) approach to integrate the outputs from these two models to classify flood severity at watershed level worldwide, and send alerts based on severity similar to the USGS PAGER (used for severity alerting and impact analysis for earthquakes) to flood impacted communities. The alerts containing flood impacts and severity information will be disseminated through the DisasterAWARE platform, operated by the Pacific Disaster Center (PDC), that provides global multi-hazard alerting and Situational Awareness information to the emergency management community and public. The current version of the MoM approach was implemented for a case study flood event that occurred during January and February of 2020 in South and Central Africa. The findings of the case study event reveal that the approach is effective in identifying potential flood impact areas and the spatio-temporal variation of flood severity, flood depth and extent at watershed level, which will be used to assess infrastructure and societal impacts using earth-observation data and for alerting.

Kar, Bandana↗

Modeling Sepsis: Establishment and Validation of a 72-Hour Swine Model of Penetrating Abdominal Trauma

Background/Objectives: Fecal peritonitis following penetrating abdominal trauma is a serious condition that often results in sepsis and organ failure. The aim of our study was to develop a novel conscious porcine model of sepsis and organ dysfunction caused by multiple penetrating injuries to the small and large intestines. Methods: Twelve female Yorkshire pigs (average weight 50.6 ± 6.5 kg) were divided into two groups: Penetrating Abdominal Trauma (PAT) (n = 8) and Control (n = 4). All surgical procedures were performed under anesthesia with adequate analgesia. In the PAT group, the small and large intestines were punctured, and feces mixed with saline were introduced into the abdominal cavity to induce peritonitis. The Control group received sham surgery with only saline solution. The animals were observed in a conscious state over a period of 72 h, vital parameters were recorded, and blood samples were taken regularly. We adapted a pig-specific SOFA score and developed pig-specific SIRS criteria and NEWS2 score to assess organ function. The model was validated by independent investigators. Results: The survival rate in the PAT group was 75%, with an average survival time of 58.5 h, while all animals in the Control group survived to euthanasia. Monitoring showed pathophysiological changes, such as tachycardia, leucopenia, and thrombocytopenia, indicative of sepsis and organ dysfunction. Blinded investigators independently confirmed the model’s validity. Conclusions: A new swine model of penetrating abdominal trauma and sepsis has been successfully developed that demonstrates significant physiological and immunologic changes comparable to human sepsis. This new model provides a realistic platform for future research into sepsis, its diagnostics, and the evaluation of therapeutic strategies.

General & Internal Medicine↗