Numerical analysis of bulk drag coefficient in dense vegetation by immersed boundary method