Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model types”

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 163 records · Page 9

Analysis of Element Size and Element Type on SS304L Material Model Performance

Finite element solutions are only valid if the information inputted into the models is a valid assumption. This is why it is important to determine if the model behaves as expected especially in the material response. Material models that consider deformation in the elastic and plastic regions are desirable. As a material deforms plastically, hardening will occur until the ultimate strength is reached due to the reduction of the cross sectional area in a uniaxial tensile test environment. A phenomenological approximation of this behavior which is known as power law hardening (plh) can be expressed as σ = σ y + K $ϵ^{n}_{p}$ where K and n are empirically determined constants, and σ y is the proportional limit of the material. The first objective of this study is to verify the behavior for a calibrated plh material model for selective deviatoric hexahedral elements when applied to a different material stock. The Sierra SM input for this calibrated model can be found in the appendix. The second objective of this study is to determine the differences in failure characteristics for considerably large shell elements versus that of fine hexahedral elements.

36 MATERIALS SCIENCE↗

P-model v1.0: an optimality-based light use efficiency model for simulating ecosystem gross primary production

Terrestrial photosynthesis is the basis for vegetation growth and drives the land carbon cycle. Accurately simulating gross primary production (GPP, ecosystem-level apparent photosynthesis) is key for satellite monitoring and Earth system model predictions under climate change. While robust models exist for describing leaf-level photosynthesis, predictions diverge due to uncertain photosynthetic traits and parameters which vary on multiple spatial and temporal scales. Here, we describe and evaluate a GPP (photosynthesis per unit ground area) model, the P-model, that combines the Farquhar–von Caemmerer–Berry model for C 3 photosynthesis with an optimality principle for the carbon assimilation–transpiration trade-off, and predicts a multi-day average light use efficiency (LUE) for any climate and C 3 vegetation type. The model builds on the theory developed in Prentice et al. (2014) and Wang et al. (2017a) and is extended to include low temperature effects on the intrinsic quantum yield and an empirical soil moisture stress factor. The model is forced with site-level data of the fraction of absorbed photosynthetically active radiation (fAPAR) and meteorological data and is evaluated against GPP estimates from a globally distributed network of ecosystem flux measurements. Although the P-model requires relatively few inputs, the R 2 for predicted versus observed GPP based on the full model setup is 0.75 (8 d mean, 126 sites) – similar to comparable satellite-data-driven GPP models but without predefined vegetation-type-specific parameters. The R 2 is reduced to 0.70 when not accounting for the reduction in quantum yield at low temperatures and effects of low soil moisture on LUE. The R 2 for the P-model-predicted LUE is 0.32 (means by site) and 0.48 (means by vegetation type). Applying this model for global-scale simulations yields a total global GPP of 106–122 Pg C yr –1 (mean of 2001–2011), depending on the fAPAR forcing data. The P-model provides a simple but powerful method for predicting – rather than prescribing – light use efficiency and simulating terrestrial photosynthesis across a wide range of conditions. The model is available as an R package (rpmodel).

54 ENVIRONMENTAL SCIENCES↗

Integrating New Arctic Plant Functional Types in a Land Surface Model Using Above- and Belowground Field Observations: Modeling Archive

