Accurate and efficient constitutive modeling remains a cornerstone issue for solid mechanics analysis. Over the years, the LAMÉ advanced material model library has grown to address this challenge by implementing models capable of describing material systems spanning soft polymers to stiff ceramics including both isotropic and anisotropic responses. Inelastic behaviors including (visco)plasticity, damage, and fracture have all incorporated for use in various analyses. This multitude of options and flexibility, however, comes at the cost of many capabilities, features, and responses and the ensuing complexity in the resulting implementation. Therefore, to enhance confidence and enable the utilization of the LAMÉ library in application, this effort seeks to document and verify the various models in the LAMÉ library. Specifically, the broader strategy, organization, and interface of the library itself is first presented. The physical theory, numerical implementation, and user guide for a large set of models is then discussed. Importantly, a number of verification tests are performed with each model to not only have confidence in the model itself but also highlight some important response characteristics and features that may be of interest to end-users. Finally, in looking ahead to the future, approaches to add material models to this library and further expand the capabilities are presented.
This report summarizes the international collaborations conducted by Sandia funded by the US Department of Energy Office (DOE) of Nuclear Energy Spent Fuel and Waste Science & Technology (SFWST) as part of the Sandia National Laboratories Salt R&D and Salt International work packages. This report satisfies the level-three milestone M3SF-22SN010303063. Several stand-alone sections make up this summary report, each completed by the participants. The sections discuss international collaborations on geomechanical benchmarking exercises (WEIMOS), granular salt reconsolidation (KOMPASS), engineered barriers (RANGERS), numerical model comparison (DECOVALEX) and an NEA Salt Club working group on the development of scenarios as part of the performance assessment development process. Finally, we summarize events related to the US/German Workshop on Repository Research, Design and Operations. The work summarized in this annual update has occurred during the COVID-19 pandemic, and little international or domestic travel has occurred. Most of the collaborations have been conducted via email or as virtual meetings, but a slow return to travel and in-person meetings has begun.
Accurate and efficient constitutive modeling remains a cornerstone issue for solid mechanics analysis. Over the years, the LAMÉ advanced material model library has grown to address this challenge by implementing models capable of describing material systems spanning soft polymers to stiff ceramics including both isotropic and anisotropic responses. Inelastic behaviors including (visco)plasticity, damage, and fracture have all incorporated for use in various analyses. This multitude of options and flexibility, however, comes at the cost of many capabilities, features, and responses and the ensuing complexity in the resulting implementation. Therefore, to enhance confidence and enable the utilization of the LAMÉ library in application, this effort seeks to document and verify the various models in the LAMÉ library. Specifically, the broader strategy, organization, and interface of the library itself is first presented. The physical theory, numerical implementation, and user guide for a large set of models is then discussed. Importantly, a number of verification tests are performed with each model to not only have confidence in the model itself but also highlight some important response characteristics and features that may be of interest to end-users. Finally, in looking ahead to the future, approaches to add material models to this library and further expand the capabilities are presented.
Accurate and efficient constitutive modeling remains a cornerstone issue for solid mechanics analysis. Over the years, the LAMÉ advanced material model library has grown to address this challenge by implementing models capable of describing material systems spanning soft polymers to stiff ceramics including both isotropic and anisotropic responses. Inelastic behaviors including (visco)plasticity, damage, and fracture have all incorporated for use in various analyses. This multitude of options and flexibility, however, comes at the cost of many capabilities, features, and responses and the ensuing complexity in the resulting implementation. Therefore, to enhance confidence and enable the utilization of the LAMÉ library in application, this effort seeks to document and verify the various models in the LAMÉ library. Specifically, the broader strategy, organization, and interface of the library itself is first presented. The physical theory, numerical implementation, and user guide for a large set of models is then discussed. Importantly, a number of verification tests are performed with each model to not only have confidence in the model itself but also highlight some important response characteristics and features that may be of interest to end-users. Finally, in looking ahead to the future, approaches to add material models to this library and further expand the capabilities are presented.
Accurate predictions of room closure are important for hazardous waste repositories in rock salt formations, such as the Waste Isolation Pilot Plant (WIPP). When Munson and co-workers simulated several room closure experiments conducted at the WIPP during the 1980's and 1990's, their simulated closure curves closely agreed with the closure measurements. A careful review of their work, however, raised concerns and prompted the reinvestigation in this paper. To begin the reinvestigation, Munson's legacy Room D closure simulation was reasonably recreated in a current-day finite element code. Next, special care was taken to obtain numerically converged results, re-introduce the anhydrite strata intermittently ignored by Munson, and calibrate the Munson–Dawson (M–D) constitutive model for salt as much as possible from laboratory test measurements. When this new model was used to simulate Room D's closure, it under-predicted the horizontal and vertical closure rates by 2.34× and 3.10×, respectively, at 5.7 years after room excavation. As a result, the M–D model was extended to capture the newly established creep behavior at low equivalent stresses (<8MPa) and replace the Tresca with the Hosford equivalent stress. Simulations using the new M–D model over-predicted the horizontal closure rate by 1.15× and under-predicted the vertical closure rate by 1.08× at 5.7 years, averaged over three room closure experiments. Although further improvements could be made, the new model has a stronger scientific foundation than Munson's legacy model and appears ready for careful engineering use.
Based on the rationale presented, nuclear criticality is improbable after salt creep causes compaction of criticality control overpacks (CCOs) disposed at the Waste Isolation Pilot Plant, an operating repository in bedded salt for the disposal of transuranic (TRU) waste from atomic energy defense activities. For most TRU waste, the possibility of post-closure criticality is exceedingly small either because the salt neutronically isolates TRU waste canisters or because closure of a disposal room from salt creep does not sufficiently compact the low mass of fissile material. The criticality potential has been updated here because of the introduction of CCOs, which may dispose up to 380 fissile gram equivalent plutonium-239 in each container. The criticality potential is evaluated through high-fidelity geomechanical modeling of a disposal room filled with CCOs during two representative conditions: (1) large salt block fall, and (2) gradual salt compaction (without brine seepage and subsequent gas generation to permit maximum room closure). Geomechanical models of rock fall demonstrate three tiers of CCOs are not greatly disrupted. Geomechanical models of gradual room closure from salt creep predict irregular arrays of closely packed CCOs after 1000 years, when room closure has asymptotically approached maximum compaction. Criticality models of spheres and cylinders of 380 fissile gram equivalent of plutonium (as oxide) at the predicted irregular spacing demonstrate that an array of CCOs is not critical when surrounded by salt and magnesium oxide, provided the amount of hydrogenous material shipped in the CCO (usually water and plastics) is controlled or boron carbide (a neutron poison) is mixed with the fissile contents.
Accurate and efficient constitutive modeling remains a cornerstone issue for solid mechanics analysis. Over the years, the LAMÉ advanced material model library has grown to address this challenge by implementing models capable of describing material systems spanning soft polymers to stiff ceramics including both isotropic and anisotropic responses. Inelastic behaviors including (visco)plasticity, damage, and fracture have all incorporated for use in various analyses. This multitude of options and flexibility, however, comes at the cost of many capabilities, features, and responses and the ensuing complexity in the resulting implementation. Therefore, to enhance confidence and enable the utilization of the LAMÉ library in application, this effort seeks to document and verify the various models in the LAMÉ library. Specifically, the broader strategy, organization, and interface of the library itself is first presented. The physical theory, numerical implementation, and user guide for a large set of models is then discussed. Importantly, a number of verification tests are performed with each model to not only have confidence in the model itself but also highlight some important response characteristics and features that may be of interest to end-users. Finally, in looking ahead to the future, approaches to add material models to this library and further expand the capabilities are presented.
This report summarizes the international collaboration work conducted by Sandia and funded by the US Department of Energy Office (DOE) of Nuclear Energy Spent Fuel and Waste Science & Technology (SFWST) as part of the Sandia National Laboratories Salt R&D and Salt International work packages. This report satisfies the level-three milestone M3SF-20SN010303062. Several stand-alone sections make up this summary report, each completed by the participants. The sections discuss international collaborations on geomechanical benchmarking exercises (WEIMOS), granular salt reconsolidation (KOMPASS), engineered barriers (RANGERS), and model comparison (DECOVALEX). Lastly, the report summarizes a newly developed working group on the development of scenarios as part of the performance assessment development process, and the activities related to the Nuclear Energy Agency (NEA) Salt club and the US/German Workshop on Repository Research, Design and Operations.
Experimental measurements of room closure in salt repositories are valuable for understanding the evolution of the underground and for validating geomechanical models. Room closure was measured during a number of experiments at the Waste Isolation Pilot Plant (WIPP) during the 1980's and 1990's. Most rooms were excavated using a multi-pass mining sequence, where each pass necessarily destroyed some of the mining sequence closure measurement points. These destroyed points were promptly reinstalled to capture the closure after the mining pass. After the room was complete, the mining sequence closure measurement stations were supplemented with remotely read closure measurement stations. Although many aspects of these experiments were thoroughly documented, the digital copies of the closure data were inadvertently destroyed, the non-trivial process of zeroing and shifting the raw closure measurements after each mining pass was not precisely described, the various closure measurements within a given room were not directly compared on the same plot, and the measurements were collected for several years longer than previously reported. Consequently, the hand-written mining sequence closure measurements for Rooms D, B, G, and Q were located in the WIPP archives, digitized, and reanalyzed for this report. The process of reconstructing the mining sequence closure histories was documented in detail and the raw data can be found in the appendices. Within the mid-section of a given room, the reconstructed closure histories were largely consistent with other mining sequence and remotely read closure histories, which builds confidence in the experiments and suggests that plane strain is an appropriate modeling assumption. The reconstructed closure histories were also reasonably consistent with previously published results, except in one notable case: the reconstructed Room Q closure histories 30 days after excavation were about 45 % less than the corresponding closures reported in Munson's 1997 capstone paper.
Accurate and efficient constitutive modeling remains a cornerstone issue for solid mechanics analysis. Over the years, the LAMÉ advanced material model library has grown to address this challenge by implementing models capable of describing material systems spanning soft polymers to stiff ceramics including both isotropic and anisotropic responses. Inelastic behaviors including (visco)plasticity, damage, and fracture have all incorporated for use in various analyses. This multitude of options and flexibility, however, comes at the cost of many capabilities, features, and responses and the ensuing complexity in the resulting implementation. Therefore, to enhance confidence and enable the utilization of the LAMÉ library in application, this effort seeks to document and verify the various models in the LAMÉ library. Specifically, the broader strategy, organization, and interface of the library itself is first presented. The physical theory, numerical implementation, and user guide for a large set of models is then discussed. Importantly, a number of verification tests are performed with each model to not only have confidence in the model itself but also highlight some important response characteristics and features that may be of interest to end-users. Finally, in looking ahead to the future, approaches to add material models to this library and further expand the capabilities are presented
Accurate and efficient constitutive modeling remains a cornerstone issue for solid mechanics analysis. Over the years, the LAMÉ advanced material model library has grown to address this challenge by implementing models capable of describing material systems spanning soft polymers to stiff ceramics including both isotropic and anisotropic responses. Inelastic behaviors including (visco)plasticity, damage, and fracture have all incorporated for use in various analyses. This multitude of options and flexibility, however, comes at the cost of many capabilities, features, and responses and the ensuing complexity in the resulting implementation. Therefore, to enhance confidence and enable the utilization of the LAMÉ library in application, this effort seeks to document and verify the various models in the LAMÉ library. Specifically, the broader strategy, organization, and interface of the library itself is first presented. The physical theory, numerical implementation, and user guide for a large set of models is then discussed. Importantly, a number of verification tests are performed with each model to not only have confidence in the model itself but also highlight some important response characteristics and features that may be of interest to end-users. Finally, in looking ahead to the future, approaches to add material models to this library and further expand the capabilities are presented.
This report is a summary of the international collaboration work conducted by Sandia and funded by the US Department of Energy Office (DOE) of Nuclear Energy Spent Fuel and Waste Science & Technology (SFWST) as part of the Sandia National Laboratories Salt R&D and Salt International work packages. This report satisfies milestone level-three milestone M3SF-205N010303062. Several stand-alone sections make up this summary report, each completed by the participants. The first two sections discuss international collaborations on geomechanical benchmarking exercises (WEIMOS), granular salt reconsolidation (KOMPASS), engineered barriers (RANGERS), and documentation of Features, Events, and Processes (FEPs).
The Waste Isolation Pilot Plant (WIPP) is a geologic repository for defense-related nuclear waste. If left undisturbed, the virtually impermeable rock salt surrounding the repository will isolate the nuclear waste from the biosphere. If humans accidentally intrude into the repository in the future, then the likelihood of a radionuclide release to the biosphere will depend significantly on the porosity and permeability of the repository itself. Room ceilings and walls at the WIPP tend to collapse over time, causing rubble piles to form on floors of empty rooms. The surrounding rock formation will gradually compact these rubble piles until they eventually become solid salt, but the length of time for a rubble pile to reach a certain porosity and permeability is unknown. This report details the first efforts to build models to predict the porosity and permeability evolution of an empty room as it closes. Conventional geomechanical numerical methods would struggle to model empty room collapse and rubble pile consolidation, so three different meshless methods, the Immersed Isogeometric Analysis Meshfree, Reproducing Kernel Particle Method (RKPM), and the Conformal Reproducing Kernel method, were assessed. First, the meshless methods and the finite element method each simulated gradual room closure, without ceiling or wall collapse. All three methods produced equivalent room closure predictions with comparable computational speed. Second, the Immersed Isogeometric Analysis Meshfree method and RKPM simulated two-dimensional empty room collapse and rubble pile consolidation. Both methods successfully simulated large viscoplastic deformations, fracture, and rubble pile rearrangement to produce qualitatively realistic results. In addition to geomechanical simulations, the flow channels in damaged salt and crushed salt were measured using micro-computed tomography, and input into a computational fluid dynamics simulation to predict the salt's permeability. Although room for improvement exists, the current simulation approaches appear promising.