In this thesis, first a one-dimensional single-layer model is presented for simulating the skin burn process resulting from contact of a skin surface with a high temperature heat source. The Pennes bioheat equation is used to simulate the thermal effects of blood perfusion in tissue. During the exposure, the heat source temperature is assumed to be constant and the thermal contact resistance between the skin surface and the heat source is neglected. The Pennes bioheat equation is solved for the steady temperature distribution in skin tissue with a convective heat transfer boundary. Then, the method of separation of variables is used to solve the Pennes bioheat equation with a constant skin surface temperature during the contact, using the obtained steady temperature distribution as the initial temperature distribution. The desired temperature function is obtained with respect to the parameters involved in the problem such as blood perfusion rate, thermophysical properties of tissue and blood, heat source temperature, arterial blood temperature, body core temperature and depth, ambient temperature and convective heat transfer coefficient. The effects of blood perfusion, ambient temperature and convective heat transfer coefficient on initial and transient temperature distributions are investigated. Comparisons between the results and previous studies show a good agreement. In the second part, a double-layer model is used instead of the single-layer one and so a higher accuracy is obtained. In this double-layer model, the skin and underlying tissue is divided into two distinguished and attached layers, a tissue layer containing blood vessels and a tissue layer containing no blood vessels. In order to simulate the thermal behavior of the tissue containing blood vessels, the Pennes bioheat equation is used. The common differential equation of heat conduction is used for the tissue containing no blood vessels. The Laplace transform is performed on the equations and the initial and boundary conditions. Then, the inversion theorem for Laplace transforms and the Cauchy residue theorem are used to obtain the desired skin temperature function with respect to the parameters involved in the problem. Applying temperature histories to a damage model, the severity and degree of tissue damage can be determined. Unlike most previous studies, the presented closed-form parametric solutions allow for theoretically predicting the extent of tissue damage in different individuals, body parts, heat source temperatures, contact durations, ambient conditions and even mental and physical conditions of the injured person. In addition, due to the high accuracy and speed of obtaining results using the closed-form analytic and parametric solutions, it is hoped that this work will be helpful in assisting physicians to identify the severity and degree of burns. Keywords: Heat transfer, Skin tissue, Thermal damage, Temperature distribution, Blood perfusion, Analytical solution