Condensation in a square minichannel: application of the VOF method

Abstract

A number of steady-state simulations of condensation of R134a at mass fluxes G=100 kg m-2s-1 and G=800 kg m-2s-1 inside a square cross section minichannel (Dh=1mm) are here proposed and compared against similar simulations in a circular cross section channel with same diameter. The VOF (Volume Of Fluid) method is used to track the vapour-liquid interface and the effects of interfacial shear stress and surface tension are both taken into account. A uniform wall temperature is fixed as boundary condition. At G=100 kg m-2s-1 the liquid film is assumed laminar and the vapour flow is turbulent; turbulence is handled by a low-Reynolds number form of the standard k-w model (Wilcox, 1998), which was modified in order to suppress the turbulent viscosity inside the liquid phase. At G=800 kg m-2s-1 a low Re form of the SST k-w model (Menter, 1994) has been used for turbulence modeling through both the liquid and vapour phases. Numerical simulations are validated against experimental data. The present paper looks at the effect of surface tension. Its influence on the shape of the vapour-liquid interface provides some heat transfer enhancement in non-circular minichannels. In circular minichannels, the overall effect of surface tension is shown to be not significant. On the contrary, the effect of surface tension in square channels provides a large enhancement at low mass flux.