In the present study, the thermal performance of multiple rock caverns for large-scale thermal energy storage (TES) was numerically investigated for different separation distances between the caverns through heat transfer analysis using a computational fluid dynamics code, FLUENT. The thermal performance of multiple caverns was assessed in terms of the thermal stratification within the caverns and the heat loss to the surroundings, and the heating characteristics of the rock around the caverns were investigated. The results of numerical simulation showed that there was little difference in thermal performance between multiple TES caverns with different separation distances when the surrounding rock was less heated and it reached thermal steady-state, which represent the thermal states of the surrounding rock at the early and long-term operational stages of the TES caverns, respectively. However, as the separation distance decreased, the rock between the caverns reached thermal steady-state more quickly, and thus the heat loss from the caverns tended to converge rapidly to the value of heat loss occurred under thermal steady-state conditions in the surrounding rock. This result implies that the operating cost of heating the surrounding rock (i.e., rock heating) can be reduced with a reduction in the separation distance between multiple caverns, and suggests that the separation distance should be determined by considering the operating cost of rock heating as well as the construction cost of the caverns.