Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SYSTEMS MODELING”

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

On the Power of Uncertainties in Microbial System Modeling: No Need To Hide Them Anymore

For decades, microbiologists have considered uncertainties as an undesired side effect of experimental protocols. As a consequence, standard microbial system modeling strives to hide uncertainties for the sake of deterministic understanding. However, recent studies have highlighted greater experimental variability than expected and emphasized uncertainties not as a weakness but as a necessary feature of complex microbial systems. We therefore advocate that biological uncertainties need to be considered foundational facets that must be incorporated in models. Not only will understanding these uncertainties improve our understanding and identification of microbial traits, it will also provide fundamental insights on microbial systems as a whole. Taking into account uncertainties within microbial models calls for new validation techniques. Formal verification already overcomes this shortcoming by proposing modeling frameworks and validation techniques dedicated to probabilistic models. However, further work remains to extract the full potential of such techniques in the context of microbial models. Herein, we demonstrate how statistical model checking can enhance the development of microbial models by building confidence in the estimation of critical parameters and through improved sensitivity analyses.

59 BASIC BIOLOGICAL SCIENCES↗

An Open Simulation System Model for Scientific Applications

A model for a generic and open environment for running multi-code or multi-application simulations - called the open Simulation System Model (OSSM) - is proposed and defined. This model attempts to meet the requirements of complex systems like the Numerical Propulsion Simulator System (NPSS). OSSM places no restrictions on the types of applications that can be integrated at any state of its evolution. This includes applications of different disciplines, fidelities, etc. An implementation strategy is proposed that starts with a basic prototype, and evolves over time to accommodate an increasing number of applications. Potential (standard) software is also identified which may aid in the design and implementation of the system.

Williams, Anthony D.↗

Recreating the California New Year's Flood Event of 1997 in a Regionally Refined Earth System Model

Abstract The 1997 New Year's flood event was the most costly in California's history. This compound extreme event was driven by a category 5 atmospheric river that led to widespread snowmelt. Extreme precipitation, snowmelt, and saturated soils produced heavy runoff causing widespread inundation in the Sacramento Valley. This study recreates the 1997 flood using the Regionally Refined Mesh capabilities of the Energy Exascale Earth System Model (RRM‐E3SM) under prescribed ocean conditions. Understanding the processes causing extreme events informs practical efforts to anticipate and prepare for such events in the future, and also provides a rich context to evaluate model skill in representing extremes. Three California‐focused RRM grids, with horizontal resolution refinement of 14 km down to 3.5 km, and six forecast lead times, 28 December 1996 at 00Z through 30 December 1996 at 12Z, are assessed for their ability to recreate the 1997 flood. Planetary to synoptic scale atmospheric circulations and integrated vapor transport are weakly influenced by horizontal resolution refinement over California. Topography and mesoscale circulations, such as the Sierra barrier jet, are better represented at finer horizontal resolutions resulting in better estimates of storm total precipitation and storm duration snowpack changes. Traditional time‐series and causal analysis frameworks are used to examine runoff sensitivities state‐wide and above major reservoirs. These frameworks show that horizontal resolution plays a more prominent role in shaping reservoir inflows, namely the magnitude and time‐series shape, than forecast lead time, 2‐to‐4 days prior to the 1997 flood onset.

54 ENVIRONMENTAL SCIENCES↗

Characterization of Extreme Hydroclimate Events in Earth System Models using ML/AI

Focal Area(s): (1) We put forward the concepts of data assimilation enabled by machine learning, AI, and advanced methods including experimental/network design/optimization and unsupervised learning applied to downscale information within Earth System Models (ESMs). (2) We discuss predictive modeling through the use of AI techniques and other tools to design a prediction system comprising of a hierarchy of models (e.g., AI-driven model/component/parameterization selection) to improve the characterization of extreme hydroclimate events in ESMs. The AI-based models will run five-six order of magnitude times faster, yet will provide similar accuracy, allowing us to provide range bounds on uncertainty faster and thus enabling faster extreme event identification. (3) Further, we consider the insight gleaned from complex data (observed/simulated) using AI, big data analytics, and other advanced methods, including explainable AI and physics- or knowledge-guided AI to improve the characterization of extreme hydroclimate events in ESMs.

