A nonlinear wave model is constructed in the current work to examine the nonlinear dynamics of extremely oblique whistler waves in radiation belts. For this purpose, numerical simulation technique has been employed. The coupled normalized equations of highly oblique whistler wave and low frequency slow Alfvén wave have been developed. The ponderomotive force of high frequency whistler wave creates density perturbations in the low frequency wave. The coupled nonlinear dynamical equations of these two waves are then modeled in the form of modified Zakharov system of equations. The temporal evolution of whistler wave and the turbulent spectrum obtained suggests the energy cascade process. Nonlinear Schrödinger equation is implemented on a parallel computing platform. The goal is to study the potential performance improvements of the algorithm. The time efficiency of the simulation on a serial and parallel computing platform is studied.