This study is aimed to simulate the bifurcation behavior occurring in a two-dimensional electrical hydrodynamic nonlinear dynamical system by using the convection-diffusion-reaction (CDR) numerical method. This nonlinear system includes the Poisson equation for the external electric field that is subject to the space charge density, incompressible Navier-Stokes equations driven by the Coulomb force, continuity equation, and the charge conservation equation which describes ions distribution in the hydrodynamic field. The dynamics in the investigated nonlinear system include the pitchfork-bifurcation, frequency-doubling, Hopf-bifurcation, and the chaotic dynamics. The simulation results are presented for the case of unipolar injection between two plane electrodes. Finally, we will use the numerical analysis to understand bifurcation behavior in the electrohydrodynamic.