Next Article in Journal
Active Learning: Encoder-Decoder-Outlayer and Vector Space Diversification Sampling
Previous Article in Journal
The Impact of Carbon Allowance Allocation Rules on Remanufacturing Decisions in a Closed-Loop Supply Chain
Previous Article in Special Issue
Amplitude and Phase Angle of Oscillatory Heat Transfer and Current Density along a Nonconducting Cylinder with Reduced Gravity and Thermal Stratification Effects
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Analysis of Entropy Generation in a Double Stage Triangular Solar Still Using CNT-Nanofluid under Double-Diffusive Natural Convection

by
Chemseddine Maatki
1,2
1
Department of Mechanical Engineering, College of Engineering, Imam Mohammad Ibn Saud Islamic University, Riyadh 11432, Saudi Arabia
2
Laboratory of Metrology and Energy System, University of Monastir, Monastir 5000, Tunisia
Mathematics 2023, 11(13), 2818; https://doi.org/10.3390/math11132818
Submission received: 23 May 2023 / Revised: 13 June 2023 / Accepted: 21 June 2023 / Published: 23 June 2023
(This article belongs to the Special Issue Advances in Computational Fluid Dynamics with Applications)

Abstract

:
The analysis of entropy generation provides valuable information for the design and optimization of thermal systems. Solar stills are used for water desalination and purification. Using renewable energies, they provide a sustainable solution for drinking water supply in remote areas and off-grid situations. This work focuses on the 3D numerical study of entropy generation in a two-stage solar still subjected to the natural double diffusion convection phenomenon in the presence of CNT nanoparticles. The effects of Rayleigh number, buoyancy ratio, and nanofluid concentration on thermal, solutal, and viscous irreversibilities and flow structure were studied. The results show that increasing the buoyancy ratio leads to an increase in thermal and solutal entropy generation. The results of this study also show that total entropy is minimal for positive volume force ratios, N, at a nanoparticle volume fraction of around 3%, and for negative N ratios, at a volume fraction of around 4%.

1. Introduction

