Linear stability of incompressible fluid flow in a cavity using finite element method

Numerical methods have been applied to theoretical studies of instability and transition to turbulence. In this study an analysis of the linear stability of incompressible flow is undertaken. By means of the finite element method the two-dimensional base flow is computed numerically over a range of Reynolds numbers and is perturbed with three-dimensional disturbances. The partial differential equations governing the evolution of perturbation are obtained from the non-linear Navier–Stokes equations with a slight compressibility by using linear stability and normal mode analysis. In terms of the finite element discretization a non-singular generalized eigenproblem is formulated from these equations whose solution gives the dispersion relation between complex growth rate and wave number. This study presents stability curves to identify the critical Reynolds number and critical wavelength of the neutral mode and discusses the mechanism of instability. The stability of lid-driven cavity flow is examined. Taylor–Goertler-like vortices in the cavity are obtained by means of reconstruction of three-dimensional flows. © 1998 John Wiley & Sons, Ltd.

[1]  A 3-D Driven Cavity Flow Simulation with N3S Code , 1992 .

[2]  Robert L. Ash,et al.  Application of spectral collocation techniques to the stability of swirling flows , 1989 .

[3]  D. Malkus Eigenproblems associated with the discrete LBB condition for incompressible finite elements , 1981 .

[4]  Y. Saad Variations on Arnoldi's method for computing eigenelements of large unsymmetric matrices , 1980 .

[5]  M. F. Webster,et al.  A Taylor–Galerkin‐based algorithm for viscous incompressible flow , 1990 .

[6]  W. Arnoldi The principle of minimized iterations in the solution of the matrix eigenvalue problem , 1951 .

[7]  Y. S. Li,et al.  One‐dimensional finite element method in hydrodynamic stability , 1981 .

[8]  M. Kawahara,et al.  A three-step finite element method for unsteady incompressible flows , 1993 .

[9]  C. P. Jackson A finite-element study of the onset of vortex shedding in flow past variously shaped bodies , 1987, Journal of Fluid Mechanics.

[10]  P. Monkewitz,et al.  LOCAL AND GLOBAL INSTABILITIES IN SPATIALLY DEVELOPING FLOWS , 1990 .

[11]  William J. Stewart,et al.  A Simultaneous Iteration Algorithm for Real Matrices , 1981, TOMS.

[12]  M. Morzynski,et al.  Numerical Investigation of Wake Instabilities , 1993 .

[13]  Leonhard Kleiser,et al.  Numerical simulation of transition in wall-bounded shear flows , 1991 .

[14]  Natarajan Ramanan,et al.  Linear stability of lid‐driven cavity flow , 1994 .

[15]  R. Kessler,et al.  Nonlinear transition in three-dimensional convection , 1987, Journal of Fluid Mechanics.

[16]  M. R. Malik,et al.  Effect of curvature on three-dimensional boundary-layer stability , 1985 .

[17]  M. Kawahara,et al.  Bifurcation Analysis of Brown Tide by Reaction-Diffusion Equation Using Finite Element Method , 1997 .

[18]  S. Orszag Accurate solution of the Orr–Sommerfeld stability equation , 1971, Journal of Fluid Mechanics.

[19]  Cyrus K. Aidun,et al.  Global stability of a lid‐driven cavity with throughflow: Flow visualization studies , 1991 .

[20]  M. Morzynski,et al.  Numerical stability analysis of a flow about a cylinder , 1991 .

[21]  M. Malik Numerical methods for hypersonic boundary layer stability , 1990 .

[22]  P. Moin,et al.  Application of a Fractional-Step Method to Incompressible Navier-Stokes Equations , 1984 .

[23]  J. Ortega,et al.  Computation of selected eigenvalues of generalized eigenvalue problems , 1993 .