×

On the solution of the Neumann Poisson problem arising from a compact differencing scheme using the full multi-grid method. (English) Zbl 1360.65263

Summary: A methodology for the numerical solution of the Neumann-Poisson problem for pressure that arises during the simulation of the incompressible Navier-Stokes equations on non-staggered grids is presented in this study. A sixth order compact differencing scheme has been used for discretizing the governing equation. A general procedure for implementing the discretized form of the integral constraint is proposed. Furthermore, different strategies for handling the corners in a rectangular domain are proposed and evaluated. Solutions for a model problem with a known analytical solution have been obtained on different grids using the full multi-grid technique with Bi-Conjugate Gradient Stabilized method as the smoother. Systematic order studies have been carried out. These bring out the fact that the overall order of the numerical solution is determined by the order of the discretization used for the boundary condition in the case of the Neumann-Poisson problem.

MSC:

65N06 Finite difference methods for boundary value problems involving PDEs
65N55 Multigrid methods; domain decomposition for boundary value problems involving PDEs

Software:

Armadillo; LSQR; MINRES
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Moin, P., Mahesh, K.: Direct numerical simulation: a tool in turbulence research. Ann. Rev. Fluid Mech. 30, 539-578 (1998) · Zbl 1398.76073
[2] Lele, S.K.: Compact finite difference schemes with spectral-like resolution. J. Comput. Phys. 103, 16-42 (1992) · Zbl 0759.65006
[3] Babu, V., Korpela, S.A.: Numerical solution of the incompressible three-dimensional navier-stokes equations. Comput. Fluids 23(5), 675-691 (1994) · Zbl 0811.76041
[4] Ghia, U., Ghia, K.N., Ramamurti, R.: Multi-grid solution of Neumann pressure problem for viscous flows using primitive variables. In: 21st Aerospace Sciences Meeting. American Institute of Aeronautics and Astronautics, Reno (1983). doi:10.2514/6.1983-557 · Zbl 0664.76032
[5] Ferziger, J.H., Peric, M.: Comput. Methods Fluid Dyn. Springer, Berlin (1992)
[6] Turkel, E., Gordon, D., Gordon, R., Tsynkov, S.: Compact 2d and 3d sixth order schemes for the helmholtz equation with variable wave number. J. Comput. Phys. 232(1), 272-287 (2013) · Zbl 1291.65273
[7] Carpenter, M.H., Gottileb, D., Abarbanel, S.: The stability of numerical boundary treatments for compact high-order finite difference schemes. J. Comput. Phys. 108, 272-295 (1993) · Zbl 0791.76052
[8] Tyliszczak, A.: A high-order compact difference algorithm for half-staggered grids for laminar and turbulent incompressible flows. J. Comput. Phys. 276, 438 - 467 (2014). ISSN 0021-9991. doi:10.1016/j.jcp.2014.07.043.http://www.sciencedirect.com/science/article/pii/S0021999114005361 · Zbl 1349.76541
[9] Wilson, R.V., Demuren, A.O., Carpenter, M.: Higher order compact schemes for numerical simulation of incompressible flows. ICSE Rep. 98(13), 1-40 (1998)
[10] Strikwerda, J.C.: High-order-accurate schemes for incompressible viscous flow. Int. J. Numer. Methods Fluids 24(7), 715-734 (1997) · Zbl 0889.76050
[11] Sanderson, C.: Armadillo: an open source c++ linear algebra library for fast prototyping and computationally intensive experiments. Technical Report NICTA, pp. 1-16 (2010). http://arma.sourceforge.net/armadillo_nicta_2010.pdf · Zbl 1157.65469
[12] Babu, V., Korpela, S.A.: On the direct solution of Poisson’s equation on a non-uniform grid. J. Comput. Phys. 104(1), 93-98 (1993) · Zbl 0765.65095
[13] Choi, S-C.: Iterative methods for singular linear equations and least-squares problems. Ph.D. thesis, Institute of Computational and Mathematical Engineering, Stanford University (2006)
[14] Paige, C.C., Saunders, M.A.: Lsqr: an algorithm for sparse linear equations and sparse least squares. ACM Trans. Math. Softw. 8(1), 43-71 (1982) · Zbl 0478.65016
[15] Van der Vorst, H.A.: Bi-cgstab: a fast and smoothly converging variant of bi-cg for the solution of nonsymmetric linear systems. SIAM J. Sci. Stat. Comput. 13(2), 631-644 (1992) · Zbl 0761.65023
[16] Wang, Y., Zhang, J.: Sixth order compact scheme combined with multigrid method and extrapolation technique for 2d Poisson equation. J. Comput. Phys. 228, 137-146 (2009) · Zbl 1157.65469
[17] Briggs, W.L., Van Henson, E., McCormick, S.F.: Multigrid Tutorial. SIAM, California (2000). doi:10.1137/1.9780898719505
[18] Oberkampf, W.L., Roy, C.J.: Verification and Validation in Scientific Computing. Cambridge University Press, New York (2010) · Zbl 1211.68499
[19] Smith, G.D.: Numerical Solution of Partial Differential Equations: Finite Difference Methods. Clarendon Press, Oxford (1985) · Zbl 0576.65089
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.