A direct numerical simulation of transition phenomena in a mixed convection channel flow

This study intends to provide an increased understanding of the laminar-turbulent transition phenomena for the buoyancy-assisted heated vertical channel flow during the early transient stage. The spectral method with weak formulation is applied in the direct numerical simulation. Initial disturbances consist of the finite-amplitude two-dimensional TS wave and a pair of three-dimensional oblique waves for the K-type disturbances. The results from the harmonic energy competitions of different wave modes show that for the buoyancy-assisted heated flow, the (k_x=1, k_z=1) or (1,1) and (1,0) modes would gain energy immediately and start to rise at almost the same rate. This phenomenon is different from that of the buoyancy-opposed flow, where the (1,1) mode decays slowly in the beginning until other modes gain enough energy and then it begins to grow quickly and overtakes the (1,0) mode after a short time period. These different transition patterns match with the experimental results that the flow transition is supercritical and subcritical for the buoyancy-assisted and -opposed flows, respectively. Buoyancy-assisted heated flow transition follows the general trend of an isothermal flow in the beginning, but the thermal-buoyant force is crucial in accelerating the instability and also causing notable differences during the subsequent transition process. All of the results for the vortex structures development, kinetic energy budget of the disturbances, flow visualization by tagged fluid particles, and the local temperature fluctuations are consistent in pointing to a clear pattern for the buoyancy-assisted heated flow transition.