In this paper, we consider a novel numerical scheme for solving incompressible flows on collocated grids. The implicit potential method utilizes an implicit potential velocity obtained from a Helmholtz decomposition for the mass conservation and employs a modified form of Bernoulli's law for the coupling of the velocity-pressure corrections. It requires the solution only of the momentum equations, does not involve the solution of additional partial differential equations for the pressure, and is applied on a collocated grid. The accuracy of the method is tested through comparison with analytical, experimental, and numerical data from the literature, and its efficiency and robustness are evaluated by solving several benchmark problems such as flow around a circular cylinder and in curved square and circular ducts.
AN IMPLICIT POTENTIAL METHOD FOR INCOMPRESSIBLE FLOWS
673The implicit potential (IPOT) scheme involves the solution only of the momentum equations, whereas the pressure field is evaluated from an algebraic relation. The fact that no partial differential equation (PDE) must be solved for the determination of the pressure field leads to a simpler algorithm compared with all other methods and enhances computational economy by decreasing the number of PDEs that need to be solved. The present scheme uses a collocated grid and requires cell vertex interpolation for the determination of the pressure derivatives.The accuracy, robustness, and computational efficiency of the proposed scheme are tested through comparison with other numerical methods and experimental measurements from the existing literature. The test cases used for the evaluation of the IPOT method include the Taylor-Green vortex problem, flow around a circular cylinder, and flow in curved ducts of square and circular crosssections. These benchmark problems are chosen so that the performance of the method is scrutinized over internal or external, steady or unsteady flow problems on Cartesian or curvilinear coordinate systems for a wide range of parameters, within the laminar flow regime. Table II. Dependence of the results for flow in curved circular pipes on the grid size. dp d´D 4 10 4 dp d´D 8 10 4 Ä D 0.05 Ä D 0.1 ÄD 0.25 Ä D 0.05 Ä D 0.1 % diff. % diff. % diff. % diff. % diff. from from from from from Grid sizeIn the Taylor-Green case study, we examine the accuracy of the IPOT method by comparing the computational results with the exact solution (28)(29)(30). The domain where we define the problem is the square OE0, 2 OE0, 2 , and we impose cyclic conditions on the boundaries. The solution procedure begins at t D 0, and it is stopped at t D 0.34, which is close to the half life of the decaying