skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Nonlinear to Linear Elastic Code Coupling in 2-D Axisymmetric Media.

Abstract

Explosions within the earth nonlinearly deform the local media, but at typical seismological observation distances, the seismic waves can be considered linear. Although nonlinear algorithms can simulate explosions in the very near field well, these codes are computationally expensive and inaccurate at propagating these signals to great distances. A linearized wave propagation code, coupled to a nonlinear code, provides an efficient mechanism to both accurately simulate the explosion itself and to propagate these signals to distant receivers. To this end we have coupled Sandia's nonlinear simulation algorithm CTH to a linearized elastic wave propagation code for 2-D axisymmetric media (axiElasti) by passing information from the nonlinear to the linear code via time-varying boundary conditions. In this report, we first develop the 2-D axisymmetric elastic wave equations in cylindrical coordinates. Next we show how we design the time-varying boundary conditions passing information from CTH to axiElasti, and finally we demonstrate the coupling code via a simple study of the elastic radius.

Authors:
 [1]
  1. Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
Publication Date:
Research Org.:
Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
Sponsoring Org.:
USDOE National Nuclear Security Administration (NNSA)
OSTI Identifier:
1376284
Report Number(s):
SAND-2017-8848
656347
DOE Contract Number:
AC04-94AL85000
Resource Type:
Technical Report
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICS AND COMPUTING

Citation Formats

Preston, Leiph. Nonlinear to Linear Elastic Code Coupling in 2-D Axisymmetric Media.. United States: N. p., 2017. Web. doi:10.2172/1376284.
Preston, Leiph. Nonlinear to Linear Elastic Code Coupling in 2-D Axisymmetric Media.. United States. doi:10.2172/1376284.
Preston, Leiph. 2017. "Nonlinear to Linear Elastic Code Coupling in 2-D Axisymmetric Media.". United States. doi:10.2172/1376284. https://www.osti.gov/servlets/purl/1376284.
@article{osti_1376284,
title = {Nonlinear to Linear Elastic Code Coupling in 2-D Axisymmetric Media.},
author = {Preston, Leiph},
abstractNote = {Explosions within the earth nonlinearly deform the local media, but at typical seismological observation distances, the seismic waves can be considered linear. Although nonlinear algorithms can simulate explosions in the very near field well, these codes are computationally expensive and inaccurate at propagating these signals to great distances. A linearized wave propagation code, coupled to a nonlinear code, provides an efficient mechanism to both accurately simulate the explosion itself and to propagate these signals to distant receivers. To this end we have coupled Sandia's nonlinear simulation algorithm CTH to a linearized elastic wave propagation code for 2-D axisymmetric media (axiElasti) by passing information from the nonlinear to the linear code via time-varying boundary conditions. In this report, we first develop the 2-D axisymmetric elastic wave equations in cylindrical coordinates. Next we show how we design the time-varying boundary conditions passing information from CTH to axiElasti, and finally we demonstrate the coupling code via a simple study of the elastic radius.},
doi = {10.2172/1376284},
journal = {},
number = ,
volume = ,
place = {United States},
year = 2017,
month = 8
}

Technical Report:

Save / Share:
  • One of the key issues related to location and design of a nuclear waste repository is the dissipation of the heat generated by the isolated waste. The Determination of Temperature (DOT) code is a general purpose finite element heat transfer code developed by Polivka and Wilson in 1976 at the University of California, Berkeley, that is appropriate for predicting thermal environments. The code deals with linear and nonlinear transient of steady-state heat conduction in two-dimensional planar or axisymmetric representations of structures. Capabilities are provided for modeling anisotropic heterogeneous materials with temperature-dependent thermal properties and time-dependent temperature, heat flux, convection andmore » radiation boundary conditions, together with time-dependent internal heat generation. This documentation describes the mathematical model, solution techniques and computer code design. It also contains a detailed user manual. The use of DOT as an integrated part of the SCEPTER technology package has been emphasized throughout the report. Links with the thermomechanical codes MATLOC, VISCOT, and UTAH2 have been identified and described in detail, as have stand-alone applications of DOT and its use with the waste package performance assessment code WAPPA. Output from DOT consists of a steady-state temperature distribution or a set of nodal temperatures at various times for a transient problem. Strengths of the DOT code include its flexibility, wide range of possible boundary conditions, nonlinear material properties, and its efficient equation solution algorithm. Limitations include the lack of a three-dimensional analysis capability, no radiative or convective internal heat transfer, and the need to maintain a constant time step in each program execution.« less
  • VISCOT is a non-linear, transient, thermal-stress finite-element code designed to determine the viscoelastic, fiscoplastic, or elastoplastic deformation of a rock mass due to mechanical and thermal loading. The numerical solution of the nonlinear incremental equilibrium equations within VISCOT is performed by using an explicit Euler time-stepping scheme. The rock mass may be modeled as a viscoplastic or viscoelastic material. The viscoplastic material model can be described by a Tresca, von Mises, Drucker-Prager or Mohr-Coulomb yield criteria (with or without strain hardening) with an associated flow rule which can be a power or an exponential law. The viscoelastic material model withinmore » VISCOT is a temperature- and stress-dependent law which has been developed specifically for salt rock masses by Pfeifle, Mellegard and Senseny in ONWI-314 topical report (1981). Site specific parameters for this creep law at the Richton, Permian, Paradox and Vacherie salt sites have been calculated and are given in ONWI-314 topical report (1981). A major application of VISCOT (in conjunction with a SCEPTER heat transfer code such as DOT) is the thermomechanical analysis of a rock mass such as salt in which significant time-dependent nonlinear deformations are expected to occur. Such problems include room- and canister-scale studies during the excavation, operation, and long-term post-closure stages in a salt repository. In Section 1.5 of this document the code custodianship and control is described along with the status of verification, validation and peer review of this report.« less
  • Variational equations for determining the deformed configuration of a pressurized axisymmetric membrane undergoing constant acceleration are derived. The general formulation is similar to that of Tielking and Feng. A Mooney-Rivlin material model is used for the membrane and large deformations and strains are accounted for. A Ritz procedure reduces the governing equations to a nonlinear algebraic form. Tensile instabilities are discussed in the context of the analytical solution available for a sphere. In a forthcoming report a computer program written to solve these equations will be described. The analysis and software will be verified by comparison with existing finite elasticitymore » solutions. Solutions to new problems will be given. Ongoing work considering the coupling between structural deformation and the thermodynamics of inflation will be described and the effects of the acceleration on the pressure loads due to the gas will be discussed.« less
  • A computer subprogram is offered as a module for coupling an axisymmetric or plane linear elastic thin shell to two-dimensional fluid or solid continuum codes. Given the current geometry and particle velocities and new interface pressures, the subprogram computes the new geometry and particle velocities. The method is based upon the unconditionally stable Newmark or trapezoidal method of implicit time integration. The shell formulation is taken from the thin shell theory of Kraus with shear deformation. The finite element formulation utilizes piecewise linear expansions of membrane and transverse displacements as well as the rotation. Utilizing the one point quadrature ofmore » Hughes, the method is accurate in the thin shell limit.« less