The diffusive transport properties in microscale convection flows are studied by using the direct simulation Monte Carlo method. The effective diffusion coefficient D* is computed from the mean square displacements of simulated molecules based on the Einstein diffusion equation D*=<delta_x^2(t)>/2t. Two typical convection flows, namely, thermal creep convection and Rayleigh–Bénard convection, are investigated. The thermal creep convection in our simulation is in the noncontinuum regime, with the characteristic scale of the vortex varying from 1 to 100 molecular mean free paths. The diffusion is shown to be enhanced only when the vortex scale exceeds a certain critical value, while the diffusion is reduced when the vortex scale is less than the critical value. The reason for phenomenon of diffusion reduction in the noncontinuum regime is that the reduction effect due to solid wall is dominant while the enhancement effect due to convection is negligible. A molecule will lose its memory of macroscopic velocity when it collides with the walls, and thus molecules are hard to diffuse away if they are confined between very close walls. The Rayleigh–Bénard convection in our simulation is in the continuum regime, with the characteristic length of 1000 molecular mean free paths. Under such condition, the effect of solid wall on diffusion is negligible. The diffusion enhancement due to convection is shown to scale as the square root of the Péclet number in the steady convection regime, which is in agreement with previous theoretical and experimental results. In the oscillation convection regime, the diffusion is more strongly enhanced because the molecules can easily advect from one roll to its neighbor due to an oscillation mechanism.
This paper was published as:
Zhang J,Fan J,Fei F. Effects of convection and solid wall on the diffusion in microscale convection flows. PHYSICS OF FLUIDS, 22(12):122005(2010)