This paper deals with the mathematical modeling of the second wave of COVID-19 and verifies the current Omicron variant pandemic data in India. We also we discussed such as uniformly bounded of the system, Equilibrium analysis and basic reproduction number $R_0$. 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$ to form numerical simulation, it helps to identify the whole system fluctuations. Finally, our proposed model validated (from real life data) the highly affected five states of COVID-19 and the Omicron variant. The algorithm guidelines are used for international arrivals, with Omicron variant cases updated by the Union Health Ministry in January 2022. Right now, the third wave is underway in India, and we conclude that it may peak by the end of May 2022.