The value function associated with an optimal control problem subject to the Navier-Stokes equations in dimension two is analyzed. Its smoothness is established around a steady state, moreover, its derivatives are shown to satisfy a Riccati equation at the order two and generalized Lyapunov equations at the higher orders. An approximation of the optimal feedback law is then derived from the Taylor expansion of the value function. A convergence rate for the resulting controls and closed-loop systems is demonstrated.