Application of computer algebra systems to the construction of the collocations and least residuals method for solving the 3D Navier--Stokes equations

V. P. Shapeev ; E. V. Vorozhtsov

Modelirovanie i analiz informacionnyh sistem, Tome 21 (2014) no. 5, pp. 131-147.

Résumé

The method of collocations and least residuals (CLR), which was proposed previously for the numerical solution of two-dimensional Navier–Stokes equations governing the stationary flows of a viscous incompressible fluid, is extended here for the three-dimensional case. The solution is sought in the implemented version of the method in the form of an expansion in the basis solenoidal functions. At all stages of the CLR method construction, a computer algebra system (CAS) is applied for the derivation and verification of the formulas of the method and for their translation into arithmetic operators of the Fortran language. For accelerating the convergence of iterations a sufficiently universal algorithm is proposed, which is simple in its implementation and is based on the use of the Krylov's subspaces. The obtained computational formulas of the CLR method were verified on the exact analytic solution of a test problem. Comparisons with the published numerical results of solving the benchmark problem of the 3D driven cubic cavity flow show that the accuracy of the results obtained by the CLR method corresponds to the known high-accuracy solutions.

Détail

Keywords: three-dimensional Navier–Stokes equations, the method of collocations and least residuals, computer algebra, cubic cavity flow, Krylov's subspaces.

Bibliographie

[1] P. Wesseling, Principles of Computational Fluid Dynamics, Springer-Verlag, Berlin, 2001

[2] J. H. Ferziger, M. Perić, Computational Methods for Fluid Dynamics, 3rd Edition, Springer-Verlag, Berlin, 2002

[3] N. G. Burago, Computational Mechanics, Moscow State Techn. University named after N.E. Bauman, Moscow, 2012

[4] Yu. A. Bystrov, S. A. Isaev, N. A. Kudryavtsev, A. I. Leontiev, Numerical Modeling of Vortex Intensification of Heat Exchange in Pipe Packages, Sudostroenie, St. Petersburg, 2005

[5] V. Girault, P. A. Raviart, Finite Element Methods for Navier–Stokes Equations: Theory and Algorithms, Springer-Verlag, London, 2011

[6] S. G. Cherny, D. V. Chirkov, V. N. Lapin, V. A. Skorospelov, Numerical Modeling of Flows in Turbomachines, Nauka, Novosibirsk, 2006

[7] A. G. Sleptsov, “Collocation-grid solution of elliptic boundary-value problems”, Modelirovanie v mekhanike, 5 (22):2 (1991), 101–126

[8] L. G. Semin, A. G. Sleptsov, V. P. Shapeev, “Collocation and least-squares method for Stokes equations”, Computat. Technologies, 1:2 (1996), 90–98

[9] V. I. Isaev, V. P. Shapeev, “High-accuracy versions of the collocations and least squares method for the numerical solution of the Navier–Stokes equations”, Computat. Math. and Math. Phys., 50:10 (2010), 1670–1681

[10] V. I. Isaev, V. P. Shapeev, “High-order accurate collocations and least squares method for solving the Navier–Stokes equations”, Dokl. Math., 85 (2012), 71–74

[11] V. P. Shapeev, E. V. Vorozhtsov, “Symbolic-numeric implementation of the method of collocations and least squares for 3D Navier–Stokes equations”, CASC 2012, LNCS, 7442, eds. V. P. Gerdt, W. Koepf, E. W. Mayr, E. V. Vorozhtsov, Springer, Heidelberg, 2012, 321–333

[12] A. N. Krylov, “On the numerical solution of the equation, which determines in technological questions the frequencies of small oscillations of material systems”, Izv. AN SSSR, 1931, no. 4, 491–539

[13] A. G. Sleptsov, “On convergence acceleration of linear iterations. II”, Modelirovanie v mekhanike, 3:5 (1989), 118–125

[14] Y. Saad, Numerical Methods for Large Eigenvalue Problems, Manchester University Press, Manchester, 1991

[15] S. Albensoeder, H. C. Kuhlmann, “Accurate three-dimensional lid-driven cavity flow”, J. Comput. Phys., 206 (2005), 536–558

[16] O. Botella, R. Peyret, “Benchmark spectral results on the lid-driven cavity flow”, Comput. Fluids, 27 (1998), 421–433