Logarithmic boundary layers have been observed in different regions in turbulence. However, how thermal plumes correlate to the log law of temperature and how the velocity profile changes with pressure gradient are no...Logarithmic boundary layers have been observed in different regions in turbulence. However, how thermal plumes correlate to the log law of temperature and how the velocity profile changes with pressure gradient are not fully understood. Here, we perform three-dimensional simulations of turbulence in a slim-box without the front and back walls with aspect ratio, width:depth:height=L:D:H=1:1/6:1width:depth:height=L:D:H=1:1/6:1 (respectively corresponding to xx, yy and zz coordinates), in the Rayleigh number Ra=[1×10^8,1×10^10]Ra=[1×10^8,1×10^10] for Prandtl number Pr=0.7Pr=0.7. To investigate the structures of the viscous and thermal boundary layers, we examine the velocity profiles in the streamwise and vertical directions (i.e. UU and WW) along with the mean temperature profile throughout the plume-impacting, plume-ejecting, and wind-shearing regions. The velocity profile is successfully quantified by a two-layer function of a stress length, e^+u=e^+0(z^+)3/2[1+(z^+/z^+sub)4]^1/4eu+=e^+0(z+)3/2[1+(z+/zsub+)4]1/4, as proposed by She et al.(J Fluid Mech, 2017), though it is neither \pb type nor logarithmic. In contrast, the temperature profile in the plume-ejecting region is logarithmic for all simulated cases, being attributed to the emission of thermal plumes. The coefficient of the temperature log-law, AA, can be described by composition of the thermal stress length ■■θ0■θ0■ and the thicknesses of thermal boundary layer z■subzsub■ and z?bufzbuf■, i.e. A■z?sub/(■■θ0z■buf3/2)A■zsub?/(■θ0■zbuf^3/2). The adverse pressure gradient responsible for turning the wind direction contributes to intensively emitting plumes and the logarithmic temperature profile at the plume-ejecting region. The Nusselt number scaling and the local heat flux in the slim box are consistent with previous results of the confined cells. Therefore, the slim-box RBC is a preferable system for investigating in-box kinetic and thermal structures of turbulent convection with the large-scale circulation in a fixed plane.展开更多
We propose a theoretical model for spatial variations of the temperature varianceσ2(z,r)(z is the dis-tance from the sample bottom and r the radial coordinate)in turbulent Rayleigh-Bénard convection(RBC).Adaptin...We propose a theoretical model for spatial variations of the temperature varianceσ2(z,r)(z is the dis-tance from the sample bottom and r the radial coordinate)in turbulent Rayleigh-Bénard convection(RBC).Adapting the“attached-eddy”modelofshearflowtothe plumesofRBC,wederivedanequationforσ2 which is based on the universal scaling of the normalized RBC temperature spectra.This equation in-cludes both logarithmic and power-law dependences on z/λth,whereλth is the thermal boundary layer thickness.The equation parameters depend on r and the Prandtl number Pr,but have only an extremelyweak dependence on the Rayleigh number Ra Thus our model provides a near-universal equation for thetemperature variance profile in turbulent RBC.展开更多
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 Project was supported by the National Natural Science Foundation of China (Grants 11452002, 11521091, and 11372362)MOST (China) 973 Project (Grant 2009CB724100).
文摘Logarithmic boundary layers have been observed in different regions in turbulence. However, how thermal plumes correlate to the log law of temperature and how the velocity profile changes with pressure gradient are not fully understood. Here, we perform three-dimensional simulations of turbulence in a slim-box without the front and back walls with aspect ratio, width:depth:height=L:D:H=1:1/6:1width:depth:height=L:D:H=1:1/6:1 (respectively corresponding to xx, yy and zz coordinates), in the Rayleigh number Ra=[1×10^8,1×10^10]Ra=[1×10^8,1×10^10] for Prandtl number Pr=0.7Pr=0.7. To investigate the structures of the viscous and thermal boundary layers, we examine the velocity profiles in the streamwise and vertical directions (i.e. UU and WW) along with the mean temperature profile throughout the plume-impacting, plume-ejecting, and wind-shearing regions. The velocity profile is successfully quantified by a two-layer function of a stress length, e^+u=e^+0(z^+)3/2[1+(z^+/z^+sub)4]^1/4eu+=e^+0(z+)3/2[1+(z+/zsub+)4]1/4, as proposed by She et al.(J Fluid Mech, 2017), though it is neither \pb type nor logarithmic. In contrast, the temperature profile in the plume-ejecting region is logarithmic for all simulated cases, being attributed to the emission of thermal plumes. The coefficient of the temperature log-law, AA, can be described by composition of the thermal stress length ■■θ0■θ0■ and the thicknesses of thermal boundary layer z■subzsub■ and z?bufzbuf■, i.e. A■z?sub/(■■θ0z■buf3/2)A■zsub?/(■θ0■zbuf^3/2). The adverse pressure gradient responsible for turning the wind direction contributes to intensively emitting plumes and the logarithmic temperature profile at the plume-ejecting region. The Nusselt number scaling and the local heat flux in the slim box are consistent with previous results of the confined cells. Therefore, the slim-box RBC is a preferable system for investigating in-box kinetic and thermal structures of turbulent convection with the large-scale circulation in a fixed plane.
基金the National Natural Science Foundation of China(Grants 11772111 and91952101)the Max Planck Partner Group.
文摘We propose a theoretical model for spatial variations of the temperature varianceσ2(z,r)(z is the dis-tance from the sample bottom and r the radial coordinate)in turbulent Rayleigh-Bénard convection(RBC).Adapting the“attached-eddy”modelofshearflowtothe plumesofRBC,wederivedanequationforσ2 which is based on the universal scaling of the normalized RBC temperature spectra.This equation in-cludes both logarithmic and power-law dependences on z/λth,whereλth is the thermal boundary layer thickness.The equation parameters depend on r and the Prandtl number Pr,but have only an extremelyweak dependence on the Rayleigh number Ra Thus our model provides a near-universal equation for thetemperature variance profile in turbulent RBC.
文摘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 ψ =π.