Time Domain Modeling of Axisymmetric Wave Propagation in Isotropic Elastic Media With Cefit — Cylindrical Elastodynamic Finite Integration Technique

Abstract

The Elastodynamic Finite Integration Technique (EFIT), originally developed by Fellinger et al., 1-3 represents a stable and efficient numerical code to model elastic wave propagation in linearly-elastic isotropic and anisotropic, homogeneous and heterogeneous as well as dissipative and nondissipative media. In previous works, the FIT discretization of the basic equations of linear elasticity, Hooke's law and Cauchy's equation of motion, was exclusively carried out in Cartesian coordinates. For problems in cylindrical geometries it is more suitable to use cylindrical coordinates. By that, axisymmetric problems can be treated in a two-dimensional staggered grid in the r, z-plane. The paper presents an EFIT version for axisymmetric problems in cylindrical coordinates called Cylindrical EFIT (CEFIT). After demonstrating the accuracy of the numerical code by a comparison between simulation results and analytical solutions, different examples of application are given. These examples include modeling of sound fields of ultrasonic transducers, thermoelastic laser sources, geophysical borehole probes, impact-echo measurements in layered media, and load simulations of the European Spallation Source (ESS) mercury target.

Laser-induced circumferential waves on hollow cylinder and their interaction with defects by finite element method

“…The whole hollow cylinder is divided by octahedron net. To guarantee the stability of the numerical calculation, the element size L e and discrete time-step must satisfy [8] ,…”

Numerical simulation and experimental detection of leaky Lamb waves induced by pulse laser at air-solid interfaces

“…In the finite element calculation, the element size and the time step are important to make the computation efficient and accurate. According to the work of Schubert et al [12] , the element size L e and the time step 't can be evaluated by e max 1 , 10…”

Gaseous Microemboli Sizing in Extracorporeal Circuits Using Ultrasound Backscatter

“…It begins with dimensionless linear equations of continuity and motion and performs an integration of these differential equations over controlled areas in the r-z plane. The controlled areas compose a staggered grid, leading to better accuracy than alternate numerical procedures (Schubert et al 2001;Peiffer et al 1997). The input parameters of the simulation include the following:…”