Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Representatives”

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 19 records

Mesoscale Convective Systems Represented in High Resolution E3SMv2 and Impact of New Cloud and Convection Parameterizations

Mesoscale convective systems (MCSs) play an important role in modulating the global hydrological cycle, general circulation, and radiative energy budget. In this study, we evaluate MCS simulations in the second version of U.S. Department of Energy (DOE) Energy Exascale Earth System Model (E3SMv2). E3SMv2 atmosphere model (EAMv2) is run at the uniform 0.25? horizontal resolution. We track MCSs consistently in the model and observations using the PyFLEXTRKR algorithm, which defines MCS based on both cloud-top brightness temperature (Tb) and surface precipitation. Results from using Tb only to define MCS, commonly used in previous studies, are also discussed. Furthermore, sensitivity experiments are performed to examine the impact of new cloud and convection parameterizations developed for EAMv3 on simulated MCSs. Our results show that EAMv2 simulated MCS precipitation is largely underestimated in the tropics and contiguous United States. This is mainly attributed to the underestimated precipitation intensity in EAMv2. In contrast, the simulated MCS frequency becomes more comparable to observations if MCSs are defined only based on cloud-top Tb. The Tb-based MCS tracking method, however, includes many cloud systems with very weak precipitation which conflicts with the MCS definition. This result illustrates the importance of accounting for precipitation in evaluating simulated MCSs. We also find that the new physics parameterizations help increase the relative contribution of convective precipitation to total precipitation in the tropics, but the simulated MCS properties are generally not improved. This suggests that simulating MCSs will remain a challenge for the next version of E3SM.

Zhang, Meng

On the Representativity of Electrode Microstructure Parameters and Their Electrochemical Response for Lithium Ion Batteries

Lithium-ion battery electrochemical models require an accurate description of the electrodes microstructures to be predictive that can be achieved through nanoscale imaging. Such observations are however limited by their field of view (FOV), as they provide only a subset of the whole electrode volume that does not necessarily represent the whole electrode microstructure heterogeneity, and therefore can bias the microstructure analysis. A microstructure scale electrochemical model was used to investigate lithium plating onset, material non-uniform utilization, and in-plane heterogeneities for an NMC-graphite full cell. To evaluate the representativeness, and thus relevance, of these model predictions, a coupled representativity analysis has been performed on the microstructure parameters and, in a novel way, on the full cell electrochemical response. Electrode microstructure parameters representativeness has been first quantified using the representative volume element (RVE) methodology. The RVE major flaw is that ultimately it can only conclude if a FOV contains representative subvolumes of the FOV, but not if the FOV itself is representative of the electrode volume. Analysis can conclude negatively ('FOV is not representative'), but not positively ('FOV is representative'). One major contribution of this work was to quantify the convergence of the RVE size with the FOV, to actually investigate the FOV representativeness and thus partly remedy this intrinsic limitation. The analysis determined that performing a standard RVE calculation, without exploring its FOV convergence, is likely to strongly underestimate the actual RVE size. The new RVE methodology has been automated in the NREL open-source Microstructure Analysis Toolbox (MATBOX) and is available to the battery community. Representativeness of microstructure parameters is however only an intermediate step, as the end-results of an electrochemical model are performances predictions. Indeed, what is the practical consequence of a given deviation for a microstructure parameter? The microstructure parameter deviation propagations to the 3D microstructure scale electrochemical response have been then quantified for different charge rates. This defines a threshold for the microstructure parameters FOV for a desired maximum deviation of the electrochemical response. Such deviation propagation analysis is analogous to error propagation analysis and is necessary to determine the relevance of microstructure scale model predictions for macroscale predictions. Electrochemical model shows cell representative section areas are increasing with C-rate, due to higher in-plane heterogeneities, indicating larger FOVs are required specifically for fast charge modeling. Therefore, we introduced the novel concept of electrochemical RVE (eRVE) that is a function of the operating conditions (thus defined as a dynamic RVE), with an increasing dependence with the C-rate. Representativity analysis of the investigated cell determined a FOV of 144.4 x 54.4 m2 is large enough to establish a convergence on the representative section areas for low to intermediate C-rate (=2.5C), but not large enough to conclude for higher rates. This work aims to emphasize the importance of representativity analysis for LIB electrode microstructures, as it is required to estimate the error, and thus the relevance, of microstructure parameters intended to be used in macroscale models. The methodology and results can help researchers to select the relevant imaging and associated FOV required to provide accurate enough microstructure parameters.