This Modeling Archive is in support of an NGEE Arctic publication "Integrating New Arctic Plant Functional Types in a Land Surface Model Using Above- and Belowground Field Observations". We used field observations of biomass and leaf traits (Salmon et al., 2019b) across a gradient of plant communities on the Seward Peninsula in northwest Alaska (Breen et al., 2020; Iversen et al., 2017) to replace the original ELM configuration with nine, arctic-specific PFTs. Original arctic PFTs in the model were (1) broadleaf deciduous boreal shrub and (2) C3 arctic grass. New PFTs were (1) lichen, (2) bryophyte, (3) evergreen dwarf shrub, (4) deciduous dwarf shrub, (5) deciduous low shrub, (6) deciduous low to tall shrub, (7) alder shrub, (8) forb, and (9) graminoid. New PFTs were parameterized and evaluated using site-level trait and biomass measurements (Salmon et al., 2019b), root C:N ratios and rooting depth patterns (Salmon et al., 2019c), and root biomass (Salmon et al., 2020). This archive contains forcing data and model output from ELM simulations conducted at the Kougarok Hillslope site (Kougarok Road Mile Marker 64). Simulations were conducted for the Kougarok Hillslope site using meteorological driving data from the Scenarios Network for Alaska and Arctic Planning (SNAP) downscaled climate projection dataset using NCAR-CCSM forcing (Bieniek et al., 2020 dataset; Walsh et al., 2018). The archive includes three versions of model simulations representing a progression from the original model configuration to an Arctic-specific configuration (see Output section below). Simulation 1 used grid cell data from a global E3SM configuration, including the fractional area of E3SM PFTs assigned to the grid cell containing the Kougarok Hillslope site in global simulations. This simulation did not distinguish between different plant communities on the Kougarok Hillslope, but instead used a single point simulation to represent the entire area. Simulation 2 used default E3SM PFT definitions combined with adjusted depth to bedrock for plant communities with shallow rocky layers to represent the role of abiotic soil factors in driving site differences (based on depth measurements of Iversen et al., 2019a), and adjusted the relative areas of the broadleaf deciduous boreal shrub and C3 arctic grass PFTs to reflect the observed spatial coverage of shrub and non-shrub PFTs across Kougarok Hillslope plant communities. Simulation 3 used the new arctic PFT definitions and parameterizations based on the vegetation types present at the site, along with the adjusted soil depths used in Simulation 2. Included are *.pdf, *.nc (NetCDF), *tar.gz, *.py, and *.bash files. We would like to thank Mary's Igloo Native Corporation for allowing us to perform this research on their land and giving us the opportunity to share this research with their community. Thanks to Shawn Serbin and Alistair Rogers for helpful discussions. Holly Vander Stel, Joanne Childs, and Stan Wullschleger helped with collection of field data. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

Search for a light pseudoscalar Higgs boson in the boosted $\mu\mu\tau\tau$ final state in proton-proton collisions at $\sqrt{s}=$ 13 TeV

A search for a light pseudoscalar Higgs boson (a) decaying from the 125 GeV (or a heavier) scalar Higgs boson (H) is performed using the 2016 LHC proton-proton collision data at $ \sqrt{s} $ = 13 TeV, corresponding to an integrated luminosity of 35.9 fb$^{−1}$, collected by the CMS experiment. The analysis considers gluon fusion and vector boson fusion production of the H, followed by the decay H → aa → μμττ, and considers pseudoscalar masses in the range 3.6 < m$_{a}$< 21 GeV. Because of the large mass difference between the H and the a bosons and the small masses of the a boson decay products, both the μμ and the ττ pairs have high Lorentz boost and are collimated. The ττ reconstruction efficiency is increased by modifying the standard technique for hadronic τ lepton decay reconstruction to account for a nearby muon. No significant signal is observed. Model-independent limits are set at 95% confidence level, as a function of m$_{a}$, on the branching fraction (ℬ) for H → aa → μμττ, down to 1.5 (2.0) × 10$^{−4}$ for m$_{H}$ = 125 (300) GeV. Model-dependent limits on ℬ(H → aa) are set within the context of two Higgs doublets plus singlet models, with the most stringent results obtained for Type-III models. These results extend current LHC searches for heavier a bosons that decay to resolved lepton pairs and provide the first such bounds for an H boson with a mass above 125 GeV.[graphic not available: see fulltext]

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

SALT2 versus SALT3: updated model surfaces and their impacts on type Ia supernova cosmology

