Thermal protection of a vertical plate using ethylene glycol film cooling flowing down on a vertical plate

This paper concerns a numerical investigation of the evaporative cooling of liquid film falling along a vertical channel. The first plate is insulated and wetted by an ethylene glycol liquid film while the second one (y = d) is dry and isothermal. A marching procedure is employed for solution of the equation of mass momentum, energy and concentration in the flow. Numerical results for air-ethylene glycol system are presented. The effects of ambient gas temperature and inlet vapor mass fraction of ethylene glycol on the heat and mass transfer and on the ethylene glycol liquid film evaporation are investigated.


Introduction
The liquid film evaporation is found in many industrial applications such as evaporator, condensing devices, air-conditioning, refrigeration technology, absorber/ generator heat exchange, nuclear technology and separation processes, heat pumps and drying technology. Cherif et al. [1] presented an experimental study of mixed convection heat and mass transfer in a vertical channel with film evaporation. The experimental results show that evaporation takes place on the majority of the surface of the two walls and, in some cases, evaporative cooling occurs especially for small heating flux and large air velocities. Ben Jabrallah et al. [2] studied numerically and experimentally the convective heat and mass transfer with the evaporation of a falling film in a cavity. They showed that the obtained results allow us to describe the thermodynamic state of the heated film by means of the liquid temperature and evaporation flow rate. Gorjaei et al. [3] presented an analysis of the effect of introduction of Al2O3 nanoparticles in the heat transfer inside a three-dimensional annulus. The convective heat transfer coefficient of the nanofluid can be improved with a lower volume concentration of CuO nanoparticles, as showed by Lazarus et al. [4]. The effects of nanoparticles dispersion on droplet evaporation were presented by Chen et al. [5]. They demonstrated an enhancement of droplet evaporation with nanoparticles adding. A study on heat transfer caused by free convection of nanofluid flowing on a vertical plate has been effected by Siddiqa et al. [6]. They demonstrated that the heat exchange was improved by the nanoparticles' dispersion. The free convective heat transfer within a porous wavy cavity with a nanofluid was numerically examined by Sheremet et al. [7]. They demonstrated that the local heat source affected the nanofluid flow and the heat transfer. Water droplets evaporation containing a low concentration of nanoparticles were examined by Askounis et al. [8]. They showed that the rate of droplet evaporation is unaffected by the dispersion of low concentration nanoparticles. Perrin et al. [9] presented a comparison between the theoretical and experimental results of nanofluid liquid drop evaporation. A numerical study by mixed convection of the falling liquid film evaporation was presented by Yan [10]. Experimental analysis of the liquid film evaporation was conducted by Huang et al. [11]. They showed that an increase in the temperature and air flow rate accelerated evaporation. Wei et al. [12] simulated the flow and heat transfer of CuO-water nanofluid in a tube using the single-and two-phase (mixture) models. Nasr and Alzahrani [13] presented a numerical study of liquid nanofilms' evaporation inside a heat exchanger by mixed convection. Nasr et al. [14] studied the evaporation of binary liquid film by mixed convection inside heated vertical channel. Nasr and AL-Ghamdi [15] studied the liquid nanofilms' condensation inside a heat exchanger by mixed convection. Nazir et al. [16] presented a thermal and mass species transportation in tri-hybridized Sisko martial with heat source over vertical heated cylinder. Nazir et al. [17] presented a significant Production of Thermal Energy in Partially Ionized Hyperbolic Tangent Material Based on Ternary Hybrid Nanomaterials. They showed that the thermal energy expression is derived by the contribution of Joule heat and viscous dissipation. They presented a model of flow equations by using the concept of boundary layer theory, which occurs in the form of a coupled system of partial differential equations (PDEs). Muhammad Sohail et al. [18] presented a finite element analysis for ternary hybrid nanoparticles on thermal enhancement in pseudo-plastic liquid through porous stretching sheet. Muhammad Sohail et al. [19] presented a study of triple-mass diffusion species and energy transfer in Carreau-Yasuda material influenced by activation energy and heat source. Nazir et al. [20] presented finite element analysis for thermal enhancement in power law hybrid nanofluid.
As regards the prior research, the numerical study of ethylene glycol liquid film evaporation is no considered. The main purpose of the present work is to analyze the effect of ambient gas temperature and inlet vapor mass fraction of ethylene glycol on the heat and mass transfer and on the ethylene glycol film evaporation.

