Effects of Interfacial Tension Gradients
From ThermalFluidsPedia
Since surface tension depends on temperature, a permanent nonuniformity of temperature or concentration (for a multicomponent system) at a liquidvapor interface causes a surface tension gradient. The interfacial area with small surface tension expands at the expense of an area with greater surface tension, which in turn establishes a steady flow pattern in the liquid; this flow caused by the surface tension gradient along the liquidvapor (gas) interface is referred to as the Marangoni effect. The surface tension of a multicomponent liquid that is in equilibrium with the vapor is a function of temperature and composition of the mixture, i.e.,

where x_{i} is the molar fraction of the i_{th} component in the liquid phase and N is the total number of components in the liquid phase. The change of surface tension can be caused by either change of temperature or composition, i.e.,

Recall from thermodynamics that as the critical temperature for a given fluid is approached, the properties of the liquid and vapor phases of the fluid become identical, i.e., σ vanishes. It therefore follows that the surface tension decreases with temperature, i.e., []. For a pure substance, surface tension is a function of temperature only. The curvefit equations for surface tension are almost linear for most fluids and can have the form
σ = C_{0} − C_{1}T 
where C0 and C1 are empirical constants that vary between substances. For water, and . The unit of the temperature T in eq. (3) is °C.
Since surface tension varies with temperature, the surface tension will not be uniform if the temperature, along the liquidvapor interface, is nonuniform. Consequently, the liquid in the lower surface tension region near the interface will be pulled toward the region with higher surface tension. Furthermore, because the surface tension usually decreases with increasing temperature, the flow driven by surface tension moves away from the interface with high temperature and towards the interface with low temperature. As noted before, this motion of a liquid caused by a surface tension gradient at the interface is referred to as the Marangoni effect. The most wellknown example of surfacetensiondriven flow is Bernard cellular flow, which occurs in a thin horizontal liquid layer heated from below. Figure 1 illustrates steady cellular flow driven by the Marangoni effect. Once a steady Marangoni flow is established, the liquid velocity is upward at point A and downward at point B. Since the liquid at surface point A comes directly from the hot surface, the temperature at point A is higher than that at point B. Consequently, the surface tension at point A is smaller than that at point B, and the fluid at point A is pulled toward point B. Although a temperature gradient exists in the vertical direction, the actual driving force is the surface tension gradient in the horizontal direction. The liquid surface loses heat to the gas phase at temperature, TG, through convection with a heat transfer coefficient of h_{δ}.
An important boundary condition at the interface explains the liquid flow field caused by surface tension variation along the interface:

Bernard cellular flow resulting from the Marangoni effect can occur only if certain conditions are met. The conditions for the onset of cellular motion can be predicted using linear stability analysis as presented by Carey (1992). The local temperature is assumed to equal the sum of the basic temperature – given by the initial linear profile – and a small sinusoidal fluctuation T' that represents a Fourier component of random disturbances, i.e.,

where ζ = (T_{w} − T_{δ}) / δ. Since all base flow velocities are zero, the velocities are


Substituting eqs. (6) – (8) into the continuity, momentum, and energy equations, and subtracting the corresponding base flow equations, the resulting equations are solved assuming the following forms for T' and v':


where α = 2π / λ is the wave number and β can be a complex number with its real part representing the amplification factor, and its imaginary part representing the temporal frequency. The resulting stationary wave solutions (for β = 0), which are not amplified or dampened with time, are therefore interpreted as stable perturbations. The conditions for which these solutions are obtained are assumed to correspond to the onset of instability, which leads to cellular convection. The case of marginal stability corresponds to

where M_{a} is the Marangoni number:

and is the thermal diffusivity of the liquid.
The other dimensionless numbers in eq. (11) are the wave number , the B_{iot} number B_{i}, the Bond number B_{o}, and the Crispation number C_{r}. Their definitions are




The liquid film is stable if the Marangoni number is below that predicted by eq. (11). However, the liquid film is not stable if its Marangoni number is above that obtained by eq. (11). The marginal stability predicted by eq. (11) for some typical combinations of parameters is illustrated in Fig. 2. Since the Fourier components of all wavelengths can be contained in a random disturbance, the system becomes unstable when it is unstable at any wavelength. For a system with Cr < 10^{ − 4} and , Fig 2 indicates that the critical Marangoni number is about 80 and the associated dimensionless wave number
The Marangoni effect can have an important influence on heat and mass transfer processes, including the evaporation of a falling film, as well as causing vapor bubbles in a liquid with a temperature gradient to move toward the high temperature region during boiling.
References
Carey, V.P., and Wemhoff, A. P., 2005, “Disjoining Pressure Effects in UltraThin Liquid Films in Micropassages – Comparison of Thermodynamic Theory with Predictions of Molecular Dynamics Simulations,” IMECE200580234, Proceedings of 2005 ASME International Mechanical Engineering Congress and Exposition, Orlando, FL.
Faghri, A., and Zhang, Y., 2006, Transport Phenomena in Multiphase Systems, Elsevier, Burlington, MA.