ABSTRACT For the past decade, SALT2 has been the most common model used to fit Type Ia supernova (SN Ia) light curves for dark energy analyses. Recently, the SALT3 model was released, which upgraded a number of model features but has not yet been used for measurements of dark energy. Here, we evaluate the impact of switching from SALT2 to SALT3 for a SN cosmology analysis. We train SALT2 and SALT3 on an identical training sample of 1083 well-calibrated Type Ia supernovae, ensuring that any differences found come from the underlying model framework. We publicly release the results of this training (the SALT ‘surfaces’). We then run a cosmology analysis on the public Dark Energy Survey 3-Yr Supernova data sample (DES-SN3YR), and on realistic simulations of those data. We provide the first estimate of the SN + CMB systematic uncertainty arising from the choice of SALT model framework (i.e. SALT2 versus SALT3), Δw = + 0.001 ± 0.005 – a negligible effect at the current level of dark energy analyses. We also find that the updated surfaces are less sensitive to photometric calibration uncertainties than previous SALT2 surfaces, with the average spectral energy density dispersion reduced by a factor of two over optical wavelengths. This offers an opportunity to reduce the contribution of calibration errors to SN cosmology uncertainty budgets.

79 ASTRONOMY AND ASTROPHYSICS↗

Changes in water vapor adsorption and water film thickness in clayey materials as a function of relative humidity

Abstract Soil water, with adsorbed water being an important component, plays an important role in fluid flow and chemical movement in the unsaturated zone. The adsorbed water consists of water film and interlayer water, and its content is related to the relative humidity (RH). This study theoretically and experimentally focuses on the variation in adsorbed water content in clayey materials as the RH changes. Based on a slit‐pore model, three types of water (water film, interlayer water, and capillary water) could be present in the water vapor adsorption process. The variation of water film is obtained from the analysis of interfacial forces, and then the total water content can be quantitatively subdivided into these three types of water for different RH values. The slit‐pore adsorption model was used with experimental measurements (water vapor adsorption and N 2 physisorption) to quantify the amounts and interrelationships of these three different types of water to six clayey materials, namely three clay minerals (kaolinite, montmorillonite, and illite–smectite mixed layer) and three clay‐rich sediments from the Jianghan Plain in China. The main results are that (a) the volumes of water vapor adsorption are much greater than those of N 2 adsorption except for kaolinite; (b) interlayer water largely dominates the growth of total water content in montmorillonite and results in a concentration of 0.090 ml g −1 at 95% RH; and (c) the minimum thickness of water film is calculated to be 0.23 nm, and the maximum value is one‐third of the pore width by considering the interfacial forces.

54 ENVIRONMENTAL SCIENCES↗

Revisiting diffusion-based moisture functions: why do they fail?

Microbial activities in soil are controlled by the balance between physical and biological mechanisms that are currently implied by soil moisture-respiration relationships, but not explicitly represented. Diffusion based moisture functions focus on physical controls, while substrate kinetics based models primarily address biological drivers. How to harmonize these two types of models for robust predictions of soil respiration and other processes has not been addressed. Here we derived a generalizable analytical solution of diffusion-limited nonlinear microbial uptake kinetics, in which microbial metabolic rates were governed by both physical and biological drivers via substrate diffusion and microbial affinity. We proposed a soil trait concept to integrate both physical and biological perspectives into one unified function, and tested it using both laboratory and field based observations from a variety of soils. This approach provides a tractable and scalable approach that incorporates soil heterogeneity and functional complexity in macroscopic model descriptions.

59 BASIC BIOLOGICAL SCIENCES↗

sf-24-103

pyzag is a library for efficiently training generic models defined with a recursive nonlinear function. We developed it primarily to provide an efficient means to train NEML2 models it can be used to train other, more general types of models.

