Numerical modeling of ground subsidence due to mining