ADVANCED PROPULSION SYSTEMS

A Fast Algorithm for Computing Zigzag Representatives

Zigzag filtrations of simplicial complexes generalize the usual filtrations by allowing simplex deletions in addition to simplex insertions. The barcodes computed from zigzag filtrations encode the evolution of homological features. Although one can locate a particular feature at any index in the filtration using existing algorithms, the resulting representatives may not be compatible with the zigzag: a representative cycle at one index may not map into a representative cycle at its neighbor. For this, one needs to compute compatible representative cycles along each bar in the barcode. It is known that the barcode for a zigzag filtration with m insertions and deletions can be computed $O(m^ω)$ in time, where $ω < 2.373$ is the matrix multiplication exponent. However, it is not known how to compute the compatible representatives so efficiently. For a non-zigzag filtration, the classical matrix-based algorithm provides representatives in $O(m^3)$ time, which can be improved to $O(m^ω)$. However, no known algorithm for zigzag filtrations computes the representatives with the $O(m^3)$ time bound. We present an $O(m^3 n)$ time algorithm for this problem, where $n ≤ m$ is the size of the largest complex in the filtration.

Persistent homology

Microstructure Scale Lithium-Ion Battery Modeling, Part IV: The Representativity of Microstructure Parameters and Electrochemical Response

Lithium-ion battery electrochemical models require an accurate description of the electrodes microstructures to be predictive, that can be achieved through nanoscale imaging. Such observations are however limited by their field of view (FOV), as they provide only a subset of the whole electrode volume that does not necessarily represent the whole electrode microstructure heterogeneity, and therefore can bias the analysis. A representativity analysis has been performed on the microstructure parameters and, in a novel way, on the full cell electrochemical response to evaluate the predictions representativeness, and thus relevance, of a microstructure scale electrochemical model. The microstructure parameter deviation propagations to the electrochemical response have been quantified for different charge rates. This defines a threshold for the microstructure parameters FOV for a desired maximum deviation of the electrochemical response. Electrochemical model shows cell representative section areas are increasing with C-rate, due to higher in-plane heterogeneities, indicating larger FOVs are required specifically for fast charge modeling. Representativity analysis determines a cell FOV of 144.4 × 154.4 μm 2 is large enough to establish a convergence on the representative section areas for low-intermediate C-rate (≤2.5 C), therefore positively concludes on the model representativeness for these rates, but is not large enough to conclude for higher rates.

25 ENERGY STORAGE

Representativity error scaling of models of the high temperature test facility

Error scaling is a critical step toward the validation of modeling capabilities for high-temperature gas-cooled reactors (HTGRs). Extensive effort is being made to bring HTGRs into the validation basis of numerous thermal-hydraulics codes. Here, this paper demonstrates how systems-level codes can be leveraged to perform error scaling analyses between an experimental facility and a plant-to-be facility. Specifically, we focus on two conduction cooldown experiments from the High Temperature Test Facility (HTTF) and the General Atomics 350 MW th modular high-temperature gas-cooled reactor (MHTGR-350). The error scaling methodology employed in this study is representativity, which in the context of this work was used to quantify how well experiments captured the physics of the plant facility by comparing sensitivity vectors between the two facilities. In addition, two different RELAP5–3D models of the HTTF were included in the comparison to determine if modeling methodologies notably impact the representativity results. The key figures of merit are the maximum block temperature and the coolant outlet temperature. The time-dependent maximum block temperature had a low representativity of below 0.1 for both models across both transients. The coolant outlet temperature had a higher representativity of around 0.6 for both models during the pressurized conduction cooldown, but it was below 0.2 for the depressurized conduction cooldown. Overall, the HTTF experiments were not representative of the transients in the MHTGR-350. However, these results can play a significant role in informing future potential experiments for HTGR systems that iterate on what the HTTF accomplished.

