The development of appropriate photothermal detection of skin diseases to meet complex clinical demands is an urgent challenge for the prevention and therapy of skin cancer. An extensive body of literature has ignored all high-order harmonics above the second order and their influences on low-order harmonics. In this paper, a new iterative numerical method is developed for solving the nonlinear thermal diffusion equation to improve nonlinear photothermal detection for the noninvasive assessment of the thickness of port-wine stain (PWS). First, based on the anatomical and structural properties of skin tissue of PWS, a nonlinear theoretical model for photothermal detection is established. Second, a corresponding nonlinear thermal diffusion equation is solved by using the new iterative numerical method and taking into account harmonics above the second-order and their effects on lower-order harmonics. Finally, the thickness and excitation light intensity of PWS samples are numerically simulated. The simulation results show that the numerical solution converges fasterand the physical meaning of the solution is clearerwith the new method than with the traditional perturbation method. The rate of change in each harmonic with the sample thickness for the new method is higher than that for the conventional perturbation method, suggesting that the proposed numerical method may provide greater detection sensitivity. The results of the study provide a theoretical basis for the clinical treatment of PWS.
Keywords: new numerical iterative method; nonlinear thermal diffusion equation; port-wine stain; sensitivity.