On the measuring of film thickness profiles and local heat transfer coefficients in falling films

https://doi.org/10.1016/j.expthermflusci.2018.07.028Get rights and content

Highlights

  • A novel procedure for measuring falling film thickness profile.

  • A novel, flexible procedure for measuring local heat transfer coefficients.

  • Film thickness resolution ±1.2 · 10−5 m at 500 Hz and accuracy of at least ±0.1 mm.

  • Accuracy of measuring heat transfer coefficients in the range 4–20%.

  • Great relevance of simultaneously studying falling film thickness and heat transfer.

Abstract

In this work we develop and report on two novel measurement procedures for elucidating important features of vertical falling films. The first one utilizes a laser triangulation scanner to continuously measure a film thickness profile along a 0.10 m long vertical line at a frequency of 500 Hz and a resolution of ±1.2 · 10−5 m. The second procedure is a flexible method to measure local time-averaged wall and bulk temperatures, making it possible to calculate local heat transfer coefficients at desired locations. The procedures have been implemented and evaluated in a falling film facility, under a broad range of operating conditions and covering different types of fluids and wetting rates. We have shown that the used scanner performs well for several reflecting fluids, that it resolves the flow pattern in high resolution and that it measures the film thickness with an accuracy of at least ±0.1 mm. Furthermore, we show that the measurements of the local heat transfer coefficient have an accuracy of ±0.01 K, allowing estimation of the heat transfer coefficient with an accuracy of 4–20%, depending on the operating conditions. Both procedures are accurate, robust and simple to use and their combination gives valuable insights into important features of falling films, especially how the hydrodynamics and heat transfer affect each other.

Introduction

Falling liquid film is a special kind of flow where a liquid film with a distinct interface flows down an inclined or vertical wall, usually in the presence of a gas layer. High heat transfer at relatively low mass flow rates and small temperature differences, together with short product residence, makes this technique a popular choice for a wide range of engineering and technological applications. For instance, falling films are used for heat-sensitive fluids such as fruit juice, sugar and dairy products, where short residence time and temperature control during the heat transfer process are essential. The heat transfer in falling film units is often predominantly influenced by the hydrodynamics [1], hence it is important to study the two simultaneously when evaluating the performance of these units.

The flow in falling film units is often described using dimensionless numbers. The most important ones are the Reynolds number (Re), the Kapitza number (Ka), the Nusselt number (Nu), and the Prandtl number (Pr), defined as follows:Re=ΓμKa=σ·ρ1/3μ4/3·g1/3Nuhkν2g1/3Prcpμkwhere Γ is the specific mass flow rate per unit circumference, g is the gravitational acceleration, h is the heat transfer coefficient, and the following properties of the working fluid are: μ is the dynamic viscosity, ρ is the density, σ is the surface tension, k is the thermal conductivity, ν is the kinematic viscosity and cp is the specific heat capacity.

The hydrodynamics of falling films have been studied by numerous researchers, e.g. [2], [3], [4], [5], [6], [7], [8], [9], [10], [11], [12], [13], [14], [15], [16] and are usually characterized by different flow regimes [17], [18]. The transitions between different flow regimes are often described as a function of the Reynolds and Kapitza numbers. Al-Sibai [19] determined conditions for the transitions under non-evaporative conditions, where the flow was mapped into five different zones, Laminar (L), Sinus-shaped waves (S), Wavy-laminar (WL), Transition region (TR) and Turbulent (T). These regimes describe how the flow is transitioning from a smooth laminar film into a fully turbulent one, see Fig. 1. Each flow regime impacts the heat transfer differently depending on the film thickness and the degree of turbulence.

The hydrodynamics of a falling film are usually further investigated by looking at the changes of the film thickness, since the emerging waves govern much of the underlying flow structure that affects the heat transfer [20]. Numerous studies have attempted to determine the local film thickness of the liquid film in different flow regimes, and a large variety of techniques exist. Charogiannis et al. [21] recently made a thorough review of different methods.

The available techniques can usually be classified as intrusive or non-intrusive. There are simple intrusive methods, such as drainage or hold-up [22], where the flow is stopped and the liquid on the surface is collected and weighted, capable of giving an estimate of the average film thickness. Hot wire anemometry [23] or needle-contact-probes [24] have also been used to measure the thickness in a single point.

As for non-intrusive methods, one can use conductance probe [25] or capacitance probe [26], [27], [28] to estimate the amount of liquid. There are also several optical methods, for example using shadow photographs, where the film thickness is estimated from images [8]. There is also laser focus displacement [29] or chromatic confocal imaging [30], where monochromatic light is continuously emitted onto the interface and the intensity of the reflective light is measured. Furthermore, there are studies using laser-induced fluorescence [31], [32], [33], [34], [35], [36], [37] or planar laser-induced fluorescence [38], [21], [39] where the target media is illuminated, usually by a laser, and the intensity of the fluorescence is captured by a detector and can be related to various medium properties. All techniques clearly have their advantages and disadvantages, but in general, they are either insufficiently accurate, have low spatial resolution, disturb the flow, need additives or require a fluorescent or transparent fluid.