22 - GENERAL STUDIES OF NUCLEAR REACTORS

Impact of representative ground motion level on seismic PSA with the boundary between overestimation and underestimation

One commonly used approach in seismic probabilistic safety assessment (PSA) is the discrete method. This method follows the standard PSA framework and can be applied to various models, such as multi-unit models, while reducing computational costs using standard software. However, due to the inability to subdivide intervals infinitely, the discrete method approximates with a finite number of subintervals. In practice, different numbers of subintervals are applied, and the representative ground motion level is selected based on expert judgment. When employing a smaller number of subintervals, it is important to take caution to prevent underestimation. This study analyzes the impact of the representative ground motion level on seismic risk. It confirms that underestimation can occur with a small number of subintervals depending on the representative ground motion level. This study also proposes a method for determining the boundary of underestimation and overestimation. The method is demonstrated through examples, providing a mathematical foundation for selecting appropriate representative ground motion levels. By avoiding underestimation, this research helps prevent the oversight of significant risk contributors and enhances the understanding of seismic risk.

99 - GENERAL AND MISCELLANEOUS

Identification of Climatological Representative Days in the Mid-Atlantic for High-Fidelity Offshore Wind Energy Modeling

The goal of reaching 30 GW of offshore wind energy by 2030 becomes more realistic with the continued approval of offshore wind energy areas by the Biden Administration. In the Mid-Atlantic, where wind energy projects are in the most advanced stages of development, there is increased research focus on the eventual interaction of these wind farms. These interactions, in the form of wakes and cluster wakes, or wakes from multiple wind farms, could have detrimental effects on power production and forecastability for downwind wind farms (Pryor et al. 2022, Golbazi et al. 2022, Rosencrans et al. 2023). To help alleviate these issues, numerical simulations in the form of numerical weather prediction (NWP) and large eddy simulations (LES) can provide insight into when cluster wake situations may occur, but running such simulations can be expensive and difficult to run for multiple years. In this study, we leverage and build upon existing techniques in the literature (Fischereit et al. 2022) to identify climatologically representative days for wind energy areas in the Mid-Atlantic where conditions would promote cluster wake situations. We select meteorological variables (wind speed, wind direction, atmospheric stability, boundary-layer height, TKE) critical to understanding wind energy production and wake propagation. We then consider two different NWP datasets of varying spatial and temporal resolution: ERA5 provides data at hourly intervals from 1940 to present at 0.25 deg (31 km) spatial resolution (Hersbach et al. 2020), and the NOW-23 dataset provides data at 5-minute resolution for 21 years at 2-km spatial resolution (Bodini et al. 2020). Our first step is to compare these two datasets for an overlapping 21-year time period. Initial results show that the required number of days to represent the long-term climate increases with each additional variable considered. In their study of the German Bight, Fischereit et al. (2022) found that they could represent the long-term wind and wave climate in a "near-perfect" way with -180 days, by reaching a Perkins Skill Score (PSS) of 0.9; our investigation of the mid-Atlantic wind resource region with ERA5 and NOW-23 data suggests that we will need -100 days to reach a PSS of 0.9. As we expand our parameter space to include multiple variables, the number of required days will likely grow. These results will ultimately be used to select case studies to best represent cluster wake conditions that apply to this region for the lifetime of likely wind farms in this mid-Atlantic region.

clusterwakes

The Freundlich isotherm equation best represents phosphate sorption across soil orders and land use types in tropical soils of Puerto Rico

Biomass production in the lowland wet tropical forest is greater than in any other biome, and it is typically limited by soil phosphorus (P) availability. However, the mechanisms involved in the P cycle remain poorly represented in Earth System Models (ESMs). Soil P sorption processes are key in the P cycle and for understanding the extent of P limitation for plant productivity. Currently, a few ESMs include isotherm equations to model these processes. Although the Langmuir equation is widely cited, other isotherm equations may better describe sorption in tropical soils. Here, we use a diverse range of soil samples from Puerto Rico to test the validity of the Langmuir, Freundlich, and Temkin equation. We found that across four soil orders (Inceptisols, Mollisols, Oxisols, Ultisols), and forested and cultivated land use types, the Freundlich equation best represented soil P sorption. Furthermore, the Langmuir and the Temkin equations poorly represent soil P adsorption, especially at low P concentrations. Specifically, the Langmuir equation underestimated soil P adsorption by 40% and the Temkin equation overestimated adsorption by 76%. We also found, as expected, that soil clay content and pH were the most important parameters explaining the variability of the Freundlich (K f ) constant. Greater clay content and lower pH, common in highly weathered Ultisols and Oxisols which are abundant in the tropics, led to greater K f values. Overall, our results suggest that a diversity of soils can prompt underestimation of P sorption when using the Langmuir isotherm, which leads to an overestimation of available P that can have repercussions on ESM predictions of the P cycle and tropical forest productivity.

