Complete anisotropic time-dependent heat equation in KTP crystal under repetitively pulsed Gaussian beams: a numerical approach.

In this work, a thorough and detailed solution for the time-dependent heat equation for a cylindrical nonlinear potassium titanyl phosphate (KTP) crystal under a repetitively pulsed pumping source is developed. The convection and radiation boundary conditions, which are usually ignored in the literature, have been taken into account, and their importance on the temperature distribution has been discussed in detail. Moreover, the temperature dependence of thermal conductivity of KTP was considered in the calculations, and its impact is discussed. It is shown that the radiation term has a negligible effect and can be dropped safely, while the temperature dependence of thermal conductivity is more influential, such that ignorance of it brings some errors into the modeling. The time evolution of the temperature while the crystal is pumping with a train of successive Gaussian pulses until reaching equilibrium is shown. To accomplish numerical calculations, we developed a homemade code written with the finite difference time domain method in Intel Fortran (ifort) and ran it with the Linux operating system.