Two-dimensional velocity-vorticity based LES for the solution of natural convection in a differentially heated enclosure by wavelet transform based BEM and FEM

• Published in: Engineering analysis with boundary elements 2006-08, Vol.30 (8), p.671-686
• Main Authors: Ravnik, J., Škerget, L., Hriberšek, M.
• Format: Article
• Language: English
• Subjects: Boundary element method; Classical statistical mechanics; Convection and heat transfer; Discrete wavelet transform; Exact sciences and technology; Fluid dynamics; Fundamental areas of phenomenology (including applications); General theory; Kinetic theory; Large eddy simulation; Marine; Natural convection; Physics; Rotational flow and vorticity; Statistical physics, thermodynamics, and nonlinear dynamical systems; Turbulent flows, convection, and heat transfer; Velocity-vorticity formulation
• ISSN: 0955-7997; 1873-197X
• DOI: 10.1016/j.enganabound.2006.02.008

A wavelet transform based boundary element method (BEM) numerical scheme is proposed for the solution of the kinematics equation of the velocity-vorticity formulation of Navier–Stokes equations. FEM is used to solve the kinetics equations. The proposed numerical approach is used to perform two-dimensional vorticity transfer based large eddy simulation on grids with 10 5 nodes. Turbulent natural convection in a differentially heated enclosure of aspect ratio 4 for Rayleigh number values Ra = 10 7 – 10 9 is simulated. Unstable boundary layer leads to the formation of eddies in the downstream parts of both vertical walls. At the lowest Rayleigh number value an oscillatory flow regime is observed, while the flow becomes increasingly irregular, non-repeating, unsymmetric and chaotic at higher Rayleigh number values. The transition to turbulence is studied with time series plots, temperature–vorticity phase diagrams and with power spectra. The enclosure is found to be only partially turbulent, what is qualitatively shown with second order statistics—Reynolds stresses, turbulent kinetic energy, turbulent heat fluxes and temperature variance. Heat transfer is studied via the average Nusselt number value, its time series and its relationship to the Rayleigh number value.