This study developed a mathematical model to simulate temperature control of a focal cooling device using a Peltier module on human skin. Experiments using a phantom validated the model, with errors within 1.16% compared to proportional control simulations. Adding integral control improved stability and reduced steady-state error to within 0.45-1.92% of experimental temperatures. Future work will refine model parameters and test new devices to enable precise temperature regulation for medical applications.