Discussions on Driven Cavity Flow

The widely studied benchmark problem, 2-D driven cavity flow problem is discussed in details in terms of physical and mathematical and also numerical aspects. A very brief literature survey on studies on the driven cavity flow is given. Based on the several numerical and experimental studies, the fact of the matter is, above moderate Reynolds numbers physically the flow in a driven cavity is not two-dimensional. However there exist numerical solutions for 2-D driven cavity flow at high Reynolds numbers

Stability Analysis in Spanwise-Periodic Double-Sided Lid-Driven Cavity Flows With Complex Cross-Sectional Profiles

Three-dimensional linear instability analyses are presented of steady two-dimensional laminar flows in the lid-driven cavity defined by [15] and further analyzed in the present volume [1], as well as in a derivative of the same geometry. It is shown that in both of the geometries considered three-dimensional BiGlobal instability leads to deviation of the flow from the two-dimensional solution; the analysis results are used to define low- and high-Reynolds number solutions by reference to the flow physics. Critical conditions for linear global instability and neutral loops are presented in both geometries

Order 10 4 speedup in global linear instability analysis using matrix formation

A unified solution framework is presented for one-, two- or three-dimensional complex non-symmetric eigenvalue problems, respectively governing linear modal instability of incompressible fluid flows in rectangular domains having two, one or no homogeneous spatial directions. The solution algorithm is based on subspace iteration in which the spatial discretization matrix is formed, stored and inverted serially. Results delivered by spectral collocation based on the Chebyshev-Gauss-Lobatto (CGL) points and a suite of high-order finite-difference methods comprising the previously employed for this type of work Dispersion-Relation-Preserving (DRP) and Padé finite-difference schemes, as well as the Summationby- parts (SBP) and the new high-order finite-difference scheme of order q (FD-q) have been compared from the point of view of accuracy and efficiency in standard validation cases of temporal local and BiGlobal linear instability. The FD-q method has been found to significantly outperform all other finite difference schemes in solving classic linear local, BiGlobal, and TriGlobal eigenvalue problems, as regards both memory and CPU time requirements. Results shown in the present study disprove the paradigm that spectral methods are superior to finite difference methods in terms of computational cost, at equal accuracy, FD-q spatial discretization delivering a speedup of ð (10 4). Consequently, accurate solutions of the three-dimensional (TriGlobal) eigenvalue problems may be solved on typical desktop computers with modest computational effort

Real-time flow simulation of indoor environments using lattice Boltzmann method

A novel lattice Boltzmann method (LBM) based 3D computational fluid dynamics (CFD) technique has been implemented on the graphics processing unit (GPU) for the purpose of simulating the indoor environment in real-time. We study the time evolution of the turbulent airflow and temperature inside a test chamber and in a simple model of a four-bed hospital room. The predicted results from LBM are compared with traditional CFD based large eddy simulations (LES). Reasonable agreement between LBM results and LES method is observed with significantly faster computational times

Finite difference and cubic interpolated profile lattice boltzmann method for prediction of two-dimensional lid-driven shallow cavity flow

In this paper, two-dimensional lid-driven cavity flow phenomena at steady state were simulated using two different scales of numerical method: the finite difference solution to the Navier–Stokes equation and the cubic interpolated pseudo-particle lattice Boltzmann method. The aspect ratio of cavity was set at 1, 2/3, 1/2 and 1/3 and the Reynolds number of 100, 400 and 1,000 for every simulation condition. The results were presented in terms of the location of the center of main vortex, the streamline plots and the velocity profiles at vertical and horizontal midsections. In this study, it is found that at the simulation of Reynolds numbers 100 and 400, both methods demonstrate a good agreement with each other; however, small discrepancies appeared for the simulation at the Reynolds number of 1,000. We also found that the number, size and formation of vortices strongly depend on the Reynolds number. The effect of the aspect ratio on the fluid flow behavior is also presented