Water scarcity appears to be an increasingly alarming prospect for the human race. The number of countries experiencing water stress is growing while the demand for water is on the rise. In light of this situation, a range of technical solutions has been developed, including seawater desalination. This practice has recently experienced significant growth due to continuous improvement of techniques and decreased costs. Desalination systems that utilize solar energy may be a revolutionary solution for regions facing water scarcity or energy shortages. Solar-powered desalination options, such as solar stills, can offer investment and operational cost savings of over 65% compared to traditional systems, as they require no energy costs or fossil fuels.
Numerous research endeavors have been undertaken to study the operational performance of solar stills, as documented in references [1,2,3,4,5,6,7,8,9,10]. Several techniques and designs have been introduced for thermal systems to improve energy efficiency [11]. The optimization of energy efficiency in thermal systems is related to the reduction of entropy generation. Within the context of solar stills, the intricate dynamics of flow, heat transfer, and mass transfer are governed by the phenomenon of double-diffusive convection, where the driving fluid is an air-vapor mixture. Consequently, the irreversibilities arising in solar stills stem from thermal, mass, and frictional influences. Various scientific investigations have explored the generation and mitigation of irreversibilities in cavities employing diverse techniques. In one such study, Dhivagar et al. [12] considered a solar still equipped with crushed gravel sand and a biomass evaporator. The researchers observed that the maximum entropy observed in this system was 61% higher than that of a conventional solar still. Ghachem et al. [13] and Borjini et al. [14] conducted comprehensive studies on three-dimensional double diffusion convection and irreversibilities production within cavities. Specifically, Ghachem et al. [13] focused their investigation on solar stills, while their subsequent study [15] explored entropy production during free convection in a cubic enclosure containing an adiabatic baffle. Maatki et al. [16] delved into examining entropy generation in the presence of discrete hot surfaces. Das et al. [17] analyzed entropy generation via free convection in various triangular cavity configurations, employing a discrete heating strategy. Their findings revealed that the configuration featuring two equally sized heating elements located in the middle and bottom halves exhibited optimal performance. Zafar et al. [18] performed an irreversibility analysis on the triple diffusive flow occurring in a porous cavity while considering the chemical reaction effect. The research demonstrated that an increase in Darcy and Lewis numbers led to a reduction in entropy generation rates. Avellaneda et al. [19] employed variational methods with weighting factors to optimize the heat transfer rate and minimize convective irreversibility production in a gas flow channel. The study showed that the optimized velocity fields led to lower entropy generations, accompanied by decreased temperatures of the heated plate. Tayebi [20] explored entropy production due to the convective flow in a porous square cavity, considering the influence of local thermal equilibrium. Chen and Jian [21] analyzed entropy generated by the convective heat transfer in a cavity involving two immiscible fluids. Surprisingly, their findings revealed that the use of two immiscible fluids generates less entropy compared to that of a single fluid. Mayeli and Sheard [22] studied the entropy generation in a shallow enclosure. The investigation showed that beyond a particular threshold of the Gay-Lussac number exceeding 0.5, the irreversibilities induced by conduction and convection manifest comparable magnitudes. Vijaybabu [23] explored the entropy production during the MDH double-diffusive convection within a porous cylinder. The author observed that the irreversibility increases with the intensity of the magnetic field when its direction aligns with the flow direction. Furthermore, the study revealed that augmenting the buoyancy ratio tends to suppress irreversibilities arising from fluid friction and the magnetic field. Kumar et al. [24] endeavored to model and numerically analyze the frictional entropy generation, as well as heat and mass flux diffusion, under the effect of a magnetic field within a saturated porous cavity containing doubly stratified fluid. Their investigation unveiled that magnetohydrodynamic forces exert a significant reduction in entropy generation resulting from viscous dissipation. Finally, Maatki et al. [25] discovered an intriguing phenomenon where intensifying the magnetic field diminishes the overall entropy generated by double-diffusive convection within a cavity.
In an intricate investigation, Munawar et al. [26] conducted a comprehensive examination of free convection occurring within a heated triangular corrugated ring, wherein a fluid was filled and augmented with a hybrid nanofluid. Astonishingly, their findings unveiled that the proposed design not only mitigates entropy production but also retains its effectiveness even when employing a high thermal conductivity hybrid nanofluid. Moreover, the study demonstrated that by incorporating corrugated walls instead of flat walls, the overall entropy is minimized, presenting a remarkable improvement in system performance. Several esteemed researchers have undertaken numerical investigations to delve into the intriguing realm of entropy generation induced by magneto-convection in the presence of nanofluids, as documented in references [27,28,29,30]. Korei et al. [28] meticulously explored a lidded cavity scenario wherein two rounded corners were filled with a nanofluid. Their study highlighted a captivating finding, i.e., that by reducing the radius of the corners, they were able to enhance heat transfer while simultaneously diminishing irreversibility. The authors further unveiled that heat transfer rate and irreversibility both exhibited a downward trend with an increase in magnetic field strength, alongside a decrease in the influence of mixed convection. Akhter et al. [29] embarked on an investigation of entropy production of a magnetoconvective flow within a porous cavity filled with a hybrid nanofluid in the presence of an external magnetic field. Their investigation disclosed intriguing correlations, revealing that the rate of heat transfers and average entropy generation intensified with an escalation in the Rayleigh number, concurrently accelerated by an augmentation in cavity porosity and nanofluid volume fraction. Sáchica et al. [30] performed meticulous numerical simulations to unravel the MHD transient mixed convection of a nanofluid in a 2D cavity having an internal hot cylinder. The authors astutely elucidated the impact of buoyancy and magnetic parameters on minimizing entropy generation within the nanofluid flow. Notably, their study highlighted the predominant influence of irreversibilities due to heat transfer in the overall entropy generation. In a distinct investigation, Aich et al. [31] meticulously explored entropy generation within a hybrid nanofluid-filled C-shaped enclosure. Chen et al. [32] undertook a comprehensive parametric study, scrutinizing entropy generation within natural double-diffusive nanofluids flow confined within a 2D cavity, specifically focusing on solar thermal systems. Their findings revealed that entropy generation intensified notably when the flow regime transitioned to turbulence. Additionally, they unveiled that the total entropy generation reached its minimum value when the buoyancy force ratio equaled 1 and exhibited a decreasing trend with an increase in nanoparticle volume fraction. Alsarraf et al. [33] conducted an insightful examination, unraveling the entropy production within turbulent forced convection of a two-phase fluid flow, accentuating the presence of various nanofluid shapes within a flat plate solar collector. Their study displayed a remarkable revelation: the case with a nanofluid concentration of ϕ = 4%, taking the shape of a blade, yielded the minimum entropy production, thus demonstrating the potential of this specific configuration. Other interesting papers related to the subject can be found in the literature [34,35,36,37].
Upon meticulous examination of the aforementioned literature review, it became apparent that a dearth of studies focusing on entropy generation within three-dimensional (3D) solar stills exists. Consequently, the present study endeavors to bridge this knowledge gap by employing the finite element method to assess entropy generation within a 3D double-stage triangular solar still. The novelty in this work lies in the study of the efficiency and entropy generation of a solar still featuring a two-stage configuration at the evaporation surface. The seawater feed zone is thus designed with two stages. The upper-level water feed plate overflows centrally, causing the water to flow from both sides towards the lower plates, thereby increasing the evaporation surface area. The surface of both feeder stages is coated with a thin film of CNT nanoparticles. Moreover, the investigation explores the influence of incorporating carbon nanotube (CNT) nanoparticles on various forms of irreversibilities. The set of equations governing the considered configuration is adeptly converted to the dimensionless form and subsequently solved using the COMSOL Multiphysics software version 5.5 (Comsol AB, Stockholm, Suède).
The outcomes of this numerical investigation are meticulously portrayed, encompassing a diverse range of variables. The depiction of the flow structure offers a detailed visualization of the intricate patterns and dynamics within the system. The temperature and concentration fields provide a profound understanding of the distribution and variation of these essential parameters throughout the domain. Furthermore, the local entropy generation fields offer valuable insights into the irreversible processes occurring within the system, shedding light on the sources and magnitude of entropy production. Additionally, the averaged entropy generations serve as crucial metrics, summarizing the overall system performance and providing a comprehensive assessment of its thermodynamic efficiency. Through the meticulous application of numerical techniques and the astute analysis of the results, a comprehensive depiction of the system’s behavior is attained. The intricate interplay between flow, temperature, concentration, and entropy generation is unraveled, offering researchers and experts a profound understanding of the system’s characteristics and enabling informed decisions and further advancements in the field.

2. Studied Configuration