Earth system models

A New Shutdown Dose Rate Benchmark Problem for Representative Fusion Applications

Here, this work introduces a new benchmark problem for calculating shutdown dose rates (SDDRs) aimed at fusion reactor applications. The model is designed to represent a simplified version of a typical ITER port plug. The responses of interest include neutron flux, gamma flux, and gamma SDDR at 12 different locations scattered throughout the port. This article outlines the geometry specifications of the problem, provides material definitions for the components, specifies the required responses to be calculated, and presents the source definition information. The need for this benchmark arises from the limited availability of publicly accessible references, with only one benchmark representing the typical dimensions and materials found in fusion systems. This existing benchmark has been cited extensively, reflecting the demand within the scientific community to test both established and novel workflows for SDDR calculations. However, since its presentation at a conference in 2011, the results have become increasingly well known. Moreover, the absence of formal publication and peer review has led to the details of this benchmark being extracted from secondary sources, such as subsequent studies that reference it. As a result, analysts are left with significant flexibility in interpreting the key parameters, which can be adjusted to account for unknown systematic errors, ultimately reproducing the already well-known responses. This new benchmark serves as an updated version of that earlier work, with the aim of providing a more reliable description of the materials and their impurities, which is crucial for assessing activation and subsequent gamma emission. Additionally, it seeks to provide a geometry that more closely represents an ITER port plug. The improvements in the problem definition will lead to a more reproducible benchmark problem, while also presenting the radiation transport community with a completely new challenge. The results will be published in a future article to allow analysts adequate time to analyze this problem independently.

Benchmark

On the bulk compaction of brittle granular materials, Part II: Representing interparticle comminution by K COMP

An applied and empirical compaction fracture resistance parameter, K COMP , was developed to represent the compression-induced interparticle-comminution response of granular brittle materials (e.g., ceramic or glass particles). The development of K COMP is an outcome from Part I of this three-paper series. The K COMP represents a continuum response of compaction fracture resistance and was adapted from established Griffith linear elastic fracture mechanics theory. From that, the K COMP relates macroscopically applied compressive stress to a corresponding inverse-square-root representative particle size, where the latter is estimated from an entire particle size distribution or specific surface area measurement. In conclusion, the K COMP value was not constant over a wide range of stresses for all examined brittle granular materials, and this is indicative of a change in the dominant mode of permanent deformation at higher compaction stresses.

brittle

Do we have globally representative data to understand soil processes?

Understanding and modeling soils and soil organic matter (SOM) are central to a variety of human needs, from food production to ecosystem management. Soil data have been collected for over a century, but the global spatial and process representativeness of soil data remains unclear. We assessed the representativeness of currently available soil data that could be used to understand a variety of SOM processes. We used 16 open-source soil databases and data from over 281,000 unique locations globally, categorizing the databases into three main data types necessary to understand SOM processes: soil carbon stocks and fluxes, mechanistic drivers of these stocks and fluxes, and soil carbon gain or loss potential. We found that stock and driver data have extensive global coverage. However, data on soil carbon gain or loss potential, particularly data describing change in soils over time such as time series data, are severely limited in their global coverage. We conclude that while significant strides have been made in measuring soil carbon stocks and fluxes, and their drivers, we are limited in global data related to changes in soils over time. Our recommendations for soil data generators are to ensure precise metadata reporting and prioritizing sampling in underrepresented areas like tropical, arctic, mountainous, wetland and arid regions. We also encourage designing revisit schemes that explicitly support change detection and reporting multi-modal datasets that can aid in model development. Targeted measurement of low coverage soil data types and regions is necessary for a range of applications including current and future biogeochemical predictions, and their management and policy implications.

