A Discrete Duality Finite Volume (DDFV) method to solve on unstructured meshes the flow problems in anisotropic nonhomogeneous porous media with full Neumann boundary conditions is proposed in the present work. We start with the derivation of the discrete problem. A result of existence and uniqueness of a solution for that problem is given thanks to the properties of its associated matrix combined with adequate assumptions on data. Their theoretical properties, namely, stability and error estimates (in discrete energy norms and L2-norm), are investigated. Numerical test is provided.