Analysis
The present work deals with a numerical analysis of evaporation of pure liquid film by mixed convection flowing along one of the channel vertical plates (Fig. 1). The studied channel is made up of two vertical and parallel plates. The first plate is subjected to a uniform heat flux q 1 and wetted by a pure ethylene glycol liquid film while the second one (y = d) is dry and isothermal. The pure liquid film flowing down with an inlet temperature T 0L and an inlet mass flow rate m L0 .The air enters the channel with a temperature T 0 , a water and ethylene-glycol vapour concentrations c 01 and c 02 and velocity u 0 .

For The liquid phase
The Continuity equation, x-momentum equation, Energy equation, diffusion equation are respectively: Continuity equation x-momentum equation

For the gaseous phase
Continuity equation x-momentum equation The heat balance at the interface implies The continuities of shear stress and local evaporated mass flux of species i give: The latent heat flux of mixture is given The sensible heat flux is given by: The total cumulated evaporation rate of mixture at the interface is given by:

Results and discussions
All the results of this study have been obtained for the case of c 02 = 0, T 0 = 293.15 K, T 0L = 293.15 K,T w = 293.15 K,m 0L = 0.01 5 k g/m .s,u 0 = 1 m/,the geometrical ratio is d/H = 0.015 and the imposed wall heat flux is q 1 = 0W/m 2 .  The evolution of the sensible heat flux of mixture along the channel for several values of inlet temperature of the gas is illustrated in Fig. 3a. This figure shows that the sensible heat flux exchange increases with the inlet temperature of the gas. This result has been explained by the fact that the interfacial temperature increase with the inlet temperature and consequently the sensible heat flux exchange increases. It is observed from Fig. 3b that the latent heat flux exchange increases with the inlet temperature of the gas. Apparently, the values of the latent heat flux are all positive, and this indicates that the direction of the latent heat flux is from the interface to the hot gas stream. It is noted that the latent heat flux increases in the flow direction.
It is apparent that from Fig. 4 that the local mass evaporation rate increases with the inlet temperature of the gas. The axial distributions of non-dimensional accumulated mass evaporation rate are shown in Fig. 5. It is clear that an increase in T 0 causes a larger Mr. This can be explained by the fact that when we increase the inlet temperature of the gas, the latent heat flux is from the interface increases and consequently the accumulated mass evaporation rate increases. Figure 6 shows that an increase in the inlet gas mass fraction of ethylene glycol induces an increase of the temperature at the interface liquid-gas. It is shown from Fig. 7 that an increase in the inlet gas mass fraction of ethylene glycol inhibits the ethylene glycol evaporation causing an increase in the temperature at the interface liquid-gas. This result has been confirmed in the Fig. 6.   (Fig. 1). The studied channel is made up of two vertical and parallel plates. The first plate is insulated and wetted by a pure liquid film (ethylene glycol) while the second one (y = d) is dry and isothermal. The liquid film flowing down with an inlet temperature T 0L and an inlet mass flow rate m L0 .The air enters the channel with a temperature T 0 , a water and ethylene-glycol vapour concentrations c 0 and velocity u 0 . It is shown that an augmentation in the gas temperature induces an increase in the temperature and ethylene glycol mass fraction at the interface and a decrease in the latent and sensible heat fluxes. The increase in the inlet gas temperature enhances the ethylene glycol film evaporation.
Acknowledgements Not Applicable.
Author contributions AN constructed the model; wrote the introduction; reviewed the paper, solved the equations system; figured out the results; reviewed the draft, wrote the discussions; wrote the conclusion; revised the English.

Funding Not applicable.
Availability of data and materials Not applicable.

Conflict of interest Not applicable.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons. org/licenses/by/4.0/.

Conclusions
We presented a numerical study of evaporation of ethylene glycol liquid film by mixed convection flowing along