Correlations that estimate the film thickness are based on the results from different measurements techniques. The film thickness for flat laminar falling films without surface evaporation can be obtained from the Nusselt film theory as [40]:δNusselt=3ν2Reg1/3

Lukach et al. proposed the following correlation for the mean film thickness under wavy-laminar conditions [41]:δLukach=1.34ν2g1/3·Re0.368and Brötz [42] proposed the following correlation for turbulent conditions (Re > 600):δBrötz=0.172ν2g1/3·Re2/3

Local heat transfer measurements in a subcooled system are typically performed by measuring the temperature difference between the wall and the bulk [43] or the interface [44], [45], [46], [47]. Usually the most difficult property to obtain is the wall temperature, which, for example, can be measured by installing a thermocouple (e.g. by soldering) inside the tube wall [48], [49]. However, this approach requires considerable effort, the number of measurement points is limited, and the position of the thermocouples must be predetermined before the start of the experimental campaign.

The heat transfer in falling films is usually described as a function of the Reynolds and Prandtl numbers. The degree of influence of each dimensionless number depends on the flow regime. For Laminar flow regime (see Fig. 1), the heat transfer takes place by pure conduction, it is rate-determined by the film thickness and it thus decreases with the Reynolds number. In Turbulent regime, however, the Nusselt number increases with the Reynolds number [50] and is also dependent on the Prandtl number. For an intermediate regime (Transition regime), a superposition is often used. Schnabel and Schlünder [51] developed a correlation for non-evaporative conditions under a constant heat flux:NuHeating=max1.43·Re-1/3(Laminar)0.042·Re0.2·Pr0.344(Transtion)0.014·Re0.4·Pr0.344(Turbulent)

The above correlation can be used to estimate the magnitude of the heat transfer coefficient but it does not give any information on how that coefficient varies in the streamwise direction or how it relates to the film thickness. Therefore, the objective of this work is to investigate two novel experimental procedures for characterizing falling films. The procedures will help us to gain a deeper insight into the hydrodynamics and the heat transfer of falling films, but also into their interaction. A quantitative method is assessed for obtaining wave characteristics and film thicknesses over a section of the surface. In the same time, we evaluate a method to measure the local wall and bulk temperatures in order to obtain the local heat transfer coefficients. Both procedures are investigated in an experimental falling film setup designed to study falling film characteristics.

Section snippets

Experimental setup and materials

The proposed measuring procedures are evaluated in an atmospheric falling film facility, see Fig. 2(Left). A copper tube, with the length of L = 0.8 m, outer diameter do = 60 mm and 5 mm thickness, is connected to an overflow distributor, illustrated in Fig. 2 (Right) to achieve an evenly distributed liquid film on the outer side of the tube. A displacement pump recirculates the liquid over the tube. An electrical heater (SAS Steatite cartridge) is mounted inside the tube and heats the

Laser triangulation scanner

The procedure we propose in this work for measuring the film thickness uses an optical triangulation scanner (light intersection method) of the model scanCONTROL 2950-100 from Micro-Epsilon [54]. The scanner uses a 20 mW power source to produce a laser beam of the wave length 658 nm that is enlarged through special lenses to form a static laser line on the target surface. A sensor matrix detects the diffused light and measures the reflected angle at 1280 points along the projected line. The

Laser triangulation scanner

Fig. 11 presents a comparison between the film thickness measurement (left) and the high-speed image (right) at two separate times (top and bottom) for the dairy product fluid at conditions Ka = 97 and the wetting rate of 0.48 kg/(m s). In the figure, it can be seen that the laser triangulation scanner is indeed able to detect the fundamental features of the liquid film flow and that it captures the wave dynamics in the same way as does the high speed imaging. The laser is capable of detecting

Conclusion

In this paper, we develop and report on two novel procedures for studying a number of fundamental features of vertical falling films. The film thickness measurement procedure is proven to capture the wave dynamics accurately along an approximately 0.1 m long vertical line, with a high spatial (±1.2 · 10−5 m) and temporal resolution (500 Hz). No additives to the liquid are required as long as the fluid investigated is reflective. There is also an advantage that the measurement equipment is

Conflict of interest

The authors declared that there is no conflict of interest.

Acknowledgment

This work was co-funded by the Swedish Energy Agency (P40550-1), Valmet AB and Tetra Pak Processing systems.

