One of the most important challenges for the research and industrial communities in the numerical solution of different problems is the time as well as the accuracy of computations in the field of computational fluid dynamics (CFD). In order to alleviate these issues, many researchers have focused on developing more accurate, efficient and robust algorithms and solvers for the solution of fluid flow equations. Among these algorithms are the pressure-based coupled algorithms using which the computational costs in solving different numerical problems are considerably reduced. In the present research, a pressure-based coupled algorithm for the solution of 2D incompressible laminar fluid flow equations is presented. This algorithm uses a collocated grid and is implemented in the OpenFOAM ® CFD toolbox. The effects of using this type of solvers are investigated in comparison with the use of common segregated solvers such as the SIMPLE algorithm. The implicit pressure-velocity coupling is accomplished by deriving a pressure equation in a procedure similar to the segregated SIMPLE algorithm using the Rhie–Chow interpolation technique and assembling the coefficients of the momentum and continuity equations into one diagonally dominant matrix. The extended systems of continuity and momentum equations are solved simultaneously and their convergence is accelerated by using an iterative linear solver. For the linear solution of the coupled system of equations, the common GMRES algorithm is used which is implemented as two internal and external solvers. The linear GMFOAM solver is implemented as an internal justify; MARGIN: 0cm 0cm 0pt" to have more options in choosing even any other arbitrary linear solver. The performance of the coupled approach in comparison with that of the segregated approach (here, we used the SIMPLE algorithm) is tested by solving four laminar flow problems. Namely, the flow in a square cavity, skew cavity, square tank with one inlet and one outlet and a channel with a T-junction. In this comparison, the computational performance and costs are considered. The number of iterations and the time needed to obtain converged solutions down to 3 different residuals (e.g. 10 -5 , 10 -8 and 10 -10 ) is obtained for each case using both algorithms for 3 linear solvers. The ratio of the required iterations between the coupled and segregated solvers and the percentage of the CPU time reduction in the coupled algorithm is also determined. Results show that the coupled approach substantially reduces the computational costs as compared with that of the segregated approach. They also show that the reduction rate increases as the grid size increases (i.e. the larger the size of the computational domain, the more the reduction in the computational costs). Keywords: Coupled Solver,Segregated Solver,Pressure-based,OpenFOAM