carbon fluxes

RAVEN Template for Dynamic Representativity Analysis of the High Temperature Test Facility

These slides present a walkthrough of the template that has been developed for using RAVEN to perform representativity analysis using models of the High Temperature Test Facility and the General Atomics Modular High Temperature Gas-cooled Reactor. The presentation provides participants in the HTTF benchmark with a walkthrough on how to use RAVEN for their sensitivity analysis and how to read results from the MHTGR-350 to perform representativity

22 - GENERAL STUDIES OF NUCLEAR REACTORS

Representing cropping systems with the MEMS 2 ecosystem model

Abstract Croplands have been the focus of substantial investigation due to their considerable potential for sequestering carbon. Understanding the potential for soil organic carbon (SOC) sequestration and necessary management strategies will be enabled with accurate process‐based models. Accurately representing crop growth and agricultural practices will be critical for realistic SOC modeling. The MEMS 2 model incorporates a current understanding of SOC formation and stabilization, measurable SOC pools, and deep SOC dynamics and is seen as a highly promising tool to inform management intervention for SOC sequestration. Thus far, MEMS 2 has been developed to represent grasslands. In this study, we further developed MEMS 2 to model annual grain crops and common agricultural practices, such as irrigation, fertilization, harvesting, and tillage. Using four Ameriflux sites, we demonstrated an accurate simulation of crop growth and development. Model performance was strong for simulating aboveground biomass (index of agreement [d] range of 0.89–0.98) and green leaf area index (dfrom 0.90 to 0.96) across corn, soybean, and winter wheat. Good agreement with observations was also achieved for net ecosystem CO 2 exchange (dfrom 0.90 to 0.96), evapotranspiration (dfrom 0.91 to 0.94), and soil temperature (dof 0.96), while discrepancy with the available soil water content data remain (dfrom 0.14 to 0.81 at four depths to 100 cm). While we will continue model testing and improvement, MEMS 2 (version 2.14) has now demonstrated its ability to effectively simulate the growth of common grain crops and practices.

Agriculture

Fully implicit crystal plasticity models representing orientations with modified Rodrigues parameters

Here, this work describes a crystal plasticity formulation combining several mathematical, numerical, and implementation choices to produce a highly efficient model. Specifically, the key choices in the implementation are (1) representing orientations with modified Rodrigues parameters, (2) implementing a fully coupled implicit time integration for the elastic stretch, the crystal orientations, and the model internal variables, (3) implementing the model in the NEML2 constitutive modeling framework, based on PyTorch, to vectorize the calculations and port the computation to GPUs and other hardware accelerators, and (4) an exact implementation of the consistent tangent matrix, even for arbitrary coupling to other field variables beyond the displacements, like temperature, neutron fluence, etc. The first two features of the model are, to our knowledge, novel. The paper considers each of these choices individually as well as the final model as a whole. This includes a full description of modified Rodrigues parameters, their advantages over other representations of orientations, the mathematical formulae and tools required to implement a model with modified Rodrigues parameters, and a detailed description of the geometry of the space of modified Rodrigues parameters (in an appendix). It also includes a description of a fully implicit time integration scheme for the orientations and the advantages in representing orientations with modified Rodrigues parameters in implementing such a model. The work then assess, via numerical examples, the advantages of fully coupled implicit time integration versus more common decoupled and explicit time integration schemes. These studies demonstrate the computational advantages of fully coupled integration versus other time integration algorithms, though the performance of the competing models depends on the complexity of the underlying single crystal model. The study concludes by demonstrating that the choice of time integration method affects the sharpness of the predicted texture, with explicit methods for integrating the orientations overestimating texture sharpness and implicit methods underestimating texture sharpness.

Crystal plasticity

Representing Anthropogenic Dust From Agricultural Sources in E3SMv1: Implementation, Evaluation, and Assessment of the Radiative Forcing

