Research article Topical Sections

Solar energy storage in evacuated tubes solar collector using nanofluid embedded in a saturated porous media in the fully developed region: Al2O3 nanofluid embedded in graphite as a saturated porous media

  • The evacuated tube solar collector is considered an efficient, convenient, and economical option used to convert solar energy into heat. In this work, enhancement of evacuated tubes solar collector performance and the potential for energy storage by using Al2O3 water-based nanofluid embedded in Graphite as a saturated porous media was presented and studied theoretically. The Governing equations derived from the principles of conservation of mass, momentum, and energy were written in a dimensionless form, and solved these equations analytically for the fully developed region, and numerically throughout the entrance region, and the heat flux data for Amman, Jordan was used in this study. The analysis of the effect of different parameters as porosity, pore diameter, nanoparticles solid volume fraction, pressure, and radius of conduit on temperature variation through the pipe and Nusselt number was done. The results show that the flow of nanofluid in porous media enhances the evacuated tube's performance compared to just using water in the same medium under the same conditions. In August, the temperature variation in the evacuated tube at the same condition for Al2O3 water-based nanofluid for 1.5% volume fraction reach 42 ℃ and 73 ℃ for 6% volume fraction. In contrast, In January for the 1.5% volume fraction the temperature variation reached 25 ℃ and 43 ℃ for 6% volume fraction. On the other hand, the stored energy approximately reached 30 kJ/kg and 45 kJ/kg for solid volume fraction 1.5% and 6% respectively at the same conditions.

    Citation: Mohannad B. Khair, Hamzeh M. Duwairi. Solar energy storage in evacuated tubes solar collector using nanofluid embedded in a saturated porous media in the fully developed region: Al2O3 nanofluid embedded in graphite as a saturated porous media[J]. AIMS Energy, 2021, 9(4): 854-881. doi: 10.3934/energy.2021040

    Related Papers:

    [1] Duong Van, Diaz Gerardo . Carbon dioxide as working fluid for medium and high-temperature concentrated solar thermal systems. AIMS Energy, 2014, 1(1): 99-115. doi: 10.3934/energy.2014.1.99
    [2] Akthem Mohi Al-Abdali, Handri Ammari . Thermal energy storage using phase-change material in evacuated-tubes solar collector. AIMS Energy, 2022, 10(3): 486-505. doi: 10.3934/energy.2022024
    [3] Eman S. Maayah, Hamzeh M. Duwairi, Banan Maayah . Analytical model of solar energy storage using non—Newtonian Fluid in a saturated porous media in fully developed region: carboxymethyl cellulose (CMC) and graphite model. AIMS Energy, 2021, 9(2): 213-237. doi: 10.3934/energy.2021012
    [4] Hayder S. Al-Madhhachi, Ahmed M. Ajeena, Nihad A. Al-Bughaebi . Dynamic simulation and energy analysis of forced circulation solar thermal system in two various climate cities in Iraq. AIMS Energy, 2021, 9(1): 138-149. doi: 10.3934/energy.2021008
    [5] Daniel Chuquin-Vasco, Cristina Calderón-Tapia, Nelson Chuquin-Vasco, María Núñez-Moreno, Diana Aguirre-Ruiz, Vanesa G. Lo-Iacono-Ferreira . Mathematical modeling of a binary ORC operated with solar collectors. Case study—Ecuador. AIMS Energy, 2023, 11(6): 1153-1178. doi: 10.3934/energy.2023053
    [6] Nassir D. Mokhlif, Muhammad Asmail Eleiwi, Tadahmun A. Yassen . Experimental evaluation of a solar water heater integrated with a corrugated absorber plate and insulated flat reflectors. AIMS Energy, 2023, 11(3): 522-539. doi: 10.3934/energy.2023027
    [7] Ruaa Al Mezrakchi . Investigation of various hybrid nanofluids to enhance the performance of a shell and tube heat exchanger. AIMS Energy, 2024, 12(1): 235-255. doi: 10.3934/energy.2024011
    [8] Van Duong Nguyen, Duy Duc Ha, Thi Hoai Thu Nguyen, Duc Huy Nguyen . Generation expansion planning toward Net Zero target considering co-firing green fuel, CCUS installation and early decommissioning in thermal power plants: A case study of Vietnam. AIMS Energy, 2025, 13(3): 540-568. doi: 10.3934/energy.2025021
    [9] Wei An, Yifan Zhang, Bo Pang, Jun Wu . Synergistic design of an integrated pv/distillation solar system based on nanofluid spectral splitting technique. AIMS Energy, 2021, 9(3): 534-557. doi: 10.3934/energy.2021026
    [10] Xian-long Meng, Cun-liang Liu, Xiao-hui Bai, De-hai Kong, Kun Du . Improvement of the performance of parabolic trough solar concentrator using freeform optics and CPV/T design. AIMS Energy, 2021, 9(2): 286-304. doi: 10.3934/energy.2021015
  • The evacuated tube solar collector is considered an efficient, convenient, and economical option used to convert solar energy into heat. In this work, enhancement of evacuated tubes solar collector performance and the potential for energy storage by using Al2O3 water-based nanofluid embedded in Graphite as a saturated porous media was presented and studied theoretically. The Governing equations derived from the principles of conservation of mass, momentum, and energy were written in a dimensionless form, and solved these equations analytically for the fully developed region, and numerically throughout the entrance region, and the heat flux data for Amman, Jordan was used in this study. The analysis of the effect of different parameters as porosity, pore diameter, nanoparticles solid volume fraction, pressure, and radius of conduit on temperature variation through the pipe and Nusselt number was done. The results show that the flow of nanofluid in porous media enhances the evacuated tube's performance compared to just using water in the same medium under the same conditions. In August, the temperature variation in the evacuated tube at the same condition for Al2O3 water-based nanofluid for 1.5% volume fraction reach 42 ℃ and 73 ℃ for 6% volume fraction. In contrast, In January for the 1.5% volume fraction the temperature variation reached 25 ℃ and 43 ℃ for 6% volume fraction. On the other hand, the stored energy approximately reached 30 kJ/kg and 45 kJ/kg for solid volume fraction 1.5% and 6% respectively at the same conditions.



    Abbreviations: Q: amount of heat stored, J; dp: average particle or fiber diameter of the porous media, m; N: avogadro number; h: convective heat transfer coefficient, W/(m2.K); D: diameter of the conduit, m; z: dimensionless length; r: dimensionless radius of the conduit; u: dimensionless velocity; keff: effective thermal conductivity of the porous medium, W/(m.K); df: equivalent diameter of a base fluid molecule, m; dnp: equivalent diameter of solid nanoparticles, m; Tf: final temperature, ℃; p: fluid pressure, Pa; q: heat flux per unit area, J/m2; Ti: initial temperature, ℃; NuD: local nsselt number; m: Mass flow rate, kg/s; m: mass of heat storage medium, kg; Num: mean nusselt number; M: molecular weight of the base fluid, kg/mol; Nu: nusselt number; R: outer radius of the conduit, m; P: perimeter of the tube, m; K: permeability of porous media, m2; Pr: prandtl number; m2; P/L: pressure gradient, Pa/m; ReD: reynolds number based on darcy velocity; Cp: specific heat at constant pressure, J/(kg.k); q: stored heat per unit mass, J/kg; T: temperature, ℃; Te: temperature at the entrance, ℃; Tm: the mean temperature, ℃; υ: the velocity at radial (r) direction, m/s; k: thermal conductivity, W/(m.K); Vt: total volume of the porous medium, m3; uD: darcy velocity, m/s; u: velocity in z-direction, m/s; Vs: volume of the solid phase in the porous media, m3; ρ: density, kg/m3; ρf0: density of the base fluid at 293 K, kg/m3; θ: dimensionless temperature; μ: dynamic viscosity, Pa.s; αeff: effective thermal diffusivity, m2/s; μeff: effective viscosity, Pa.s; : nanofluid volume fraction; ε: porosity of porous media

    Renewable energy is considered an important source of clean energy. The negative impact of conventional energy sources on the environment and the fact that these resources will deplete, made most of the researches turn for the non-conventional renewable energy like solar, wind, geothermal, and hydropower.

    Solar energy is the most available, environmentally friendly renewable energy source, a possible utilization by converting it to useful heat energy by using solar collectors. The evacuated tube solar collector is an efficient, convenient, and economical option for doing this [1].

    The most straightforward, best promising technique for enhancing the performance of the evacuated tube solar collectors is using nanofluid as a working fluid, thus enhancing the thermal performance of the heat transfer fluid [2].

    Nanofluids are an innovative, promising way to increase the heat transfer rate in the solar thermal collectors, by using them as working fluids instead of the conventional base fluids used before. Choi first proposed nanofluid in 1995 [3], Nanofluids formed by dispersing nano-solid particles from metal or metallic oxide into the base fluid, such as oil, water, and ethylene glycol. The nanofluids are engineered colloids made of solid nanoparticles with a size typically of 1–100 nm, which are suspended in liquids [4].

    The use of proper techniques for heat transfer enhancement is vital to optimizing energy devices, and nanofluid utilization is one of the relatively new solutions to accommodate this purpose. Since nanofluids have high thermal conductivity compared to the base fluids used before, they can improve the heat transfer rate in solar thermal collectors as one of the engineering systems that may be used with it. Currently, nanofluid is a hot topic, and many researchers around the world are working on it [5].

    From the need to store the solar energy during its availability and reuse it at the time of requirement and maximize the efficiency of the energy conversion from solar to thermal energy, porous materials were widely used in solar systems, because it can improve the heat transfer characteristics in solar systems by offering a high ratio of surface area to volume. It can act as an energy storage medium since it represents a packed bed storage system at a macroscopic level [6].

    An attractive technique to enhance the heat transfer is the usage of porous media, such as the application of metal-based porous materials such as copper foams in channels and heat exchangers. Porous media consists of a solid matrix containing pores (voids), which are typically filled with fluid. Porous media are rigid, open-cell, and saturated, which means that the pores are connected and filled with the fluid completely so that the fluid can flow through the voids. Recently, the technique of using both nanofluid and porous media has received considerable attention and has led to extensive investigations in this field [7].

    Porous media has the potential to increase the contact surface area between the fluid and solid surface. On the other hand, nanoparticles dispersed in nanofluid enhance the effective thermal conductivity. Therefore, it seems that using both porous media and nanofluid can augment the efficiency of typical thermal systems dramatically [5].

    One of the techniques to enhance the convective heat transfer is to employ solid materials, such as nanoparticles or porous media, with higher thermal conductivity concerning the traditional base fluids. The usage of such materials also allows a significant increase of heat transfer contact area inside the fluids, and the heat transfer passes from a surface exchange to a bulk or mass exchange. The coupling of nanoparticles and porous media with high thermal conductivity is an interesting and promising solution. Many research activities have been developed in the convective heat transfer of nanofluids in porous media [8].

    Muhammad [9] reviewed different models to enhance the thermal performance of flat plate and evacuated tube solar collectors by using nanofluids.

    Gautam and Saini [6] considered the packed bed as a porous media at the macroscopic level, and they did a review on technical, applications, and economic aspects of packed bed solar thermal energy storage systems. The description of the packed bed storage system and its working principle and design aspects were discussed. Various factors affecting the heat transfer involved within a packed bed were described. The packed bed storage system (PBSS) has three major components, which are the storage medium, energy transfer mechanism, and contaminated systems. The storage medium keeps the sensible heat inside it in the form of sensible heat, latent heat, or thermo-chemical form. The energy transfer mechanism used to extract or supply energy to the storage medium by the flow of the heat transfer fluid. The job for the contaminated system is to hold storage medium, energy transfer mechanism, and different accessories needed for the PBSS.

    Rashidi [10] reviewed the applications of porous materials in solar systems, their advantages, and their limitations. Different energy conversion and storage systems that use porous materials were presented. It was concluded that the use of porous media inside a solar system is considered as an easy and simple modification that can enhance and improve the efficiency of that system. The effect of the porous media is positive for heat transfer. In comparison, the opposite happens with the pressure drop. Thus, it required more pumping power.

    Ahmed and Mohammed [11] studied the performance of the double pass PV/thermal collector air collector. They concluded that the heat transfer area is increased due to the use of the porous medium. Thus, it reflects directly to increase in the temperature of the exit air from the hybrid PV/thermal collector and the thermal efficiency of the system was increased with the presence of the porous media almost 29%.

    Saryazdi [4] studied numerically the forced convection flow and heat transfer of nanofluid flowing inside a straight circular pipe filled with a saturated porous medium, the finite volume method used to analyze the heat transfer coefficient of Al2O3-water nanofluid. It is approved that the local Nusselt number is increased with the increase of the nanoparticles volume fraction. Still, the increase is not significant and almost not clearly observed. That is due to the effective thermal conductivity of porous media that has great value, and this originally comes from the high thermal conductivity of the porous media material used.

    Tajik Jamal-Abad [12] experimented with the effect of partially metal foam inside the absorber of parabolic trough solar collectors. The copper foam was used, and the results showed that the absorber filled with metal foam enhanced and improved the efficiency because the enhancement happens in the thermal conductivity.

    AlMasa'deh and Duwairi [13] modeled analytically and numerically the fluid flow and heat transfer inside a saturated porous conduit at constant surface heat flux. The energy, momentum, and continuity equations were solved by using a theoretical model with constant heat flux boundary conditions through the fully developed and entrance region. The results indicate that the porous media has very good potential to enhance the heat transfer rate.

    Saedodin [14] carried out experimental and numerical investigations on the performance of a flat plate solar collector filled with porous copper metal foam. The numerical and experimental results showed that the use of porous media enhances the flat plate solar collector's performance, where the Nusselt number and the maximum thermal efficiency are enhanced by 82% and 18.5%, respectively.

    Esfe [15] reviewed the heat transfer of nanofluids in porous media for different nanofluids and porous media in a forced, natural, and mixed convection systems. It is found that porous media can improve the performance and heat transfer in various applications. While the nanofluids considerably enhance the thermal conductivity of the working fluid. Consequently, using both porous media and nanofluids can significantly improve the thermal efficiency in the thermal systems.

    Alihosseini and Jafari [16] studied the effect of porous medium configuration on nanofluid heat transfer using the computational fluid dynamics model. Al2O3 water-based nanofluid was used as a working heat transfer fluid with aluminum foam as a porous medium. Different configurations were used and compared with each regarding the enhancement parameter and the Nusselt number. It was concluded that, the fully saturated porous cylinder gave the best results among the other configurations.

    Esmaeili [17] investigated numerically the performance of a direct absorption solar collector using CuO water-based nanofluid and copper oxide porous foam by using the finite volume technique to solve the mathematical model suggested. The best performance and efficiency obtained when the collector was fully filled with porous metal oxide foam. Moreover, the results proved that using nanofluid and porous media enhances efficiency by 26.8% compared to the flow of water without porous media.

    Tay [18] reviewed different techniques for phase change material that moved in the thermal energy storage systems. The results showed that the use of transporting phase change material gave an efficient economic low cost heat transfer and a high compactness factor for the thermal energy storage. Another work was done by Liu [19] were they reviewed several models to enhance the thermal conductivity of phase change material for thermal energy storage. Where the importance of enhance the conductivity of the phase change material is essential due to reach a high energy charging/discharging rates.

    The enhancement of the finned heat pipe performance using nano-enhanced H2O as phase change material was studied by Lohrasbi [20], the study was based on a latent heat thermal energy storage system in a finned heat pipe impeded into the nano-enhanced phase change material. The results were obtained numerically and showed that it can be achieve the efficient design of the latent heat thermal energy storage system by using the maximum energy storage capacity as an evaluation parameter.

    Another model for energy storage was studied by Bazri [21], the latent heat storage tank was used in combination with a solar heater system using evacuated heat pipe solar collectors connected to a tank full with paraffin was as a phase change material. The results showed that there is a promising enhancement in the efficiency and reached to 58% compared to the baseline.

    Naghavi [22] studied an on-demand dynamic performance of a thermal battery in tankless domestic solar water heating, that done by collecting and storing the solar energy in a latent heat storage tank that used as a thermal battery, this solar heat used in a heat exchanger to provide a hot water for the domestic use. The results showed a promising results with a 2 m2 solar collector size that proved it can provide at least 112 liters per day in the worst weather condition.

    This study deal with the forced convection problem inside a saturated porous media with a constant surface heat flux, the Darcy model is taken for this case.

    A theoretical model for nanofluid flow in a saturated porous media was build and solved analytically, with studying the effect of some factors (porosity, pressure.etc) on the temperature variation to illustrate the difference of using water and nanofluid in energy storage. The study shows the energy storage in August (summer) and January (winter) in Amman city.

    For understanding the nanofluid thermophysical properties, the nanofluid volume fraction should be defined. The nanofluid volume fraction () is the constituent volume divided by all constituent volumes of the mixture before mixing, as shown in Eq 1. It is found that increasing the nanoparticles volume fraction increases thermal conductivity [23].

    =VnpVtotal (1)

    A weighted average process with the volume concentration between the nanoparticles and the base fluid is accomplished to find the nanofluid's thermal properties. For the density, a simple mass balance is applied, as shown in Eq 2.

    ρnf=ρnp+(1)ρf (2)

    For the specific heat, the simple weighted average with the concentration as in Eq 3 will not give an accurate result.

    Cpnf=(Cp)np+(1)(Cp)f (3)

    It is better to use an analytical model by assuming thermal equilibrium between nanoparticles and the base fluid as in Eq 4 [24].

    Cpnf=(ρCp)np+(1)(ρCp)fρnf (4)

    The Corcione Eq [25] will be used in this work to simulate the dynamic viscosity

    μnf=μf1134.87(dnpdf)0.31.03

    where df is the equivalent diameter of a base fluid molecule, given by

    df=0.1(6MNπρf0)1/3 (5)

    In which M is the molecular weight of the base fluid, N is the Avogadro number, and ρf0 is the density of the base fluid calculated at temperature T0=293K.

    For this work, Aluminum Oxide nanoparticles (Al2O3) with water as a base fluid were used as the working heat transfer fluid. The thermophysical properties of the Al2O3 and water are given in Table 1 [26].

    Table 1.  Thermo-physical properties of the water and Al2O3.
    ρ (kg/m3) Cp (J/kg K) k (W/m K) dnp (nm)
    Pure water 997.1 4179 0.613 -
    Al2O3 3970 765 40 25

     | Show Table
    DownLoad: CSV

    The properties of Graphite were used in calculations of this study are shown in Table 2 [27,28].

    Table 2.  Properties of graphite are used in analytical calculation.
    Thermal conductivity (W/m K) 25–470
    Heat capacity (J/g K) 0.7
    Average particles or pore diameter (mm) 350
    Density (kg/m3) 2260

     | Show Table
    DownLoad: CSV

    If total volume of the media, and the solid volume are known, literally porosity could be defined as follow.

    ε=1VsVt (6)

    where ε is the porosity of the medium, Vs and Vt are respectively the volume of the solid phase and total volume of the medium.

    Permeability is a geometrical characteristic of the porous media. It is important because it has a direct effect on the directional movement and the flow rate of the working fluids. For particle beds and fibers, permeability can be defined utilizing the Carman-Kozeny theory as follows [7,29].

    K=ε2dp2150(1ε)2 (7)

    where dp is the average particle or fiber diameter.

    Graphite has an important advantage that it does not react with any of the organic material. Thus, it considered having chemically inert characteristics at the room temperature. Graphite is being used in various applications such as electrical transfer devices, thermal transfer devices, and electrochemical super capacitors. Graphite is well known for the high thermal conductivity that it has. Thus, it gave it the capability to serve as a good porous media, where the high thermal conductivity will lead to enhance the heat transfer rate [30,31].

    The Darcy velocity of the nanofluid (uD) is given by the following Eq,

    uD=KμeffPL (8)

    where μeff is the effective viscosity, μ is the nanofluid viscosity and K is the permeability of the porous media. The Einstein's formula for dilute suspensions can be used [32].

    μeff=μnf[1+2.5(1ε)] (9)

    Reynold number can be calculated for nanofluid based on pores diameter or based on permeability. The studies found that calculating Reynolds number based on pore diameter is more accurate since permeability and inertial resistance are dependent on pore size [23].

    Reynolds number based on pores diameter can be found using the below Eq.

    ReD=ρnfudpμnf (10)

    And Reynolds number based on permeability is:

    ReK=ρnfuKμnf (11)

    The suggested model examines enhancement of the performance of evacuated tube by nanofluid flowing in a saturated porous media, the tube is taken to be evacuated to restrict convective and conductive heat losses to the ambient. The governing equations were written and solved in the two-dimensional system (r, z). The continuity, momentum, and energy Eqs were solved simultaneously for intention to Darcy effects on fluid flow and heat transfer inside a saturated porous conduit.

    As a governing principle in fluids, the fluid flow is controlled and governed by the conservation laws of momentum, energy and mass.

    In the present study, the suggested model was developed using the above conservation equations. Taking into consideration that the momentum equations for fluid flow with the presence of porous media differ from the absence of it.

    A two-dimensional forced convection incompressible nanofluid flow in a circular tube saturated with porous media, as illustrated in Figure 1.

    Figure 1.  Physical model and coordinate system.

    For the mathematical model in this study, the following assumptions were taken:

    1.The nanofluid flow is steady, incompressible, laminar, forced fluid flow in two-dimensions.

    2.The temperature of the nanofluid is below the boiling point.

    3.The nanofluid and solid matrix are everywhere in thermodynamics equilibrium.

    4. All the physical properties of the fluid are isotropic and homogeneous

    The dimensional governing equations that describe the convection flow can be written as follows:

    ● The continuity Eq (the conservation of mass)

    The continuity equation is based on the principle of conservation of mass, since the total mass flow in will equal the total mass flow out from a specified control volume if there is no mass storage or loss in the system.

    1r(υr)r+uz=0 (12)

    ● The z-momentum Eq (conservation of momentum):

    Since the velocity is pure Darcy velocity, the below equation can be written:

    u=uD=Kμpz (13)

    ● The r-momentum Eq (conservation of momentum):

    υ=Kμpr (14)

    ● The energy Eq (conservation of energy):

    uTz+υTr=αeff(2Tz2+2Tr2+1r(T)r) (15)

    where K, p, μ, T, and αeff is permeability, fluid pressure, dynamic viscosity, fluid temperature, and the effective thermal diffusivity (m2/s) respectively.

    where αeff can be calculated using the following Eq:

    αeff=keff(ρCp)nf (16)

    where

    keff=εknf+(1ε)ks (17)

    Is the effective thermal conductivity of the porous medium. knf and ks are the thermal conductivity of the nanofluid, and the solid porous media respectively, and ε is porosity of the porous media.

    The governing boundary conditions for the model suggested is:

    ● The velocity in the r-direction can be neglected (υ = 0).

    ● The heat transfer in the 𝑟 direction at the center of the pipe is equal to zero (𝜕T/r=0atr=0).

    Thus, the boundary conditions above were used in this work as the following formula

    ● At the outer surface:

    Atr=Rqw=+kTr=h(TwTc)

    ● At the center:

    Atr=0Tr=0

    Since the velocity in the r-direction was neglected (υ = 0), and there is no heat transfer in the radial direction at the center (T/r=0atr=0). The heat diffused in the z-direction were neglected compared to the heat diffused in the radial direction (r).

    Energy Eq can be reduced to:

    uTz=αeff(2Tz2+2Tr2+1rTr) (18)

    The initial point of the conduit has no heat transfer so the temperature is uniform and constant and can be presented as =Te, at z = 0.

    In the fully developed region, the relative shape of the temperature profile remains the same but Tm always vary in the Z-direction i.e., TmZ0.

    The temperature profile in the fully developed region can be found by the following Eq [33].

    z[Tw(z)T(z)Tw(z)Tc(z)]=0 (19)

    The conduction heat transfer rate in the axial direction along the pipe will be negligible compared to the one that is in the radial direction, so 2Tz2 was ignored compared to 2Tr2 [33].

    Thus, the energy Eq becomes:

    uTz=αeff(2Tr2+1rTr) (20)

    By arranging Eq 19 and substitute it into Eq 20, the energy Eq becomes

    udTwdz=αeff(2Tr2+1rTr) (21)

    Substitute temperature profile and Darcy velocity into energy Eq:

    (Kμpz)dTwdz=αeff(2Tr2+1rTr) (22)

    The left term of the above Eq is a constant, the symbol will be used instead of the left term in Eq 22 to simplify the solution.

    By substitute the left term in Eq 22 by and multiplying the Eq by (r):

    r(×αeff=r×(2Tr2+Tr)
    rαeff=r(rTr) (23)

    By integrating Eq 23 twice and applying the boundary conditions and because the model is subjected by constant heat flux we get

    T=Tw+q2Rk(r2R2) (24)

    Now to find Nusselt number we should find mean temperature Tm to substitute it in the heat transfer coefficient. So, from Newton's Law of Cooling heat transfer coefficient is:

    h=qTwTm (25)

    The mean temperature Tm can be found as the following [33]:

    Tm=2umR2R0uTrdr (26)

    Substitution Eq 24 into Eq 26 and the result is integrated over the interval [0,R]. Finally, mean temperature represented as follows:

    Tm=Twq4kR (27)

    Substitute Eq 27 into 25 to get:

    h=4kR (28)

    Now the Nusselt number is founded by substitute the heat transfer coefficient found in Eq 28 in Nusselt number Eq (hDk), where D=2R. Since the heat flux concentrated on a one half of the tube, the Nusselt number will be as:

    Nu=4

    In the entrance region or developing region where the boundary layer develops and reaches the tube center, the developing in velocity profile is then essentially fully developed before the heat transfer take place, and because of the presence of porous media, the velocity profile is not changing with z along the conduit which means there is only a thermal entrance region. The following dimensionless variables are used in case of constant heat flux at the wall:

    θ=TTeqD/ku=uuD
    r=rDz=zD (29)

    where θ is the dimensionless temperature, Te is the temperature at the entrance point, uis dimensionless velocity, uD is the Darcy velocity in the conduit, D is the diameter of the conduit (D=2R), z is the dimensionless length and r is the dimensionless radius of the conduit.

    By substituting dimensionless variables in the energy Eq, the energy Eq will be as the following after substitute the effective thermal diffusivity:

    uθz=1PrRe(2θr2+1rθr) (30)

    The dimensionless velocity u can be eliminated, since the velocity through porous media is constant and here is Darcy velocity which means u=uD so u=1.

    The following is the dimensionless boundary conditions for the above Eq.

    The initial condition:

    θ=0asz=0

    The boundary conditions

    θr=0asr=0and
    θr=1asr=0.5 (31)

    The linear second-order differential Eq 30 concerning boundary conditions Eq 31 was solved by finite difference method using MATLAB program.

    Thus, using the FDM approach we get the following Eq:

    θ(z,r)θ(zz,r)z=1PrRe(θ(z,r+r)2θ(z,r)+θ(z,rr)r2+1rθ(z,r)θ(z,rr)r) (32)

    We can write the dimensionless temperature at constant heat flux at the wall as:

    θ|r=0.5=TwTeqD/k (33)

    So, we can calculate local Nusselt number by the following formula.

    NuD=1θ|r=0.5 (34)

    where NuD is the local Nusselt number.

    The Nusselt number is defined in terms of the difference between the wall and the mean temperatures (TwTm).

    The mean temperature Tm can be found as the following [33]:

    Tm=2umR2R0uTrdr (35)

    Temperature T can be obtained from dimensionless form Eq 31:

    T=θ(qDk)+Te (36)

    Substitution Eq 34 into Eq 33 and u=um we get:

    TmTeqD/k=2R2R0θrdr (37)

    A right side represent a non-dimensional mean temperature θm.

    After substituting the dimensionless form from Eq 31 the result will be:

    θm=8r0θrdr (38)

    Equation 38 was solved using MATLAB, and the results used to find the Nusselt number using the following Eq.

    NuDm=12(θwθm) (39)

    Figures 2 and 3 show the results of these Eqs.

    Figure 2.  Mean Nusselt number along the conduit at different velocities for Al2O3 nanofluid.
    Figure 3.  Mean Nusselt number along the conduit at different Prandtl number for Al2O3 nanofluid.

    The enhancement of the convection heat transfer relative to the conduction through the same fluid layer is expressed by the Nusselt number. Where the higher the Nusselt number is, the more the convection heat transfer is effective.

    The Nusselt number shows the enhancement of heat transfer through a fluid boundary as a result of convection relative to conduction across the same fluid layer. The larger the Nusselt number, the more effective the convection.

    In the entrance region, the mean Nusselt number was found and it can be presented in an exponential distribution. The Nusselt number starts with a high value then dramatically decrease to reach the fully developed region. In the fully developed region, the Nusselt number that was found analytically remains constant. Obviously, the Nusselt number in the entrance region is decreased to reach the value of the fully developed region. The value of ∆z and ∆r in the numerical solution is the largest value when does not change the solution in order to the staked solution.

    Figures 4 and 5 present the dimensionless temperature θ with axial coordinate z* at different fluid velocities and Prandtl numbers (Pr) respectively for water and Al2O3 water-based nanofluid. It is noted that the temperature increases in the axial direction of the flow, thus because the continuous heating of the fluid in the entrance region. From the mentioned figures before, it is concluded that the increase in the fluid velocity and Prandtl number cause to decrease in temperature profile and this because of increasing in the fluid velocities means increasing in Reynolds numbers and thus increasing in inertia forces and increasing in heat transfer coefficient. Also, the increase in Prandtl number caused a decrease in fluid temperatures, that enhance heat transfer rates. In the same conditions, nanofluid reach higher dimensionless temperature than water did, and especially when the fluid velocities were low the dimensionless temperature difference was higher.

    Figure 4.  Dimensionless temperature θ at different axial coordinate z* at different velocities for water and Al2O3 nanofluid.
    Figure 5.  Dimensionless temperature θ at different axial coordinate z* at different Prandtl number for water and Al2O3 nanofluid.

    Along the pipe, the temperature is varied in the z-direction. The mean temperature variation in the axial direction can be calculated by applying the energy balance Eq as presented in Figure 6.

    Figure 6.  Energy and mass balance through conduit.

    A differential control volume is taken, and the energy conservation is applied at constant heat flux surface, that leads for the following:

    dTmdz=qwP˙mcp (40)

    where P is the perimeter of the tube, andcp is the specific heat.

    Since the temperature gradient is independent of the axial position in the fully developed flow, the shape of the temperature profile will be the same without any change along the tube in that region.

    By Integration of Eq 40 at (Tm=Tmi as z=0) and substituting into it P=2πR, and ˙m=ρumAc=ρumπR2 leads to:

    Tm=2qw(ρcp)effumR+Tmi (41)

    where

    (ρcp)eff=ε(ρcp)nf+(1ε)(ρcp)s (42)

    Is the heat capacity per unit volume and um is the mean velocity of the fluid which is equivalent to the Darcy velocity.

    Heat flux at the wall in Eq 41 is the irradiance (W/m2), the irradiance for Amman city in Jordan was taken in our work, and the data of the irradiance were taken from "Photovoltaic Geographical Information System" (PVGIS). The data for Amman irradiance is given in Table 3.

    Table 3.  Irradiance data for Amman. **From JRC Photovoltaic Geographical Information System (PVGIS).
    Time (hour)
    Month
    Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec
    5 0 0 0 0 0 0 0 0 0 0 0 0
    6 0 0 5.12 20.77 38.7 55.42 45.92 34.1 32.63 15.87 0.04 0
    7 41.56 51.47 92.5 130.71 179.5 219.23 213.7 195.69 191.11 167.3 113.67 66.21
    8 175.83 178.93 245.29 295.9 352.95 420.84 414.36 396.95 399.63 358.03 279.59 219.62
    9 315.54 318.47 405.65 477.34 519.7 592.1 596.66 589.21 596.79 535.94 440.39 357.06
    10 431.58 435.37 542.23 601.52 650.64 730.05 733.83 736.3 737.72 660.28 517.48 474.67
    11 485.98 509.91 603.29 669.99 706.8 758.03 814.76 818.19 805.38 715.28 560.88 520.97
    12 485.71 516.23 608.75 677.95 708.72 808.7 821.81 825.98 804.31 715.79 576.95 498.86
    13 435 470.95 574.22 620.6 669.94 747.64 762.42 765.26 704.74 618.68 495.33 440.45
    14 321.76 362.39 458.05 499.04 540.38 617.7 643.16 630.5 505.29 384.95 344.23 306.76
    15 198.71 237.23 319.87 353.13 371.14 450.3 469.56 455.88 381.99 278.13 170.96 170.01
    16 59.63 91 156.07 182.88 202.44 264.33 282.3 260.5 202.21 90.57 37.51 19.28
    17 0 4 25.21 38.68 51.73 84.69 97.69 78.55 26.61 0 0 0
    18 0 0 0 0 0 0 0 0 0 0 0 0
    **https://ec.europa.eu/jrc/en/pvgis

     | Show Table
    DownLoad: CSV

    The stored energy can be calculated by using the following Eq.

    Q=TfTimCpdT=mCp(TfTi)=mCpT (43)

    The sensible heat storage for the material depends on its heat capacity which is the energy density, and the thermal diffusivity which is the rate that the heat extracted and released.

    Where Q is the amount of heat stored, Ti and Tf are respectively, the initial temperature, and the final temperature, m is the mass of heat storage medium, and Cp is the specific heat.

    The stored energy per unit of mass can be calculated using the following Eq.

    Qm=ɋ=CpT (44)

    The enhancement of the convection heat transfer relative to the conduction through the same fluid layer is expressed by the Nusselt number. Where the higher the Nusselt number is, the more the convection heat transfer is effective.

    The Nusselt number shows the enhancement of heat transfer through a fluid boundary as a result of convection relative to conduction across the same fluid layer. The larger the Nusselt number, the more effective the convection.

    The local Nusselt number at different nanofluid velocities and Prandtl numbers were illustrated respectively, in Figures 7 and 8. The local Nusselt number was found by taking the difference between the outer surface and the entrance temperatures. The local Nusselt number profile was facing a sharp decrease until it hit the zero, that is because of the temperature difference decrease between the outer surface temperature and the entrance temperature along the tube. It is noted that the higher the nanofluid velocities and Prandtl numbers are, the longer the entrance region length is, and that is unpreferable. For that reason, it is recommended for the better performance to keep the nanofluid velocity and the Prandtl numbers as low as possible.

    Figure 7.  Local Nusselt number along the conduit at different velocities for Al2O3 nanofluid.
    Figure 8.  Local Nusselt number along the conduit at different Prandtl number for Al2O3 nanofluid.

    Figures 9 and 10 present the mean Nusselt number with different nanofluid velocities and Prandtl numbers. The sharp decreasing profile of the Nusselt number, that it starts with a high value then it faced a sharp sudden exponential decrease in the entrance region until it reached a constant value, which remains constant along the fully developed region. Noting that, friction factor and heat transfer coefficient stay constant in the fully developed region without any effect from the nanofluid velocity or the Prandtl number.

    Figure 9.  Mean Nusselt number along the conduit at different velocities for Al2O3 nanofluid.
    Figure 10.  Mean Nusselt number along the conduit at different Prandtl number for Al2O3 nanofluid.

    In the fully developed region, Figures 11 to 14 show a comparison between nanofluid with different solid volume fractions and water. The mean temperature variation with axial coordinate for nanofluid and water at different heat flux are shown in the figures. The performance of the nanofluid is better than water at the same heat flux. It is concluded that the solid volume fraction of the nanoparticles in the nanofluid has a significant effect on the temperature variation, where the increase in the fraction lead for increase in the temperature variation.

    Figure 11.  Temperature difference with Z coordinate at 1.5% volume fraction and different q for Al2O3 nanofluid compared with water.
    Figure 12.  Temperature difference with Z coordinate at 3% volume fraction and different q for Al2O3 nanofluid compared with water.
    Figure 13.  Temperature difference with Z coordinate at 4.5% volume fraction and different q for Al2O3 nanofluid compared with water.
    Figure 14.  Temperature difference with Z coordinate at 6% volume fraction and different q for Al2O3 nanofluid compared with water.

    Figures from 15 to 17 show the effect of pressure, pore diameter, porosity, and tube diameter at a constant heat flux (500 W/m2) on the temperature variation along the tube.

    Figure 15.  Temperature difference with Z coordinate at different volume fraction and pressure for Al2O3 nanofluid at 500 W/m2.
    Figure 16.  Temperature difference with Z coordinate at different volume fraction and porosity for Al2O3 nanofluid at 500 W/m2.
    Figure 17.  Temperature difference with Z coordinate at different volume fraction and particles for Al2O3 nanofluid diameter at 500 W/m2.

    Figure 15 presents how the temperature affected by the pressure at different φ for the nanofluid. It is noted that, as the pressure increase, the temperature decrease, that is because the increasing in pressure leads to increasing in the nanofluid velocity. Figure 16 shows the effect of porosity on temperature variation. The porosity has a significant inverse effect with the temperature, where the slight increase in porosity, lead for a significant decrease in the temperature variation.

    Figure 17 shows the effect of solid particle diameters of the porous media on the temperature variation. The same effect as porosity, particle diameters has also a significant inverse effect with the temperature, where the slight increase in the particle diameter, lead for a significant decrease in the temperature variation.

    Figure 18.  Temperature difference with Z coordinate at different volume fraction and tube radius (r) for Al2O3 nanofluid at 500 W/m2.

    The previous calculations were based on the tube radius of 3 cm, the effect of the tube radius is illustrated in Figure 18, the increase in the tube radius caused a decrease in the temperature variation as shown in the figure.

    Figures 19 and 20 present the temperature variation at different heat fluxes in Amman for different for the nanofluid. The data of irradiances were taken from the from "Photovoltaic Geographical Information System" (PVGIS) and the temperature variation was calculated depending on that data. The high values of temperature were found in summer months, July and August. But the lowest were in the winter in December and January. Thus, due to the irradiance be maximum in summer and minimum in winter. The daily peak for the curve is different from month to month, that is due to the difference of solar noon between the months.

    Figure 19.  Temperature variation T(℃) during the day time(hr) at different heat flux for different months in Amman for Al2O3 with 6% solid volume fraction.
    Figure 20.  Temperature variation T(℃) during the day time(hr) at different heat flux for different months in Amman for Al2O3 with 1.5% solid volume fraction.

    Figures 21 and 22 present the storage capacity in winter (January) and summer (August) at specified conditions and parameters shown on the figures. The conditions and the parameters can be changed to meet the required temperature. It is noted that the increase in the nanoparticles solid volume fraction will increase the temperature that can be achieved.

    Figure 21.  Temperature variation T(℃) during the day time(hr) in Amman for Al2O3 in January.
    Figure 22.  Temperature variation T(℃) during the day time(hr) in Amman for Al2O3 in August.

    In the present work, a novel model of Al2O3 water-based nanofluid embedded in a saturated porous media was investigated theoretically. The governing equations were solved numerically in the entrance region using the finite difference method and analytically in the fully developed region. From the obtained results, it is concluded that:

    ⅰ. The porous media enhance the performance of the energy storage since it has a high thermal conductivity.

    ⅱ. The use of nanofluid gave a better result than pure water in the presence of porous media, where the temperature difference for Al2O3 nanofluid with a solid volume fraction of 4.5% almost reached 40 ℃ with heat flux of 500 Watt/m2, compared to a temperature difference around 25 ℃ when using water under the same conditions.

    ⅲ. The analysis of the behave of the Nusselt number shows that through the entrance region, it starts with a maximum high value, then it sharply decreases to reach a constant value. While the value of the Nusselt number remains constant in the fully developed region, and it is calculated and found to be equal to 4.

    ⅳ. The effect of different parameters as the heat flux, the porosity of the porous media, solid particle diameter of the porous media, nanoparticles volume fraction in the nanofluid, pressure, and tube diameter was investigated and studied.

    ⅴ. It is found that the pressure, porosity, tube radius, and solid particles diameter of the porous media had an inverse relation with the temperature variation along the tube, where the increase of any them, holding the other parameters not changed, leads for a decrease in the temperature and the heat energy stored. On the other hand, the increase of heat flux, and the nanoparticles enhance the performance and increase the temperature and the stored energy.

    ⅵ. The results prove that increasing the nanoparticles solid volume fraction in the nanofluid will enhance the system's performance. For instance, 1.5% solid volume fraction for the nanofluid with heat flux of 500 Watt/m2 gave a temperature difference around 26 ℃, while it was almost 45 ℃ when the solid volume fraction for the nanofluid was 6% under the same conditions.

    ⅶ. The best results were obtained at low velocities and Prandtl number, where the low velocity means low Reynolds number.

    ⅷ. The decrease in the Prandtl number and nanofluid velocity leads to a decrease in entrance length, and thus the fully developed region achieved faster.

    ⅸ. The results for a natural heat flux was illustrated by taking Amman in Jordan as an example; it is clear that the required temperature variation can be achieved through changing the conditions and parameters such as the pressure, nanoparticles solid volume fraction in the nanofluid, the porosity of the porous media, solid particles diameter of the porous media, and the tube radius.

    As a future work, the following are suggested:

    ⅰ. Study the effect of changing the porous media, base fluid, or nanoparticles material.

    ⅱ. Study the suggested model experimentally and compare the outcomes and the results with this work that is done theoretically.

    The authors declare no conflict of interest.



    [1] Sabiha MA, Saidur R, Mekhilef S, et al. (2015) Progress and latest developments of evacuated tube solar collectors. Renewable Sustainable Energy Rev 51: 1038-1054. doi: 10.1016/j.rser.2015.07.016
    [2] Chopra K, Tyagi VV, Pandey AK, et al. (2018) Global advancement on experimental and thermal analysis of evacuated tube collector with and without heat pipe systems and possible applications. Appl Energy 228: 351-389. doi: 10.1016/j.apenergy.2018.06.067
    [3] Choi SUS (1995) Enhancing thermal conductivity of fluids with nanoparticles. American Society of Mechanical Engineers, Fluids Eng Div 231: 99-105.
    [4] Baqaie Saryazdi A, Talebi F, Armaghani T, et al. (2016) Numerical study of forced convection flow and heat transfer of a nanofluid flowing inside a straight circular pipe filled with a saturated porous medium. Eur Phys J Plus 131.
    [5] Kasaeian A, Daneshazarian R, Mahian O, et al. (2017) Nanofluid flow and heat transfer in porous media: a review of the latest developments. Int J Heat Mass Transfer 107: 778-791. doi: 10.1016/j.ijheatmasstransfer.2016.11.074
    [6] Gautam A, Saini RP (2020) A review on technical, applications and economic aspect of packed bed solar thermal energy storage system. J Energy Storage 27: 101046. doi: 10.1016/j.est.2019.101046
    [7] Nield, DA, Bejan A (2017) Convection in Porous Media. Available from: https://doi.org/10.1007/978-3-319-49562-0.
    [8] Mahmoudi Y, Hooman K, Vafai K (2019) Convective Heat Transfer in Porous Media.
    [9] Muhammad MJ, Muhammad IA, Sidik NAC, et al. (2016) Thermal performance enhancement of flat-plate and evacuated tube solar collectors using nanofluid: a review. Int Commun Heat Mass Transfer 76: 6-15. doi: 10.1016/j.icheatmasstransfer.2016.05.009
    [10] Rashidi S, Esfahani JA, Rashidi A (2017) A review on the applications of porous materials in solar energy systems. Renewable Sustainable Energy Rev 73: 1198-1210. doi: 10.1016/j.rser.2017.02.028
    [11] Ahmed OK, Mohammed ZA (2017) Influence of porous media on the performance of hybrid PV/Thermal collector. Renewable Energy 112: 378-387. doi: 10.1016/j.renene.2017.05.061
    [12] Tajik Jamal-Abad M (2017) Experimental investigation on the effect of partially metal foam inside the absorber of parabolic trough solar collector. Int J Eng 30: 281-287.
    [13] AlMasa'deh HA, Duwairi HM. Modeling of Fluid Flow and Heat Transfer inside a Saturated Porous Conduit at Constant Surface Heat Flux. Available from: https://doi.org/10.18280/mmc_b.860309.
    [14] Saedodin SAHZS, Zamzamian SAH, Nimvari ME, et al. (2017) Performance evaluation of a flat-plate solar collector filled with porous metal foam: Experimental and numerical analysis. Energy Convers Manage 153: 278-287. doi: 10.1016/j.enconman.2017.09.072
    [15] Esfe MH, Bahiraei M, Hajbarati H, et al. (2020) A comprehensive review on convective heat transfer of nanofluids in porous media: Energy-related and thermohydraulic characteristics. Appl Therm Eng, 115487. doi: 10.1016/j.applthermaleng.2020.115487
    [16] Alihosseini S, Jafari A (2020) The effect of porous medium configuration on nanofluid heat transfer. Appl Nanosci 10: 895-906. doi: 10.1007/s13204-019-01192-1
    [17] Esmaeili M, Karami M, Delfani S (2020) Performance enhancement of a direct absorption solar collector using copper oxide porous foam and nanofluid. Int J Energy Res 44: 5527-5544. doi: 10.1002/er.5305
    [18] Tay NHS, Liu M, Belusko M, et al, (2017) Review on transportable phase change material in thermal energy storage systems. Renewable Sustainable Energy Rev 75: 264-277. doi: 10.1016/j.rser.2016.10.069
    [19] Liu L, Su D, Tang Y, et al. (2016) Thermal conductivity enhancement of phase change materials for thermal energy storage: A review. Renewable Sustainable Energy Rev 62: 305-317. doi: 10.1016/j.rser.2016.04.057
    [20] Lohrasbi S, Miry SZ, Gorji-Bandpy M, et al. (2017) Performance enhancement of finned heat pipe assisted latent heat thermal energy storage system in the presence of nano-enhanced H2O as phase change material. Int J Hydrogen Energy 42: 6526-6546. doi: 10.1016/j.ijhydene.2017.01.045
    [21] Bazri S, Badruddin IA, Naghavi MS, et al. (2019) An analytical and comparative study of the charging and discharging processes in a latent heat thermal storage tank for solar water heater system. Sol Energy 185: 424-438. doi: 10.1016/j.solener.2019.04.046
    [22] Naghavi MS, Ang BC, Rahmanian B, et al. (2020) On-demand dynamic performance of a thermal battery in tankless domestic solar water heating in the tropical region. Appl Therm Eng 167: 114790. doi: 10.1016/j.applthermaleng.2019.114790
    [23] Mahdi RA, Mohammed HA, Munisamy KM, et al. (2015) Review of convection heat transfer and fluid flow in porous media with nanofluid. Renewable Sustainable Energy Rev 41: 715-734. doi: 10.1016/j.rser.2014.08.040
    [24] Khanafer K, Tavakkoli F, Vafai K, et al. (2015) A critical investigation of the anomalous behavior of molten salt-based nanofluids. Int Commun Heat Mass Transfer 69: 51-58. doi: 10.1016/j.icheatmasstransfer.2015.10.002
    [25] Corcione M (2011) Empirical correlating equations for predicting the effective thermal conductivity and dynamic viscosity of nanofluids. Energy Convers Manage 52: 789-793. doi: 10.1016/j.enconman.2010.06.072
    [26] Kandelousi MS (2014) KKL correlation for simulation of nanofluid flow and heat transfer in a permeable channel. Phys Lett A 378: 3331-3339. doi: 10.1016/j.physleta.2014.09.046
    [27] Doerr M, Frommherz M (2002) Graphite (C) Classifications, Properties & Applications.
    [28] Sanchez-Coronado J, Chung DDL (2003) Thermomechanical behavior of a graphite foam. Carbon 41: 1175-1180. doi: 10.1016/S0008-6223(03)00025-3
    [29] Vafai K (2015) Handbook of porous media. Crc Press.
    [30] Karaipekli A, Biçer A, Sarı A, et al. (2017) Thermal characteristics of expanded perlite/paraffin composite phase change material with enhanced thermal conductivity using carbon nanotubes. Energy Convers Manage 134: 373-381. doi: 10.1016/j.enconman.2016.12.053
    [31] Kim JH, Jeong E, Lee YS (2015) Preparation and characterization of graphite foams. J Ind Eng Chem 32: 21-33. doi: 10.1016/j.jiec.2015.09.003
    [32] Bianco V, Manca O, Nardini S, et al. (2015) Heat Transfer Enhancement with Nanofluids.
    [33] Oosthuizen P, Naylor D (1999) Introduction to Convective Heat Transfer Analysis. WCB/McGrawHill.
  • This article has been cited by:

    1. Aslı Akyol İnada, Samaneh Arman, Babak Safaei, A novel review on the efficiency of nanomaterials for solar energy storage systems, 2022, 55, 2352152X, 105661, 10.1016/j.est.2022.105661
    2. Muneesh Sethi, R.K. Tripathi, Birajashis Pattnaik, Sushil Kumar, Rohit Khargotra, Satish Chand, Abhishek Thakur, Recent developments in design of evacuated tube solar collectors integrated with thermal energy storage: A review, 2022, 52, 22147853, 1689, 10.1016/j.matpr.2021.11.324
    3. Berrin YILMAZ, Evren OSMA, FARKLI DUYULUR ISI DEPOLAMA MADDELERİNİN DEŞARJ PERİYODUNA ETKİLERİNİN DENEYSEL ARAŞTIRILMASI, 2023, 24, 2147-0308, 39, 10.59314/tujes.1378892
    4. Sally F. Naji, Hosham S. Anead, Khalid F. Sultan, 2024, 3092, 0094-243X, 050008, 10.1063/5.0199869
    5. Lukman Ahmed Omeiza, Muhammad Abid, Yathavan Subramanian, Anitha Dhanasekaran, Saifullah Abu Bakar, Abul Kalam Azad, Challenges, limitations, and applications of nanofluids in solar thermal collectors—a comprehensive review, 2023, 1614-7499, 10.1007/s11356-023-30656-9
  • Reader Comments
  • © 2021 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(2672) PDF downloads(130) Cited by(5)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog