Horizontal convection driven by nonuniform radiative heating in liquids with different surface behavior

https://doi.org/10.1016/j.ijheatmasstransfer.2018.06.036Get rights and content

Highlights

  • Convection in liquids with free surface heated with infrared radiation is studied.

  • Instantaneous temperature and velocity fields are measured using BOS and PIV.

  • Results for water and ethanol are compared to analyze influence of surface behavior.

  • Marangoni convection in ethanol leads to quick dispersion of heat over surface layer.

  • Surface film in water corresponds to no-slip or zero surface divergence condition.

Abstract

Horizontal convection development in liquids with free surface nonuniformly heated by infrared radiation is investigated. Experimental measurements of velocity and temperature fields are performed by Particle Image Velocimetry (PIV) and Background Oriented Schlieren (BOS). Absence of Marangoni convection in distilled water due to presence of surface film leads to significant difference between the results obtained in distilled water and in ethanol. Strong Marangoni flow in ethanol results in quick dispersion of heat over the surface layer and acceleration of convective flow. In distilled water heat propagation along the surface is inhibited, because only weak secondary flow associated with 3D effects is observed at the surface. Experimental results are compared to numerical simulations with different boundary conditions for horizontal velocity at the free surface: no-slip condition for distilled water, representing stagnant surface film, and tangential-stress condition, taking into account Marangoni effect, for ethanol. Good agreement is obtained, demonstrating the importance of prescribing appropriate boundary condition, which corresponds to surface behavior of the considered liquid.

Introduction

When nonuniform heating or cooling is applied along one of the horizontal boundaries of the liquid layer, natural convection circulation develops, which is called horizontal convection [1]. It received some attention, mainly in the context of geophysical flows. Meridional overturning circulation in oceans, which advects warm water from subtropical regions to high latitudes, can be partly explained by horizontal convection, driven by solar radiation heat flux difference near the equator and near the poles. At smaller scale, horizontal convection, associated with nonuniform radiation and evaporation, can provide mixing of water in lakes. Also, it is important for industrial processes, which involve potentially nonuniform heating of liquid from above, in particular for glass production. High-temperature glass melt is processed in glass furnace with several burners heating the melt surface from above in order to provide melting of the floating raw material and to avoid inhomogeneous crystallization. Horizontal convection promotes mixing of the glass melt and thus improves the homogeneity of the produced glass sheets.

In most papers devoted to horizontal convection nonuniform heating was applied along the base of liquid tank. In these studies [2], [3], [4] numerical simulations of the flow were performed for various ranges of Rayleigh number and tank aspect ratio and scaling of Nusselt number, boundary layer thickness and circulation intensity was obtained. Boundary layer stability was also analyzed [5]. Experimental investigations, including temperature profiles measurements using thermistors, qualitative schlieren visualizations, dye visualizations of the flow and measurements of heat flux distribution at the tank bottom, were conducted by Mullarney et al. [2] and Sanmiguel Vila et al. [6]. Horizontal convection in liquid nonuniformly heated from above was mostly studied numerically for large Prandtl numbers typical for glass production problems [7], [8], [9]. Note that in all these studies nonuniform temperature profile was prescribed at the top boundary and no-slip conditions were imposed at all the boundaries of liquid volume. Simulations for no-slip and stress-free conditions at the top boundary were performed by Chiu-Webster et al. [10], and it was shown that in the limit of infinite Prandtl number different boundary conditions yield the same scaling, though Nusselt number obtained in stress-free case is about 60% larger. Both configurations with heating applied along the top and bottom rigid boundaries were considered in experimental study by Wang and Huang [11]. Velocity fields for horizontal convection steady state were obtained using PIV.

