Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Predictive modelling”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

AlphaBuilding - Synthetic Buildings Operation Dataset

This is a synthetic building operation dataset which includes HVAC, lighting, miscellaneous electric loads (MELs) system operating conditions, occupant counts, environmental parameters, end-use and whole-building energy consumptions at 10-minute intervals. The data is created with 1395 annual simulations using the U.S. DOE detailed medium-sized reference office building, and 30 years' historical weather data in three typical climates including Miami, San Francisco, and Chicago. Three energy efficiency levels of the building and systems are considered. Assumptions regarding occupant movements, occupants' diverse temperature preferences, lighting, and MELs are adopted to reflect realistic building operations. A semantic building metadata schema - BRICK, is used to store the building metadata. The dataset is saved in a 1.2 TB of compressed HDF5 file. This dataset can be used in various applications, including building energy and load shape benchmarking, energy model calibration, evaluation of occupant and weather variability and their influences on building performance, algorithm development and testing for thermal and energy load prediction, model predictive control, policy development for reinforcement learning based building controls.

AlphaBuilding↗
DOE OSTI2022

Using aircraft measurements to characterize subgrid-scale variability of aerosol properties near the Atmospheric Radiation Measurement Southern Great Plains site

Abstract. Complex distributions of aerosol properties evolve in space and time as a function of emissions, new particle formation, coagulation, condensational growth, chemical transformation, phase changes, turbulent mixing and transport, removal processes, and ambient meteorological conditions. The ability of chemical transport models to represent the multi-scale processes affecting the life cycle of aerosols depends on their spatial resolution since aerosol properties are assumed to be constant within a grid cell. Subgrid-scale-dependent processes that affect aerosol populations could have a significant impact on the formation of particles, their growth to cloud condensation nuclei (CCN) sizes, aerosol–cloud interactions, dry deposition and rainout and hence their burdens, lifetimes, and radiative forcing. To address this issue, we characterize subgrid-scale variability in terms of measured aerosol number, size, composition, hygroscopicity, and CCN concentrations made by repeated aircraft flight paths over the Atmospheric Radiation Measurement (ARM) program's Southern Great Plains (SGP) site during the Holistic Interactions of Shallow Clouds, Aerosols and Land Ecosystem (HI-SCALE) campaign. Subgrid variability is quantified in terms of both normalized frequency distributions and percentage difference percentiles using grid spacings of 3, 9, 27, and 81 km that represent those typically used by cloud-system-resolving models as well as the current and next-generation climate models. Even though the SGP site is a rural location, surprisingly large horizontal gradients in aerosol properties were frequently observed. For example, 90 % of the 3, 9, and 27 km cell mean organic matter concentrations differed from the 81 km cell around the SGP site by as much as ∼ 46 %, large spatial variability in aerosol number concentrations and size distributions were found during new particle formation events, and consequently 90 % of the 3, 9, and 27 km cell mean CCN number concentrations differed from the 81 km cell mean by as much as ∼ 38 %. The spatial variability varied seasonally for some aerosol properties, with some having larger spatial variability during the spring and others having larger variability during the late summer. While measurements at a single surface site cannot reflect the surrounding variability of aerosol properties at a given time, aircraft measurements that are averaged within an 81 km cell were found to be similar to many, but not all, aerosol properties measured at the ground SGP site. This analysis suggests that it is reasonable to directly compare most ground SGP site aerosol measurements with coarse global climate model predictions. In addition, the variability quantified by the aircraft can be used as an uncertainty range when comparing the surface point measurements with model predictions that use coarse grid spacings.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

3D Characterization of the Columnar-to-Equiaxed Transition in Additively Manufactured Inconel 718

Additive manufacturing (AM) provides enormous processing flexibility, enabling novel part geometries and optimized designs. Access to a local heat source further permits the potential for local microstructure control on the scale of individual melt pools, which can enable local control of part properties. In order to design tailored processing strategies for target microstructures, models predicting the columnar-to-equiaxed transition must be extended to the high solidification velocities and complex thermal histories present in AM. Here, we combine 3D characterization with advanced modeling techniques to develop a more complete understanding of the solidification process and evolution of microstructure during electron beam melting (EBM) of Inconel 718. Full calibration of existing microstructure prediction models demonstrates the differences between AM processes and more conventional welding techniques, underlying the need for accurate determination of key parameters that can only be measured directly in 3D. The ability to combine multisensor data in a consistent 3D framework via data fusion algorithms is essential to fully leverage these advanced characterization approaches. Thermal modeling provides insight on microstructure development within isolated solidification events and demonstrates the role of Marangoni effects on controlling solidification behavior.

