Premium
Computational bifurcation and stability studies of the 8: 1 thermal cavity problem
Author(s) -
Salinger Andrew G.,
Lehoucq Richard B.,
Pawlowski Roger P.,
Shadid John N.
Publication year - 2002
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.392
Subject(s) - mathematics , upwind scheme , discretization , finite element method , galerkin method , steady state (chemistry) , hopf bifurcation , bifurcation , method of mean weighted residuals , rayleigh number , mathematical analysis , nonlinear system , mechanics , physics , convection , natural convection , chemistry , quantum mechanics , thermodynamics
Stability analysis algorithms coupled with a robust Newton–Krylov steady‐state iterative solver are used to understand the behavior of the 2D model problem of thermal convection in a 8: 1 differentially heated cavity. Parameter continuation methods along with bifurcation and linear stability analysis are used to study transition from steady to transient flow as a function of Rayleigh number. To carry out this study the steady‐state form of the governing PDEs is discretized using a Galerkin/least‐squares finite element formulation, and solved on parallel computers using a fully coupled Newton method and preconditioned Krylov iterative linear solvers. Linear stability analysis employing a large‐scale eigenvalue capability is used to determine the stability of the steady solutions. The boundary between steady and time‐dependent flows is determined by a Hopf bifurcation tracking capability that is used to directly track the instability with respect to the aspect ratio of the system and with respect to mesh resolution. The effect of upwinding stabilization terms in the finite element formulation on the computed value of critical Rayleigh number is investigated. The Hopf bifurcation signaling the onset of flow is determined to occur at a critical Rayleigh number of Ra =3.0604×10 5 . Copyright © 2002 John Wiley & Sons, Ltd.