Studies of horizontal convection in liquids with free surface heated by radiation are less numerous, though this problem formulation is more relevant both to ocean circulation and to glass production. Steady-state flow patterns and temperature fields were visualized by Kurosaki et al. [12] using holographic interferometry and tracer particles. Silicon oil was used as working liquid, and surface flow associated with Marangoni convection was observed. In contrast, PIV measurements for horizontal convection in water, carried out by Wȧhlin et al. [13], showed that surface velocity is close to zero. The authors of [13] explained this by the influence of surface film, as experiments were performed in large water tanks for several days and it was not possible to avoid presence of surfactants in water. Though they did not perform measurements in other liquids, they predicted that ’the flow may be very sensitive to the boundary conditions for the horizontal velocity at the free surface’. Horizontal convection driven by radiative heating of the free surface was also investigated by Shmyrov et al. [14] for bidistilled water with addition of insoluble surfactant. Surfactant transport by the surface flow resulted in formation of two flow regions, separated by the surface stagnation point. Marangoni convection was observed in region with clean surface, whereas in region, covered with surfactant, surface layer was stagnant and convective vortex was located below. Note that water surface was cleaned by aspiration prior to addition of surfactant, which was required for Marangoni convection to be observed.

In the present study we investigate the effect of different boundary conditions for horizontal velocity at the free surface upon the flow structure and heat transfer during the development of horizontal convection in liquid heated by infrared radiation from above. Instant velocity and temperature fields are measured using PIV and BOS for convection in distilled water and ethanol, which exhibit different surface behavior. Earlier study of convective plume generated by a horizontal heated wire and impacting the free surface [15] in these two liquids showed that velocity of heat wave propagation along the surface is several times larger in ethanol due to presence of Marangoni flow, absent in distilled water. Using Marangoni flow, produced by local heating of the liquid surface with laser radiation, for surface cleaning and microparticles manipulation was proposed in [16], [17], [18], which implies that presence or absence of Marangoni flow can significantly alter horizontal convection flow, at least in shallow tanks. As in [15], we directly compare experimental data to results of numerical simulations, performed with different boundary conditions for horizontal velocity, representing surface behavior of water and ethanol. The paper is organized as follows. Section 2 describes the employed experimental techniques and equipment. In Section 3 mathematical problem is formulated and numerical method is outlined. The results of experimental measurements and numerical simulations are presented and discussed in Section 4. Conclusions are drawn in Section 5.

Section snippets

Experimental techniques and equipment

Experiments are performed in rectangular tank with internal dimensions 141×72×80 mm, made of 3-mm window glass (Fig. 1). Radiation from a ceramic heater 240×60 mm passes through rectangular aperture with dimensions 72×40 mm, cut in a special screen, and heats the liquid surface. In order to minimize heat transfer through the screen, the screen is made from a 40-mm-high cardboard box with all the walls covered with aluminium foil, reflecting radiation. Thus, nearly uniform irradiation of

Numerical modeling

Liquid flow in vertical plane XY is simulated by solving unsteady 2D Navier-Stokes equations in low-Mach approximationρt+(ρVi)Xi=0ρVit+VjViXj=-pXi-23XiηVjXj+XjηViXj+VjXi+ρgiρcpTt+ViTXi=XiλTXi+-23ηVkXkδij+ηViXj+VjXiViXjwith empirical correlations for temperature dependencies of liquid density, viscosity, thermal conductivity and specific heat given in Appendix A.

Bottom and side walls of the tank are assumed adiabatic. Third-kind boundary condition is imposed

Results and discussion

Heating the liquid surface results in pressure imbalance and free surface deformation. A convex crest is formed in irradiated region, which corresponds to positive pressure perturbation. Warm liquid below the surface is forced to flow towards the side walls and two counter-rotating vortices are formed below the edges of irradiated region (Fig. 5). This is buoyancy-driven horizontal convection, which is observed both in distilled water and ethanol. However, in liquids, which exhibit Marangoni

Conclusions

In this study, liquid temperature and velocity fields during the development of horizontal convection driven by nonuniform heating of free surface have been measured using BOS and PIV. Experiments have been performed in distilled water and ethanol, which are known to exhibit different surface behavior. In both liquids ’partial penetration’ is observed – the circulation flow is mostly restricted to a shallow surface layer. However, comparison of the results, obtained in different liquids, shows