Polonsky, Andrew↗
DOE OSTI2021

A comparative study of machine learning models for predicting the state of reactive mixing

Mixing phenomena are important mechanisms controlling flow, species transport, and reaction processes in fluids and porous media. Accurate predictions of reactive mixing are critical for many Earth and environmental science problems such as contaminant fate and remediation, macroalgae growth, and plankton biomass evolution. Here, to investigate the evolution of mixing dynamics under different scenarios (e.g., anisotropy, fluctuating velocity fields), a finite-element-based numerical model was built to solve the fast, irreversible bimolecular reaction-diffusion equations to simulate a range of reactive-mixing scenarios. A total of 2,315 simulations were performed using different sets of model input parameters comprising various spatial scales of vortex structures in the velocity field, time-scales associated with velocity oscillations, the perturbation parameter for the vortex-based velocity, anisotropic dispersion contrast (i.e., ratio of longitudinal-to-transverse dispersion), and molecular diffusion. The outputs comprised concentration profiles of reactants and products. The inputs to and outputs from these simulations were concatenated into feature and label matrices, respectively, to train 20 different machine learning (ML) models intended to emulate system behavior. These 20 ML emulators, based on linear methods, Bayesian methods, ensemble learning methods, and multilayer perceptrons (MLPs), were trained to classify the state of mixing and predict three quantities of interest (QoIs) characterizing species production, decay (i.e., average concentration, square of average concentration), and degree of mixing (i.e., variances of species concentration). Unsurprisingly, linear classifiers and regressors failed to reproduce the QoIs; however, ensemble methods (classifiers and regressors) and the MLP model accurately classified the state of reactive mixing and the QoIs. Among ensemble methods, random forest and decision-tree-based AdaBoost faithfully predicted the QoIs. At run time, trained ML emulators produced results times faster than the finite-element simulations. Due to their low computational expense and high accuracy, ensemble and MLP models are excellent emulators for these numerical simulations and great utilities in uncertainty quantification exercises, which can require 1,000s of forward model runs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Simulation and Experimental Validation of an Integrated Heat Pump – Thermal Energy Storage Using a Room-Temperature Phase Change Material

As the dependence on electrical heat pumps (HPs) and intermittent renewables increases, grid strains are expected to grow. This necessitates an energy storage system to reduce the mismatch between energy supply and demand. Thus, a proposed dual-mode commercially available 14.1 kW HP was integrated with a single 22°C phase change material (PCM) thermal storage system (TES) to load-shift both cooling and heating loads. The HP-TES system was manufactured and experimentally tested using a novel test matrix based on AHRI 210/240 psychrometric conditions. Furthermore, transient dual-mode system-level HP-TES models were developed in Modelica and validated using the experimental test conditions. Base HP cooling and heating experimental tests at ambient temperatures of 35°C and −8.3°C show that the modified HP-TES maintained the rated system capacity and performance. The HP-TES discharge provided approximately 30% and 50% reductions in cooling and heating demand, respectively. The transient HP-TES models predicted system capacity and total power input for discharge and recharge operating modes within ±4% mean percentage error, and recharge power input within ±2%, with maximum errors occurring at the equipment startup. During system operation, the sources of model deviations are first-order polynomial fits of the PCM digital scanning calorimetry (DSC) data and unaccounted supercooling in the PCM during solidification. Nonetheless, the model predictions agree with the experimental tests, demonstrating the availability of robust, accurate, and validated transient models that can be used for further validation and the development of system controls.

25 ENERGY STORAGE↗
DOE OSTI2023

Measurement of the total and differential Higgs boson production cross-sections at $ \sqrt{s} $ = 13 TeV with the ATLAS detector by combining the H → ZZ* → 4ℓ and H → γγ decay channels

