arXiv:2610.21721v1 Announce Type: new Abstract: Steady solutions at high Reynolds numbers are computed using an iterative method that linearizes the nonlinear convective terms. This linearization reduces the original nonlinear boundary-value problem to a sequence of linear problems, one solved at each iteration, thereby improving the numerical robustness of the algorithm and enabling convergence over a wide range of governing parameters. The stability of steady incompressible fluid flows is studied within the framework of a linear problem for perturbations. The original steady velocity and pressure fields are subjected to small perturbations, after which the equations are linearized with respect to the perturbations. The resulting linear problem allows one to analyze the evolution of small perturbations and to establish whether the steady solutions under consideration are stable with respect to infinitesimal perturbations. In this work, the computational implementation is based on spatial discretization by the finite element method on triangular meshes. The stability of the obtained solutions is analyzed by solving a generalized algebraic eigenvalue problem formed after discretization of the linear problem for perturbations. The presence of eigenvalues with positive real part indicates instability, whereas their absence indicates stability. Calculations were performed on a sequence of refined meshes at various Reynolds numbers for the flow in a square cavity.