Numerical simulation of linear and nonlinear waves in hypoelastic solids by the CESE method

Abstract In this paper, we present a comprehensive model for the numerical simulation of linear and nonlinear elastic waves in slender rods. The mathematical model, based on conservation of mass, conservation of linear momentum, and a hypoelastic constitutive equation, consists of a set of three first-order, fully-coupled, nonlinear, hyperbolic partial differential equations. Three forms of the model equations are presented: the conservative form, the non-conservative form, and the diagonal form. The conservative form is solved numerically using the Conservation Element and Solution Element (CESE) method, a highly-accurate explicit space–time finite-volume scheme for solving nonlinear hyperbolic systems. To demonstrate this numerical approach, two elastic wave problems are directly calculated: (i) resonant standing waves arising from a time-harmonic external axial load and (ii) propagating compression waves arising from a bi-material collinear impact. Simulations of the resonance problem illustrate a linear-to-nonlinear evolution of the resonant vibrations, the emergence of super-harmonics of the forcing frequency, and the distribution of wave energy among multiple modes. For the bi-material impact problem, the CESE method successfully captures the sharp traveling wavefronts, wave reflection and transmission, and the different wave propagation speeds in each material.

[1]  Randall J. LeVeque,et al.  Finite‐volume methods for non‐linear elasticity in heterogeneous media , 2002 .

[2]  René Marklein,et al.  Numerical modeling of elastic wave propagation and scattering with EFIT — elastodynamic finite integration technique , 1995 .

[3]  Zhang Jianfeng,et al.  Elastic wave modelling in 3D heterogeneous media: 3D grid method , 2002 .

[4]  S.-T. John Yu,et al.  Direct Calculation of Wave Implosion for Detonation Initiation , 2005 .

[5]  Minghao Cai.  Acousto-Plastic deformation of metals by nonlinear stress waves , 2006 .

[6]  Moujin Zhang,et al.  Solving the MHD equations by the space-time conservation element and solution element method , 2006, J. Comput. Phys..

[7]  Sin-Chung Chang.  The Method of Space-Time Conservation Element and Solution Element-A New Approach for Solving the Navier-Stokes and Euler Equations , 1995 .

[8]  Moujin Zhang,et al.  Solving Magnetohydrodynamic Equations Without Special Treatment for Divergence-Free Magnetic Field , 2004 .

[9]  Martin Käser,et al.  Numerical simulation of 2D wave propagation on unstructured grids using explicit differential operators , 2001 .

[10]  Emmanuel Dormy,et al.  Numerical simulation of elastic wave propagation using a finite volume method , 1995 .

[11]  Koen Van Den Abeele,et al.  Two-dimensional modeling of wave propagation in materials with hysteretic nonlinearity. , 2007, The Journal of the Acoustical Society of America.

[12]  Glaucio H. Paulino,et al.  Wave propagation and dynamic analysis of smoothly graded heterogeneous continua using graded finite elements , 2007 .

[13]  Hrvoje Jasak,et al.  Application of the finite volume method and unstructured meshes to linear elasticity , 2000 .

[14]  Chang-Kee Kim,et al.  Cavity Flows in a Scramjet Engine by the Space-Time Conservation and Solution Element Method , 2004 .

[15]  A. Kulikovskii,et al.  Mathematical Aspects of Numerical Solution of Hyperbolic Systems , 1998, physics/9807053.

[16]  Sin-Chung Chang,et al.  Regular Article: The Space-Time Conservation Element and Solution Element Method: A New High-Resolution and Genuinely Multidimensional Paradigm for Solving Conservation Laws , 1999 .

[17]  Frank Schubert,et al.  Numerical time-domain modeling of linear and nonlinear ultrasonic wave propagation using finite integration techniques--theory and applications. , 2004, Ultrasonics.

[18]  Wai Ming To,et al.  A brief description of a new numerical framework for solving conservation laws — The method of space-time conservation element and solution element , 1991 .

[19]  Sin-Chung Chang,et al.  Robust and Simple Non-Reflecting Boundary Conditions for the Euler Equations - A New Approach based on the Space-Time CE/SE Method , 2003 .

[20]  Wai Ming To,et al.  Application of the space-time conservation element and solution element method to one-dimensional convection-diffusion problems , 2000 .

[21]  D. Komatitsch,et al.  Introduction to the spectral element method for three-dimensional seismic wave propagation , 1999 .

[22]  Bernd Köhler,et al.  The elastodynamic finite integration technique for waves in cylindrical geometries , 1998 .

[23]  Sin-Chung Chang,et al.  A space-time conservation element and solution element method for solving the two- and three-dimensional unsteady euler equations using quadrilateral and hexahedral meshes , 2002 .

[24]  Mrinal K. Sen,et al.  Finite‐difference modelling of S‐wave splitting in anisotropic media , 2008 .

[25]  Arben Pitarka,et al.  3D Elastic Finite-Difference Modeling of Seismic Motion Using Staggered Grids with Nonuniform Spacing , 1999 .

[26]  Susumu Kase,et al.  Studies on melt spinning. I. Fundamental equations on the dynamics of melt spinning , 1965 .