Dust emissions related to anthropogenic activities (i.e., anthropogenic dust) is not represented in most global climate models and its radiative impact remains unassessed. In this study, we develop a new and physically based method to parameterize anthropogenic dust emission from agricultural sources (i.e., agricultural dust, or AD) based on the DOE's Energy Exascale Earth System Model version 1 (E3SMv1). This method relates AD emission to the crop land use fraction in the E3SMv1 land component. Major AD sources simulated by our parameterization include those over Central America, the Sahel, North India, and North China. The annual averaged AD emission with diameter less than 10 μm is 567 Tg yr −1 in present-day (year 2000), which contributes to 13.3% of total dust emission. Model evaluation against satellite and ground-based observations shows that the new parameterization can represent AD emissions and global dust cycle reasonably well. We find that the total dust emission increases by 13.1% (495 Tg yr −1 ) from 1850 to 2000 due to the cropland land use fraction changes without considering the impact from climate changes. This induces a net dust direct effective radiative forcing of −0.041 W m −2 at top of the atmosphere. This dust-induced cooling exceeds 10% of the total anthropogenic aerosol direct effective radiative forcing from 1750 to 2014 estimated by the Intergovernmental Panel on Climate Change Sixth Assessment Report. Our findings indicate an important role of anthropogenic dust in global climate change, which should be included in future climate change assessments.

agricultural dust

An ontology-based knowledge graph for representing interactions involving RNA molecules

The "RNA world" represents a novel frontier for the study of fundamental biological processes and human diseases and is paving the way for the development of new drugs tailored to each patient's biomolecular characteristics. Although scientific data about coding and non-coding RNA molecules are constantly produced and available from public repositories, they are scattered across different databases and a centralized, uniform, and semantically consistent representation of the "RNA world" is still lacking. We propose RNA-KG, a knowledge graph (KG) encompassing biological knowledge about RNAs gathered from more than 60 public databases, integrating functional relationships with genes, proteins, and chemicals and ontologically grounded biomedical concepts. To develop RNA-KG, we first identified, pre-processed, and characterized each data source; next, we built a meta-graph that provides an ontological description of the KG by representing all the bio-molecular entities and medical concepts of interest in this domain, as well as the types of interactions connecting them. Finally, we leveraged an instance-based semantically abstracted knowledge model to specify the ontological alignment according to which RNA-KG was generated. RNA-KG can be downloaded in different formats and also queried by a SPARQL endpoint. A thorough topological analysis of the resulting heterogeneous graph provides further insights into the characteristics of the "RNA world". RNA-KG can be both directly explored and visualized, and/or analyzed by applying computational methods to infer bio-medical knowledge from its heterogeneous nodes and edges. The resource can be easily updated with new experimental data, and specific views of the overall KG can be extracted according to the bio-medical problem to be studied.

59 BASIC BIOLOGICAL SCIENCES

Generalized representative structures for atomistic systems

A new method is presented to generate atomic structures that reproduce the essential characteristics of arbitrary material systems, phases, or ensembles. Previous methods allow one to reproduce the essential characteristics (e.g. the chemical disorder) of a large random alloy within a small crystal structure. The ability to generate small representations of random alloys, along with the restriction to crystal systems, results from using the fixed-lattice cluster correlations to describe structural characteristics. A more general description of the structural characteristics of atomic systems is obtained using complete sets of atomic environment descriptors. These are used within for generating representative atomic structures without restriction to fixed lattices. A general data-driven approach is provided here utilizing the atomic cluster expansion (ACE) basis. The N-body ACE descriptors are a complete set of atomic environment descriptors that span both chemical and spatial degrees of freedom and are used within for describing atomic structures. The generalized representative structure (GRS) method presented within generates small atomic structures that reproduce ACE descriptor distributions corresponding to arbitrary structural and chemical complexity. It is shown that systematically improvable representations of crystalline systems on fixed parent lattices, amorphous materials, liquids, and ensembles of atomic structures may be produced efficiently through optimization algorithms. With the GRS method, we highlight reduced representations of atomistic machine-learning training datasets that contain similar amounts of information and small 40–72 atom representations of liquid phases. The ability to use GRS methodology as a driver for informed novel structure generation is also demonstrated. The advantages over other data-driven methods and state-of-the-art methods restricted to high-symmetry systems are highlighted.

atomic cluster expansion