Abstarct: Aabstract Hyperthermia methods are one of the conventional approaches in cancer treatment. In this method, the temperature of the cancerous tissue is raised for a specified duration to around 42 to 46 degrees Celsius. One of the existing problems with this treatment is the damage to healthy tissue due to the increase in temperature. To prevent thermal damage to healthy tissue, the heat flux or heat source within the unhealthy tissue must be such that only the temperature of the targeted areas is increased. The use of the inverse heat transfer method can resolve this issue. In this method, control parameters such as heat flux or the power of the heat source within the tissue are estimated in a way that minimizes damage to healthy tissues. In conventional heat transfer models that use Fourier’s model, the speed of thermal wave propagation is assumed to be infinite. This means that as soon as the laser hits the tissue, heat penetrates into the depth of the tissue without delay, which is often not baxsed on reality. To address the issue of the instantaneous response of changes to the applied heat flux, a non-Fourier heat transfer model is used. In the non-Fourier model, the heat transfer equations are transformed into wave equations (hyperbolic), and two time constants or delay times appear in the equations. In this thesis, a non-Fourier heat transfer model with dual phase delay is used to estimate the boundary heat flux at the tissue surface. The proposed inverse method is the conjugate gradient method along with an auxiliary problem, and the results obtained indicate a very good accuracy of the proposed method in controlling the process of hyperthermic treatment of tumors and preventing damage to healthy tissue.