When solar energy is absorbed by the solar still, temperature and concentration gradients create heat and mass fluxes that promote the evaporation of saltwater and the condensation of pure water vapor to obtain fresh water. The region of the solar still bounded by the first formed vapor layer (at the evaporation surface) and the last layer that disappears through condensation at the cold surface is characterized by an air-vapor mixture. At the evaporation surface, the vapor concentration is maximum, while at the condensation surface, the vapor concentration is minimum. This is how the phenomenon of double-diffusive convection of the air-vapor mixture appears in the solar still.
Increasing the heat exchange surface enhances the performance of the solar still. This led to the introduction of a triangular geometric design, which provides a larger condensation surface than the rectangular solar still. Additionally, a new design for the evaporation surface was achieved by adding a second stage in the center, allowing the flow of saltwater on both sides of this stage towards the primary stage, resulting in an increased evaporation surface.
The studied configuration (Figure 1) consists of a 3D triangular solar still equipped with a double stage at the bottom and filled with humid air (air-vapor mixture). A cold temperature (Tc) and low concentration (Cl) are imposed at the left and right inclined walls, where condensation occurs. On the double-stage surface, a thin film of nanofluid (CNT-water) exists (gray domain in Figure 1), and it is maintained at a hot temperature (Th) and high concentration (Ch) where the evaporation occurs. The remaining walls are considered thermally insulated and impermeable. The properties of CNT and water are presented in Table 1.

3. Mathematical Formulation

3.1. Assumptions and Governing Equations

This study delves into the intricate dynamics of the laminar, Newtonian, and incompressible flow of moist air, considering the interplay between thermal and solutal forces. Furthermore, it conjectures that the buoyancy forces operate in a negative direction along the y-axis. The properties of the humid air and the nanofluid remain constant throughout the analysis, except for the density term, which is subject to variation. To incorporate solutal and thermal buoyancy forces, the density term is approximated using the Boussinesq approximation.
The Soret and Dufour effects, while pertinent in other scenarios, are deliberately neglected in this analysis. Within the context of solar still operation, a continuous feed flow is assumed to be present to replenish evaporated water, sustain the volume fraction of the nanofluid, and curtail the potential risks associated with sedimentation.
The region occupied by the CNT water nanofluid is commonly mentioned as the “CNT water nanofluid zone” in this study. This region is characterized by carbon nanotubes (CNTs) suspended in water, which leads to unique thermophysical properties that differ from those of pure water. The CNT water nanofluid is utilized in the operation of the solar still, and its concentration is continuously maintained through a feed flow to prevent sedimentation and ensure the system’s efficient functioning.
The mathematical model employed in this study comprises five equations, namely the equation of continuity, the equation of momentum, the equation of energy governing the region of humid air, the equation of energy governing the region of nanofluid, and the equation of diffusion governing the species.
× V = 0
V t + V × V = 1 ρ P + ν Δ V + [ 1 β t T T 0 β c C C 0 ] g
T t + V × T = α 2 T
T t = α n f 2 T I n   C N T w a t e r   z o n e
C t + V × C = D 2 C
The equations are subject to further transformation to get the dimensionless form based on specific parameters. These parameters are crucial in the analysis of the system’s behavior and its physical properties. The following transformations are used to get the dimensionless governing equations, where the dimensionless numbers (Pr, Le, Ra, and N) appear.
t = t l 2 α ; x , y , z = x , y , z l ; V x , V y , V z = V x , V y , V z α / l ; P = P ρ ( α l ) ; C = C C l C h C l   and   T = T T c T h T c
× V = 0
V t + V × V = P + P r Δ V + 0 R a P r T + N C 0
T t + V × T = 2 T
T t = α n f α 2 T I n   C N T w a t e r   z o n e
C t + V × C = 1 L e 2 C
The mathematical expressions presented in Equations (6)–(10) have been transformed into dimensionless forms by utilizing certain dimensionless parameters. These equations are characterized by the appearance of four dimensionless numbers:
P r = v α , R a = g β t ( T h T c ) l 3 v α , N = β c ( C h C l ) β t ( T h T c ) a n d   L e = α D
In this study, the effective properties of the nanofluid, namely density, heat capacity, and thermal conductivity, are evaluated using Equations (12)–(14).
The effective density is given by Equation (12):
ρ n f = 1 ϕ ρ f + ϕ ρ s
The heat capacity is expressed by Equation (13):
ρ C p n f = 1 ϕ ρ C p f + ϕ ρ C p s
The effective thermal conductivity is determined by Equation (14), which is based on the Xue model.
k n f k f = 1 ϕ + 2 ϕ k s k s k f l n k s + k f 2 k f 1 ϕ + 2 ϕ k f k s k f l n k s + k f 2 k f .

3.2. Entropy Generation Equations

The thermal and solutal gradients near the active cavity walls also cause entropy generation in the system. The local entropy generation in a three-dimensional flow is given by:
S g e n = k f T 0 2 T x 2 + T y 2 + T z 2 + µ f T 0 2 V x x 2 + V y y 2 + V z z 2 + V y x + V x y 2 + V z y + V y z 2 + V x z + V z x 2 + R D C 0 C x 2 + C y 2 + C z 2 + R D T 0 T x C x + T y C y + T z C z
where C0 and T0 are the reference concentration and temperature, respectively. R is the universal gas constant.
The dimensionless local entropy generation can be expressed as:
N s = T x 2 + T y 2 + T z 2   + I s 1 2 V x x 2 + V y y 2 + V z z 2 + V y x + V x y 2 + V z y + V y z 2 + V x z + V z x 2   + I s 2 C x 2 + C y 2 + C z 2 + I s 3 T x C x + T y C y + T z C z
The dimensionless ratios I s 1 , I s 2 , and I s 3 are associated with irreversibility distribution ratios for velocity gradients, concentration gradients, and the product of concentration and temperature gradients, respectively.
The irreversibility distribution ratios are given respectively by:
I s 1 = µ f α 2 k l 2 Δ T T 0 ,   I s 2 = R D T 0 k C 0 Δ C Δ T 2 ,   I s 3 = R D k ( Δ C Δ T )

3.3. Boundary Conditions

