Liposuction continues to be one of the most popular procedures performed in cosmetic surgery. As the public's demand for body contouring continues, laser lipolysis has been proposed to improve results, minimize risk, optimize patient comfort, and reduce the recovery period. Mathematical modeling of laser lipolysis could provide a better understanding of the laser lipolysis process and could determine the optimal dosage as a function of fat volume to be removed.An Optical-Thermal-Damage Model was formulated using finite-element modeling software (Femlab 3.1, Comsol Inc). The general model simulated light distribution using the diffusion approximation of the transport theory, temperature rise using the bioheat equation and laser-induced injury using the Arrhenius damage model. Biological tissue was represented by two homogenous regions (dermis and fat layer) with a nonlinear air-tissue boundary condition including free convection. Video recordings were used to gain a better understanding of the back and forth movement of the cannula during laser lipolysis in order to consider them in our mathematical model. Infrared video recordings were also performed in order to compare the actual surface temperatures to our calculations. The reduction in fat volume was determined as a function of the total applied energy and subsequently compared to clinical data reported in the literature.In patients, when using cooled tumescent anesthesia, 1064 nm Nd:YAG laser or 980 nm diode laser: (6 W, back and forth motion: 100 mm/s) give similar skin surface temperature (max: 41 degrees C). These measurements are in accordance with those obtained by mathematical modeling performed with a 1 mm cannula inserted inside the hypodermis layer at 0.8 cm below the surface. Similarly, the fat volume reduction observed in patients at 6-month follow up can be determined by mathematical modeling. This fat reduction depends on the applied energy, typically 5 cm3 for 3000 J. At last, skin retraction was observed in patients at 6-month follow up. This observation can be easily explained by mathematical modeling showing that the temperature increase inside the lower dermis is sufficient (48-50 degrees C) to induce skin tighteningLaser lipolysis can be described by a theoretical model. Fat volume reduction observed in patients is in accordance with model calculations. Due to heat diffusion, temperature elevation is also produced inside the lower reticular dermis. This interesting observation can explain remodeling of the collagenous tissue, with clinically evident skin tightening. In conclusion, while the heat generated by interstitial laser irradiation provides stimulate lipolysis of the fat cells, the collagen and elastin are also stimulated resulting in a tightening in the skin. This mathematical model should serve as a useful tool to simulate and better understand the mechanism of action of the laser lipolysis.