The total and differential Higgs boson production cross-sections are measured through a combined statistical analysis of the H → ZZ* → 4ℓ and H → γγ decay channels. The results are based on a dataset of 139 fb –1 of proton–proton collisions at a centre-of-mass energy of 13 TeV, recorded by the ATLAS detector at the Large Hadron Collider. The measured total Higgs boson production cross-section is $55.5$$^{+4.0}_{–3.8}$ pb, consistent with the Standard Model prediction of 55.6 ± 2.5 pb. All results from the two decay channels are compatible with each other, and their combination agrees with the Standard Model predictions. A combined statistical interpretation of the measured fiducial cross-sections as a function of the Higgs boson transverse momentum is performed in order to probe the Yukawa couplings to the bottom and charm quarks. A similar interpretation is performed by including also the constraints from the measurements of Higgs boson production in association with a W or Z boson in the H → bb¯ and cc¯ decay channels.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Muon g-2 Experiment and SM

The anomalous magnetic moment of muon has been used as an indication of physics beyond the Standard Model. Fermilab Muon g-2 experiment has published the most precise measurement of the anomalous magnetic moment of muon with an uncertainty of 460 part-per-billion which showed 4.2 sigma discrepancy between the experiment average and the 2020 g-2 Theory Initiative Standard Model prediction. This talk will show the recent results from the Fermilab Muon g-2 experiment while also briefly cover the status of the Standard Model prediction on calculating the anomalous magnetic moment of muon.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Deformation behavior of CP-titanium under strain path changes: Experiment and crystal plasticity modeling

The deformation behavior of commercially pure rolled titanium subjected to strain path changes is studied using experiments and a crystal plasticity model. Four different loading combinations are performed at room temperature to study the activation of slip, twinning, de-twinning and double-twinning in hexagonal closed packed titanium. The strain paths considered are: rolling direction compression (RDC) followed by normal direction compression (NDC), RDC followed by transverse direction compression (TDC), NDC followed by RDC, and NDC followed by TDC. An EBSD-based analysis of the distribution of boundary misorientation angles before and after reload was developed to analyze the evolution of {$10\bar{1}2$} tensile and {$11\bar{2}2$} compression twins. This analysis supports the model results concerning the treatment of twin reorientation. A de-twinning and double-twinning model accounting for back stress effects, an important feature of strain path changes, is implemented within the framework of the visco-plastic self-consistent (VPSC) model along with a dislocation density (DD) based hardening scheme. In the model, plasticity is accommodated by prismatic $\langle$a$\rangle$, basal $\langle$a$\rangle$ and pyramidal $\langle$c+a$\rangle$ slip modes, and {$10\bar{1}2$} tensile and {$11\bar{2}2$} compression twinning modes. The VPSC model predicts the evolution of twinning, de-twinning and double-twinning processes for both tensile and compression twinning modes under strain path change. The model predicts macroscopic stress-strain response, texture evolution, and twin volume fraction that are in agreement with experimental observations. The evolution of texture is investigated in detail by separately analyzing the twinned domains, rather than the evolution of the global texture.

36 MATERIALS SCIENCE↗
DOE OSTI2022

High Temperature Irradiation Resistant Thermocouple Drift Model Application to Commercial Thermocouples

The Advance Sensors and Instrumentation fiscal year report on the direct application of the HTIR-TC drift model on other commercial thermocouples. The Type N and K thermocouples are useful in irradiation environments up to 1100°C (or 1290°C for short bursts and/or single use), but over temperature and high fluence can cause the TCs to drift. The HTIR-TC drift model predicts the end use environment of temperature and irradiation and the total time scale to show the amount of drift a thermocouple will undergo. The model predicts that HTIR-TCs at full temperature range and fluence of 10^21 N/cm^2 will exhibit a drift of under -1%. The type N and K TCs will see a drift of -0.5% at the same fluence, but temperatures below 1100°C. Above this temperature the Type N and K TCs quickly de-calibrate and will give erroneous temperatures measurements.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Dynamic domain kinematic modelling for predicting interflow over leaky impeding layers

