To assess the groundwater flow near high-level radioactive waste repositories, it is important to understand the effect of coupling among thermal, hydraulic, and mechanical effects. In this paper, detailed mathematical approach to model the groundwater flow near the waste form surrounded by buffer, influenced by decay heat of radioactive waste along with stress change is developed. Two cases(1) before the full expansion of buffer and (2) after the full expansion of buffer are modelled. Based on the mathematical models in this paper, detailed numerical study shall be pursued later.