Acupuncture is one of the most extensively used complementary and alternative medicine therapies worldwide. In this study, we explore the use of near-infrared light-emitting diodes (LEDs) to provide acupuncture-like physical stimulus to the skin tissue, but in a completely non-invasive way. A computational modeling framework has been developed to investigate the light-tissue interaction within a three-dimensional multi-layer model of skin tissue. Finite element-based analysis has been conducted, to obtain the spatiotemporal temperature distribution within the skin tissue, by solving Pennes’ bioheat transfer equation, coupled with the Beer-Lambert law. The irradiation profile of the LED has been experimentally characterized and imposed in the numerical model. The experimental validation of the developed model has been conducted through comparing the numerical model predictions with those obtained experimentally on the agar phantom. The effects of the LED power, treatment duration, LED distance from the skin surface, and usage of multiple LEDs on the temperature distribution attained within the skin tissue have been systematically investigated, highlighting the safe operating power of the selected LEDs. The presented information about the spatiotemporal temperature distribution, and critical factors affecting it, would assist in better optimizing the desired thermal dosage, thereby enabling a safe and effective LED-based photothermal therapy.