Pol Grasland-Mongrain edited Simu disp maps.tex  over 8 years ago

Commit id: ffce19df4eee67d2faee6a7091df1e4d6c0b9e8d

deletions | additions      

       

  This physical phenomenon was then modeled numerically. The vaporization was modeled as a point force directed along Z direction with a depth of 50 $\mu$m and increasing linearly from -2.5 to 0 mm and decreasing symmetrically from 0 to 2.5 mm, to simulate an approximate Gaussian shape. Propagation as a shear wave was calculated using Green operators $G_y$ and $G_z$ as calculated by Aki Richards \cite{aki1980quantitative}:  \begin{equation}  G_y (x,y,z)= (r,\theta,z)=  \frac{\cos \beta \sin \beta}{4\pi \theta}{4\pi  \rho c_p^2 r} \delta_P + \frac{-\sin \beta \theta  \cos \beta}{4\pi \theta}{4\pi  \rho c_s^2 r} \delta_S + \frac{3\cos \beta \theta  \sin \beta}{4\pi \rho r^3} \int_{r/\alpha}^{r_\beta}{\tau \int_{r/c_p}^{r/c_s}{\tau  \delta_{NF}}\\ G_z (x,y,z)= (r,\theta,z)=  \frac{\cos^2 \beta}{4\pi \theta}{4\pi  \rho c_p^2 r} \delta_P + \frac{\sin^2 \beta}{4\pi \theta}{4\pi  \rho c_s^2 r} \delta_S + \frac{3\cos^2 \beta-1}{4\pi \theta-1}{4\pi  \rho r^3} \int_{r/\alpha}^{r_\beta}{\tau \int_{r/c_p}^{r/c_s}{\tau  \delta_{NF}} \label{eq:akirichards}  \end{equation}  where $\beta$ $\theta$  is the angle between the applied force and thedirection of the  considered point (x,y,z), (r,$\theta$,z),  $\rho$ the medium density, $c_p$ and $c_s$ the compression and shear wave speed respectively, $\delta_S$ and $\delta_P$ Dirac distribution indicating the position of the compression and shear waves along space and time, $\tau$ the time and $\delta_{NF}$ representing near-field effects. The three terms correspond respectively to the far-field compression wave, the far-field shear wave and the near-field component. Displacement can then be computed by convoluting $G_y$ and $G_z$ with time and spatial extent of the force:  \begin{equation}