54 ENVIRONMENTAL SCIENCES↗

Clay Generic Disposal System Model - Sensitivity Analysis for 32 PWR Assembly Canisters (+2 associated model files).

The Used Fuel Disposition Campaign (UFDC), as part of the DOE Office of Nuclear Energy’s (DOE-NE) Fuel Cycle Technology program (FCT) is investigating the disposal of high level radioactive waste (HLW) and spent nuclear fuela (SNF) in a variety of geologic media. The feasibility of disposing SNF and HLW in clay media has been investigated and has been shown to be promising [Ref. 1]. In addition the disposal of these wastes in clay media is being investigated in Belgium, France, and Switzerland. Thus, Argillaceous media is one of the environments being considered by UFDC. As identified by researchers at Sandia National Laboratory, potentially suitable formations that may exist in the U.S. include mudstone, clay, shale, and argillite formations [Ref. 1]. These formations encompass a broad range of material properties. In this report, reference to clay media is intended to cover the full range of material properties. This report presents the status of the development of a simulation model for evaluating the performance of generic clay media. The clay Generic Disposal System Model (GDSM) repository performance simulation tool has been developed with the flexibility to evaluate not only different properties, but different waste streams/forms and different repository designs and engineered barrier configurations/ materials that could be used to dispose of these wastes.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Structural dynamics system model reduction

Loads analysis for structural dynamic systems is usually performed by finite element models. Because of the complexity of the structural system, the model contains large number of degree-of-freedom. The large model is necessary since details of the stress, loads and responses due to mission environments are computed. However, a simplified model is needed for other tasks such as pre-test analysis for modal testing, and control-structural interaction studies. A systematic method of model reduction for modal test analysis is presented. Perhaps it will be of some help in developing a simplified model for the control studies.

Chen, J. C.↗

Joint emulation of Earth System Model temperature-precipitation realizations with internal variability and space-time and cross-variable correlation: fldgen v2.0 software description

