Department of physic, Faculty of sciences and technology, University Cheikh Anta Diop of Dakar, Dakar, Senegal.
World Journal of Advanced Research and Reviews, 2025, 28(03), 945-953
Article DOI: 10.30574/wjarr.2025.28.3.4145
Received on 05 November 2025; revised on 12 December 2025; accepted on 15 December 2025
In this paper, we seek to determine the critical Reynolds number of a viscoelastic fluid flowing in a cylindrical pipe with a horizontal axis. The problem obtained is a generalized eigenvalue problem 
. A Gauss-Lobatto-Tchebyshev method was adopted to discretize this equation and the QZ algorithm combined with the Newton-Raphson method was used to determine this critical value of the Reynolds number. It is obtained by searching for two successive and very close values for which correspond two eigenvalues whose maximum real parts are respectively negative and positive. In other words, the critical value is the smallest value of the Reynolds number for which instability occurs. The code for performing this calculation was written in FORTRAN.
The flow is stable if all the real parts of the eigenvalues obtained are negative and unstable if only one of these values is positive.
Viscoelastic fluids; Linear instability; Petrov-Galerkin; Generalized eigenvalue problem; Algorithm QZ;Critical Reynolds’ number
Get Your e Certificate of Publication using below link
Preview Article PDF
Ibrahima Kama, Mamadou Yacine Ba, Alpha Malick Ndiaye and Cheikh Mbow. Theoretical and numerical study of the critical threshold of linear stability for the flow of a weakly viscoelastic fluid in a cylindrical pipe with a horizontal axis. World Journal of Advanced Research and Reviews, 2025, 28(03), 945-953. Article DOI: https://doi.org/10.30574/wjarr.2025.28.3.4145.
Copyright © 2025 Author(s) retain the copyright of this article. This article is published under the terms of the Creative Commons Attribution Liscense 4.0