Numerical study of divertor plasma transport with thermal force due to temperature gradient