Numerical modeling of thermal fracture development in unconventional reservoirs