The following are the definitions of the dimensionless boundary conditions:
Temperature: at y = 0   T = 1 , at inclined wall T = 0 ; T n = 0 on remaining walls.
Concentration: at y = 0.1   C = 1 , at inclined wall C = 0 ; C n = 0 on remaining walls.
Velocity: V x = V y = V z = 0 on all walls.
At the humid air-nanofluid interface, the continuity condition for both temperature and heat flux is expressed as follows:
T n f = R c T n i
where R c = k n f / k .
The local Nusselt and Sherwood numbers are:
N u = T y | y = 0 a n d S h = C y | y = 0.1
The averaged Nusselt and Sherwood numbers are as follows:
N u a v = 0 1 0 1 N u . y z a n d S h a v = 0 1 0 1 S h . y z

4. Numerical Method and Validation

A computational model was developed using COMSOL Multiphysics 5.5 (Tegnérgatan 23, 111 40 Stockholm, Suède) to investigate the flow structure and entropy generation in a double-stage solar still. COMSOL Multiphysics uses the finite element method. In keeping with the software’s Multiphysics concept, four coupled “application modules” are used to model the phenomenon of three-dimensional diffusive double convection: fluid flow, heat transfer, mass transfer, and chemical engineering. Each of these models corresponds to a physical phenomenon, such as fluid flow, heat or energy transfer, and vapor transport in the air-humid mixture, and also has a differential equation (Equations (1)–(3) and (5)). We specify the model region and boundary conditions, as described in Section 3.3. We have defined the lower region corresponding to the two evaporation stages as being formed by water, and a thin film of nanofluid is applied; heat transfer is assumed to be by conduction in this zone. The upper region is formed by the air-vapor mixture and is where double-diffusive convection takes place. The simulation employed a second-order upwind discretization technique and implemented a species transport model for the air-vapor mixture. We have opted to use the default quadratic Lagrange shape functions on a mesh composed of triangular finite elements for all variables as our discretization method.
The simulation converged on the chosen mesh for most parameter variations, presented below, using the solver’s default options. We also opted for transient simulations to check the stability of the steady-state results. When the steady state is established, the entropy generation related to each irreversibility cause: friction, thermal, and solute is calculated. The total entropy production of all causes of irreversibility is then also calculated.
The convergence of the solution was determined by the scaled residuals for the energy equation, which were required to be less than 10−6, and 10−5 for the other equations. This rigorous approach ensures the accuracy and reliability of the simulation results.
Figure 2 compares the 3D flow structure with the results of Kolsi et al. [34] for Pr = 0.7, Le = 0.85, N = −1, and Ra = 105. The authors studied the effect of conductive fins on the double-diffusive convection in a trapezoidal solar still, where the bottom wall is maintained at a hot temperature and the top inclined wall is kept cold. A very good concordance between the results ensures the validity of the present numerical code.

5. Results and Discussions

