We present numerical simulation results of the nonlinear evolution of the (1,1) internal kink mode in the presence of various kinds of equilibrium plasma flows. The present studies are carried out in the framework of a two fluid model to extend our past investigations done with a reduced magneto-hydrodynamic (RMHD) model. Two-fluid effects are found to significantly influence the mode dynamics in a number of different ways. In the linear regime diamagnetic effects in combination with flows provide a synergistic stabilizing influence that also carries over to the nonlinear regime. In addition one observes novel symmetry breaking phenomena in the linear growth rates as well as in the nonlinear saturated states of the mode. Our study also explores the influence of strong viscosity on the mode evolution and finds interesting modifications in the real frequency of the mode in the linear regime.