Messner, Mark [Argonne National Laboratory (ANL), ↗

Pretest Modeling A Spent Nuclear Fuel Seismic Shake Test

The U.S. Department of Energy Spent Fuel and Waste Science and Technology (SFWST) program is planning to conduct a series of full-scale shake table tests to simulate hypothetical earthquake conditions and record the response of surrogate spent nuclear fuel (SNF) assemblies in a dry canister storage system mockup. The shake table motions will represent a range of hypothetical earthquake conditions at hypothetical locations in the continental U.S. to generally define the range of mechanical loads that SNF can be expected to experience during extended dry storage periods. This paper describes the pretest predictions made with LS-DYNA models of the mockup storage systems. The test will use two dry storage system configurations, a mockup vertical concrete cask (VCC) and a mockup horizontal storage module (HSM). The test will use a production-quality canister and basket. Within the canister will be four instrumented fuel assemblies with fuel rods containing surrogate mass and 28 instrumented dummy assemblies that are intended to match the mass and outer dimension of a fuel assembly. The finite element models include models of the VCC and HSM on the shake table to calculate the system level dynamic responses and separate single fuel assembly models to calculate stress and strain in fuel assembly components. Both types of models include nonlinear behavior like rod-to-rod contact and the ability for VCC’s to rock and slide. This paper presents the expected response of the VCC, HSM, and fuel assemblies to the shake table testing that is planned to start in April of 2024. The earthquake conditions represent seismic hazards in the 2,000-to-20,000-year return period range. The test data is expected to confirm the expectation that fuel rod cladding will remain intact, fuel assembly structural components like guide tubes will remain intact, and no significant VCC sliding or tipping will occur in the range of conditions to be tested with the shake table.

Klymyshyn, Nicholas A.↗

Optimization of Canister Loading Patterns in Dual Purpose Canisters for Criticality Suppression

This report documents work performed in support of the US Department of Energy Office of Nuclear Energy (NE) Spent Fuel and Waste Disposition, Spent Fuel and Waste Science and Technology, under work breakdown structure element 1.08.01.03.05, “Direct Disposal of Dual Purpose Canisters.” In particular, this report fulfills milestone M3SF-22PN010305094, “Application of AI/ML techniques to optimize DPC loading,” within work package SF-22PN01030509, “Direct Disposal of Dual Purpose Canisters - PNNL.” This report continues the process of examining the potential for using loading optimization as a disposal criticality suppression technique. This latest update to the report added: 1. A validation of the artificial neural network (ANN) reactivity prediction tool against as-loaded dual-purpose canisters from the UNF-ST&DARDS database 2. Refinement of the algorithm to incorporate an additional DPC design to improve the predictions 3. Investigation of a storage and transportation focused loading optimization algorithm. A key result from this year’s work is that the ANN performance was significantly improved by including the additional model for the NUHOMS canisters and it is apparent that fuel type specific modeling considerations should be incorporated into the ANN in future work. The GRASP-enabled adaptive multi-objective memetic algorithm with partial clustering (GAMMA-PC) algorithm was evaluated as a tool for storage and transportation oriented optimization. The GAMMA-PC optimization routine considers decay heat at loading, the time between loading and when the canister is eligible for transportation and minimizes the number of casks loaded. Future work will develop a dose minimization optimization routine.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

An Electromagnetic Transient and Three-Phase Phasor Co-Simulation Platform for Studying Distribution System Transients with High Pentration of DERs

Inverter vendors usually have electro-magnetic transient (EMT) models of distributed energy resources (DERs) for validation purposes and these types of models are usually black-boxed considering intellectual property (IP) protection. Limited by small simulation time step for the entire power system network, studying the system-level responses of DERs in EMT simulators is more time-consuming as the network size increases. Meanwhile, vendors usually do not have inverter models that can be used in three-phase phasor-domain simulators. This paper proposes a platform that can dynamically co-simulate the EMT inverter models and three-phase phasor-domain distribution networks considering high penetration of inverter-based DERs. This platform is very useful for evaluating the impacts of vendor-provided black-box DER model on distribution system dynamics. The co-simulation platform is validated against the full EMT model using a simple two-inverter test system and tested using IEEE 34-node test feeder.

Liu, Yuan↗

Probing lepton number violation and Majorana nature of neutrinos at the LHC

Abstract Observation of lepton number (L) violation by two units at colliders would provide evidence for the Majorana nature of neutrinos. We study signals ofL-violation in the context of two popular models of neutrino masses, the type-II seesaw model and the Zee model, wherein small neutrino masses arise at the tree-level and one-loop level, respectively. We focus onL-violation signals at the LHC arising through the processpp→ℓ ± ℓ ′± + jets within these frameworks. We obtain sensitivity toL-violation in the type-II seesaw model for triplet scalar masses up to 700 GeV and in the Zee model for charged scalar masses up to 4.8 TeV at the high-luminosity LHC with an integrated luminosity of 3 ab −1 .

Physics↗

Using Parameter Sweep in WaterTAP to Analyze New Water Treatment Technologies

We describe a powerful and generalized parameter sweep tool in this report that was originally developed to analyze the performance of existing and novel water treatment models being developed in WaterTAP. Since WaterTAP is built upon IDAES and Pyomo, the parameter sweep tool can be used to systematically explore and debug the behavior of most Pyomo and IDAES numerical models. In order to enable meaningful analyses, the parameter sweep tool has been designed with the following features: 1) Model flexibility: The parameter sweep tool does not enforce any restrictions on the types of models that can be used with it. As long as a Pyomo model can be solved and the parameter is active and mutable, the tool only needs functions that describe how to run the model, the sweep parameters, and the output quantities of interest. 2) Flexible sampling: The parameter sweep tool has inbuilt functions to generate samples from a random distribution or a multidimensional Euclidean space. Furthermore, the users have to ability to supply samples generated from a tool of their choice. 3) Multiple sweep types: A user can choose from one of 3 types of parameter sweeps depending on their needs. 4) Detailed outputs: Outputs generated by the parameter sweep tool can be stored in detailed H5 file or user-friendly CSV files for post processing. 5) Parallel computing: The parameter sweep supports shared and distributed memory parallel computing to enable the use of high performance computers (HPC) for large-scale analyses. 6) Modular: The parameter sweep tool is self-contained and can easily be integrated within an outer-loop analysis or as desired by the user. 7) Ease of use: The tool is well documented and a simple sweep can be easily executed by following the online documentation in a few lines of code. We demonstrate the use of the parameter sweep tool on a simple water treatment system from the WaterTAP repository and show its parallel scaling performance on an Apple laptop and NREL's Eagle HPC. The parameter sweep tool is actively being used with models currently being developed within WaterTAP and we expect its use to grow beyond it to other IDAES and Pyomo models.

