Acne vulgaris, mostly known as pimple, is common among humans in various ages. Many treatments are available to treat acnes. With the current most effective treatment of ace, isotretinoin, having side effects, photothermolysis therapy, a process that uses absorption of light to produce thermal damages, offers a new solution for acne or skin scar removals. Many different types of lasers are available for this photothermolysis procedure. For the simulation of this study, we choose to use the 1450nm laser to model the heat transfer process for the skin layers and the acne region and examine the effects of the whole procedure to the acne region and the surrounding healthy skin tissues. The model of the acne treatment involves a total of five steps of cooling process using low temperature refrigerant and four steps of heating process using the1450nm laser. To simplify the complexity of the skin tissues for this model, we use 2-D axis symmetry for the skin layers, including epidermis, sebaceous glands (acne region), and dermis, and each layer?s thickness is followed as closely as possible to human skin surfaces to make the model more realistic. With the diffusion approximation finite element analysis, we are able to evaluate the temperature profiles in different skin layers and the target acne region and gain insights on the potential acne and healthy skin thermal damages with the 1450nm laser.