Premium
Bifurcation points and bifurcated branches in fluids mechanics by high‐order mesh‐free geometric progression algorithms
Author(s) -
Rammane Mohammed,
Mesmoudi Said,
Tri Abdeljalil,
Braikat Bouazza,
Damil Noureddine
Publication year - 2021
Publication title -
international journal for numerical methods in fluids
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.938
H-Index - 112
eISSN - 1097-0363
pISSN - 0271-2091
DOI - 10.1002/fld.4910
Subject(s) - mathematics , bifurcation , reynolds number , galerkin method , finite element method , moving least squares , nonlinear system , algorithm , flow (mathematics) , computational fluid dynamics , fluid mechanics , mathematical analysis , geometry , physics , mechanics , quantum mechanics , turbulence , thermodynamics
In this article, we propose to investigate numerically the steady bifurcation points and bifurcated branches in fluid mechanics by employing high‐order mesh‐free geometric progression algorithms. These algorithms are based on the use of the geometric progression (GP) in a high‐order mesh‐free approach. The first proposed algorithm is applied on a strong formulation using the moving least squares (MLS) approximation coupled with GP (HO‐MLS‐GPM). While the second proposed algorithm is applied on a weak formulation using the element‐free Galerkin (EFG) coupled also with GP (HO‐EFG‐GPM). The incompressibility condition is taken by introducing the penalty technique to transform the stationary Navier–Stokes equations verified by the pressure and velocity into ones verified by only the velocity. The high‐order mesh‐free algorithm permits to transform this nonlinear equations into a succession of linear ones. The GP allows to detect with precision the bifurcation points and the Lyapunov–Schmidt reduction is coupled with HO‐MLS‐GPM and HO‐EFG‐GPM as a continuation procedure to follow the many bifurcated branches. The aim of this resolution strategy concerns the treatment of the bifurcation phenomena for a fluid flow through an expansion in several geometries, where the steady flow becomes unstable after a critical Reynolds value. The obtained results are compared with those presented in literature and with those computed using the high‐order finite element algorithm coupled with GP.