Theoretical study of lower ionospheric response to solar flares: sluggishness of D-region and peak time delay