Earth System Models (ESMs) are excellent tools for quantifying many aspects of future climate dynamics but are too computationally expensive to produce large collections of scenarios for downstream users of ESM data. In particular, many researchers focused on the impacts of climate change require large collections of ESM runs to rigorously study the impacts to both human and natural systems of low-frequency high-importance events, such as multi-year droughts. Climate model emulators provide an effective mechanism for filling this gap, reproducing many aspects of ESMs rapidly but with lower precision. The fldgen v1.0 R package quickly generates thousands of realizations of gridded temperature fields by randomizing the residuals of pattern scaling temperature output from any single ESM, retaining the spatial and temporal variance and covariance structures of the input data at a low computational cost. The fldgen v2.0 R package described here extends this capability to produce joint realizations of multiple variables, with a focus on temperature and precipitation in an open source software package available for community use (https://github.com/jgcri/fldgen). This substantially improves the fldgen package by removing the requirement that the ESM variables be normally distributed, and will enable researchers to quickly generate covarying temperature and precipitation data that are synthetic but faithful to the characteristics of the original ESM.

54 ENVIRONMENTAL SCIENCES↗

Surrogate multi-fidelity data and model fusion for scientific discovery and uncertainty quantification in Earth System Models

This whitepaper addresses the Earth and Environmental Systems Sciences Division (EESSD)’s predictability challenges in modeling the integrated water cycle and data-model integration. Specifically, it focuses on reducing and characterizing the uncertainty in the representation of process models for unresolved physics, either due to model resolution or limited by the physical under standing or computational efficiency, and the use of observational data for in-situ process parameter optimization within ESM. The described methods may also be used to determine the nature of responses (e.g. strength and direction), and hence to identify critical processes that drive the overall ESM responses to perturbation in the forcing

54 ENVIRONMENTAL SCIENCES↗

Application of Generic Disposal System Models

This report describes specific GDSA activities in fiscal year 2015 (FY2015) toward the development of the enhanced disposal system modeling and analysis capability for geologic disposal of nuclear waste. The GDSA framework employs the PFLOTRAN thermal-hydrologic-chemical multi-physics code (Hammond et al., 2011) and the Dakota uncertainty sampling and propagation code (Adams et al., 2013). Each code is designed for massively-parallel processing in a high-performance computing (HPC) environment. Multi-physics representations in PFLOTRAN are used to simulate various coupled processes including heat flow, fluid flow, waste dissolution, radionuclide release, radionuclide decay and ingrowth, precipitation and dissolution of secondary phases, and radionuclide transport through the engineered barriers and natural geologic barriers to a well location in an overlying or underlying aquifer. Dakota is used to generate sets of representative realizations and to analyze parameter sensitivity.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Surface analysis of model systems: From a metal-graphite interface to an intermetallic catalyst

This thesis summarizes research completed on two different model systems. In the first system, we investigate the deposition of the elemental metal dysprosium on highly-oriented pyrolytic graphite (HOPG) and its resulting nucleation and growth. The goal of this research is to better understand the metal-carbon interactions that occur on HOPG and to apply those to an array of other carbon surfaces. This insight may prove beneficial to developing and using new materials for electronic applications, magnetic applications and catalysis.

36 MATERIALS SCIENCE↗

Surface analysis of model systems: From a metal-graphite interface to an intermetallic catalyst

This thesis summarizes research completed on two different model systems. In the first system, we investigate the deposition of the elemental metal dysprosium on highlyoriented pyrolytic graphite (HOPG) and its resulting nucleation and growth. The goal of this research is to better understand the metal-carbon interactions that occur on HOPG and to apply those to an array of other carbon surfaces. This insight may prove beneficial to developing and using new materials for electronic applications, magnetic applications and catalysis. In the second system, we investigate the intermetallic single crystal NaAu2. We conducted a characterization of the clean (111) surface along with a determination of the effects of gas exposures on the surface. No surfaces of bulk NaAu2 have been previously characterized, although the closely related system of Na films on Au(111) has been investigated. NaAu2 is an active catalyst for CO oxidation. To better understand the catalytic reaction occurring, we studied the NaAu2 surface after exposure to O 2 ,g, CO g , CO 2,g ,and H 2 O g . The goal of this research is to find a possible alternative for supported gold catalysts by utilizing an Au-rich intermetallic compound and to determine how its catalytic reaction proceeds. This insight may be valuable in determining other catalytic alternatives as well as helping understand catalytic reactions that occur with other intermetallic compounds. In order to characterize and examine these systems, we utilized several ultrahigh vacuum (UHV) techniques, including scanning tunneling microscopy (STM), X-ray photoelectron spectroscopy (XPS) and low-energy electron diffraction (LEED).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evaluation of SO2, SO4(2−) and an Updated SO2 Dry Deposition Parameterization in the United Kingdom Earth System Model

In this study we evaluate simulated surface SO2 and sulphate (SO4(2-)) concentrations from the United Kingdom Earth System Model (UKESM1) against observations from ground based measurement networks in the USA and Europe for the period 1987 to 2014. We find that UKESM1 captures the historical trend for decreasing concentrations of atmospheric SO2 and SO4(2-) in both Europe and the USA over the period 1987 to 2014. However, in the polluted regions of the eastern USA and Europe, UKESM1 over-predicts surface SO2 concentrations by a factor of 3, while under-predicting surface SO4(2-) concentrations by 25-35%. In the cleaner western USA, the model over-predicts both surface SO2 and SO4(2-) concentrations by a factor of 12 and 1.5 respectively. We find that UKESM1’s bias in surface SO2 and SO4(2-) concentrations is variable according to region and season. We also evaluate UKESM1 against total column SO2 from the Ozone Monitoring Instrument (OMI) using an updated data product. This comparison provides information about the model’s global performance, finding that UKESM1 over predicts total column SO2 over much of the globe, including the large source regions of India, China, the USA and Europe as well as over outflow regions. Finally, we assess the impact of a more realistic treatment of the model’s SO2 dry deposition parameterization. This change increases SO2 dry deposition to the land and ocean surfaces, thus reducing the atmospheric loading of SO2 and SO(sup 2-)(sub 4). In comparison with the ground-based and satellite observations, we find that the modified parameterization reduces the models over prediction of surface SO2 concentrations and total column SO2. Relative to the ground-based observations the simulated surface SO4(2-) concentrations are also reduced, while the simulated SO2 dry deposition fluxes increase.

SO2↗

Parallel, iterative solution of sparse linear systems - Models and architectures

Solving large, sparse, linear systems of equations is a fundamental problem in large scale scientific and engineering computation. A model of a general class of asynchronous, iterative solution methods for linear systems is developed. In the model, the system is solved by creating several cooperating tasks that each compute a portion of the solution vector. A data transfer model predicting both the probability that data must be transferred between two tasks and the amount of data to be transferred is presented. This model is used to derive an execution time model for predicting parallel execution time and an optimal number of tasks given the dimension and sparsity of the coefficient matrix and the costs of computation, synchronization, and communication. The suitability of different parallel architectures for solving randomly sparse linear systems is discussed. Based on the complexity of task scheduling, one parallel architecture, based on a broadcast bus, is presented and analyzed.

Reed, D. A.↗

Reconciling and Improving Formulations for Thermodynamics and Conservation Principles in Earth System Models (ESMs)

Abstract This paper provides a comprehensive derivation of the total energy equations for the atmospheric components of Earth System Models (ESMs). The assumptions and approximations made in this derivation are motivated and discussed. In particular, it is emphasized that closing the energy budget is conceptually challenging and hard to achieve in practice without resorting to ad hoc fixers. As a concrete example, the energy budget terms are diagnosed in a realistic climate simulation using a global atmosphere model. The largest total energy errors in this example are spurious dynamical core energy dissipation, thermodynamic inconsistencies (e.g., coupling parameterizations with the host model) and missing processes/terms associated with falling precipitation and evaporation (e.g., enthalpy flux between components). The latter two errors are not, in general, reduced by increasing horizontal resolution. They are due to incomplete thermodynamic and dynamic formulations. Future research directions are proposed to reconcile and improve thermodynamics formulations and conservation principles.

54 ENVIRONMENTAL SCIENCES↗

Characterizing the Variation and Covariation of Cloud Microphysical Properties and Implications for Simulation of Subgrid-scale Warm-Rain Processes in Earth System Models (Final DOE-ASR Report)

Warm marine boundary layer (MBL) clouds constitute an important component in the global climate system, and precipitation plays a central role in controlling the water budget, radiative effects, and lifetime of these MBL clouds. Unfortunately, because of the relatively coarse effective grid resolution of the current generation of Earth system models (ESMs), the variety of cloud microphysical processes occurring inside an ESM grid cell are often oversimplified or unconstrained by observations. For example, the warm rain processes (e.g., autoconversion and accretion) are usually parameterized as nonlinear functions of grid-mean cloud properties. Because of the nonlinear nature of these functions, neglecting variability within the ESM grid volume can lead to substantial biases in precipitation production, cloud cover, and surface radiative fluxes. In state-of-the art ESMs, the influence of subgrid-scale variability is represented as an enhancement factor (EF) coefficient to the autoconversion, and accretion rates calculated from the model variables. However, EF is typically taken to be a constant or even used as a knob to tune model cloud properties to match observations, an ad hoc approach that may yield a desired cloud outcome yet introduce compensating errors. In this project, we used the combination of in situ cloud microphysics measurements from the ACE-ENA field campaign and large-eddy simulations (LES) to characterize and understand subgrid-scale variations and co-variations of cloud microphysical properties and use the results to evaluate and improve the representation of subgrid warm-rain processes in ESMs, in particular the EF used to tune the autoconversion and accretion processes. In this final report, we summarize our research activities and main findings in Section 2, provide a list of publications (Section 3) and presentations (Section 4) resulted from our research, and briefly discuss the student activities supported by this project.

54 ENVIRONMENTAL SCIENCES↗

Improved methodologies for Earth system modelling of atmospheric soluble iron and observation comparisons using the Mechanism of Intermediate complexity for Modelling Iron (MIMI v1.0)

Abstract. Herein, we present a description of the Mechanism of Intermediate complexity for Modelling Iron (MIMI v1.0). This iron processing module was developed for use within Earth system models and has been updated within a modal aerosol framework from the original implementation in a bulk aerosol model. MIMI simulates the emission and atmospheric processing of two main sources of iron in aerosol prior to deposition: mineral dust and combustion processes. Atmospheric dissolution of insoluble to soluble iron is parameterized by an acidic interstitial aerosol reaction and a separate in-cloud aerosol reaction scheme based on observations of enhanced aerosol iron solubility in the presence of oxalate. Updates include a more comprehensive treatment of combustion iron emissions, improvements to the iron dissolution scheme, and an improved physical dust mobilization scheme. An extensive dataset consisting predominantly of cruise-based observations was compiled to compare to the model. The annual mean modelled concentration of surface-level total iron compared well with observations but less so in the soluble fraction (iron solubility) for which observations are much more variable in space and time. Comparing model and observational data is sensitive to the definition of the average as well as the temporal and spatial range over which it is calculated. Through statistical analysis and examples, we show that a median or log-normal distribution is preferred when comparing with soluble iron observations. The iron solubility calculated at each model time step versus that calculated based on a ratio of the monthly mean values, which is routinely presented in aerosol studies and used in ocean biogeochemistry models, is on average globally one-third (34 %) higher. We redefined ocean deposition regions based on dominant iron emission sources and found that the daily variability in soluble iron simulated by MIMI was larger than that of previous model simulations. MIMI simulated a general increase in soluble iron deposition to Southern Hemisphere oceans by a factor of 2 to 4 compared with the previous version, which has implications for our understanding of the ocean biogeochemistry of these predominantly iron-limited ocean regions.

54 ENVIRONMENTAL SCIENCES↗

Systematic Underestimation of Canopy Conductance Sensitivity to Drought by Earth System Models

The response of vegetation canopy conductance (g c ) to changes in moisture availability ($γ^{m}_{gc}$) is a major source of uncertainty in climate projections. While vegetation typically reduces stomatal conductance during drought, accurately modeling how and to what degree stomata respond to changes in moisture availability at global scales is particularly challenging, because no global scale g c observations exist. Here, we leverage a collection of satellite, reanalysis and station-based near-surface air and surface temperature estimates, which are physically and statistically linked to $γ^{m}_{gc}$ due to the local cooling effect of g c through transpiration, to develop a novel emergent constraint of $γ^{m}_{gc}$ in an ensemble of Earth System Models (ESMs). We find that ESMs systematically underestimate $γ^{m}_{gc}$ by ~33%, particularly in grasslands, croplands, and savannas in semi-arid and bordering regions of the Central United States, Central Europe, Southeastern South America, Southern Africa, Eastern Australia, and parts of East Asia. We show that this underestimation occurs because ESMs inadequately reduce g c when soil moisture decreases. As g c controls carbon, water and energy fluxes, the misrepresentation of modeled $γ^{m}_{gc}$ contributes to biases in ESM projections of gross primary production, transpiration, and temperature during droughts. Our results suggest that the severity and duration of droughts may be misrepresented in ESMs due to the impact of sustained g c on both soil moisture dynamics and the biosphere-atmosphere feedbacks that affect local temperatures and regional weather patterns.

54 ENVIRONMENTAL SCIENCES↗