A method for computation of discontinuous wave propagation in heterogeneous solids: basic algorithm description and application to one-dimensional problems

K.C. Park (1, 2) , S.J. Lim (3) , H. Huh (3)

1 Division of Ocean Systems Engineering [Daejeon]
2 Department of Aerospace Engineering Sciences
3 KAIST - Korea Advanced Institute of Science and Technology

An explicit integration algorithm for computations of discontinuous wave propagation in heterogeneous solids is presented, which is aimed at minimizing spurious oscillations when the wave fronts pass through several zones of different wave speeds. The essence of the present method is a combination of two wave capturing characteristics: a new integration formula that is obtained by pushforward–pullback operations in time designed to filter post-shock oscillations, and the central difference method that intrinsically filters front-shock oscillations. It is shown that a judicious combination of these two characteristics substantially reduces both spurious front-shock and post-shock oscillations. The performance of the new method is demonstrated as applied to wave propagation through a uniform bar with varying courant numbers, then to heterogeneous bars.

Mots clés: explicit integrator; minimal spurious oscillations; heterogeneous bar; propagation of stress waves

Domaines: Mécanique [physics.med-ph]

K.C. Park, S.J. Lim, H. Huh. A method for computation of discontinuous wave propagation in heterogeneous solids: basic algorithm description and application to one-dimensional problems. International Journal for Numerical Methods in Engineering, 2012, 91, pp.622 - 643. ⟨10.1002/nme.4285⟩. ⟨hal-01443190⟩