References (56)

  • W. Ambrosini et al.

    Statistical characteristics of a water film falling down a flat plate at different inclinations and temperatures

    Int. J. Multiph. Flow

    (2002)
  • T. Takamasa et al.

    Measuring interfacial waves on film flowing down a vertical plate wall in the entry region using laser focus displacement meters

    Int. J. Heat Mass Transf.

    (2000)
  • A.V. Cherdantsev et al.

    Study of gas-sheared liquid film in horizontal rectangular duct using high-speed LIF technique: three-dimensional wavy structure and its relation to liquid entrainment

    Int. J. Multiph. Flow

    (2014)
  • R. Mathie et al.

    Heat transfer augmentation in unsteady conjugate thermal systems – Part II: applications

    Int. J. Heat Mass Transf.

    (2013)
  • I. Zadrazil et al.

    An experimental characterization of downwards gas–liquid annular flow by laser-induced fluorescence: flow regimes and film statistics

    Int. J. Multiph. Flow

    (2014)
  • I. Zadrazil et al.

    An experimental characterization of liquid films in downwards co-current gas–liquid annular flow by particle image and tracking velocimetry

    Int. J. Multiph. Flow

    (2014)
  • F. Al-Sibai et al.

    Local and instantaneous distribution of heat transfer rates through wavy films

    Int. J. Therm. Sci.

    (2002)
  • C.N. Markides et al.

    An experimental study of spatiotemporally resolved heat transfer in thin liquid-film flows falling over an inclined heated foil

    Int. J. Heat Mass Transf.

    (2016)
  • A. Charogiannis et al.

    Thermographic particle velocimetry (TPV) for simultaneous interfacial temperature and velocity measurements

    Int. J. Heat Mass Transf.

    (2016)
  • M. Johansson et al.

    Heat transfer in evaporating black liquor falling film

    Int. J. Heat Mass Transf.

    (2009)
  • M. Madoumier et al.

    Modelling the properties of liquid foods for use of process flowsheeting simulators: application to milk concentration

    J. Food Eng.

    (2015)
  • A. Åkesjö et al.

    Experimental and numerical study of heat transfer in a large-scale vertical falling film pilot unit

    Int. J. Heat Mass Transf.

    (2018)
  • A.E. Dukler

    The role of waves in two-phase flow: some new understanding

    Chem. Eng. Educ.

    (1976)
  • W. Nußelt

    Die Abhängigkeit der Wärmeübergangszahl von der Rohrlänge

    Z. Ver. Dtsch. Ing.

    (1910)
  • W. Nußelt

    Die oberflachenkondensation des wasserdampfes

    Z. Ver. Dtsch. Ing.

    (1916)
  • P.L. Kapitza et al.

    Wave flow of thin layers of a viscous fluid

  • S.V. Alekseenko et al.

    Wave Flow of Liquid Films

    (1994)
  • A. Charogiannis et al.

    Detailed hydrodynamic characterization of harmonically excited falling-film flows: a combined experimental and computational study

    Phys. Rev. Fluids

    (2017)
  • Cited by (16)

    • Heat transfer mechanism and characteristics of fluid film on multi-faceted surface under constant cross-flow conditions

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

      Variation in Pr is achieved by varying the feedwater inlet temperature [32]. In addition, laminar fluid film flow is not a function of Pr [33]. Hence, the key observations are made by varying the liquid load, which allows us to analyze the surface effect on the heat transfer mechanism at the gas-liquid interface based on the hydrodynamics and kinematics of the fluid film.

    • Hydrodynamics numerical simulation of a vertical falling film evaporator for ionic liquid systems

      2021, Chemical Engineering Science
      Citation Excerpt :

      Numerous studies have attempted to determine the local film thickness of the liquid film in different flow regimes. Hence, a series of models for the prediction of film thickness were proposed (Akesjo et al., 2018; Min and Park, 2011). For example, Gao et al. exhibited a relative accurate prediction of the film thickness in laminar flow (Gao et al., 2018), while Lukach et al. proposed a correlation for the mean film thickness under wavy-laminar conditions (Lukach et al., 1972).

    • Evaporation heat transfer characteristics of falling film in small diameter fabricated tubes of absorption refrigeration system: An experimental investigation

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

      The heat transfer coefficient was found to be independent of salt concentration, flow rate, and tube bundle geometries. Akesjo et al. [12] investigated the heat transfer and hydrodynamic characteristics of falling film and concluded that the heat transfer and hydrodynamics affect each other at different rates; hence, it is necessary to measure both of them individually to fully understand the mechanism of the falling liquid film. The purpose of this research is to investigate the falling-film evaporation heat transfer characteristics of newly developed enhanced tubes of an absorption refrigeration system.

    View all citing articles on Scopus
    View full text