Fractional thermo-viscoelastic Rayleigh wave modeling in laser irradiated skin tissue with nonlocal effects


Khan M. A., Mahmoud E. E., Jahangir A., Abualnaja K. M., Riaz U., YAYLACI M.

International Journal of Thermal Sciences, cilt.227, 2026 (SCI-Expanded, Scopus) identifier

  • Yayın Türü: Makale / Tam Makale
  • Cilt numarası: 227
  • Basım Tarihi: 2026
  • Doi Numarası: 10.1016/j.ijthermalsci.2026.110905
  • Dergi Adı: International Journal of Thermal Sciences
  • Derginin Tarandığı İndeksler: Science Citation Index Expanded (SCI-EXPANDED), Scopus, Compendex, INSPEC
  • Anahtar Kelimeler: Bioheat transfer, Fractional thermoelasticity, Hydrostatic stress, Laser pulse, Rayleigh waves, Sensitivity analysis, Skin tissue
  • Recep Tayyip Erdoğan Üniversitesi Adresli: Evet

Özet

Rayleigh wave propagation in skin tissue under pulsed laser heating is at the heart of advancing biomedical applications such as therapeutic ultrasonics, elastography, and laser-assisted surgery. This paper presents an integrated fractional-order thermo-viscoelastic model that combines Caputo-type three-phase-lag (TPL) heat conduction, Eringen's nonlocal elasticity, hydrostatic prestress, rotational motions, and pulsed laser thermal loading a combination that has not been addressed in the existing literature. The coupled partial differential equations are solved analytically for a homogeneous, isotropic, viscoelastic half-space, resulting in a sixth-degree characteristic equation whose admissible roots satisfy the surface wave attenuation conditions. Variance-based Sobol Global Sensitivity Analysis (GSA) with N=104 Latin Hypercube samples systematically assess the relative importance of nine model parameters for four quantities of interest: phase velocity, attenuation coefficient, penetration depth, and specific heat loss. The parameter ranges are anchored to experimental data for skin tissue where available. The results indicate that the laser pulse duration (tp) is the dominant factor for penetration depth (93% variance contribution) and phase velocity (84%), while the fractional order (α) is the dominant factor for attenuation (78%). Mechanical nonlocality (ϵ1) is always more influential than thermal nonlocality (ϵ2). The model correctly captures the classical Rayleigh wave speed in the two limiting scenarios, thus validating the model. A quantitative validation table is provided, and trends are compared with experimental elastography and photoacoustic studies. The low-sensitivity parameters are fixed to achieve model reduction with less than 3% loss of accuracy. The results offer a quantitative basis for optimizing laser pulse duration and fractional tissue properties to improve thermal dose control, diagnostic accuracy in elastography, and safety margins in laser-based medical interventions. The paper recognizes the following limitations: the model assumes homogeneous tissue, and there is no direct experimental validation.