Laser Schlieren device was established to observe Rayleigh-Bénard-Marangoni effect produced in thick liquid depth by the diffusion of binary liquid mixtures.A recorder and a camera separately recorded dynamic and...Laser Schlieren device was established to observe Rayleigh-Bénard-Marangoni effect produced in thick liquid depth by the diffusion of binary liquid mixtures.A recorder and a camera separately recorded dynamic and static convective flow patterns during experiment.According to different signs of Ma number and Ra number,different binary liquid mixtures were selected to investigate the RBM effect induced by different driven mechanism including density gradient or surface tension gradient or combination of both.Qualitative analysis to experiment results was presented.展开更多
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.展开更多
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 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 ψ =π.展开更多
The fully nonlinear equations of gas dynamics are solved in the framework of a numerical approach in order to study the stability of the steady mode of Rayleigh-Bénard convection in compressible,viscous and heat-...The fully nonlinear equations of gas dynamics are solved in the framework of a numerical approach in order to study the stability of the steady mode of Rayleigh-Bénard convection in compressible,viscous and heat-conducting gases encapsulated in containers with no-slip boundaries and isothermal top and bottom walls.An initial linear temperature profile is assumed.A map of the possible convective modes is presented assuming the height of the region and the value of the temperature gradient as influential parameters.For a relatively small height,isobaric convection is found to take place,which is taken over by an adiabatic mode when the height exceeds the critical value,or by a super-adiabatic mode in case of a relatively high temperature gradient.In the adiabatic mode,convective flow develops due to adiabatic processes given a stable initial stratification.An analytic formula for the critical height of the region is derived taking into account and neglecting the dependence of the gas viscosity on the temperature.Moreover,an analytic formula is obtained for the upper boundary of the region of applicability of the Boussinesq approximation for incompressible gases.These models for compressible gases are relevant to practical situations such as the study of convective flows in spatially extended gas mixtures when dealing with safety issues related to hydrocarbons stored in gas stations.A dangerous situation arises when the tank is almost empty but some hydrocarbon is left at the bottom of the tank.In the presence of convective flows,the vaporized fuel is mixed with the oxidizer(air)forming a gas-vapor medium.However,if the volumetric concentration of fuel vapor(hydrocarbon)is in the interval between the lower and upper concentration limits of ignition,then the gas-vapor mixture becomes explosive and any accidental spark is sufficient to cause an emergency.展开更多
文摘Laser Schlieren device was established to observe Rayleigh-Bénard-Marangoni effect produced in thick liquid depth by the diffusion of binary liquid mixtures.A recorder and a camera separately recorded dynamic and static convective flow patterns during experiment.According to different signs of Ma number and Ra number,different binary liquid mixtures were selected to investigate the RBM effect induced by different driven mechanism including density gradient or surface tension gradient or combination of both.Qualitative analysis to experiment results was presented.
文摘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.
文摘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 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 ψ =π.
文摘The fully nonlinear equations of gas dynamics are solved in the framework of a numerical approach in order to study the stability of the steady mode of Rayleigh-Bénard convection in compressible,viscous and heat-conducting gases encapsulated in containers with no-slip boundaries and isothermal top and bottom walls.An initial linear temperature profile is assumed.A map of the possible convective modes is presented assuming the height of the region and the value of the temperature gradient as influential parameters.For a relatively small height,isobaric convection is found to take place,which is taken over by an adiabatic mode when the height exceeds the critical value,or by a super-adiabatic mode in case of a relatively high temperature gradient.In the adiabatic mode,convective flow develops due to adiabatic processes given a stable initial stratification.An analytic formula for the critical height of the region is derived taking into account and neglecting the dependence of the gas viscosity on the temperature.Moreover,an analytic formula is obtained for the upper boundary of the region of applicability of the Boussinesq approximation for incompressible gases.These models for compressible gases are relevant to practical situations such as the study of convective flows in spatially extended gas mixtures when dealing with safety issues related to hydrocarbons stored in gas stations.A dangerous situation arises when the tank is almost empty but some hydrocarbon is left at the bottom of the tank.In the presence of convective flows,the vaporized fuel is mixed with the oxidizer(air)forming a gas-vapor medium.However,if the volumetric concentration of fuel vapor(hydrocarbon)is in the interval between the lower and upper concentration limits of ignition,then the gas-vapor mixture becomes explosive and any accidental spark is sufficient to cause an emergency.