Traditional Boussinesq or kinematic simulations of interflow (i.e., lateral subsurface flow) assume no leakage through the impeding layer and require a no-flow boundary condition at the ridge top. However, recent analyses of many interflow-producing landscapes indicate that leaky impeding layers are common, that most interflow percolates well before reaching the toe slope, and therefore, the downslope contributing length is shorter than the hillslope length. In watersheds characterised by perched interflow over a low conductivity layer through permeable topsoil, interflow with percolation may be modelled with a kinematic wave model using a mobile upslope boundary condition defining the hillslope portion contributing interflow to valleys. Here, we developed and applied a dynamic interflow model to simulate interflow using a downslope travel distance concept such that only the active contributing length is modelled at any time. The model defines a variable active area based on the depth of the perched layer, the topographic slope and the ratio of the hydraulic conductivity of topsoil to that of the impeding layer. It incorporates a two-layer soil moisture accounting water balance analysis, a pedo-transfer function, and percolation and evaporation routines to predict interflow rates in continuous and event-based scenarios. We tested the modelling concept on two sets of data (2-year dataset of rainfall observations for the continuous simulation and a multi-day irrigation experiment for the event simulation) from a 121-m-long open interflow collection trench on an experimental hillslope at the Savannah River Site, South Carolina. The continuous model simulation partially represented the observed interflow hydrograph and perched water depth in the experimental hillslope with correlation coefficients of 0.85 and 0.35, respectively. Model performance improved significantly at event-scale analysis. The modelling approach realistically represents interflow dynamics in hillslopes with leaky impeding layers and can be integrated into catchment-scale hydrology models for more detailed hillslope process modelling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Predicting metal–metal interactions. I. The influence of strain on nanoparticle and metal adlayer stabilities

Strain-engineering of bimetallic nanomaterials is an important design strategy for developing new catalysts. In this paper, we introduce an approach for including strain effects into a recently introduced, density functional theory (DFT)-based alloy stability model. The model predicts adsorption site stabilities in nanoparticles and connects these site stabilities with catalytic reactivity and selectivity. Strain-based dependencies will increase the model’s accuracy for nanoparticles affected by finite-size effects. In addition to the stability of small nanoparticles, strain also influences the heat of adsorption of epitaxially grown metal-on-metal adlayers. In this respect, we successfully benchmark the strain-including alloy stability model with previous experimentally determined trends in the heats of adsorption of Au and Cu adlayers on Pt (111). For these systems, our model predicts stronger bimetallic interactions in the first monolayer than monometallic interactions in the second monolayer. We explicitly quantify the interplay between destabilizing strain effects and the energy gained by forming new metal–metal bonds. While tensile strain in the first Cu monolayer significantly destabilizes the adsorption strength, compressive strain in the first Au monolayer has a minimal impact on the heat of adsorption. Hence, this study introduces and, by comparison with previous experiments, validates an efficient DFT-based approach for strain-engineering the stability, and, in turn, the catalytic performance, of active sites in bimetallic alloys with atomic level resolution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Damped Ly α absorbers and atomic hydrogen in galaxies: the view of the GAEA model

ABSTRACT Using the GAEA semi-analytic model, we analyse the connection between Damped Ly α systems (DLAs) and H i in galaxies. Our state-of-the-art semi-analytic model is tuned to reproduce the local galaxy H i mass function, and that also reproduces other important galaxy properties, including the galaxy mass–gas metallicity relation. To produce catalogues of simulated DLAs we throw 105 random lines of sight in a composite simulated volume: dark matter haloes with log$(\frac{M_{200}}{ {\rm M}_{\odot }}) \ge 11.5$ are extracted from the Millennium Simulation, while for $9.2 \le \log (\frac{M_{200}}{ \mathrm{M}_{\odot }})\lt 11.5$ we use the Millennium II, and for $8 \le \log (\frac{M_{200}}{\mathrm{M}_{\odot }}) \lt 9.2$ a halo occupation distribution model. At 2 < z < 3, where observational data are more accurate, our fiducial model predicts the correct shape of the column density distribution function, but its normalization falls short of the observations, with the discrepancy increasing at higher redshift. The agreement with observations is significantly improved increasing both the H i masses and the disc radii of model galaxies by a factor of 2, as implemented ‘a posteriori’ in our 2M−2R model. In the redshift range of interest, haloes with $M_{200} \ge {10}^{11} \, \mathrm{M}_{\odot }$ give the major contribution to ΩDLA, and the typical DLA host halo mass is $\sim \!{10}^{11} \, \mathrm{M}_{\odot }$. The simulated DLA metallicity distribution is in relatively good agreement with observations, but our model predicts an excess of DLAs at low metallicities. Our results suggest possible improvements for the adopted modelling of the filtering mass and metal ejection in low-mass haloes.

Di Gioia, Serafina↗
DOE OSTI2021

