The finite element method is used to visualize the development of flow patern and temperature distribution in two-dimensional liquid-dominated geothernal system. This system is assumed as a homogenous porous host rock that is fully saturated by pure water. FEMLAB software is used to solve the differential equations of flow and transport simultaneously. The different flow pattern and temperature distribution appear with increasing simulation time.