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

Commit id: 7434655fe48a15527d8ace8209ba6afa2dcf9179

deletions | additions      

       

Let's describe now the other regime. Solving equation \ref{eq:eqChaleurApprox} with same experimental parameters as before but a laser energy of 200 mJ, we find a maximum increase of temperature of 60 K, giving a maximum medium temperature of about 298+60 = 360 K. While slightly below to the vaporization point of our medium, supposed close to 373K (water vaporization temperature), it can be sufficient to vaporize the medium, as it has been demonstrated that the presence of small particles like the graphite particles acts as nucleation sites and facilitate thus the vaporization of the medium at lower temperature \cite{Alimpiev_1995}. Vaporization of the particles leads by reaction to displacements inside the medium, which generate then shear waves: this constitutes the \textit{ablative regime}.  To describe physically this regime, we supposed the medium as homogeneous and isotropic, and we adopted again a 1-D model. estimated displacements in the middle of the laser beam.  Similarly to the previous section, we calculated the stress $\sigma_{zz}$, which is the sum of the axial strain component and a term given by the second law of motion due to the reaction to the particles ejected outside the medium when they reach vaporization point \cite{ready1971effects}: \begin{equation}   \sigma_{zz} = (\lambda + 2 \mu) \frac{\partial u_z}{\partial z} - \frac{1}{\rho}\frac{I^2}{(L+C(T_V-T_0))^2}  \label{eq:stressAbla} 

u_z = \frac{\zeta}{\rho \lambda}\frac{I^2}{(L+C(T_V-T_0))^2}  \label{eq:deplAblaApprox}  \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$ approximately equal to 2.9 $\mu$m. While slightly higher, this value is in good agreement with experimental displacement (2.5 $\mu$m). It is directed inside the medium, as illustrated by the white circle 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 $G_{zz}$ \cite{aki1980quantitative}:  \begin{equation}