The results presented in this study pertain to examining intricate interactions and patterns within an air-vapor mixture that arise when simultaneous gradients of temperature and concentration exist in a double-stage water-feed solar still. The surface of both feeder stages is coated with a thin film of CNT nanoparticles. Iso-concentration and iso-temperature profiles are employed to delineate the dispersion characteristics of the vapor within the moist air mixture. Furthermore, examining local entropy production fields offers invaluable insights into the irreversible processes transpiring within the system, thereby shedding light on the origins and magnitude of entropy generation. The outcomes were derived within specified ranges of guiding parameters, encompassing Rayleigh numbers spanning from 104 to 106, buoyancy values ranging from −2 to 2, and nanofluid concentrations varying between 0 and 0.45. The Prandtl number was rigorously set to Pr = 0.7, while the Lewis number was obtained as Le = 1.2.
Figure 3 shows the particle trajectories (left side) and the irreversibilities generated by viscous effects for Ra = 105, ϕ = 0%, and different values of buoyancy ratio, N = 0, N = 1, and N = −1. The particle trajectory defining the flow structure in the solar still is characterized by four cells for the three values of buoyancy forces. When N = 0, the cells are governed by the thermal volume forces, with the cells on the right rotating in the clockwise direction and those on the left rotating counterclockwise. When N = 1, the thermal and solutal volume forces cooperate to increase the intensity of the particle velocity in the solar still. The maximum dimensionless velocity increases from 7.28 for N = 0 to 22.6 for N = 1. For N = −1, there is an equitable competition between the thermal and solute volume forces, resulting in a reversal of the direction of rotation of the cells in the still compared to the N = 0 and N = 1 cases. The particle velocity in the solar still is therefore reduced and weakened. The viscous entropy generation is related to the velocity component’s gradients. For N = 1, the maximum values of viscous entropy generation occur at the corners of the solar still due to the important friction, especially between the fluid and the wall at these regions. For N = 1, a similar field of the viscous entropy generation is encountered but with higher values. In this case, the increase of the viscous entropy is justified by the cooperative behavior of the thermal and solutal forces that lead to higher velocities and, thus, more friction. For N = −1, due to the competitive interplay between the solutal and thermal forces, the flow intensity is reduced, and lower viscous entropy is generated.
Figure 4 illustrates the iso-surfaces of temperature (left side) and iso-surfaces of the thermal entropy generation (right side) for Ra = 105, ϕ = 0%, and various buoyancy ratio values. For all the considered values of N, the iso-surfaces of temperature are parallel to the active wall with distortions at the interface separating the humid air and the nanofluid film. These distortions are due to the difference in the thermal conductivity between the humid air and the water film. For N = −1, the temperature field is characterized by a vertical stratification at the core region of the cavity, indicating the dominance of the conductive heat transfer mode. This is due to the opposition between the thermal and solutal forces that reduce the flow intensity, as discussed in Figure 3. For higher values of N, some distortions of the iso-surface of temperature are encountered at the center of the cavity due to the increase of the fluid velocities. It is also to be mentioned that at the left and right corners of the solar still, the iso-surfaces of temperature are piled up due to the existence of the collector that connects the hot and cold actives walls. The iso-surfaces of entropy generation are mainly in the central and lower regions of the cavity. This is due to the important temperature gradients in these regions. In fact, the top region is mainly delimited by two cold walls, and thus, the temperature gradients are very low. It is also to be mentioned that the maximum values of the thermal entropy generation are at the left and right corners. In fact, these regions exist close to the connections between the hot and cold walls, where the density variations are very important and cause higher temperature gradients.
Figure 5 depicts the effects of the buoyancy ratio on the iso-surfaces of concentration and solutal entropy generation for Ra = 105 and ϕ = 0%. Due to the imposed constant concentration boundary condition, the iso-surfaces of concentration are parallel to the active wall. In the central region of the cavity, the concentration field is characterized by vertical stratification. The concentration gradients are more important close to the cold and hot wall walls, especially close to the corner of the cavity. This fact leads to the generation of higher solutal entropy generation values at these regions, as presented in the right column of the figure. It is also to be mentioned that the concentration gradients are more important for N = 1 and N = 0 compared to N = −1, and thus, lower solutal entropy generation occurs for N = −1. This is due to the equilibrium between the solutal and thermal buoyancy forces at this specific value of the buoyancy ratio.
Figure 6 illustrates the iso-surfaces of total entropy generation for Ra = 105, ϕ = 0%, and different values of buoyancy ratio. As described in Figure 4 and Figure 5, the solutal and thermal entropy generations are concentrated close to the active walls, and their lowest values occur for N = −1. Obviously, the total entropy generation follows the same behavior.
Figure 7 presents the effect of the nanoparticles’ volume fraction on the flow structures (top rows) and frictional entropy generation (bottom row) for Ra = 106 and different values of the buoyancy ratio. For N = 1, the flow structure is characterized by four vortexes, and compared to Ra = 105 (Figure 3), it is noticed that due to the intensification of the flow, the particle trajectories are no more symmetrical. The increase of the nanoparticles’ volume fraction does not cause a considerable change in the flow structure and leads only to an increase in the velocity magnitude. This increase in the velocity causes an increase in the viscous effect entropy generation, as presented in the second row of the figure. It also noticed that compared to Ra = 105 (Figure 3), the viscous entropy generation becomes distributed on the whole domain, but still, the maximum values occur at the corners. For N = −1, the flow structure is again symmetrical due to the reduction of the flow intensity caused by the concurrence between the solutal and thermal forces. A considerable reduction in the velocity magnitude is encountered compared to N = 1 cases, and by increasing the nanoparticles’ volume fraction, it becomes lower. This reduction in the velocity causes an obvious reduction in the frictionally generated entropy generation.
Figure 8 is plotted to present the effect of the nanoparticles’ volume fraction on the temperature field and thermal entropy generation for Ra = 106 and different values of the buoyancy ratio. Similarly, Figure 9 depicts the same effect on the concentration field and solutal entropy generation. Due to the fixed Lewis number (Le = 0.85) corresponding to the air vapor mixture, which is close to 1, the temperature and concentration fields are similar and the same for the thermal and solutal entropy generations. For N = −1, the temperature (concentration) field is characterized by a central vertical stratification. For N = 1, the iso-surfaces of temperature (concentration) at the central region become distorted due to the intensification of the convective flow, as can be noticed for the velocity magnitude in Figure 7. This intensification is due to cooperative solutal and thermal buoyancy forces for positive N values. The addition of nanoparticles enhances the thermos-physical properties of the fluid and leads to more intense flow, and thus, the iso-surface distortions become more pronounced. The thermal and solutal entropy generation is mainly concentrated close to the active walls for N = −1. For N = 1, due to the flow intensification, the hot and highly concentrated fluid is carried away toward the central region. Thus, the thermal and solutal entropy generations become distributed in the whole domain. The increase in the nanoparticle volume fraction leads to an increase in the thermal and solutal entropy generations.
Figure 10 presents the iso-surfaces of total entropy generation for Ra = 106. Similar to the viscous and thermal entropy generations, the iso-surfaces of total entropy generation are mainly concentrated close to the active walls. It is also to be mentioned that for N = −1, the values of the produced entropy generation are the lowest due to the equilibrium between the thermal and solutal forces. In addition, the increase in the CNT volume fraction leads to an increase in the total entropy generation.
Figure 11 presents the variations of the total entropy generation versus the nanoparticles volume fraction for various buoyancy ratios and Rayleigh number values. It is obvious that the increase of Ra increases the produced irreversibilities. For Ra = 104 (Figure 11a), the heat transfer regime is mainly conductive, and the buoyancy ratio is not very effective; thus, the curves are grouped into two groups; positive N values and negative N values (including N = 0). For Ra = 105 (Figure 11b), the convective effects become more important, and the curves of total entropy generation become no longer superposed. It is also to be mentioned that for negative N values, the maximum entropy generation occurs for ϕ = 0.03, and the minimum occurs for ϕ = 0.04. For positive N values, the maximum entropy generation occurs for ϕ = 0.02, and the minimum occurs for ϕ = 0.03. This result is very interesting when the aim is to enhance the heat and mass transfer while producing less irreversibility leading to higher performances and lower energy losses.

6. Conclusions

