In the current study, a very simple preconditioned Navier-Stokes solver based on Roe type numerical fluxes is developed to investigate the unsteady aerodynamics around the flapping wings. A modified Osher-Chakravarthy (MOC) upwind finite-volume scheme is