Stability analysis of the steady solution for incompressible lid-driven cavity flow
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.
Related
- Sharp stability and instability of stratified steady states for the incompressible porous media equation
- Pressure estimates and classifications for steady H1loc-solutions to two-dimensional incompressible Euler equations
- The unique solvability of strong solution to the multi-dimensional nonhomogeneous incompressible two-phase magnetohydrodynamic model
- Navier-Stokes flow past a translating body: attainability of the steady state as limit of an Lp weak solution
- Global nonlinear stability of the plasma-vacuum free boundary problem for the 2D ideal incompressible MHD
BibTeX
@misc{lomasov2026stability,
title = {{Stability analysis of the steady solution for incompressible lid-driven cavity flow}},
author = {Dmitri Lomasov and Petr Vabishchevich},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.21721},
eprint = {2610.21721},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.21721},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Stability analysis of the steady solution for incompressible lid-driven cavity flow
AU - Dmitri Lomasov
AU - Petr Vabishchevich
PY - 2026
DA - 2026/10/15
PB - arXiv
DO - 10.48550/arxiv.2610.21721
UR - https://arxiv.org/abs/2610.21721
AB - 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.
ER -