We present a thermohydraulic EGS modeling study for the hottest geo-thermal field in Türkiye, the Bozköy geothermal field. The transient 3D model, which is 300 m × 50 m × 300 m in size, was parameterized using laboratory measurements on outcrop samples, i.e., density, porosity, specific heat capacity and thermal conductivity were measured on samples. The entire modeled domain was assumed as porous medium and freshwater as working medium. We performed simulations for flow rates of 50, 75, 100, 125 and 150 l/s for a production period of 20 years, always with an injection temperature of 60 °C. The simulated production temperatures for flow rates of 150 and 50 l/s after 20 years of production are 187 and 283 °C, respectively. Thus, lower flow rates result in higher production temperatures. To determine the optimum flow rate, we plan further techno-economic analyses.