RESUMO
A model of thermal ablation with application to multi-pulsed laser lithotripsy is presented. The approach is based on a one-sided Stefan-Signorini model for thermal ablation, and relies on a level-set function to represent the moving interface between the solid phase and a fictitious gas phase (representing the ablated material). The model is discretized with an embedded finite element method, wherein the interface geometry can be arbitrarily located relative to the background mesh. Nitsche's method is adopted to impose the Signorini condition on the moving interface. A bound constraint is also imposed to deal with thermal shocks that can arise during representative simulations of pulsed ablation with high-power lasers. We report simulation results based on experiments for pulsed laser ablation of wet BegoStone samples treated in air, where Begostone has been used as a phantom material for kidney stone. The model is calibrated against experimental measurements by adjusting the percentage of incoming laser energy absorbed at the surface of the stone sample. Simulation results are then validated against experimental observations for the crater area, volume, and geometry as a function of laser pulse energy and duration. Our studies illustrate how the spreading of the laser beam from the laser fiber tip with concomitantly reduced incident laser irradiance on the damaged crater surface explains trends in both the experimental observations and the model-based simulation results.
RESUMO
Holmium:yttrium-aluminum-garnet (Ho:YAG) laser lithotripsy (LL) has been the treatment of choice for kidney stone disease for more than two decades, yet the mechanisms of action are not completely clear. Besides photothermal ablation, recent evidence suggests that cavitation bubble collapse is pivotal in kidney stone dusting when the Ho:YAG laser operates at low pulse energy (Ep) and high frequency (F). In this work, we perform a comprehensive series of experiments and model-based simulations to dissect the complex physical processes in LL. Under clinically relevant dusting settings (Ep = 0.2 J, F = 20 Hz), our results suggest that majority of the irradiated laser energy (>90 %) is dissipated by heat generation in the fluid surrounding the fiber tip and the irradiated stone surface, while only about 1 % may be consumed for photothermal ablation, and less than 0.7 % is converted into the potential energy at the maximum bubble expansion. We reveal that photothermal ablation is confined locally to the laser irradiation spot, whereas cavitation erosion is most pronounced at a fiber tip-stone surface distance about 0.5 mm where multi foci ring-like damage outside the thermal ablation zone is observed. The cavitation erosion is caused by the progressively intensified collapse of jet-induced toroidal bubble near the stone surface (<100 µm), as a result of Raleigh-Taylor and Richtmyer-Meshkov instabilities. The ensuing shock wave-stone interaction and resultant leaky Rayleigh waves on the stone surface may lead to dynamic fatigue and superficial material removal under repeated bombardments of toroidal bubble collapses during dusting procedures in LL.