In this paper, a solution-adaptive algorithm is presented for the simulation of incompressible viscous flows. The framework of this method consists of an adaptive local stencil refinement algorithm and 3-points central difference discretization. The adaptive local stencil refinement is designed in such a manner that 5-points symmetric stencil is guaranteed at each interior node, so that convent...