Skip to content
Licensed Unlicensed Requires Authentication Published by De Gruyter March 15, 2021

Collocated finite-volume method for the incompressible Navier–Stokes problem

  • Kirill M. Terekhov EMAIL logo


A collocated finite-volume method for the incompressible Navier–Stokes problem is introduced. The method applies to general polyhedral grids and demonstrates higher than the first order of convergence. The velocity components and the pressure are approximated by piecewise-linear continuous and piecewise-constant fields, respectively. The method does not require artificial boundary conditions for pressure but requires stabilization term to suppress the error introduced by piecewise-constant pressure for convection-dominated problems. Both the momentum and continuity equations are approximated in a flux-conservative fashion, i.e., the conservation for both quantities is discretely exact. The attractive side of the method is a simple flux-based finite-volume construction of the scheme. Applicability of the method is demonstrated on several numerical tests using general polyhedral grids.

JEL Classification: 76M12

Funding statement: This work was supported by the Russian Science Foundation through the grant No. 19-71-10094.


[1] L. Agélas, R. Eymard, and R. Herbin, A nine-point finite volume scheme for the simulation of diffusion in heterogeneous media, Comptes Rendus Mathematique, Vol. 347, Elsevier France, Paris, 2009, pp. 673–676.10.1016/j.crma.2009.03.013Search in Google Scholar

[2] A. Arakawa, Computational design for long-term numerical integration of the equations of fluid motion: Two-dimensional incompressible flow. Part I, J. Comp. Phys., 1 (1966), No. 1, 119–143.10.1016/0021-9991(66)90015-5Search in Google Scholar

[3] A. Arakawa and V. R. Lamb, A potential enstrophy and energy conserving scheme for the shallow water equations, Monthly Weather Review, 109 (1981), No. 1, 18–36.10.1175/1520-0493(1981)109<0018:APEAEC>2.0.CO;2Search in Google Scholar

[4] E. Barragy and G. F. Carey, Stream function-vorticity driven cavity solution using pfinite elements, Comput. & Fluids, 26 (1997), No. 5, 453–468.10.1016/S0045-7930(97)00004-2Search in Google Scholar

[5] M. Bollhöfer, A robust ILU with pivoting based on monitoring the growth of the inverse factors, Linear Algebra Appl., 338 (2001), No. 1, 201–218.10.1016/S0024-3795(01)00385-8Search in Google Scholar

[6] A. Bouchnita, K. Terekhov, P. Nony, Yu. Vassilevski, and V. Volpert, A mathematical model to quantify the effects of platelet count, shear rate, and injury size on the initiation of blood coagulation under venous flow conditions, PLoS ONE, 15 (2020), No. 7, e0235392.10.1371/journal.pone.0235392Search in Google Scholar

[7] M. Braack and T. Richter, Solutions of 3D Navier–Stokes benchmark problems with adaptive finite elements, Comput. & Fluids, 35 (2006), No. 4, 372–392.10.1016/j.compfluid.2005.02.001Search in Google Scholar

[8] S. Charnyi, T. Heister, M. A. Olshanskii, and L. G. Rebholz, On conservation laws of Navier–Stokes Galerkin discretizations, J. Comp. Phys., 337 (2017), 289–308.10.1016/ in Google Scholar

[9] A. J. Chorin, The numerical solution of the Navier–Stokes equations for an incompressible fluid, Bull. Amer. Math. Soc., 73 (1967), No. 6, 928–931.10.1090/S0002-9904-1967-11853-6Search in Google Scholar

[10] C. R. Ethier and D. A. Steinman, Exact fully 3D Navier–Stokes solutions for benchmarking, Int. J. Numer. Methods Fluids, 19 (1994), No. 5, 369–375.10.1002/fld.1650190502Search in Google Scholar

[11] C. Geuzaine and J.-F. Remacle, Gmsh: A 3-D finite element mesh generator with built-in pre-and post-processing facilities, Int. J. Numer. Methods Engrg., 79 (2009), No. 11, 1309–1331.10.1002/nme.2579Search in Google Scholar

[12] U. Ghia, K. N. Ghia, and C. T. Shin, High-Re solutions for incompressible flow using the Navier–Stokes equations and a multigrid method, J. Comp. Phys., 48 (1982), No. 3, 387–411.10.1016/0021-9991(82)90058-4Search in Google Scholar

[13] P. M. Gresho, On the theory of semi-implicit projection methods for viscous incompressible flow and its implementation via a finite element method that also introduces a nearly consistent mass matrix. Part 1: Theory, Int. J. Numer. Methods Fluids, 11 (1990), No. 5, 587–620.10.1002/fld.1650110510Search in Google Scholar

[14] F. H. Harlow and J. E. Welch, Numerical calculation of time-dependent viscous incompressible flow of fluid with free surface, Phys. Fluids, 8 (1965), No. 12, 2182–2189.10.1063/1.1761178Search in Google Scholar