The present study deals with the investigation of the effects of using CNT-nanofluid on the entropy generation in a 3D double-stage solar still. A numerical model based on the FEM is established using COMSOL Multiphysics 5.5 (Comsol AB, Stockholm, Suède) to solve the dimensionless governing equations. The effects of the governing parameters such as Rayleigh number, buoyancy ratio, and CNT volume fraction on the flow structure, temperature and concentration field, and entropy generations are investigated. The main findings can be summarized as follows:
-
The increase of the buoyancy ratio, Rayleigh number, and CNT volume fraction leads to the intensification of the flow;
-
Due to the equal competitively between the solutal and thermal forces, the lower fluid velocity magnitudes and viscous entropy generation occur for N = −1;
-
The maximum values of viscous entropy generation occur close to the corners of the solar still;
-
The increase of the buoyancy ratio leads to the increase of the thermal and solutal entropy generations;
-
The volume fraction of CNT added in a film bonded to the two water supply stages of a Solar still improves heat transfer performance at the evaporation surface, which can be used to optimize the entropy generated;
-
Depending on the values of Ra and N, an optimum value occurs on all the variations of the entropy generations versus the volume fraction. For our two-stage water supply design, the total entropy is minimal for positive volume force ratios, N, at a nanoparticle volume fraction of around 3%, and for negative N ratios, at a volume fraction of around 4%.

Funding

This research was supported by the Deanship of Scientific Research, Imam Mohammad Ibn Saud Islamic University (IMSIU), Saudi Arabia, Grant No. (221414002).

Data Availability Statement

Not applicable.

Acknowledgments

The authors extend their appreciation to the Deanship of Scientific Research, Imam Mohammad Ibn Saud Islamic University (IMSIU), Saudi Arabia, for funding this research work through Grant No. (221414002).

Conflicts of Interest

The author declares no conflict of interest.

Nomenclature

CDimensionless concentration
DMass diffusivity (m2/s)
GrGrashof number
kThermal conductivity (W/m·K)
LeLewis number
NBuoyancy ratio
NsDimensionless local entropy generation
NuNusselt number
PrPrandtl number
RaRayleigh number
S’genLocal entropy generation
ShSherwood number
tDimensionless time (=)
TDimensionless temperature
VDimensionless velocity vector (=)
Greek symbols
α Thermal diffusivity (m2/s)
β t Coefficient of thermal expansion (K−1)
β c Coefficient of solutal expansion (K−1)
µ Dynamic viscosity (kg/m.s)
ϕ Nanoparticles volume fraction
Subscripts
x, y, zCartesian coordinates
hhot, high
cCold
lLow
avAverage
nfNanofluid
sNanoparticle
0Reference