Conflict of interest

The authors declare that there is no conflict of interest.

References (37)

  • G.O. Hughes et al.

    Horizontal convection

    Annu. Rev. Fluid Mech.

    (2008)
  • J.C. Mullarney et al.

    Convection driven by differential heating at a horizontal boundary

    J. Fluid Mech.

    (2004)
  • O. Shishkina

    Mean flow structure in horizontal convection

    J. Fluid Mech.

    (2017)
  • I.E. Sarris et al.

    Natural convection in a 2D enclosure with sinusoidal upper wall temperature

    Numer. Heat Transf. A

    (2002)
  • S. Chiu-Webster et al.

    Very viscous horizontal convection

    J. Fluid Mech.

    (2008)
  • W. Wang et al.

    An experimental study on thermal circulation driven by horizontal differential heating

    J. Fluid Mech.

    (2005)
  • A.K. Wȧhlin et al.

    Horizontal convection in water heated by infrared radiation and cooled by evaporation: scaling analysis and experimental results

    Tellus

    (2010)
  • A. Shmyrov et al.

    On the extent of surface stagnation produced jointly by insoluble surfactant and thermocapillary flow

    Adv. Colloid Interface Sci.

    (2017)
  • Cited by (21)

    • Transitional fluctuations of thermal boundary layer in horizontal convection

      2023, International Journal of Heat and Mass Transfer
      Citation Excerpt :

      HC also occurs in engineering process such as glass melting and furnaces, which results from temperature gradient on the top surface and is vital for the quality control in manufacturing [6–8]. A great number of investigations for various boundary conditions have been done in theories [1,9–14], experiments [15–22] and simulations [23–32], which are dedicated to a better understanding of transport mechanism and its connection with the mean field structures. The remaining part of the paper is organised as follows.

    • Obtaining velocity and pressure distributions in natural convection flows using experimental temperature fields

      2022, Applied Thermal Engineering
      Citation Excerpt :

      Uniform cartesian grid, formed by the centers of interrogation windows in BOS image processing, is employed in hybrid simulations to avoid spatial interpolation of experimental data. More detailed descriptions of the numerical simulations for the considered flows, including the specified boundary conditions, code verification and grid convergence studies, can be found in [36,37]. In this study new data assimilation technique is proposed for instantaneous velocity and pressure fields reconstruction from a sequence of experimental temperature fields.

    • Effect of radiative heat transfer on thermal-solutal Marangoni convection in a shallow rectangular cavity with mutually perpendicular temperature and concentration gradients

      2022, International Journal of Heat and Mass Transfer
      Citation Excerpt :

      In addition, it was verified that the spoke patterns were resulting from the thermal Marangoni instability due to the reverse temperature gradient in the thin thermal boundary layer. Vinnichenko et al. [15] examined the effect of local radiative heating on the characteristics of horizontal convection in a rectangular tank filled with ethanol. They found that, compared with buoyancy-driven convection, the thermal Marangoni convection plays a predominant role in the thin thermal boundary layer and enhances heat transport to the whole free surface.

    • A method to measure vapor concentration of droplet evaporation based on background oriented Schlieren

      2021, International Journal of Heat and Mass Transfer
      Citation Excerpt :

      As a non-invasive optical method, the Schlieren technology can convert the density change of the fluid that is not visible to the human eye into light intensity information, and has led to great progress in qualitative measurements in recent years [25–27]. To obtain quantitative information, several derivative forms of the Schlieren method have been developed including the rainbow Schlieren method [28,29] and the background oriented Schlieren method (BOS) [30–33]. The BOS method combines the basic principles of the traditional Schlieren technology and the particle tracking technology, which can be used to measure the density gradient in two directions.

    View all citing articles on Scopus
    View full text