Investigation of heat transport using directly driven gold spheres

Recently, heat transport was investigated using a directly driven beryllium sphere [Farmer et al., Phys. Plasmas 27, 082701 (2020)]. Models that overly restrict heat transport were rejected. This paper extends work to directly driven gold spheres where radiation loss is more important. Here, gold coated spheres are directly driven at the OMEGA laser facility at intensities of 5×10 14 W/cm 2 . Plasma conditions, laser coupling, and x-ray flux are all measured. Additionally, comparisons to 2D radiation-hydrodynamic simulations are performed. Simulations use three common heat transport models: local transport with flux limiters of f = 0.15 and f = 0.03, and the nonlocal Schurtz–Nicolai–Busquet (SNB) model. It is shown that both the SNB model and f = 0.15 match the measured plasma conditions with the SNB model better capturing the temporal evolution of electron temperature. The f = 0.03 model predicts too low of an electron density and too hot of a temperature. The measured scattered light is roughly 6% of the incident energy, the f = 0.15 and SNB models predict 0.5% uncoupled light, and f = 0.03, 38% uncoupled light. The x-ray fluxes in the f = 0.15 and SNB simulations rise too quickly and are just outside the measurement's error, while the x-ray flux in the f = 0.03 simulation is low by a factor of two-three. For these reasons, the f = 0.03 model is rejected.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

RELAP5-3D Modeling of High Temperature Test Facility (HTTF) Test PG-26

The High Temperature Test Facility (HTTF) at Oregon State University (OSU) is a scaled integral effects experiment designed to investigate transient behavior in high-temperature gas-cooled nuclear reactors with prismatic fuel and reflector blocks. Several tests have been completed, and more are still planned to at the HTTF, including depressurized conduction cooldown (DCC) and pressurized conduction cooldown (PCC) transients. This report analyses test PG-26, a progression of the Double Ended Inlet-Outlet Crossover Duct Break transient that is referred to as a DCC. PG-26 has been performed at the HTTF between May 30 and June 30, 2019. Core initial conditions (i.e., before the DCC started) have been met using low power (<100 kW) and two of ten available electric heaters. The DCC transient was initiated during the 50 th hour of the test. The break valves were opened, and hot helium from the core and cold helium from the reactor cavity simulation tank (RCST) started mixing. The gases flowed in a countercurrent fashion, where the top half of the hot duct contained hot helium that flowed in one direction and cold helium that flowed in the other direction in the bottom half of the duct. After the pressure and density reached equilibrium, the event entered a diffusion mode. The onset of a reverse natural circulation was not observed during the DCC period of the test. Version 4.4.2ie of the RELAP5-3D computer code has been used to model the HTTF PG-26 test, and results have been compared to available high-quality measured data. The model used in this study is the quality-controlled HTTF RELAP5-3D model (HTTF base 2018-04-19 QA), originally developed by P. Bayless. The report includes RELAP5-3D results of the “base calculations” as well as some sensitivities to important uncertain model inputs, such as primary helium mass flow rate, core ceramic thermal properties, as well as heat evacuation and loop friction models. Using the base RELAP5-3D model predicts a countercurrent helium flow in the hot duct observed at the beginning of the DCC, but instead of going into a molecular diffusion mode, the model predicts the onset of natural convection. Increasing friction in the core and hot duct prevents the natural convection from happening in some of the simulations. Although some temperatures are well predicted (and even overpredicted), the general tendency is to underpredict the ceramic and helium temperatures and heat removal rates during the DCC, resulting in many of the assessment findings being in minimal or insufficient agreement with the data. It is worth noting that the described discrepancies between measured data and RELAP5-3D predictions are not RELAP5-3D code limitations. More so, they reflect limitations in boundary condition and thermal property knowledge. While the RELAP5-3D calculations of the test provide some insights into what happens during the transient, and point to missing or potentially uncertain data to which the experimenters can direct their attention, the principal conclusion is that the PG-26 test data are insufficient for a system code assessment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Special Observing Period (SOP) data for the Year of Polar Prediction site Model Intercomparison Project (YOPPsiteMIP)

