An efficient code to compute non-parallel steady flows and their linear stability
An algorithm which is not only simple, but fast as well as efficient has been presented. It is used to compute steady non-parallel flows and their linear stability in parameter space. The pseudo-arclength continuation method is used to trace branches of steady states as one of the parameters is varied. In order to determine the linear stability of each state computed, a generalized eigenvalue problem of large order is to be solved. Only a prescribed number of eigenvalues, those closest to the imaginary axis, are calculated by a combination of a complex mapping and the simultaneous iteration technique. The underlying linear system are solved with preconditioned Bi-CGSTAB. The method is applied efficiently to problems up to \(O(10^5)\) degrees of freedom. The authors use the algorithm to compute the bifurcation structure of steady two-dimensional Rayleigh-Bénard flows in large rectangular containers up to aspect ratio 20. The results compare qualitatively well with those of weakly nonlinear theory. The authors' hypothesis about the selection induced by saddle node bifurcation fails in this particular case. The bifurcation results do not clarify the experimentally observed increase of wavelength of pattern with increasing Rayleigh numbers. The numerical method is deficient in that the test functions used to detect simple bifurcation points are not robust and have to be reconsidered for each application. Further, the linear solver works very well for linear systems which are diagonally dominant, but for systems which are not, the amount of fill-in in the preconditioning matrix and thereby the required CPU time and use of memory increase. Notwithstanding these limitations, the numerical approach can be used to study new problems like stability of large scale ocean circulation etc..
- Large-scale eigenvalue calculations for stability analysis of steady flows on massively parallel computers
- Computational bifurcation and stability studies of the 8: 1 thermal cavity problem
- A continuation and bifurcation technique for Navier-Stokes flows
- Parallel pseudo-transient Newton-Krylov-Schwarz continuation algorithms for bifurcation analysis of incompressible sudden expansion flows
- scientific article; zbMATH DE number 1569062
- A finite element stability analysis for the Marangoni problem in a rectangular container with rigid sidewalls
- A numerical and experimental study of anomalous modes in the Taylor experiment
- A Simultaneous Iteration Algorithm for Real Matrices
- An efficient method for computing leading eigenvalues and eigenvectors of large asymmetric matrices
- Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems
- Finding leading modes of a viscous free surface flow: An asymmetric generalized eigenproblem
- scientific article; zbMATH DE number 420978 (Why is no real title available?)
- scientific article; zbMATH DE number 3841164 (Why is no real title available?)
- scientific article; zbMATH DE number 4156922 (Why is no real title available?)
- scientific article; zbMATH DE number 3930251 (Why is no real title available?)
- scientific article; zbMATH DE number 3551785 (Why is no real title available?)
- scientific article; zbMATH DE number 665522 (Why is no real title available?)
- scientific article; zbMATH DE number 3806869 (Why is no real title available?)
- scientific article; zbMATH DE number 3892457 (Why is no real title available?)
- Krylov Subspace Methods on Supercomputers
- Modal exchange mechanisms in Lapwood convection
- Numerical calculations of two-cell and single-cell Taylor flows
- Numerical computational of branch points in nonlinear equations
- On the structure of cellular solutions in Rayleigh–Bénard–Marangoni flows in small-aspect-ratio containers
- Phase-winding solutions in a finite container above the convective threshold
- The effect of distant side-walls on the evolution and stability of finite-amplitude Rayleigh-Bénard convection
- The effect of rigid sidewalls on nonlinear two-dimensional Bénard convection
- The rate of convergence of conjugate gradients
- The structure of (linearly) stable double diffusive flow patterns in a laterally heated stratified liquid
- Variations on Arnoldi's method for computing eigenelements of large unsymmetric matrices
- A parallel computer implementation of the asymptotic numerical method to study thermal convection instabilities
- Pseudospectra for matrix pencils and stability of equilibria
- Computing eigenvalues occurring in continuation methods with the Jacobi-Davidson QZ method
- Different modes of Rayleigh-Bénard instability in two- and three-dimensional rectangular enclosures
- Solution of the eigenvalue problems resulting from global non-parallel flow stability analysis
- Linear stability analysis for the differentially heated rotating annulus
- BIFURCATION TRACKING ALGORITHMS AND SOFTWARE FOR LARGE SCALE APPLICATIONS
- A pseudospectral mapping theorem for operator pencils
- A fully implicit model of the three-dimensional thermohaline ocean circulation
- Most likely noise-induced tipping of the overturning circulation in a two-dimensional Boussinesq fluid model
- Parallel pseudo-transient Newton-Krylov-Schwarz continuation algorithms for bifurcation analysis of incompressible sudden expansion flows
This page was built for publication: An efficient code to compute non-parallel steady flows and their linear stability
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q1905635)
