Hot dry rock geothermal energy extraction experiments in artificial reservoirs created by hydraulic stimulation in naturally fractured crystalline rocks have been undertaken in several countries over the last 20 years. The experiments have had mixed results in terms of fluid recovery, system impedance, and heat extraction. Numerical models have not yet delivered a generally agreed understanding of the processes and sensitivities involved in reservoir creation and circulation. In this paper a two-dimensional fracture network model is described, which attempts to address the problems of both reservoir creation and circulation using rock mass characterization and in situ stress data as the primary inputs with a view to encapsulating our present understanding of how such systems work. The model is a prototype for a three-dimensional version currently under development and is intended mainly for engineering sensitivity studies. The basis of the model lies in approximations of fracture mechanical behavior drawn from the rock mechanics literature, a very simplified analysis of the operative physical processes, and mapping of the connectivity of fracture networks to a fine resolution regular grid. Taken together, these permit the approximate resolution of what is normally a supercomputer problem on a personal computer. The model is applied to field data gathered at Fenton Hill, New Mexico; Hijiori, Yamagata Prefecture, Japan; and Rosemanowes, Cornwall, England. ¿ American Geophysical Union 1996