The rapid changes occurring in the polar regions require an improved understanding of the processes that are driving these changes. At the same time, increased human activities such as marine navigation, resource exploitation, aviation, commercial fishing, and tourism require reliable and relevant weather information. One of the primary goals of the World Meteorological Organization's Year of Polar Prediction (YOPP) project is to improve the accuracy of numerical weather prediction (NWP) at high latitudes. During YOPP, two Canadian “supersites” were commissioned and equipped with new ground-based instruments for enhanced meteorological and system process observations. Additional pre-existing supersites in Canada, the United States, Norway, Finland, and Russia also provided data from ongoing long-term observing programs. These supersites collected a wealth of observations that are well suited to address YOPP objectives. In order to increase data useability and station interoperability, novel Merged Observatory Data Files (MODFs) were created for the seven supersites over two Special Observing Periods (February to March 2018 and July to September 2018). All observations collected at the supersites were compiled into this standardized NetCDF MODF format, simplifying the process of conducting pan-Arctic NWP verification and process evaluation studies. This paper describes the seven Arctic YOPP supersites, their instrumentation, data collection and processing methods, the novel MODF format, and examples of the observations contained therein. MODFs comprise the observational contribution to the model intercomparison effort, termed YOPP site Model Intercomparison Project (YOPPsiteMIP). All YOPPsiteMIP MODFs are publicly accessible via the YOPP Data Portal (Whitehorse: https://doi.org/10.21343/a33e-j150, Huang et al., 2023a; Iqaluit: https://doi.org/10.21343/yrnf-ck57, Huang et al., 2023b; Sodankylä: https://doi.org/10.21343/m16p-pq17, O'Connor, 2023; Utqiagvik: https://doi.org/10.21343/a2dx-nq55, Akish and Morris, 2023c; Tiksi: https://doi.org/10.21343/5bwn-w881, Akish and Morris, 2023b; Ny-Ålesund: https://doi.org/10.21343/y89m-6393, Holt, 2023; and Eureka: https://doi.org/10.21343/r85j-tc61, Akish and Morris, 2023a), which is hosted by MET Norway, with corresponding output from NWP models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A microstructure-based modeling approach to predict the mechanical properties of Zr alloy with hydride precipitates

In nuclear reactors, hydrides can form in fuel cladding due to hydrogen absorption in ZIrcaloy and cause embrittlement. This work presents a microstructure-based finite element model to predict the stress-strain response of Zircaloy containing hydrides. Quantitative microstructural details extracted from scanning electron microscopy (SEM) images were used to generate heterogeneous microstructures including the morphology and spatial distribution of hydrides. The constitutive material model for zircaloy in this study is based on crystal plasticity theory which considers the hexagonal close-packed (HCP) atomic structure of Zircaloy material. The hydrides were modeled as brittle material along with a damage model. Hydride formation inside the zircoloy matrix results in residual stress. This phenomenon is also captured in this model. A parametric study has been conducted to understand the effect of volume fraction, orientation, and lamellae thickness of the hydride phase on the mechanical properties of the overall material.

Kulkarni, Shank S.↗
DOE OSTI2023

Improved Compaction Experiments and Modeling of Waste Isolation Pilot Plant Standard, Non-degraded, Waste Containers

A credible simulation of disposal room porosity at the Waste Isolation Pilot Plant (WIPP) requires a tenable compaction model for the 55-gallon waste containers within the room. A review of the legacy waste material model, however, revealed several out-of-date and untested assumptions that could affect the model’s compaction behavior. For example, the legacy model predicted non-physical tensile out-of-plane stresses under plane strain compression. (Plane strain compression is similar to waste compaction in the middle of a long drift.) Consequently, a suite of new compaction experiments were performed on containers filled with surrogate, non-degraded, waste. The new experiments involved uniaxial, triaxial, and hydrostatic compaction tests on quarter-scale and full-scale containers. Special effort was made to measure the volume strain during uniaxial and triaxial tests, so that the lateral strain could be inferred from the axial and volume strain. These experimental measurements were then used to calibrate a pressure dependent, viscoplastic, constitutive model for the homogenized compaction behavior of the waste containers. This new waste material model’s predictions agreed far better with the experimental measurements than the legacy model’s predictions, especially under triaxial and hydrostatic conditions. Under plane strain compression, the new model predicted reasonable compressive out-of-plane stresses, instead of tensile stresses. Moreover, the new model’s plane strain behavior was substantially weaker for the same strain, yet substantially stronger for the same porosity, than the legacy model’s behavior. Although room for improvement exists, the new model appears ready for prudent engineering use.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