This paper deals the mathematical modeling of second wave COVID19 pandemic in India, also we discussed such as uniformly bounded of the system, Equilibrium analysis and basic reproduction number R0. We calculated the analytic solutions by HPM (Homotopy Perturbation Method) and used Mathematica 12 software for numerical analysis up to 8th order approximation. It checked the error values of the approximation while the system has residual error, absolute error and h curve initial derivation of square error at up to 8th order approximation. The basic reproduction number ranges between 0.8454 and 2.0317 form numerical simulation, it helps to identify the whole system fluctuations. Finally, our proposed model validated from real life data for highly affected 5 states