97 MATHEMATICS AND COMPUTING↗

Upscaling Wetland Methane Emissions From the FLUXNET–CH4 Eddy Covariance Network (UpCH4 v1.0): Model Development, Network Assessment, and Budget Comparison

Wetlands are responsible for 20%–31% of global methane (CH 4 ) emissions and account for a large source of uncertainty in the global CH 4 budget. Data-driven upscaling of CH 4 fluxes from eddy covariance measurements can provide new and independent bottom-up estimates of wetland CH 4 emissions. Here, we develop a six-predictor random forest upscaling model (UpCH4), trained on 119 site-years of eddy covariance CH 4 flux data from 43 freshwater wetland sites in the FLUXNET-CH4 Community Product. Network patterns in site-level annual means and mean seasonal cycles of CH 4 fluxes were reproduced accurately in tundra, boreal, and temperate regions (Nash-Sutcliffe Efficiency ~0.52–0.63 and 0.53). UpCH4 estimated annual global wetland CH 4 emissions of 146 ± 43 TgCH 4 y –1 for 2001–2018 which agrees closely with current bottom-up land surface models (102–181 TgCH 4 y –1 ) and overlaps with top-down atmospheric inversion models (155–200 TgCH 4 y –1 ). However, UpCH4 diverged from both types of models in the spatial pattern and seasonal dynamics of tropical wetland emissions. We conclude that upscaling of eddy covariance CH 4 fluxes has the potential to produce realistic extra-tropical wetland CH 4 emissions estimates which will improve with more flux data. To reduce uncertainty in upscaled estimates, researchers could prioritize new wetland flux sites along humid-to-arid tropical climate gradients, from major rainforest basins (Congo, Amazon, and SE Asia), into monsoon (Bangladesh and India) and savannah regions (African Sahel) and be paired with improved knowledge of wetland extent seasonal dynamics in these regions.

54 ENVIRONMENTAL SCIENCES↗

