This paper conducts a Large Eddy Simulation (LES) of Rayleigh Bénard convection in a cubic cavity based on the WMLES S-Omega subgrid-scale model. For a cubic cavity with a vertical temperature difference of 6.7...This paper conducts a Large Eddy Simulation (LES) of Rayleigh Bénard convection in a cubic cavity based on the WMLES S-Omega subgrid-scale model. For a cubic cavity with a vertical temperature difference of 6.7°C and 20°C, the velocity pulsation profiles and the mean velocity profiles of the vertical section in the middle of the cubic cavity were simulated, respectively. And they are consistent with the experiment results. Furthermore, the mean velocity field of the vertical cross-section in the middle of the cavity was calculated. Structures of the mean velocity field in the two cases are similar. A counterclockwise large vortex is found to occupy the cavity, and there are two small clockwise vortices in the lower left and upper right corners, and the mean velocity fields at two different temperature differences are consistent with the experimental results. The two-dimensional instantaneous temperature field and mean temperature field with different cross-sections in the z-direction, as well as the three-dimensional instantaneous isothermal surface structure, indicate that the large-scale circulation motion within the cubic cavity is moving diagonally. In addition, the structure of the mean streamline also illustrates this viewpoint. For the reverse vortex formed at two corners in the mean streamline structure, we used the Q criterion to identify and obtain two vortex structures similar to boomerangs. The basic turbulent structure in RB thermal convection includes the rising and falling plumes generated by buoyancy effects.展开更多
The effects of result from the substitution of the classical Fourier law by the non-classical Maxwell-Cattaneo law on the Rayleigh-Bénard Magneto-convection in an electrically conducting micropolar fluid is studi...The effects of result from the substitution of the classical Fourier law by the non-classical Maxwell-Cattaneo law on the Rayleigh-Bénard Magneto-convection in an electrically conducting micropolar fluid is studied using the Galerkin technique. The eigenvalue is obtained for free-free, rigid-free and rigid-rigid velocity boundary combinations with isothermal or adiabatic temperature on the spin-vanishing boundaries. The influences of various micropolar fluid parameters are analyzed on the onset of convection. The classical approach predicts an infinite speed for the propagation of heat. The present non-classical theory involves a wave type heat transport (SECOND SOUND) and does not suffer from the physically unacceptable drawback of infinite heat propagation speed. It is found that the results are noteworthy at short times and the critical eigenvalues are less than the classical ones.展开更多
The effects of electric field and non-uniform basic temperature gradient on the onset of Rayleigh-Bénard-Marangoni convection in a micropolar fluid are studied using the Galerkin technique. The eigenvalues are ob...The effects of electric field and non-uniform basic temperature gradient on the onset of Rayleigh-Bénard-Marangoni convection in a micropolar fluid are studied using the Galerkin technique. The eigenvalues are obtained for an upper free/adiabatic and lower rigid/isothermal boundaries. The microrotation is assumed to vanish at the boundaries. A linear stability analysis is performed. The influence of various micropolar fluid parameters and electric Rayleigh number on the onset of convection has been analysed. Six different non-uniform temperature profiles are considered and their comparative influence on onset is discussed.展开更多
The onset of Rayleigh-Bnard convection in a fluid layer dispersed with phase-change-material particles was studied numerically by using the linear stability theory.The dimensionless fluctuation of specific heat Q wi...The onset of Rayleigh-Bnard convection in a fluid layer dispersed with phase-change-material particles was studied numerically by using the linear stability theory.The dimensionless fluctuation of specific heat Q with dimensionless temperature T was given as a form of sine-function Q =1+ b sin( ψT ).Two kinds of numerical methods were used separately in the calculation of critical Rayleigh number Ra _ cr and wave number k _ cr .One was the numerical integration method using Simpson 1/3 rule,and the other was the numerical difference method of Runge-Kutta with Newton-Raphson iteration. Both methods showed the same calculation results that the critical Rayleigh number Ra _ cr decreased monotonically with increase in the amplitude b of the sine-function,however,the critical wave number k _ cr did not show much difference with the amplitude b of the sine-function while ψ =π/2,but exponentially increased while ψ =π.展开更多
文摘This paper conducts a Large Eddy Simulation (LES) of Rayleigh Bénard convection in a cubic cavity based on the WMLES S-Omega subgrid-scale model. For a cubic cavity with a vertical temperature difference of 6.7°C and 20°C, the velocity pulsation profiles and the mean velocity profiles of the vertical section in the middle of the cubic cavity were simulated, respectively. And they are consistent with the experiment results. Furthermore, the mean velocity field of the vertical cross-section in the middle of the cavity was calculated. Structures of the mean velocity field in the two cases are similar. A counterclockwise large vortex is found to occupy the cavity, and there are two small clockwise vortices in the lower left and upper right corners, and the mean velocity fields at two different temperature differences are consistent with the experimental results. The two-dimensional instantaneous temperature field and mean temperature field with different cross-sections in the z-direction, as well as the three-dimensional instantaneous isothermal surface structure, indicate that the large-scale circulation motion within the cubic cavity is moving diagonally. In addition, the structure of the mean streamline also illustrates this viewpoint. For the reverse vortex formed at two corners in the mean streamline structure, we used the Q criterion to identify and obtain two vortex structures similar to boomerangs. The basic turbulent structure in RB thermal convection includes the rising and falling plumes generated by buoyancy effects.
文摘The effects of result from the substitution of the classical Fourier law by the non-classical Maxwell-Cattaneo law on the Rayleigh-Bénard Magneto-convection in an electrically conducting micropolar fluid is studied using the Galerkin technique. The eigenvalue is obtained for free-free, rigid-free and rigid-rigid velocity boundary combinations with isothermal or adiabatic temperature on the spin-vanishing boundaries. The influences of various micropolar fluid parameters are analyzed on the onset of convection. The classical approach predicts an infinite speed for the propagation of heat. The present non-classical theory involves a wave type heat transport (SECOND SOUND) and does not suffer from the physically unacceptable drawback of infinite heat propagation speed. It is found that the results are noteworthy at short times and the critical eigenvalues are less than the classical ones.
文摘The effects of electric field and non-uniform basic temperature gradient on the onset of Rayleigh-Bénard-Marangoni convection in a micropolar fluid are studied using the Galerkin technique. The eigenvalues are obtained for an upper free/adiabatic and lower rigid/isothermal boundaries. The microrotation is assumed to vanish at the boundaries. A linear stability analysis is performed. The influence of various micropolar fluid parameters and electric Rayleigh number on the onset of convection has been analysed. Six different non-uniform temperature profiles are considered and their comparative influence on onset is discussed.
文摘The onset of Rayleigh-Bnard convection in a fluid layer dispersed with phase-change-material particles was studied numerically by using the linear stability theory.The dimensionless fluctuation of specific heat Q with dimensionless temperature T was given as a form of sine-function Q =1+ b sin( ψT ).Two kinds of numerical methods were used separately in the calculation of critical Rayleigh number Ra _ cr and wave number k _ cr .One was the numerical integration method using Simpson 1/3 rule,and the other was the numerical difference method of Runge-Kutta with Newton-Raphson iteration. Both methods showed the same calculation results that the critical Rayleigh number Ra _ cr decreased monotonically with increase in the amplitude b of the sine-function,however,the critical wave number k _ cr did not show much difference with the amplitude b of the sine-function while ψ =π/2,but exponentially increased while ψ =π.