References

  1. Rahbar, N.; Asadi, A.; Fotouhi-Bafghi, E. Performance evaluation of two solar stills of different geometries: Tubular versus triangular: Experimental study, numerical simulation, and second law analysis. Desalination 2018, 443, 44–55. [Google Scholar] [CrossRef]
  2. Al-Madhhachi, H.; Smaisim, G.F. Experimental and numerical investigations with environmental impacts of affordable square pyramid solar still. Sol. Energy 2021, 216, 303–314. [Google Scholar] [CrossRef]
  3. Moussa, R.R.; Hatem, T.M. Estimating the potential of desalinate seawater using solar glass pyramid (SGP) in hot arid zones. Clean. Eng. Technol. 2021, 4, 100189. [Google Scholar] [CrossRef]
  4. Dev, R.; Tiwari, G. Characteristic equation of a passive solar still. Desalination 2009, 245, 246–265. [Google Scholar] [CrossRef]
  5. Kabeel, A.; Abdelgaied, M. Enhancement of pyramid-shaped solar stills performance using a high thermal conductivity absorber plate and cooling the glass cover. Renew. Energy 2020, 146, 769–775. [Google Scholar] [CrossRef]
  6. Xiong, J.; Xie, G.; Zheng, H. Experimental and numerical study on a new multi-effect solar still with enhanced condensation surface. Energy Convers. Manag. 2013, 73, 176–185. [Google Scholar] [CrossRef]
  7. Alvarado-Juáreza, R.; Álvarez, G.; Xamána, J.; Hernández-Lópe, I. Numerical study of conjugate heat and mass transfer in a solar still device. Desalination 2013, 325, 84–94. [Google Scholar] [CrossRef]
  8. Alvarado-Juárez, R.; Xamána, J.; Álvarez, G.; Hernández-López, I. Numerical study of heat and mass transfer in a solar still device: Effect of the glass cover. Desalination 2015, 359, 200–211. [Google Scholar] [CrossRef]
  9. Rahman, M.; Öztop, H.F.; Ahsan, A.; Kalam, M.; Varol, Y. Double-diffusive natural convection in a triangular solar collector. Int. Commun. Heat Mass Transf. 2012, 39, 264–269. [Google Scholar] [CrossRef]
  10. Maatki, C. Heat transfer enhancement using CNT-water nanofluids and two stages of seawater supply in the triangular solar still. Case Stud. Therm. Eng. 2022, 30, 101753. [Google Scholar] [CrossRef]
  11. Delouei, A.A.; Sajjadi, H.; Ahmadi, G. Ultrasonic Vibration Technology to Improve the Thermal Performance of CPU Water-Cooling Systems: Experimental Investigation. Water 2022, 14, 4000. [Google Scholar] [CrossRef]
  12. Dhivagar, R.; Deepanraj, B.; Mohanraj, M.; Ong, H.C. Second law based thermodynamic analysis of crushed gravel sand and biomass evaporator assisted solar still. Sustain. Energy Technol. Assess. 2022, 52, 102160. [Google Scholar] [CrossRef]
  13. Ghachem, K.; Kolsi, L.; Mâatki, C.; Hussein, A.K.; Borjini, M.N. Numerical simulation of three-dimensional double diffusive free convection flow and irreversibility studies in a solar distiller. Int. Commun. Heat Mass Transf. 2012, 39, 869–876. [Google Scholar] [CrossRef]
  14. Borjini, N.; Kolsi, L.; Ghachem, K.; Maatki, C.; Aissia, H.B. Entropy generation of double diffusive natural convection in a three dimensional differentially heated enclosure. Int. J. Eng. 2014, 27, 215–226. [Google Scholar]
  15. Ghachem, K.; Hassen, W.; Maatki, C.; Kolsi, L.; Al-Rashed, A.; Borjini, M. Numerical simulation of 3D natural convection and entropy generation in a cubic cavity equipped with an adiabatic baffle. Int. J. Heat Technol. 2018, 36, 1047–1054. [Google Scholar] [CrossRef] [Green Version]
  16. Maatki, C.; Ghachem, K.; Almeshaal, M.A.; Ben Khedher, N.; Kolsi, L. Numerical Investigation of Rayleigh–Benard Natural Convection and Entropy Generation in a Cubic Cavity With Discrete Heat Sources. J. Therm. Sci. Eng. Appl. 2021, 13, 051023. [Google Scholar] [CrossRef]
  17. Das, D.; Lukose, L.; Basak, T. Role of multiple discrete heaters to minimize entropy generation during natural convection in fluid filled square and triangular enclosures. Int. J. Heat Mass Transf. 2018, 127, 1290–1312. [Google Scholar] [CrossRef]
  18. Khan, Z.H.; Khan, W.A.; Tang, J.; Sheremet, M.A. Entropy generation analysis of triple diffusive flow past a horizontal plate in porous medium. Chem. Eng. Sci. 2020, 228, 115980. [Google Scholar] [CrossRef]
  19. Avellaneda, J.; Bataille, F.; Toutant, A.; Flamant, G. Variational entropy generation minimization of a channel flow: Convective heat transfer in a gas flow. Int. J. Heat Mass Transf. 2020, 160, 120168. [Google Scholar] [CrossRef]
  20. Tayebi, T. Analysis of the local non-equilibria on the heat transfer and entropy generation during thermal natural convection in a non-Darcy porous medium. Int. Commun. Heat Mass Transf. 2022, 135, 106133. [Google Scholar] [CrossRef]
  21. Chen, X.; Jian, Y. Entropy generation minimization analysis of two immiscible fluids. Int. J. Therm. Sci. 2022, 171, 107210. [Google Scholar] [CrossRef]
  22. Mayeli, P.; Sheard, G.J. An entropy generation analysis of horizontal convection under the centrifugal buoyancy approximation. Int. Commun. Heat Mass Transf. 2022, 133, 105923. [Google Scholar] [CrossRef]
  23. Vijaybabu, T.R. Influence of porous circular cylinder on MHD double-diffusive natural convection and entropy generation. Int. J. Mech. Sci. 2021, 206, 106625. [Google Scholar] [CrossRef]
  24. Kumar, V.; Murthy, S.K.; Kumar, B.R. Multi-force effect on fluid flow, heat and mass transfer, and entropy generation in a stratified fluid-saturated porous enclosure. Math. Comput. Simul. 2023, 203, 328–367. [Google Scholar] [CrossRef]
  25. Maatki, C.; Hassen, W.; Kolsi, L.; Alshammari, N.; Borjini, M.N.; Ben Aissia, H. 3-D Numerical Study of Hydromagnetic Double Diffusive Natural Convection and Entropy Generation in Cubic Cavity. J. Appl. Fluid Mech. 2016, 9, 1915–1925. [Google Scholar] [CrossRef]
  26. Munawar, S.; Saleem, N.; Khan, W.A. Entropy generation minimization E.G.M. analysis of free convective hybrid nanofluid flow in a corrugated triangular annulus with a central triangular heater. Chin. J. Phys. 2022, 75, 38–54. [Google Scholar] [CrossRef]
  27. Al-Rashed, A.A.; Kolsi, L.; Kalidasan, K.; Maatki, C.; Borjini, M.N.; Aichouni, M.; Kanna, P.R. Effect of magnetic field inclination on magneto-convective induced irreversibilities in a CNT-water nanofluid filled cubic cavity. Front. Heat Mass Transf. FHMT 2017, 8, 31. [Google Scholar] [CrossRef] [Green Version]
  28. Korei, Z.; Benissaad, S.; Berrahil, F.; Filali, A. MHD mixed convection and irreversibility analysis of hybrid nanofluids in a partially heated lid-driven cavity chamfered from the bottom side. Int. Commun. Heat Mass Transf. 2022, 132, 105895. [Google Scholar] [CrossRef]
  29. Akhter, R.; Ali, M.M.; Alim, M.A. Entropy generation due to hydromagnetic buoyancy-driven hybrid-nanofluid flow in partially heated porous cavity containing heat conductive obstacle. Alex. Eng. J. 2023, 62, 17–45. [Google Scholar] [CrossRef]
  30. Vijaybabu, T.R.; Karthik Babu, N.B. Magnetohydrodynamic mixed convection and entropy generation analysis of Al2O3-water nanofluid past a confined circular cylinder. Int. J. Mech. Sci. 2022, 230, 107542. [Google Scholar]
  31. Aich, W.; Chaabane, R.; Öztop, H.; Almeshaal, M.A.; Maatki, C.; Kahouli, O.; Kolsi, L. Numerical simulation of buoyancy-induced heat transfer and entropy generation in 3D C-shaped cavity filled with CNT–Al2O3/water hybrid nanofluid. Int. J. Nonlinear Sci. Numer. Simul. 2022. [Google Scholar] [CrossRef]
  32. Chen, S.; Yang, B.; Xiao, X.; Zheng, C. Analysis of entropy generation in double-diffusive natural convection of nanofluid. Int. J. Heat Mass Transf. 2015, 87, 447–463. [Google Scholar] [CrossRef]
  33. Alsarraf, J.; Shahsavar, A.; Mahani, R.B.; Talebizadehsardari, P. Turbulent forced convection and entropy production of a nanofluid in a solar collector considering various shapes for nanoparticles. Int. Commun. Heat Mass Transf. 2020, 117, 104804. [Google Scholar] [CrossRef]
  34. Kolsi, L.; Ghachem, K.; Larguech, S.; AlNemer, G. Numerical Investigation of the Double Diffusive Convection in 3D Trapezoidal Solar Still Equipped with Conductive Fins. Mathematics 2022, 10, 2115. [Google Scholar] [CrossRef]
  35. Kolsi, L.; Hussain, S.; Ghachem, K.; Jamal, M.; Maatki, C. Double Diffusive Natural Convection in a Square Cavity Filled with a Porous Media and a Power Law Fluid Separated by a Wavy Interface. Mathematics 2022, 10, 1060. [Google Scholar] [CrossRef]
  36. Saleem, K.B.; Ghachem, K.; Koufi, L.; Kolsi, L. Analysis of Double-diffusive natural convection in a solar distiller embedded with PCM and cooled with external water stream. J. Taiwan Inst. Chem. Eng. 2021, 126, 67–79. [Google Scholar] [CrossRef]
  37. Saleem, K.B.; Koufi, L.; Alshara, A.K.; Kolsi, L. Double-diffusive natural convection in a solar distiller with external fluid stream cooling. Int. J. Mech. Sci. 2020, 181, 105728. [Google Scholar] [CrossRef]