[15] INMOST – a toolkit for distributed mathematical modelling. URL: http://www.inmost.orgSearch in Google Scholar

[16] I. E. Kaporin, High quality preconditioning of a general symmetric positive definite matrix based on its UT U + UT R + RT U-decomposition, Numer. Lin. Alg. Appl., 5 (1998), No. 6, 483–509.10.1002/(SICI)1099-1506(199811/12)5:6<483::AID-NLA156>3.0.CO;2-7Search in Google Scholar

[17] O. A. Ladyzhenskaya, The mathematical theory of viscous incompressible flow, 12, Gordon & Breach, New York, 1969.Search in Google Scholar

[18] V. I. Lebedev, Difference analogues of orthogonal decompositions, basic differential operators and some boundary problems of mathematical physics. I, USSR Computational Mathematics and Mathematical Physics4 (1964), No. 3, 69–92.10.1016/0041-5553(64)90240-XSearch in Google Scholar

[19] N. Li, Y. Saad, and E. Chow, Crout versions of ILU for general sparse matrices, SIAM J. Sci. Comp., 25 (2003), No. 2, 716–728.10.1137/S1064827502405094Search in Google Scholar

[20] M. Olschowka and A. Neumaier, A new pivoting strategy for Gaussian elimination, Linear Algebra Appl., 240 (1996), 131–151.10.1016/0024-3795(94)00192-8Search in Google Scholar

[21] M. A. Olshanskii, K. M. Terekhov, and Yu.V. Vassilevski, An octree-based solver for the incompressible Navier–Stokes equations with enhanced stability and low dissipation, Comput. & Fluids, 84 (2013), 231–246.10.1016/j.compfluid.2013.04.027Search in Google Scholar

[22] B. Perot, Conservation properties of unstructured staggered mesh schemes, J. Comp. Phys., 159 (2000), No. 1, 58–89.10.1006/jcph.2000.6424Search in Google Scholar

[23] C. M. Rhie and W. L. Chow, Numerical study of the turbulent flow past an airfoil with trailing edge separation, AIAA J. 21 (1983), No. 11, 1525–1532.10.2514/6.1982-998Search in Google Scholar

[24] M. Schäfer, S. Turek, F. Durst, E. Krause, and R. Rannacher, Benchmark computations of laminar flow around a cylinder. In: Flow simulation with high-performance computers II, Springer, 1996, pp. 547–566.10.1007/978-3-322-89849-4_39Search in Google Scholar

[25] J. Schöberl, NETGEN An advancing front 2D/3D-mesh generator based on abstract rules, Comput. Visual. Sci., 1 (1997), No. 1, 41–52.10.1007/s007910050004Search in Google Scholar

[26] G. L. G. Sleijpen, H. A. Van der Vorst, and D. R. Fokkema, BiCGstab(l) and other hybrid Bi-CG methods, Numer. Algorithms, 7 (1994), No. 1, 75–109.10.1007/BF02141261Search in Google Scholar

[27] K. Terekhov, Cell-centered finite-volume method for heterogeneous anisotropic poromechanics problem, J. Comp. Appl. Math., 365 (2020), 112357.10.1016/ in Google Scholar

[28] K. Terekhov, B. Mallison, and H. Tchelepi, Cell-centered nonlinear finite-volume methods for the heterogeneous anisotropic diffusion problem, J. Comp. Phys., 330 (2017), 245–267.10.1016/ in Google Scholar

[29] K. Terekhov and H. Tchelepi, Cell-centered finite-volume method for elastic deformation of heterogeneous media with full-tensor properties, J. Comp. Appl. Math., 364 (2020), 112331.10.1016/ in Google Scholar

[30] K. Terekhov and Yu. Vassilevski, INMOST parallel platform for mathematical modeling and applications, In: Supercomputing. RuSCDays 2018. (Eds. V. Voevodin and S. Sobolev), Communications in Computer and Information Science, Vol. 965, Springer, Cham, 2018, pp. 230–241.10.1007/978-3-030-05807-4_20Search in Google Scholar

[31] K. Terekhov and Yu. Vassilevski, Finite volume method for coupled subsurface flow problems, I: Darcy problem, J. Comp. Phys., 395 (2019), 298–306.10.1016/ in Google Scholar

[32] K. Terekhov and Yu. Vassilevski, Mesh modification and adaptation within INMOST programming platform, In: Numerical Geometry, Grid Generation and Scientific Computing, Springer, 2019, pp. 243–255.10.1007/978-3-030-23436-2_18Search in Google Scholar

Received: 2020-02-10
Revised: 2020-07-06
Accepted: 2020-07-27
Published Online: 2021-03-15
Published in Print: 2021-03-26

© 2021 Walter de Gruyter GmbH, Berlin/Boston

Downloaded on 4.6.2023 from
Scroll to top button