Numerical Simulation of Focused Shock Shear Waves in Soft Solids and a Two-Dimensional Nonlinear Homogeneous Model of the Brain

Computation of Viscoelastic Shear Shock Waves Using Finite Volume Schemes With Artificial Compressibility

The formation of shear shock waves in the brain has been proposed as one of the plausible explanations for deep intracranial injuries. In fact, such singular solutions emerge naturally in soft viscoelastic tissues under dynamic loading conditions. To improve our understanding of the mechanical processes at hand, the development of dedicated computational models is needed. The present study concerns three-dimensional numerical models of incompressible viscoelastic solids whose motion is analysed by means of shock-capturing finite volume methods. More specifically, we focus on the use of the artificial compressibility method, a technique that has been frequently employed in computational fluid dynamics. The material behaviour is deduced from the Fung-Simo quasi-linear viscoelasiticity (QLV) theory where the elastic response is of Yeoh type. We analyse the accuracy of the method and demonstrate its applicability for the study of nonlinear wave propagation in soft solids. The numerical results cover accuracy tests, shock formation and wave focusing.

Generalization of the Zabolotskaya equation to all incompressible isotropic elastic solids

We study elastic shear waves of small but finite amplitude, composed of an anti-plane shear motion and a general in-plane motion. We use a multiple scales expansion to derive an asymptotic system of coupled nonlinear equations describing their propagation in all isotropic incompressible nonlinear elastic solids, generalizing the scalar Zabolotskaya equation of compressible nonlinear elasticity. We show that for a general isotropic incompressible solid, the coupling between anti-plane and in-plane motions cannot be undone and thus conclude that linear polarization is impossible for general nonlinear two-dimensional shear waves. We then use the equations to study the evolution of a nonlinear Gaussian beam in a soft solid: we show that a pure (linearly polarized) shear beam source generates only odd harmonics, but that introducing a slight in-plane noise in the source signal leads to a second harmonic, of the same magnitude as the fifth harmonic, a phenomenon recently observed experimentally. Finally, we present examples of some special shear motions with linear polarization.

Piecewise parabolic method for propagation of shear shock waves in relaxing soft solids: One-dimensional case

Shear shock waves can be generated spontaneously deep within the brain during a traumatic injury. This recently observed behavior could be a primary mechanism for the generation of traumatic brain injuries. However, shear shock wave physics and its numerical modeling are relatively unstudied. Existing numerical solvers used in biomechanics are not designed for the extremely large Mach numbers (greater than 1) observed in the brain. Furthermore, soft solids, such as the brain, have a complex nonclassical viscoleastic response, which must be accurately modeled to capture the nonlinear wave behavior. Here, we develop a 1D inviscid velocity-stress-like system to model the propagation of shear shock waves in a homogeneous medium. Then a generalized Maxwell body is used to model a relaxing medium that can describe experimentally determined attenuation laws. Finally, the resulting system is solved numerically with the piecewise parabolic method, a high-order finite volume method. The nonlinear and the relaxing components of this method are validated with theoretical predictions. Comparisons between numerical solutions obtained for the proposed model and the experiments of plane shear shock wave propagation based on high frame-rate ultrasound imaging and tracking are shown to be in excellent agreement.

Simulation of nonlinear transient elastography: finite element model for the propagation of shear waves in homogeneous soft tissues

In this study, visco-hyperelastic Landau's model, which is widely used in acoustical physic field, is introduced into a finite element formulation. It is designed to model the nonlinear behaviour of finite amplitude shear waves in soft solids, typically, in biological tissues. This law is used in finite element models based on elastography, experiments reported in Jacob et al, the simulations results show a good agreement with the experimental study: It is observed in both that a plane shear wave generates only odd harmonics and a nonplane wave generates both odd and even harmonics in the spectral domain. In the second part, a parametric study is performed to analyse the influence of different factors on the generation of odd harmonics of plane wave. A quantitative relation is fitted between the odd harmonic amplitudes and the non-linear elastic parameter of Landau's model, which provides a practical guideline to identify the non-linearity of homogeneous tissues using elastography experiment.