Solution of Partial Differential Equations (PDEs) in some region R of the space of independent variables is a function, which has all the derivatives that appear on the equation, and satisfies the equation everywhere in the region R. Some linear and most nonlinear differential equations are virtually impossible to solve using exact solutions, so it is often possible to find numerical or approximate solutions for such type of problems. Therefore, numerical methods are used to approximate the solution of such type of partial differential equation to the exact solution of partial differential equation. The finite-volume method is a method for representing and evaluating partial differential equations in the form of algebraic equations [LeVeque, 2002; Toro, 1999]. In the finite volume method, volume integrals in a partial differential equation that contain a divergence term are converted to surface integrals, using the divergence theorem. These terms are then evaluated as fluxes at the surfaces of each finite volume. Because the flux entering a given volume is identical to that leaving the adjacent volume, these methods are conservative. Another advantage of the finite volume method is that it is easily formulated to allow for unstructured meshes. The method is used in many computational fluid dynamics packages.