Tools for analysis of optimization models

The software is code for analyzing, debugging, and solving optimization models. The code implements several algorithms and provides convenient APIs to apply these algorithms to optimization models. The algorithms are primarily based on graph theory. They compute well-known partitions of graphs, and use these partitions to provide the user information about their optimization model, including diagnosing certain types of structural modeling errors. Some algorithms may be used as a subroutine to solve an optimization problem, and may call an optimization solver to facilitate this. Well-known linear algebra routines are called to provide further diagnostics. The software contains no data and no models other than toy models used for testing.

Parker, Robert↗

PARETO 0.8.0 Release

PARETO 0.8.0 Release. Highlights: Model Updates - Applied unified sets for pipeline and trucking arcs in strategic model - Apply unified sets for pipeline and trucking arcs in operational model - Added new config argument for removal efficiency calculation method - Standardized bidirectional capacity constraint - Added dependencies removed in IDAES 2.1 - Created bounding functions & utilities - Added Hydraulics module to the strategic model - Add additional arc types to strategic model Documentation and Tutorial Updates - Improved PARETO treatment document - Introduced general tutorial and treatment module Jupyter notebooks for Strategic Model - Update docs with correct support email list address - Consolidate and deduplicate Getting Started and resources for developers - Enable Black formatting for Jupyter notebooks - Add Binder configuration files and README Bug Fixes - Fix strategic model documentation typos - Removed duplicated units from output file header UI Updates - Added view for comparing different scenarios - Added functionality for manually overriding PARETO decisions

PARETO,PARETO-UI,PSE,Process Systems Engineering,P↗

Model: Year-Long Reforecast, Coldstart, type ICBC, Domain 02

The primary purpose of WFIP2 Model Development Team is to improve existing numerical weather prediction models in a manner that leads to improved wind forecasts in regions of complex terrain. Improvements in the models will come through better understanding of the physics associated with the wind flow in and around the wind plant across a range of temporal and spatial scales, which will be gained through WFIP2’s observational field study and analysis.

17 WIND ENERGY↗

Emulation of neutron-irradiated microstructure of austenitic 21Cr32Ni model alloy using dual-ion irradiation

Here, in this study, the capability of heavy-ion irradiation to emulate neutron irradiation was demonstrated on an austenitic 21Cr32Ni type ternary model alloy. The model alloy used in this study is chemically analogous but compositionally simpler than of alloy 800H, which is a candidate austenitic Ni alloy which has been proposed for use in Generation IV reactors. The microstructure of the 21Cr32Ni model alloy irradiated in the BOR-60 fast reactor to 17.1 dpa and 35 dpa at ~380°C was characterized using transmission electron microscopy (TEM). The 17.1 dpa BOR-60 irradiated microstructure was then compared with the microstructure of the same material developed under dual-ion (DI) irradiation using various He/dpa ratios between 0.1 and 16.6 appm He/dpa in the temperature range of 430°C-500°C. The results showed that both neutron and DI irradiation of 21Cr32Ni model alloy produced dislocations in the form of a dislocation network as well as {111}-type faulted dislocation loops, cavities, and radiation-induced Ni enrichment at radiation-induced sinks. When the dose and the He/dpa ratio were kept similar to those in neutron irradiation, DI irradiation of the 21Cr32Ni model alloy at 460°C resulted in over-nucleation of small cavities and in a high density of faulted dislocation loops compared to those observed in the fast-neutron irradiated alloy of the same heat irradiated at ~380°C. The optimal condition for reproducing the neutron-irradiated microstructure was DI irradiation at 460°C and 0.1 appm He/dpa. In that case, the faulted loop and cavity size distributions in the BOR-60 irradiated 21Cr32Ni model alloy samples closely matched with those measured in the DI irradiated 21Cr32Ni model alloy sample. The fact that the He/dpa is an order of magnitude smaller than the helium generation rate for fast neutron irradiation, stops over nucleation and allows for the development of a similar microstructure as for neutron irradiation.

21Cr32Ni model alloy↗