The existing approximations of diffusion flux in unstructured cell-centered finite volume methods are examined in detail with each other and clarified to have indefinite expressions in several respects. A new numerical approximation of diffusion flux at cell face center is then proposed, which is second-order accurate even on irregular grids and may be easily implemented in CFD code using cell-centered finite volume method with unstructured grids composed of arbitrary convex polyhedral shape.