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

Commit id: 9bd20c6f334cba5862ce57d87a2fbc953bc17f93

deletions | additions      

       

\end{equation}  Using high-energy experimental parameter, $\zeta \approx \gamma^{-1} = 40 \mu$m (average depth of absorption), $\lambda$ = 2 GPa (first Lamé's coefficient of water), $L$ = 2.2 MJ.kg$^{-1}$ (vaporization latent heat of water) and $T_V-T_0$ = 373-298 = 75 K, we obtain a displacement $u_z$ of 2.9 $\mu$m. This value is again in good agreement with experimental displacement (2.5 $\mu$m). Both theoretical and experimental displacements are directed inside the medium (see white circle arrow in the Figure \ref{figElastoPVA}-(B)).  To calculate the propagation of the displacement along space and time, we modeled the ablation regime as a point force directed along Z direction with a depth of 40 $\mu$m and of constant value from -2.5 to 2.5 mm. The magnitude of the force is stored in a matrix $H_z^{abla}(y,z,t)$. Propagation as a shear wave was calculated using Green operators Displacements along Z are again equal to the convolution between $H_z^{abla}$ and  $G_{zz}$ \cite{aki1980quantitative}: \begin{equation}  G_{zz}(r,\theta,t) = \frac{\cos^2 \theta}{4\pi \rho c_p^2 r} \delta(t-\frac{r}{c_p}) + \frac{\sin^2 \theta}{4\pi \rho c_s^2 r} \delta(t-\frac{r}{c_s})  \label{eq:Gzz} 

\end{equation}  with same notations as in equation \ref{eq:Gyz}.  Displacements can again be computed by convoluting the applied force 4-D matrix $H_z^{abla}(y,z,t)$ with $G_{zz}$.  Using same values for the physical quantities as previously, results are illustrated in Figure \ref{figGreen}-(B) which represents displacement maps along Z axis, 1.0, 1.5, 2.0, 2.5 and 3.0 ms after force application. The displacement maps present many similarities with the experimental results of the Figure \ref{figElastoPVA}-(B), with initial displacement directed inside the medium and propagation of only two half cycles.  Ablation regime was also confirmed visually: at high power, a disk of paler color of the same size as the beam diameter appears at the impact location of the laser on the phantom, which is consistent with the theory of a partial vaporization of the medium. This was not observed in the low-energy experiments.