Abstract
The heat and mass transfer process of the laminar falling liquid film along a vertical heated plate with constant heat flux is studied. Through thermal equilibrium analysis, the energy conservation equation is established and solved by the parametric variation method. The analytical solution considering evaporation is in good agreement with the experimental results. Based on the model, the temperature distribution of liquid film with or without evaporation is obtained. It is demonstrated that the thickness and average temperature of liquid film show decreasing and increasing trends linearly in the flow direction. Along with the increase in the evaporation rate, the average temperature decreases and the share of evaporative heat dissipation in the total heat flux increases, which demonstrates that evaporation is an important physical factor for heat dissipation when the falling liquid film flows on a hot plate.
Introduction
Vertical falling film heat dissipation refers to the process in which the film flows from top to bottom under the action of gravity, viscous shear, and wall adhesion, and the heat is transferred from the wall to the film and accompanied by evaporation. It is widely used as an enhanced cooling method in the fields of chemical process, nuclear reactor safety, and fire protection engineering (Yang and Yuen, 2018; Turkyilmazoglu, 2021; Choi and Cho, 2022). In particular, the mixed ferromagnetic nanofluid film and microchannel fluid film show a good application prospect for heat transfer (Tlili et al., 2020; Tlili et al., 2022). For large military platforms, the falling film can be used to reduce the infrared contrast between the warm vertical wall surface and the environment (Peng et al., 2012a; Peng et al., 2012b). In the present work, the authors mainly discuss the related heat and mass transfer problems in the application of liquid film to control or balance the temperature of the large vertical wall to reach or approach that of the background environment.
In recent years, many researchers have studied the flow and heat transfer process of the falling film under different conditions by theoretical, numerical, and experimental methods (Nasr et al., 2017; Malika et al., 2017; Turkyilmazoglu, 2021; Yue et al., 2021). Because of the complexity of the problem, numerical and experimental studies tend to prevail. However, analytical solutions play an irreplaceable role in universality. The following findings are of great significance. Nusselt (1916) made basic research on this topic and gave a series of equations describing the velocity and thickness distribution of the film. Saouli and Aïboud-Saouli (2004) established a theoretical model of laminar flow falling film along an inclined plate by ignoring the inlet effect and gained the velocity and temperature distributions by separating variables. Peng et al. (2012a) and Peng et al. (2012b) derived the solution of liquid film energy equation with constant wall temperature and heat flux using the integral analysis method. Meyer (2014) and Meyer (2015) introduced Laplace transform to solve the problem with constant film velocity for an isothermal or adiabatic wall boundary condition. Gu et al. (2018) analytically treated the energy and absorbate equations of the liquid film by means of integral transforms and obtained the exact series solutions for temperature and mass fraction fields under the assumptions of both the insulated and isothermal wall boundary conditions.
Another problem is the effect of evaporation on falling film heat transfer. According to the best knowledge of the authors, some scholars in this field have not considered the evaporation effect in the study of the above physical processes, which may be a factor that cannot be ignored. In a few relevant studies, the evaporation of binary liquid film was numerically analyzed by Nasr et al. (2011). Assad and Lampnen (2002), Leu et al. (2006), and Leu et al. (2009) established a mathematical model for the evaporation process on a constant temperature vertical plate between liquid film and wall surface. Cherif et al. (2010) conducted experimental studies on natural and forced convection, studying the evaporation of thin liquid films falling on the inner surface of plates in vertical channels. In addition, some studies on film flow have been based on constant film thickness (He et al., 1998; Song et al., 1999; Wang et al., 2000; Cherif et al., 2010; Nasr et al., 2011), which simplifies the problem but deviates from the physical facts to some extent.
Based on the above situation, the purpose of this work is to study the heat and mass transfer process of falling film, considering the change of water evaporation and film thickness, to determine the temperature distribution of the film under constant heat flux boundary conditions. In particular, the effect of film evaporation rate on heat transfer will be analyzed. In terms of application background, as mentioned above, the research results of this study will serve to eliminate infrared contrast between target and background, but this is not the content of this study.
Mathematical formulation
The physical model for a laminar falling liquid film due to gravity along a vertical plate is illustrated in Figure 1. The plate is subjected to a uniformed heat flux as the boundary condition. The liquid film flows down with an inlet temperature and an inlet mass flow rate .
FIGURE 1
Assumptions
For the mathematical formulation of the problem, the flowing simplified assumptions should be considered:
1. The flow is two dimensional and steady state with constant thermal-physical properties because the temperature range is not very large.
2. The flow is laminar, and the surface of the liquid film is non-oscillatory. This assumption is reasonable by virtue of the low velocities associated with the film.
3. The inertial force between liquid film and air is ignored because of the stillness of the air and the low speed of the liquid film.
4. The heat is transferred across the film by conduction only due to the laminar flow.
Liquid film flow model
According to the above assumptions and Figure 1, the Navier–Stokes equations describing the laminar falling film can be reduced to (Nusselt, 1916; Yang and Yuen, 2018):where is the velocity, is the dynamic viscosity, is the water density, and is gravity acceleration.
The momentum equation adopts the boundary conditions of no slip-on wall surface and no stress on film surface according to the assumptions mentioned above:where Δ is the thickness of the liquid film.
Combining Eqs 1–3, the velocity , thickness , mass flow rate , and mean velocity of the liquid falling film can be induced as follows (Nusselt, 1916; Yang and Yuen, 2018):where is the kinematic viscosity of water.
Energy conservation equation with evaporation considered
Considering the limit layer approximation, the energy equation iswhere is the thermal diffusion coefficient of water.
The total heat load of the wall has two dissipation ways: one is convective heat transfer between the liquid film and the wall, and the other is the evaporation of the film surface. The thermal balance in the liquid film is shown in Figure 2 and is expressed as
FIGURE 2
Eq. 9 can be simplified by removing the second term on the left and the first term on the right:where is the evaporation latent heat and is the specific heat capacity of water. Since the thermal characteristics of the fluid remain unchanged, the average temperature of the liquid film can be defined as
By substituting Eqs 6, 11 into Eq. 10, the first and second terms on the right side of Eq. 10 may be derived as the function of and :
Therefore, Eq. 10 can be written as
Furthermore, Eq. 14 is rearranged as
with the boundary conditions:
.
Here, the evaporation rate of the liquid film is supposed to be uniform. Therefore, the thickness reduction of the liquid film is proportional to the flowing distance x. Eq. 15 is a nonhomogeneous linear differential equation and can be solved by utilizing the parametric variation method:
For the fully developed region of the liquid film with constant heat flux boundary conditions, there is a relationship of (Incropera et al., 1996). By introducing the velocity distribution Eq. 4 and a newly defined dimensionless quantity, , Eq. 8 can be changed into
The boundary conditions of Eq. 18 are
By integrating Y in Eq. 18 and using boundary conditions (Eq. 19, Eq. 20), it can be obtained:
Then, the relationship between bulk and wall temperature is yielded based on Eqs. 4, 21:
In addition, the other form of temperature distribution is gained by substituting Eq. 22 into Eq. 21:in which and are obtained from Eq. 16 and Eq. 17, respectively.
Energy conservation equation without evaporation considered
For the case without considering the evaporation process, the liquid film thickness remains unchanged, which indicates that the evaporation rate and . Then, in Eq. 16 and in Eq. 17 can be expressed as
Furthermore, Eq. 23 is reduced with Eqs 24, 25 intowhich describes the temperature distribution of liquid film without evaporation considered.
Evaporation rate of liquid film
As we know, there exists an analogy between mass transfer and heat transfer (Incropera et al., 1996):where , , , and are dimensionless Nusselt, Prandtl, Sherwood, and Schmidt numbers, respectively. Among the four dimensionless variables, is the characteristic length of the liquid film (here is the length of the liquid film), and are the convection heat transfer coefficient and mass transfer coefficient, respectively, between the liquid film surface and air, is the binary mass diffusion coefficient. In addition, n is a constant number, equal to 1/3 for most applications.
Expanding the four dimensionless numbers in Eq. 27, the following relationship can be obtainedwhere is dimensionless Lewis Number and and are the density and specific heat capacity of air.
The evaporation quality of the liquid film iswhere and are the density of water vapor at water vapor saturation and infinity, respectively and and are the width and length of the liquid film.
Supposing the thickness reduction of the liquid film is proportional to the flowing distance x, can be derived from Eq. 29 as
By substituting Eq. 28 into Eq. 30, the evaporation rate of the liquid film is
Convection heat transfer coefficient between liquid film and wall
The heat transfer coefficient is defined as
If the evaporation effect were not considered, will be 0 and . Eq. 22 can be simplified by replacing with Eq. 17 as
Therefore, the convection heat transfer coefficient without evaporation considered is
If the evaporation effect were considered, the convection heat transfer coefficient can be derived from Eqs 22, 32:where is obtained from Eq. 17.
Results and discussion
In order to verify the validity, the analytical solutions with and without evaporation were compared with the experimental data of Wang et al. (2000) in Table 1. Tin in Table 1 represents the inlet temperature. Exp-Tout, Evap-Tout, and Non-evap-Tout represent the outlet temperature of in reference experiment, evaporation solution, and non-evaporation solution, respectively. The errors in Table 1 represent the relative deviation of outlet temperature between the evaporation solution and the experiment.
TABLE 1
| No. | Γ(kg/m s) | qw(W/m2) | Tin(°C) | Exp-Tout(°C) | Evap-Tout(°C) | Non-evap-Tout(°C) | Error |
|---|---|---|---|---|---|---|---|
| 1 | 0.307 | 13,240 | 17.9 | 20.8 | 21.9 | 23.1 | 5.02% |
| 2 | 0.375 | 17,420 | 19.1 | 22.1 | 23.7 | 24.7 | 6.75% |
| 3 | 0.375 | 17,220 | 18.8 | 22.1 | 23.4 | 24.3 | 5.56% |
| 4 | 0.375 | 17,630 | 18.8 | 22.2 | 23.5 | 24.5 | 5.53% |
| 5 | 0.43 | 17,850 | 16.7 | 18.7 | 20.8 | 21.7 | 10.10% |
| 6 | 0.43 | 19,460 | 16.7 | 19.1 | 21.3 | 22.1 | 10.33% |
| 7 | 0.45 | 24,860 | 17.9 | 21.9 | 23.7 | 24.5 | 7.59% |
| 8 | 0.47 | 23,800 | 17.5 | 23.0 | 22.8 | 23.6 | 0.88% |
| 9 | 0.47 | 25,465 | 17.7 | 23.7 | 23.5 | 24.2 | −0.85% |
| 10 | 0.47 | 27,200 | 17.8 | 24.1 | 24.0 | 24.8 | −0.42% |
| 11 | 0.47 | 13,950 | 17.6 | 20.7 | 20.4 | 21.2 | −1.47% |
| 12 | 0.47 | 20,140 | 17.9 | 22.5 | 22.3 | 23.1 | −0.90% |
| 13 | 0.47 | 21,480 | 18.1 | 22.9 | 22.8 | 23.6 | −0.44% |
| 14 | 0.47 | 20,800 | 18.1 | 23.0 | 22.7 | 23.4 | −1.32% |
| 15 | 0.47 | 18,870 | 17.5 | 21.8 | 21.6 | 22.3 | −0.93% |
| 16 | 0.474 | 23,500 | 17.6 | 21.3 | 22.8 | 23.6 | 6.58% |
| 17 | 0.641 | 21,710 | 18.4 | 20.6 | 21.9 | 22.5 | 5.94% |
Comparison of Tout under different conditions.
Table 1 shows that the evaporation solutions established in this article are in good agreement with the experimental values, and the error values are all less than 10%. This emphasizes that, in the field of this study, the evaporative heat transfer of liquid film cannot be ignored unless the requirements for solving accuracy are relaxed.
The model results are compared with the experiment in the data graph form, as shown in Figures 3, 4. These two graphs represent the relationship between Tout and qw. When Tin is almost the same, the outlet temperature Evap-Tout of evaporation solution has a good consistency with Exp-Tout of the experiment, whereas the Non-evap-Tout of the non-evaporation solution has a relatively poor consistency with experiment. In detail, the outlet temperature got by the analytical solution increases strictly linearly with the heat load, whereas that got by the experiment shows an approximate trend like the former. This phenomenon is easy to understand because there are always some measurement errors in the test. For example, in Figure 4, in the two experimental cases with the same inlet temperature of 18.1°C and different heat loads of 20,800 and 21,480 W/m2, the seemingly contradictory outlet temperature intensity is shown, which is 23°C and 22.9°C, respectively.
FIGURE 3
FIGURE 4
For the application of falling liquid film presented in this article, the heat flux loaded onto the film is set at 20,000 w/m2, and the mass flow rate is 0.3 kg m−1 s−1. Under the same initial and boundary conditions, the final falling film depth will depend on the evaporation rate. By calculating Eq. 31, the evaporation rate dΔ/dx of the liquid film with different air heat transfer coefficients can be obtained. When the values of heat transfer coefficients are 0, 8, 12, and 20 W/(m2°C), dΔ/dx are 0, −3.88e-7, −5.82e-7, and −9.70e-7, respectively. Figure 5 shows that the thickness of liquid film decreases with the increase in x. This figure shows that if the evaporation were not considered, the thickness of the falling film would keep the same; otherwise, that will decrease along with the flow distance. Note that the existence of evaporation does not change the thickness of the liquid film a lot, approximately 0.03%–0.04%.
FIGURE 5
Figure 6 shows the bulk temperature of the liquid film with different dΔ/dx. Along with the flow distance of film, the bulk temperature increases linearly for both evaporation and non-evaporation solutions. Nevertheless, the temperature decreases rapidly with the increase in dΔ/dx. At the outlet, the temperature with dΔ/dx = 9.7e-7 is 26.03°C, about 2°C less than that of the non-evaporation solution. This confirms the general notion that heat transfer is more efficient in the presence of phase transitions.
FIGURE 6
Eq. 10 represents the dissipation of wall energy. For the two terms on the right, the first represents the energy taken away by the flow of the liquid film, and the second term represents that by evaporation. Figure 7 shows the relationship between the percentage of energy carried by evaporation and the evaporation rate (dΔ/dx). This figure shows that the share of evaporation energy increases with the increase in dΔ/dx. When the natural convection heat transfer coefficient is 12W/(m2°C), dΔ/dx is −5.82e-7, and the evaporation rate accounts for 14.7% of the total energy.
FIGURE 7
Conclusion
The energy conservation equation of falling liquid film is established by thermal equilibrium analysis, and the analytical solutions with and without evaporation effect considered are obtained by means of parameter variational and integral methods. The validity of the analytical solution is verified by comparing it with the experimental data. It is proved that the evaporative heat transfer of the liquid film cannot be ignored unless the requirement of solving precision is relaxed.
In order to obtain the analytical solution, the physical problem is simplified appropriately. For example, the flow is assumed to be laminar flow and is quasi-two-dimensional. Fully turbulent flow and three-dimensional problems will be considered in future studies. In addition, this study applies to the case where there is no liquid film boiling. Therefore, in subsequent studies, the problem of high heat flux needs to be focused.
Statements
Data availability statement
The original contributions presented in the study are included in the article/supplementary material. Further inquiries can be directed to the corresponding author.
Author contributions
All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.
Funding
This study was funded by the National Natural Science Foundation of China (11504426) and the Independent Research Projects of Naval University of Engineering (No.20190008).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Abbreviations
thermal diffusion coefficient (m2 s−1)
width of plate or film (m)
specific heat for air (J kg−1 K−1)
specific heat for water (J kg−1 K−1)
binary mass diffusion coefficient (m2 s−1)
gravitational acceleration (m s−2)
heat transfer coefficient (W m−2 K−1)
heat transfer coefficient between film surface and air (W m−2 K−1)
evaporation heat transfer coefficient (W m−2 K−1)
non-evaporation heat transfer coefficient (W m−2 K−1)
evaporation latent heat (J kg−1)
mass transfer coefficient (m s−1)
length of plate (m)
Lewis number
evaporation quality (kg)
Nusselt number
Prandtl number
heat flux of wall (W m−2)
Schmidt number
Sherwood number
temperature (oC)
the bulk temperature (oC)
the inlet bulk temperature (oC)
velocity of the liquid film (m s−1)
mean velocity of the liquid film (m s−1)
coordinate in the flowing direction (m)
coordinate in the transverse direction (m)
Nomenclature
Abbreviations
thermal diffusion coefficient (m2 s−1)
width of plate or film (m)
specific heat for air (J kg−1 K−1)
specific heat for water (J kg−1 K−1)
binary mass diffusion coefficient (m2 s−1)
gravitational acceleration (m s−2)
heat transfer coefficient (W m−2 K−1)
heat transfer coefficient between film surface and air (W m−2 K−1)
evaporation heat transfer coefficient (W m−2 K−1)
non-evaporation heat transfer coefficient (W m−2 K−1)
evaporation latent heat (J kg−1)
mass transfer coefficient (m s−1)
length of plate (m)
Lewis number
evaporation quality (kg)
Nusselt number
Prandtl number
heat flux of wall (W m−2)
Schmidt number
Sherwood number
temperature (oC)
the bulk temperature (oC)
the inlet bulk temperature (oC)
velocity of the liquid film (m s−1)
mean velocity of the liquid film (m s−1)
coordinate in the flowing direction (m)
coordinate in the transverse direction (m)
Greek symbols
mass flow rate per unit width (kg m−1 s−1)
liquid film thickness (m)
thermal conductivity (W m−1 K−1)
dynamic viscosity of the fluid (kg m−1 s−1)
kinematic viscosity of the fluid (m2 s−1)
density of the water (kg m−3)
density of the air (kg m−3)
density of the water vapor (kg m−3)
Subscripts
component of binary mixture
air
evaporation
fluid
inlet condition
mean; mass
non-evaporation
wall; water
References
1
AssadM. H.LampnenM. J. (2002). Mathematical modeling of falling liquid film evaporation process. Int. J. Refrig.25, 985–991. 10.1016/s0140-7007(01)00064-0
2
CherifA. S.JabrallahS. B.CorriouJ. P.BelghithA. (2010). Intensification of the liquid film evaporation in a vertical channel. Desalination250, 433–437. 10.1016/j.desal.2009.09.071
3
ChoiC. -J.ChoH. K. (2022). Numerical investigations of liquid film offtake by transverse gas flow in a downcomer annulus geometry. Front. Energy Res.10, 847458. 10.3389/fenrg.2022.847458
4
GuL. D.MinJ. C.TangY. C. (2018). Effects of mass transfer on heat and mass transfer characteristics between water surface and airstream. Int. J. Heat Mass Transf.122, 1093–1102. 10.1016/j.ijheatmasstransfer.2018.02.061
5
HeS.AnP.LiJ.JacksonJ. (1998). Combined heat and mass transfer in a uniformly heated vertical tube with water film cooling. Int. J. Heat Fluid Flow19, 401–417. 10.1016/s0142-727x(98)10021-8
6
IncroperaF. P.DeWittD. P.BergmanT. L. (1996). Fundamentals of heat and mass transfer. New York: John Wiley & Sons.
7
LeuJ.JangJ.ChouY. (2006). Heat and mass transfer for liquid film evaporation along a vertical plate covered with a thin porous layer. Int. J. Heat Mass Transf.49, 1937–1945. 10.1016/j.ijheatmasstransfer.2005.11.004
8
LeuJ. S.JangJ. Y.ChouW. C. (2009). Convection heat and mass transfer along a vertical heated plate with film evaporation in a non-Darcian porous medium. Int. J. Heat Mass Transf.52, 5447–5450. 10.1016/j.ijheatmasstransfer.2009.06.033
9
MalikaI.MalekT. B.AberkaneS.MouderesM.SpitemP.GhezalA. (2017). Analytical and numerical study of the evaporation on mixed convection in vertical rectangular cavity. Fluid Dyn. Material Process13 (2), 85–105. 10.3970/fdmp.2017.013.085
10
MeyerT. (2015). Analytical solution for combined heat and mass transfer in laminar falling film absorption with uniform film velocity - diabatic wall boundary. Int. J. Heat Mass Transf.80, 802–811. 10.1016/j.ijheatmasstransfer.2014.09.049
11
MeyerT. (2014). Analytical solution for combined heat and mass transfer in laminar falling film absorption with uniform film velocity – isothermal and adiabatic wall. Int. J. Refrig.48, 74–86. 10.1016/j.ijrefrig.2014.08.005
12
NasrA.Al-GhamdiA. S.FeddaouiM’barek.AlsoufiM. S.Ben NasrallahS. (2017). Numerical study of falling binary liquid film evaporation: Liquid film thickness. Desalination Water Treat.74, 35–43. 10.5004/dwt.2017.20606
13
NasrA.Debbissi HfaiedhC.Ben NasrallahS. (2011). Numerical study of evaporation by mixed convection of a binary liquid film. Energy36, 2316–2327. 10.1016/j.energy.2011.02.039
14
NusseltW. (1916). Die Oberflächenkondensation des Wasserdampfes. Z. Ver. Deutch. Ing.60, 541–546.
15
PengY. S.YangL.DuY. C. (2012). Integral analysis of heat transfer on falling laminar liquid Film with constant heat flux. Adv. Mat. Res.516-517, 30–35. 10.4028/www.scientific.net/amr.516-517.30
16
PengY. S.YangL.DuY. C. (2012). Study on application of subcooling films in temperature field and infrared stealth technology. Infrared Laser Eng. (China)41 (10), 2572–2577.
17
SaouliS.Aïboud-SaouliS. (2004). Second law analysis of laminar falling liquid film along an inclined heated plate. Int. Commun. Heat Mass Transf.31 (6), 879–886. 10.1016/s0735-1933(04)00074-0
18
SongB.InabaH.HoribeA.OzakiK. (1999). Heat, mass and momentum transfer of a water film flowing down a tilted plate exposed to solar irradiation. Rev. Gen. Therm.38, 384–397. 10.1016/s0035-3159(99)80028-2
19
TliliI.MustafaM. T.KumarK. A.SandeepN. (2020). Effect of asymmetrical heat rise/fall on the film flow of magnetohydrodynamic hybrid ferrofluid. Sci. Rep.10, 6677. 10.1038/s41598-020-63708-y
20
TliliI.SajadiS. M.BaleanuD.GhaemiF. (2022). Flat sheet direct contact membrane distillation study to decrease the energy demand for solar desalination purposes. Sustain. Energy Technol. Assessments52, 102100. 10.1016/j.seta.2022.102100
21
TurkyilmazogluM. (2021). Heat absorption due to falling film with imposed uniform mass fraction at the wall. Int. J. Heat Mass Transf.177, 121585. 10.1016/j.ijheatmasstransfer.2021.121585
22
WangB. X.ZhangJ. T.PengX. F. (2000). Experimental study on the dryout heat flux of falling liquid film. Int. J. Heat Mass Transf.43, 1897–1903. 10.1016/s0017-9310(99)00279-3
23
YangJ. C.YuenW. W. (2018). Heat transfer in a falling laminar liquid film with in-depth radiation absorption. Int. Commun. Heat Mass Transf.94, 47–52. 10.1016/j.icheatmasstransfer.2018.03.009
24
YueY.YangJ.LiX.SongY.ZhangY.ZhangZ. (2021). Experimental research on falling film flow and heat transfer characteristics outside the vertical tube. Appl. Therm. Eng.199, 117592. 10.1016/j.applthermaleng.2021.117592
Summary
Keywords
liquid film, evaporation, heat and mass transfer, laminar flow, analytical solution
Citation
Du Y, Wang S, Li D, Peng Y and Yang L (2022) Effect of evaporation rate for the heat and mass transfer of the laminar liquid falling film in still air. Front. Energy Res. 10:899023. doi: 10.3389/fenrg.2022.899023
Received
18 March 2022
Accepted
27 July 2022
Published
30 August 2022
Volume
10 - 2022
Edited by
Kim Choon Ng, King Abdullah University of Science and Technology, Saudi Arabia
Reviewed by
Iskander Tlili, Monastir, Tunisia
Kasra Ghasemi, University of Guelph, Canada
Updates
Copyright
© 2022 Du, Wang, Li, Peng and Yang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yongcheng Du, dycheng@yeah.net
This article was submitted to Process and Energy Systems Engineering, a section of the journal Frontiers in Energy Research
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.