Abstract
A fully implicit method for the solution of time-dependent Navier-Stokes equations written in terms of vorticity and velocity is presented. The field equations are discretized using second-order finite-difference technique and cast into a large sparse linear system. The numerical solution of the algebraic problem generated is obtained using a preconditioned Bi-CGSTAB method. The numerical method, developed mainly for the solution of true transient problems, also shows good stability behavior at high Rayleigh numbers. Sample results are presented for the window cavity at 104 ≤ Ra ≤ 108 and an oscillatory flow in liquid metal.
Notes
Address correspondence to Dr. F. Stella, Dipartimento di Meccanica e Aeronautica, La Sapienza, Universitá degli Studi di Roma, Via Eudossiana 18, 00814 Rome, Italy. E-mail: [email protected]