  Understanding the dynamics of naturally fractured systems and fractured
reservoirs in terms of flow, heat transport and fracture stability (e.g.
induced seismicity) is important for a range of applications associated with
waste water injection, renewable energy (e.g. geothermal systems), and
greenhouse gas mitigation (e.g. geological sequestration of CO2). Here we
present the implementation and validation of an open source MATLAB package for
efficient numerical simulations of the coupled processes in fractured systems.
We take advantage of the embedded discrete fracture model that efficiently
accounts for discrete fractures. We perform a series of numerical benchmark
experiments to validate the implemented approach against analytical solutions
and established numerical methods. Finally, we use a simplified geomechanical
model and an integrated fracture stability analysis that allows estimating the
potential for shear stimulation, and thus a mechanistic assessment of induced
seismic risk during stimulation. The open source distribution of the source
code and results can be used as a blue print for the re-implementation of the
method in a high performance computing (HPC) framework or as a standalone
simulation package for investigating TH(m) problems in geothermal reservoirs.