Figure 1. Studied configuration: (a) 3D geometry and (b) central section with boundary conditions.
Figure 1. Studied configuration: (a) 3D geometry and (b) central section with boundary conditions.
Mathematics 11 02818 g001
Figure 2. Comparison of the 3D flow structure with the results of Kolsi et al. [34] for Pr = 0.7, Le = 0.85, N = −1, and Ra = 105.
Figure 2. Comparison of the 3D flow structure with the results of Kolsi et al. [34] for Pr = 0.7, Le = 0.85, N = −1, and Ra = 105.
Mathematics 11 02818 g002
Figure 3. Particle trajectories and Viscous effects irreversibility for Ra = 105 and ϕ = 0%.
Figure 3. Particle trajectories and Viscous effects irreversibility for Ra = 105 and ϕ = 0%.
Mathematics 11 02818 g003
Figure 4. Temperature Iso-surfaces and thermal effects irreversibility for Ra = 105 and ϕ = 0%.
Figure 4. Temperature Iso-surfaces and thermal effects irreversibility for Ra = 105 and ϕ = 0%.
Mathematics 11 02818 g004
Figure 5. Concentration Iso-surfaces and solutal irreversibility for Ra = 105 and ϕ = 0%.
Figure 5. Concentration Iso-surfaces and solutal irreversibility for Ra = 105 and ϕ = 0%.
Mathematics 11 02818 g005
Figure 6. Local total entropy generation, Ra = 105 and ϕ = 0%.
Figure 6. Local total entropy generation, Ra = 105 and ϕ = 0%.
Mathematics 11 02818 g006
Figure 7. Particle trajectories and Viscous irreversibility, Ra = 106.
Figure 7. Particle trajectories and Viscous irreversibility, Ra = 106.
Mathematics 11 02818 g007
Figure 8. Temperature Iso-surfaces (Top) and thermal effects irreversibility (Bottom) for Ra = 106.
Figure 8. Temperature Iso-surfaces (Top) and thermal effects irreversibility (Bottom) for Ra = 106.
Mathematics 11 02818 g008
Figure 9. Concentration Iso-surfaces and solutal effects irreversibility, Ra = 106.
Figure 9. Concentration Iso-surfaces and solutal effects irreversibility, Ra = 106.
Mathematics 11 02818 g009
Figure 10. Local total entropy generation, Ra = 106.
Figure 10. Local total entropy generation, Ra = 106.
Mathematics 11 02818 g010
Figure 11. Effect of the nanoparticles volume fraction on the total entropy generation for various N and Ra values.
Figure 11. Effect of the nanoparticles volume fraction on the total entropy generation for various N and Ra values.
Mathematics 11 02818 g011
Table 1. Thermophysical Properties of CNT and water [31].
Table 1. Thermophysical Properties of CNT and water [31].
Physical PropertiesCNT-NanoparticlesWater
ρ (kg/m3)2600997.1
Cp (kJ/kg·K)0.4254.179
k (kW/m·K)6.6000.613 × 10−3
β (K−1)1.6 × 10−621 × 10−5
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Maatki, C. Numerical Analysis of Entropy Generation in a Double Stage Triangular Solar Still Using CNT-Nanofluid under Double-Diffusive Natural Convection. Mathematics 2023, 11, 2818. https://doi.org/10.3390/math11132818

AMA Style

Maatki C. Numerical Analysis of Entropy Generation in a Double Stage Triangular Solar Still Using CNT-Nanofluid under Double-Diffusive Natural Convection. Mathematics. 2023; 11(13):2818. https://doi.org/10.3390/math11132818

Chicago/Turabian Style

Maatki, Chemseddine. 2023. "Numerical Analysis of Entropy Generation in a Double Stage Triangular Solar Still Using CNT-Nanofluid under Double-Diffusive Natural Convection" Mathematics 11, no. 13: 2818. https://doi.org/10.3390/math11132818

APA Style

Maatki, C. (2023). Numerical Analysis of Entropy Generation in a Double Stage Triangular Solar Still Using CNT-Nanofluid under Double-Diffusive Natural Convection. Mathematics, 11(13), 2818. https://doi.org/10.3390/math11132818

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop