Next Article in Journal
Integrated Decision-Making Method for Heterogeneous Attributes Based on Probabilistic Linguistic Cross-Entropy and Priority Relations
Next Article in Special Issue
Impacts of Uniform Magnetic Field and Internal Heated Vertical Plate on Ferrofluid Free Convection and Entropy Generation in a Square Chamber
Previous Article in Journal
Analysis of Controlled Rabi Flopping in a Double Rephasing Photon Echo Scheme for Quantum Memories
Previous Article in Special Issue
Aspects of Chemical Entropy Generation in Flow of Casson Nanofluid between Radiative Stretching Disks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Investigation of Forced Convection Enhancement and Entropy Generation of Nanofluid Flow through a Corrugated Minichannel Filled with a Porous Media

1
School of Mechanical Engineering, Iran University of Science and Technology, Tehran 16846-13114, Iran
2
Department of Microelectronics, Delft University of Technology, 2628 CD Delft, The Netherlands
*
Author to whom correspondence should be addressed.
Submission received: 19 August 2020 / Revised: 5 September 2020 / Accepted: 7 September 2020 / Published: 9 September 2020
(This article belongs to the Special Issue Thermal Radiation and Entropy Analysis)

Abstract

:
Corrugating channel wall is considered to be an efficient procedure for achieving improved heat transfer. Further enhancement can be obtained through the utilization of nanofluids and porous media with high thermal conductivity. This paper presents the effect of geometrical parameters for the determination of an appropriate configuration. Furthermore, the optimization of forced convective heat transfer and fluid/nanofluid flow through a sinusoidal wavy-channel inside a porous medium is performed through the optimization of entropy generation. The fluid flow in porous media is considered to be laminar and Darcy–Brinkman–Forchheimer model has been utilized. The obtained results were compared with the corresponding numerical data in order to ensure the accuracy and reliability of the numerical procedure. As a result, increasing the Darcy number leads to the increased portion of thermal entropy generation as well as the decreased portion of frictional entropy generation in all configurations. Moreover, configuration with wavelength of 10 mm, amplitude of 0.5 mm and phase shift of 60° was selected as an optimum geometry for further investigations on the addition of nanoparticles. Additionally, increasing trend of average Nusselt number and friction factor, besides the decreasing trend of performance evaluation criteria (PEC) index, were inferred by increasing the volume fraction of the nanofluid (Al2O3 and CuO).

1. Introduction

Achieving improved heat transfer in tubes and channels are of great significance, especially in industrial and technological applications from the energy saving viewpoint. Hence, numerous studies have been performed on evaluating the flow field and heat transfer in diverse configurations [1]. Among the proposed procedures for ameliorating the thermal performance in various systems, the utilization of wavy-wall channels is considered to be an efficient approach with the application in heat exchangers, heat sinks, and solar collectors. This issue has been addressed in several investigations [2,3,4,5,6,7]. Additionally, the heat transfer in wavy-wall channels have been studied extensively considering their various applications [8,9,10,11,12,13,14,15,16]. In an experimental study conducted by Rush et al. [17], the flow field and heat transfer characteristics were evaluated in corrugated wavy-wall minichannels. They introduced geometrical features of the channel and the Reynolds numbers as the most influential parameters on the local Nusselt number and, consequently, the local heat transfer. In a numerical investigation, Wang and Chen [18] investigated the heat transfer and forced convection of nanofluid flow (in laminar flow regime) within sinusoidal wavy-wall minichannels. They stated the effectiveness of the wavelength, amplitude ratio, and Reynolds number in improving the heat transfer in these channels. Zhang and Che examined the influence of corrugation profile on flow field and heat transfer in the cross-corrugated walls [19]. Their results declared the higher Nusselt number and friction factor in trapezoidal channels compared to the elliptical ones. A numerical study on the one-sided corrugated plate (in turbulent flow regime) was performed by Naphon and Kornkumjayrit [20]. According to their results, the superiority of the studied geometry in improving the heat transfer was concluded. In another numerical investigation, Khoshvaght-Aliabadi [21] studied the fluid flow and heat transfer in sinusoidal-corrugated channels with the addition of nanofluid. Regarding their results, Nusselt number and friction factor high were affected by channel height and wave amplitude, respectively.
The use of porous media in industrial equipment has been considered to be a practical procedure for ameliorating the heat transfer, which has attracted the attentions of many researchers [22,23,24,25,26,27,28]. The partial employment of the porous media in the core of a sinusoidal channel was explored by Akbarzadeh and Maghrebi [29]. They assessed the hydrothermal performance of different geometries and metal foam characteristics on local thermal non-equilibrium method. The heat transfer in a counter flow sinusoidal parallel-plate heat exchanger was numerically examined by Arasteh et al. [30]. The divergent sections of the desired heat exchanger were filled partially with the metal foam porous media. They reported an enhanced heat transfer rate, efficiency and overall heat transfer coefficient up to 19.2%. In recent years, the simultaneous application of porous media and nanofluids has proved to be effective in improving the heat transfer performance. In numerical and analytical investigations, Hajipour and Dehkordi [31] studied the fully-developed mixed-convective heat transfer of nanofluid flow inside the parallel-plate channel. The studied channel was filled partially with porous medium. In addition, the momentum and thermal energy equations were characterized by the corresponding Darcy–Brinkman–Forchheimer and different viscous dissipation models. The simultaneous application of nanofluid and porous medium in throttle area of the sinusoidal channel was explored by Nazari and Toghraie [32]. They utilized local thermal equilibrium model for the evaluation of fluid flow and heat transfer in the desired system. Their outcomes declared the higher improvement of heat transfer parameters in convex areas, rather than the concave areas of the throttle.
The addition of nanofluid to the fluid flow of channels has been the subject of several investigations [33,34,35,36,37,38,39]. In a numerical survey by Heidary and Kermani [40], the heat transfer and nanofluids flow within a corrugated wavy-wall channel were investigated. They stated the enhanced Nusselt number by increasing the Reynolds number, nanoparticle volume fraction, and wave amplitude. The influence of groove angel in a corrugated duct with S i O 2 /water nanofluid on heat transfer and flow properties was assessed by Sadripour [41]. Their results revealed the high dependency of hydraulic and thermal efficiencies of the desired exchanger to the variation of grooves angles. In addition, the best performance evaluation criteria (PEC) index corresponded to the configuration with the duct angle of 30°, volume fraction of 4%, and nanoparticle size of 40   nm . Ho et al. [42] carried out a study on the forced convection heat transfer of Al2O3-water nanofluid. They stated the higher effectiveness of the nanofluid with volumetric concentration of 1%, rather than 2%, according to the higher variation of dynamic viscosity by temperature changes. They also reported 70% enhancement in convective heat transfer coefficient by the utilization of Al2O3-water nanofluid with the concentration of 1%. The numerical study of CuO nanofluid in a micro-tube was performed by Akbari et al. [43] by the use of a two-phase model. As their results depicted, the addition of nanoparticles would result in a reduction of Nusselt number and friction factor. The effect of several variables, including mass flow rate of the nanofluid, concentration of the nanoparticle, as well as geometrical parameters in wavy minichannels were examined by Aliabadi and Sahamiyan [44]. In their study, they evaluated the heat transfer and pumping power for the Al2O3-water flow under different conditions.
The other advantageous parameter that is helpful for achieving more comprehensive understanding of the quality of available energy and optimal thermodynamic conditions is known as entropy generation [45,46,47,48,49,50,51]. The minimization of entropy generation is proposed as a solution in several researches in order to optimize the thermodynamic systems [52,53,54,55,56,57]. Similar to the conducted researches on the field of evaluating the heat transfer in wavy-wall channels, the entropy generation can be investigated in both natural and forced convections. Sheremet et al. examined the entropy generation of natural convection heat transfer [58]. For this purpose, they utilized computational fluid dynamics (CFD) procedure. Their studied system consisted of Cu/water nanofluid in a cavity with corrugated perpendicular wall. In a numerical investigation by Cho et al. [59], the finite volume method was utilized for the simulation of entropy generation and natural convective heat transfer in a horizontal enclosure with wavy walls. They also investigated the effect of different metal oxide nanofluids. Different indexes, such as heat transfer, pressure drop, and entropy generation in the heat exchanger with sinusoidal wavy-wall and a porous insert with the nanofluid flow were examined by Akbarzadeh et al. [60] The calculation of entropy generation in a wavy heat exchanger with nanofluid flow was performed by Esfahani et al. [61]. They used two-dimensional (2D) simulation by ANSYS-FLUENT software for this aim. In a study by Dormohammadi et al. [62], entropy generation minimization technique was used for optimizing the mixed convective heat transfer in a wavy channel with nanofluid flow. Siavashi et al. performed the numerical investigation of the flow characteristics, heat transfer, and entropy generation in the annular pipe that were filled fully or partially with the porous media and nanofluid flow [63]. They utilized two-phase mixture model for the investigations and declared the dependency of the performance and entropy generation to the configuration parameters, nanoparticles concentration, and Reynolds number. Moreover, they disclosed the existence of an optimum thickness of the porous media for each nanofluid flow in a porous medium at a specific Reynolds number according to the thermodynamics second law.
Based on the available literature, the impact of simultaneous application of wavy-wall and porous medium, as well as their individual usage on heat transfer and fluid flow have been investigated in several studies. However, to the best of author’s knowledge, only a few studies concerning the laminar forced convective heat transfer of nanofluids flow as well as the calculation of entropy generation in a wavy-wall minichannel with porous media has been conducted. Therefore, in the current study, it is aimed to conduct a numerical investigation to evaluate the heat transfer, nanofluid flow, and entropy generation in a sinusoidal-wavy minichannel with porous medium. Moreover, the dependency of the average Nusselt number, friction factor, PEC index, and thermal/frictional entropy generation on geometrical characteristics and porosity would be evaluated. Furthermore, the impact of adding nanofluids (Al2O3 and CuO) to the pure fluid and volume fraction is also performed.

2. Materials and Methods

2.1. Geometrical Specifications of Physical Model

The schematic representation of the studied arrangement is depicted in Figure 1, which shows the system consisting of three main sections; entrance, test and exit sections. The entrance section with the length of 150 mm was utilized in order to ensure the fully developed fluid flow in the test section. This section was followed by the test section with the length of 60 mm, which consisted of a two-dimensional minichannel with the opening height of D = 5   mm and two sinusoidal-wavy walls. Finally, the exit section with the length of 24 mm was implemented in order to avoid the development of adverse pressure according to the subsequent influence of the heat transfer and fluid flow properties in computational domain. In this study, the effect of wave amplitude, wave length, and phase shift were evaluated. For investigating the influence of wave amplitude, three constant values of a = 1.5 ,   1 , and 0.5 mm, as well as two variable values (with decreasing and increasing trend) were considered. In addition, the similar procedure was considered for the wavelength; three constant values of λ = 10 , 12, and 15 mm and two variable values (with decreasing and increasing trend). Moreover, the effect of phase shift was evaluated in five values of φ = 0 ° , 30 ° , 60 ° , 90 ° , and 180 ° among lower and upper wavy-walls. It should be noted that the outlet boundary conditions are pressure outlet.

2.2. Governing Equations

The mathematical formulations include mass, momentum, and energy conservation equations in a 2D steady-state flow, as below [62]:
u x + v y = 0
1 ε 2 ( u u x + v u y ) = 1 ρ f p x + ϑ m ε ( 2 u x 2 + 2 u y 2 ) ϑ f   u K ϑ m C d | u | K u
1 ε 2 ( u v x + v v y ) = 1 ρ f p y + ϑ fm ε ( 2 v x 2 + 2 v y 2 ) ϑ m   v K ϑ m C d | v | K v
u T x + v T y = α m ( 2 T x 2 + 2 T y 2 )
For dimensionless form of equations, some variables are presented, as below.
  X = x D h ,   Y = y D h ,   U = u U in ,   V = v U in ,   θ = T T in T w T in α m = k m ( ρ C p ) m ,   P = p ρ f U i n 2 ,   C d = 1.75 150 ε 3 2 ,   T b = 0 D h u T d y 0 D h u d y
Re = ρ f U in D h μ f ,   Nu = D h ( T y ) wall T wall T b ,   Da = K D h 2 ,   Pr = ϑ f α f
In Equation (6), Re ,   Nu ,   Da , and Pr represent, Reynolds, Nusselt, Darcy, and the Prandtl numbers, respectively. In order to calculate the effective thermal conductivity of the applied fluid (based on the volumetric linear models) in porous media, the following equations were utilized [64]:
k m = ( 1 φ ) k f + φ k np
k e f f = ( 1 ε ) k p + ε k m
Additionally, the hydraulic diameter of corrugated wavy-wall has been presented as:
D h = 2 D + 2 a
By substituting non-dimensional parameters into Equations (1)–(4), non-dimensional form of equations will be achieved.
U X + V Y = 0
1 ε 2 ( U U X + V U Y ) = P X + 1 Re . ε ( 2 U X 2 + 2 U Y 2 ) U Re   Da ε C d | U | Da U
1 ε 2 ( U V X + V V Y ) = P Y + 1 Re . ε ( 2 V X 2 + 2 V Y 2 ) V Re   Da ε C d | V | Da V
( U θ X + V θ Y ) = 1 Re . Pr ( ρ c p ) f ( ρ c p ) m ( 2 θ X 2 + 2 θ Y 2 )

2.3. Thermo-Physical Properties of Nanofluids

Extensive investigations have been conducted for the determination of the physical characteristics of nanofluids. In this study, the desired nanofluid was prepared by the dispersion of copper and aluminum oxide particles in water (base fluid). The physical characteristics of the desired nanofluid can be calculated through the following equations:
ρ m = ( 1 ϕ ) ρ f + ϕ ρ np
( ρ β ) m = ( 1 ϕ ) ( ρ β ) f + ϕ ( ρ β ) np
( ρ c p ) m = ( 1 ϕ ) ( ρ c p ) f + ϕ ( ρ c p ) np
In Equations (14)–(16), the indices of m, f, and np are corresponding to the respective characteristics of mixture, pure fluid, and nanoparticles. There are several models regarding dynamic viscosity that can be used, depending on the problem solving conditions. In this work, after detailed reviewing the proposed models, it is found that some researchers have been used the following experimental model for the volume fraction in the selected range (based on [65]). The nanofluid viscosity can be defined through the following equations [65]:
μ np = μ f ( 0.188 + 537.42 ϕ )
μ m = ( 1 ϕ ) μ f + ϕ μ np

2.4. Entropy Generation Analysis

The entropy generation can be attributed to the irreversibility as a consequence of heat transfer (thermal part) and fluid flow friction (frictional part). Several investigations have been dedicated to minimizing the entropy generation through the second law of thermodynamics. The entropy generation rate can be determined regarding the procedure introduced in Reference [66], as follows:
S gen = S gen , f + S gen , T
The following relation is used for the calculation of S gen , T in porous media:
S gen , T = k m T in 2 [ ( T x ) 2 + ( T y ) 2 ]
Additionally, the friction-induced generation of entropy in porous media is determined by:
S gen , f = μ m T in [ 2 ( u x ) 2 + 2 ( v y ) 2 + ( u y + v x ) 2 ] + μ m KT in [ u 2 + v 2 ]
S gen = S gen dV
In the above equations, S gen , T and S gen , f represent the entropy generation induced by heat transfer and fluid friction, respectively.

2.5. Boundary Conditions

The steady laminar flow was considered for the entrance flow to the minichannel, as depicted in Figure 1. The Reynolds number of this flow can be calculated through Equation (6) ( R e = 500 ) . In addition, constant temperatures of T w = 400   K were considered for both wavy walls, while the inlet flow temperature was assumed to be T i n = 300   K . In addition, steady-state parameters were considered for the heat transfer and fluid flow. Furthermore, fully saturated condition ( ε = 0.9 ) with Newtonian and incompressible fluid were presumed for the porous medium. The thermo-physical properties of operating fluid and porous medium were maintained constant. The local thermal equilibrium between the operating fluid and porous medium were also established. The thermo-physical properties of the used water/nanoparticles and porous medium are provided in Table 1 and Table 2, respectively.

2.6. Numerical Simulation

For the 2D numerical simulations, the commercial software of ANSYS FLUENT was utilized. The equations of energy and momentum were discretized through the second-order upwind method. The line-by-line procedure was implemented for solving the obtained equations iteratively. Moreover, velocity and pressure fields were related to each other by the use of SIMPLE algorithm (semi-implicit method for pressure-linked equations). The considered divergence criteria for the conservation equations of mass, energy and momentum was equal to 10 5 . The under-relaxation factors for energy, momentum and pressure equations were considered to be the respective values of 1, 0.7, and 0.3 for maintaining the stability of the solution.

2.7. Governing Equations

The generating grids is regarded as the preliminary and important steps in numerical simulation while considering the subsequent effect on computational time, convergence and solution results. The generated grid for the current problem is depicted in Figure 2. Better performance would be achieved through the use of regular gridding as compared to the irregular counterpart. Finer and regular gridding with radially incrementing elements should be implemented in the elements adjacent to the wall for evaluation of steep gradient of physical features.
The grid independency of the results was also evaluated and provided in Figure 3. As can be seen in this figure, grid 4 revealed approximately 0.5% deviation for the Nusselt number and pressure drop. The results for diverse grid sizes are also provided in Figure 3. According to this figure, the grid with 119,000 nodes would be efficient when considering both accuracy and computational time.

2.8. Model Validation

The present results were compared with the numerical data of Maerefat et al. [68] for the same operating conditions to validate the simulation procedure. The comparison results for Nusselt number in sinusoidal minichannel revealed an acceptable agreement, as provided in Figure 4. The same operating conditions reported by Maerefat et al. [68] are adopted here for the model validation. It is necessary to mention that, parameters of D, x, and Nu in Figure 4 are diameter, length of the pipe, and local Nusselt number along the length of the pipe.

3. Results and Discussion

3.1. Effect of Geometrical Properties on Average Nusselt Number

The diagram of average Nusselt number against different Darcy numbers are plotted in Figure 5, Figure 6, Figure 7, Figure 8 and Figure 9 in terms of wave amplitude, wavelengths, and phase shifts. Based on Figure 5, Figure 6, Figure 7, Figure 8 and Figure 9, the enhanced average Nusselt number would be achieved by the decrement of Darcy number from 0.1 to 0.0001. In Figure 5, for the desired diagram was plotted for sinusoidal minichannel with the wavelength of λ = 12   mm and phase shift of φ = 0 ° . The ameliorated average Nusselt number approximately up to 400 as a result of increasing flow velocity is evident from this figure. As it is clear, the higher average Nusselt number was achieved for all of the studied configurations of corrugated wavy minichannel regarding the introduction of superior disturbance and formation of thinner boundary and the subsequent higher temperature gradients. The highest heat transfer was in accordance with the configuration with wave amplitude of a = 1.5   mm .
The diagram of the average Nusselt number against Darcy number at wave-amplitude of a = 0.5   mm for diverse wavelengths is plotted in Figure 6. As it is obvious from this figure, the reduction of Darcy number was followed by an increment in the average Nusselt number. The smaller Darcy number represents the lower medium permeability, which resulted in an improved heat transfer and increased average Nusselt number. In addition, the highest average Nusselt number was belonged to the wavy minichannel with the wavelength of λ = 10   mm . Moreover, Figure 7 represents the average Nusselt number against Darcy number for the decreasing trend of wave amplitude. Regarding this figure, the wavy tube with the wavelength of λ = 10   mm attained the highest average Nusselt number value when compared to other configurations.
Figure 8 represent the variation of average Nusselt number against Darcy number in five different phase shifts of wavy minichannel. The amplitude and the wavelength were considered to be the constant values of a = 0.5   mm and = 12   mm , respectively. In addition, the variation of average Nusselt number against Darcy numbers for five phase shifts of wavy minichannel with the amplitude of a = 0.5   mm and augmenting wavelength is depicted in Figure 9. Regarding these two figures, the maximum average Nusselt number was achieved in the phase shift of φ = 0 ° , which was followed by the phase shifts of φ = 30 ° ,   60 ° ,   90 ° , and 180°, respectively.

3.2. Effect of Geometrical Properties on Friction Factor

The friction factor for fully developed flow can be calculated from Equation (23):
f = 2 ( L D h ) Δ p ρ U in 2
In the above relation, Δ p , L, and D h are the respective values for pressure drop, minichannel length, and hydraulic diameter. The effects of Darcy number on friction factor for different parameters are represented in Figure 10, Figure 11, Figure 12, Figure 13 and Figure 14. The augmenting trend of friction factor with the reduction of Darcy numbers is obvious from Figure 10, Figure 11, Figure 12, Figure 13 and Figure 14, while it showed a reducing behavior by increasing the porosity from D a = 0.0001 to D a = 0.1 at a constant wave amplitude. Various constant and variable wavy amplitudes in the diagram of friction factor versus Darcy number are analyzed in Figure 10. Regarding this figure, the highest friction factor was attributed to the wavy amplitude of α = 1.5   mm .
The friction factor coefficient variation against Darcy number for different wavelength are plotted in Figure 11. As it is obvious from this figure, the wavelength of λ = 10   mm was corresponding to the maximum value of friction factor, while the wavelength of λ = 15   mm attained the minimum friction factor coefficient. In addition, the variation of friction factor against Darcy numbers for various wavelengths and decreasing wave amplitude is depicted in Figure 12. This figure reveals the minimum friction factor for the configuration with the wavelength of λ = 15   mm . From Figure 11 and Figure 12, it can be inferred that, in all of the studied configurations, the friction factor showed the reducing behavior by increasing Darcy number from 0.0001 to 0.1.
Figure 13 depicts the variation of friction factor in the desired minichannel for different phase shifts and wavelength of λ = 10   mm . The maximum friction factor value was associated with the phase shift of φ = 180 ° as it is obvious from this figure. However, the configuration with the phase shift of φ = 90 ° had attained the minimum friction factor value. Moreover, in Figure 14, the variations of the friction factor in corrugated minichannel for various phase shifts and increasing wavelength are represented. The highest and lowest friction factor values were related to the phase shift of φ = 180 ° and 90°, respectively.

3.3. Effect of Geometrical Properties on the PEC Index

The PEC index has been introduced for the calculation of fluid dynamic and thermal performances. This parameter is beneficial for evaluating the heat transfer performance and pumping power in corrugated wavy-wall minichannel with variable amplitudes. This factor can be determined by the use of estimated friction factor and average Nusselt numbers coefficient, as follows [13]:
PEC = ( Nu av . w Nu av , s ) ( f w f s ) 1 / 3
In which, Nu av . w and Nu av , s as well as f w and f s represent the estimated average Nusselt number and friction factor for wavy and smooth minichannel, respectively. The variations of PEC index against Darcy number for different geometrical properties are provided in Figure 15, Figure 16, Figure 17, Figure 18 and Figure 19. Regarding these figures for all of the studied configurations, PEC index was ameliorated by increasing the Darcy number as a consequence of heat transfer coefficient effect on the flow. The highest PEC index was attributed to D a = 0.1 , as can be seen. The variation of PEC index versus Darcy numbers for different wave-amplitudes is plotted in Figure 15. From this figure, it can be realized that wavy amplitude of a = 1   mm gained the maximum PEC index among all other configurations. Moreover, the application of wave-amplitude of a = 0.5   mm , resulting in a higher Nusselt number.
The PEC index was determined by average Nusselt numbers and friction factor coefficients, as can be seen in Figure 16. Regarding this figure, a declining trend of PEC index by the increasing of Darcy number can be inferred. In addition, the determination of the PEC index from friction factor coefficients and average Nusselt numbers can be also inferred from Figure 17. The best configuration from the PEC index perspective belonged to the sinusoidal-wavy minichannel in porous media with the wavelength of λ = 10   mm .
Similar to the previous cases, according to Figure 18, the PEC was determined for various phase shifts, a = 0.5   mm and = 10   mm . This figure states the enhancement of the PEC index in minichannel by increasing Darcy number. The highest PEC index (approximately 16.95) in the sinusoidal corrugated minichannels was associated with the configuration with the phase shift of φ = 60 ° , wave amplitude of a = 0.5   mm , and wavelength of λ = 10   mm at D a = 0.1 . Hence, this configuration was selected as the optimal configuration from the PEC index viewpoint and used in further investigation of different nanofluids. In addition, the variation of PEC index against Darcy numbers for various phase shifts and increasing wavelength in the range of 9   mm λ 15   mm and a = 0.5   mm is plotted in Figure 19. As can be inferred from Figure 19, the most efficient configuration of sinusoidal-wavy minichannel in porous medium was determined as the configuration with the constant wave amplitude of a = 0.5   mm , augmenting wavelength and phase shift of φ = 60 ° .

3.4. Effect of Geometrical Properties and Porosity on Entropy Generation

The variations of generation of entropy according to the heat transfer (thermal entropy generation) and fluid flow (frictional entropy generation) are depicted in Figure 20, Figure 21, Figure 22, Figure 23 and Figure 24, respectively. These figures clarify the increasing and decreasing contributions of thermal and frictional entropy generations as a result of increased Darcy number, respectively. These figures also provide the profile of frictional entropy generation, which declares the higher contribution of frictional entropy generation at D a = 0.00001 in comparison of the thermal entropy generation.
The variation of entropy generations for different wavelengths and wave amplitudes are depicted in Figure 21 and Figure 22, respectively. The geometrical properties of different configurations for the cases with constant wave amplitude of a = 0.5   mm and reducing wave amplitude are depicted in Figure 21 and Figure 22, respectively. The approximately constant contribution of heat transfer in entropy generation by increasing wavelengths can be concluded regarding Figure 21. The variation of entropy generations against Darcy numbers for various phase shifts and decreasing wave amplitude is plotted in Figure 22. From this figure, it can be found that the use of porous medium and increasing Darcy number at all cases would result in the increased share of entropy generation (thermal portion) besides a decreased share of entropy generation (frictional portion) in the overall generated entropy.
The respective thermal and frictional entropy generations in various phase shifts, wave amplitude of a = 0.5   mm and wavelengths of λ = 12   mm , as well as the increasing wavelengths are presented in the respective diagrams of Figure 23 and Figure 24. According to the obtained results, thermal entropy generations share was increased, while the frictional share was decreased by increasing the Darcy number. As it is obvious from these figures, frictional entropy generation attained the higher contribution in total entropy generation at D a = 0.0001 and D a = 0.00001 . The highest thermal and frictional entropy generations for the desired minichannel was associated with the phase shift of φ = 0 ° and φ = 180 ° , wave amplitude of a = 0.5   mm , and wavelength of λ = 10   mm . Moreover, regarding Figure 24, the sinusoidal-wavy minichannel in porous medium with constant wave-amplitude of a = 0.5   mm , augmenting wavelength and phase shift of φ = 0 ° and φ = 180 ° , are the most efficient configuration for thermal and frictional entropy generations in comparison of other cases, respectively.

3.5. Effect of Utilizing Nanoparticles in Base Fluid and Volume Fraction

The corrugated minichannels with a = 0.5   mm , λ = 10   mm , and φ = 60 ° introduced as the most efficient configuration in terms of Nusselt number, friction factor, and PEC index among all studied ones, according to the obtained results in previous sections. Hence, this configuration was selected for conducting the further investigations on the addition of nanoparticles (Al2O3 and CuO) to the pure fluid. The diagrams of average Nusselt number, friction factor, and PEC index in terms of the volume fraction of two-phase nanofluid for two different nanofluids (Al2O3 and CuO nanofluid) are represented in Figure 25, Figure 26 and Figure 27, respectively. Regarding these figures, by augmenting the volume fraction of the nanoparticles, an increased average Nusselt number, and friction factor besides the decreased PEC index can be realized. The highest amelioration in Nusselt number was reported for the case of using Al2O3 nanofluid regarding the thermo-physical properties.

4. Conclusions

In this study, the influence of geometrical variables on the forced convective heat transfer and entropy generation of the minichannels with the simultaneous use of corrugation and porous media was investigated numerically. For this purpose, a 2D CFD scheme was proposed for modeling incompressible, steady, and laminar flow of nanofluid flow through the corrugated minichannel. The fluid flow in the porous media was characterized by Darcy–Brinkman–Forchheimer model. In addition, the effect of different parameters, including amplitude, wavelength, shift phase, and porosity of the wavy-wall were analyzed. The main findings of the current work can be summarized, as follows:
  • Enhanced average Nusselt number was achieved by decreasing the Darcy number from 0.1 to 0.0001. The maximum average Nusselt number value was obtained as 430.91.
  • Reducing behavior of friction factor with increasing the Darcy number from 0.0001 to 0.1 was observed at constant wavelength, amplitude and phase shift. Moreover, of the phase shifts of φ = 180 ° and 90° corresponded to the maximum and the minimum friction factor, respectively.
  • The PEC index of sinusoidal corrugated minichannel was maximized (maximum value of 16.95) for the phase shift of φ = 60 ° , wave amplitude of a = 0.5   mm , and wavelength of λ = 10   mm at D a = 0.1 .
  • The application of porous medium as well as increasing Darcy number at all cases, resulted in the increase and decrease in the respective thermal and frictional entropy generations.
  • According to the obtained results, the geometrical properties of the most efficient configuration were λ = 10   mm , a = 0.5   mm , and φ = 60 ° . This configuration was used for the investigations concerning the addition of nanoparticles (Al2O3 and CuO). Furthermore, the highest Nusselt number augmentation was stated for Al2O3 nanofluid, because of its thermo-physical characteristics.

Author Contributions

E.A., H.M.; Conceptualization, Methodology, Software, Writing—Original Draft Preparation, Writing—Review & Editing, H.S.; Project Administration, Supervision, A.M.G.; resources, Writing—Review & Editing, Project Administration. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

aWavy amplitude (m)
CSpecific Heat ( J   kg 1   K 1 )
DDiameter ( m )
KPermeability of Porous Medium ( m 2 )
kThermal Conductivity ( W   m 1   K 1 )
LLength (m)
PDimensionless Pressure
pPressure (Pa)
STotal Entropy
S g e n , T Thermal Entropy Generation ( W   K 1 )
S g e n , f Friction Entropy Generation ( W   K 1 )
TTemperature ( K )
U , V Dimensionless Velocity
u , v Velocity Components ( m   s 1 )
X , Y Dimensionless Coordinates
x , y Dimensional Coordinates (m)
Greek symbols
αThermal Diffusivity ( m 2   s 1 )
λWavelength (m)
βThermal Expansion Coefficient ( K 1 )
εPorosity
θDimensionless Temperature
ϑKinematic Viscosity ( m 2   s 1 )
μDynamic Viscosity ( kg   m 1 s 1 )
ρDensity ( kg   m 3 )
φ wavy-wall shift phase/degrees
ϕ Volume Fraction
Subscripts
a v Averaged
fPure Fluid
sSmooth
wWavy
g e n Generation
hHydraulic
mMixture
n p Nanoparticle
Abbreviations
CFDComputational Fluid Dynamics
D a Darcy Number
N u Nusselt Number
P r Prandtl Number
PECPerformance Evaluation Criteria
R e Reynolds Number
SIMPLESemi-Implicit Method for Pressure-Linked Equations

References

  1. Ho, C.; Liu, Y.-C.; Ghalambaz, M.; Yan, W.-M. Forced convection heat transfer of Nano-Encapsulated Phase Change Material (NEPCM) suspension in a mini-channel heatsink. Int. J. Heat Mass Transf. 2020, 155, 119858. [Google Scholar] [CrossRef]
  2. Hatami, M.; Jing, D. Evaluation of wavy direct absorption solar collector (DASC) performance using different nanofluids. J. Mol. Liq. 2017, 229, 203–211. [Google Scholar] [CrossRef]
  3. Khoshvaght-Aliabadi, M.; Tatari, M.; Salami, M. Analysis on Al2O3/water nanofluid flow in a channel by inserting corrugated/perforated fins for solar heating heat exchangers. Renew. Energy 2018, 115, 1099–1108. [Google Scholar] [CrossRef]
  4. Khoshvaght-Aliabadi, M.; Hassani, S.; Mazloumi, S.H. Performance enhancement of straight and wavy miniature heat sinks using pin-fin interruptions and nanofluids. Chem. Eng. Process. Process Intensif. 2017, 122, 90–108. [Google Scholar] [CrossRef]
  5. Ma, T.; Du, L.-X.; Sun, N.; Zeng, M.; Sundén, B.; Wang, Q.-W. Experimental and numerical study on heat transfer and pressure drop performance of Cross-Wavy primary surface channel. Energy Convers. Manag. 2016, 125, 80–90. [Google Scholar] [CrossRef]
  6. Rostami, J.; Abbassi, A.; Harting, J. Heat transfer by nanofluids in wavy microchannels. Adv. Powder Technol. 2018, 29, 925–933. [Google Scholar] [CrossRef] [Green Version]
  7. Zhu, X.W.; Fu, Y.H.; Zhao, J.Q. A novel wavy-tape insert configuration for pipe heat transfer augmentation. Energy Convers. Manag. 2016, 127, 140–148. [Google Scholar] [CrossRef]
  8. Pati, S.; Mehta, S.K.; Borah, A. Numerical investigation of thermo-hydraulic transport characteristics in wavy channels: Comparison between raccoon and serpentine channels. Int. Commun. Heat Mass Transf. 2017, 88, 171–176. [Google Scholar] [CrossRef]
  9. Aslan, E.; Taymaz, I.; Islamoglu, Y. Finite volume simulation for convective heat transfer in wavy channels. Heat Mass Transf. 2016, 52, 483–497. [Google Scholar] [CrossRef]
  10. Valinataj-Bahnemiri, P.; Ramiar, A.; Manavi, S.; Mozaffari, A. Heat transfer optimization of two phase modeling of nanofluid in a sinusoidal wavy channel using Artificial Bee Colony technique. Eng. Sci. Technol. Int. J. 2015, 18, 727–737. [Google Scholar] [CrossRef] [Green Version]
  11. Shehzad, N.; Zeeshan, A.; Ellahi, R.; Vafai, K. Convective heat transfer of nanofluid in a wavy channel: Buongiorno’s mathematical model. J. Mol. Liq. 2016, 222, 446–455. [Google Scholar] [CrossRef]
  12. Mohammed, K.A.; Talib, A.A.; Nuraini, A.; Ahmed, K. Review of forced convection nanofluids through corrugated facing step. Renew. Sustain. Energy Rev. 2017, 75, 234–241. [Google Scholar] [CrossRef]
  13. Yang, L.; Du, K.; Zhang, Z. Heat transfer and flow optimization of a novel sinusoidal minitube filled with non-Newtonian SiC/EG-water nanofluids. Int. J. Mech. Sci. 2020, 168, 105310. [Google Scholar] [CrossRef]
  14. Mohammed, H.; Gunnasegaran, P.; Shuaib, N. Numerical simulation of heat transfer enhancement in wavy microchannel heat sink. Int. Commun. Heat Mass Transf. 2011, 38, 63–68. [Google Scholar] [CrossRef]
  15. Chang, K.; Hughes, T.J.R.; Calo, V.M. Isogeometric variational multiscale large-eddy simulation of fully-developed turbulent flow over a wavy wall. Comput. Fluids 2012, 68, 94–104. [Google Scholar] [CrossRef]
  16. Lecrivain, G.; Rayan, R.; Hurtado, A.; Hampel, U. Using quasi-DNS to investigate the deposition of elongated aerosol particles in a wavy channel flow. Comput. Fluids 2016, 124, 78–85. [Google Scholar] [CrossRef]
  17. Rush, T.; Newell, T.; Jacobi, A. An experimental study of flow and heat transfer in sinusoidal wavy passages. Int. J. Heat Mass Transf. 1999, 42, 1541–1553. [Google Scholar] [CrossRef]
  18. Wang, C.-C.; Chen, C.-K. Forced convection in a wavy-wall channel. Int. J. Heat Mass Transf. 2002, 45, 2587–2595. [Google Scholar] [CrossRef]
  19. Zhang, L.; Che, D. Influence of corrugation profile on the thermalhydraulic performance of cross-corrugated plates. Numer. Heat Transf. Part A Appl. 2011, 59, 267–296. [Google Scholar] [CrossRef]
  20. Naphon, P.; Kornkumjayrit, K. Numerical analysis on the fluid flow and heat transfer in the channel with V-shaped wavy lower plate. Int. Commun. Heat Mass Transf. 2008, 35, 839–843. [Google Scholar] [CrossRef]
  21. Khoshvaght-Aliabadi, M. Influence of different design parameters and Al2O3-water nanofluid flow on heat transfer and flow characteristics of sinusoidal-corrugated channels. Energy Convers. Manag. 2014, 88, 96–105. [Google Scholar] [CrossRef]
  22. Aminian, E.; Moghadasi, H.; Saffari, H. Magnetic field effects on forced convection flow of a hybrid nanofluid in a cylinder filled with porous media: A numerical study. J. Therm. Anal. Calorim. 2020, 141, 2019–2031. [Google Scholar] [CrossRef]
  23. Siavashi, M.; Bahrami, H.R.T.; Aminian, E. Optimization of heat transfer enhancement and pumping power of a heat exchanger tube using nanofluid with gradient and multi-layered porous foams. Appl. Therm. Eng. 2018, 138, 465–474. [Google Scholar] [CrossRef]
  24. Siavashi, M.; Bahrami, H.R.T.; Aminian, E.; Saffari, H. Numerical analysis on forced convection enhancement in an annulus using porous ribs and nanoparticle addition to base fluid. J. Cent. South Univ. 2019, 26, 1089–1098. [Google Scholar] [CrossRef]
  25. Pekmen, B.; Tezer-Sezgin, M. MHD flow and heat transfer in a lid-driven porous enclosure. Comput. Fluids 2014, 89, 191–199. [Google Scholar] [CrossRef]
  26. Bahrami, H.R.T.; Aminian, E.; Saffari, H. Energy transfer enhancement inside an annulus using gradient porous ribs and nanofluids. J. Energy Resour. Technol. 2020, 142, 122102. [Google Scholar] [CrossRef]
  27. Mehryan, S.; Ghalambaz, M.; Chamkha, A.J.; Izadi, M. Numerical study on natural convection of Ag–MgO hybrid/water nanofluid inside a porous enclosure: A local thermal non-equilibrium model. Powder Technol. 2020, 367, 443–455. [Google Scholar] [CrossRef]
  28. Moghadasi, H.; Aminian, E.; Saffari, H.; Mahjoorghani, M.; Emamifar, A. Numerical analysis on Laminar Forced Convection improvement of Hybrid Nanofluid within a U-Bend Pipe in Porous Media. Int. J. Mech. Sci. 2020, 179, 105659. [Google Scholar] [CrossRef]
  29. Akbarzadeh, M.; Maghrebi, M.J. Combined effects of corrugated walls and porous inserts on performance improvement in a heat exchanger channel. Int. J. Therm. Sci. 2018, 127, 266–276. [Google Scholar] [CrossRef]
  30. Arasteh, H.; Mashayekhi, R.; Ghaneifar, M.; Toghraie, D.; Afrand, M. Heat transfer enhancement in a counter-flow sinusoidal parallel-plate heat exchanger partially filled with porous media using metal foam in the channels’ divergent sections. J. Therm. Anal. Calorim. 2019, 141, 1669–1685. [Google Scholar] [CrossRef]
  31. Hajipour, M.; Dehkordi, A.M. Analysis of nanofluid heat transfer in parallel-plate vertical channels partially filled with porous medium. Int. J. Therm. Sci. 2012, 55, 103–113. [Google Scholar] [CrossRef]
  32. Nazari, S.; Toghraie, D. Numerical simulation of heat transfer and fluid flow of Water-CuO Nanofluid in a sinusoidal channel with a porous medium. Phys. E Low-Dimens. Syst. Nanostruct. 2017, 87, 134–140. [Google Scholar] [CrossRef]
  33. Abu-Nada, E. Application of nanofluids for heat transfer enhancement of separated flows encountered in a backward facing step. Int. J. Heat Fluid Flow 2008, 29, 242–249. [Google Scholar] [CrossRef]
  34. Ahmadi, A.A.; Khodabandeh, E.; Moghadasi, H.; Malekian, N.; Akbari, O.A.; Bahiraei, M. Numerical study of flow and heat transfer of water-Al2O3 nanofluid inside a channel with an inner cylinder using Eulerian–Lagrangian approach. J. Therm. Anal. Calorim. 2018, 132, 651–665. [Google Scholar] [CrossRef]
  35. Farhangmehr, V.; Moghadasi, H.; Asiaei, S. A nanofluid MHD flow with heat and mass transfers over a sheet by nonlinear boundary conditions: Heat and mass transfers enhancement. J. Cent. South Univ. 2019, 26, 1205–1217. [Google Scholar] [CrossRef]
  36. Malekian, S.; Fathi, E.; Malekian, N.; Moghadasi, H.; Moghimi, M. Analytical and numerical investigations of unsteady graphene oxide nanofluid flow between two parallel plates. Int. J. Thermophys. 2018, 39, 100. [Google Scholar] [CrossRef]
  37. Das, K. Slip flow and convective heat transfer of nanofluids over a permeable stretching surface. Comput. Fluids 2012, 64, 34–42. [Google Scholar] [CrossRef]
  38. Ghalambaz, M.; Doostani, A.; Izadpanahi, E.; Chamkha, A.J. Conjugate natural convection flow of Ag–MgO/water hybrid nanofluid in a square cavity. J. Therm. Anal. Calorim. 2020, 139, 2321–2336. [Google Scholar] [CrossRef]
  39. Takabi, B.; Gheitaghy, A.M.; Tazraei, P. Hybrid water-based suspension of Al2O3 and Cu nanoparticles on laminar convection effectiveness. J. Thermophys. Heat Transf. 2016, 30, 523–532. [Google Scholar] [CrossRef]
  40. Heidary, H.; Kermani, M. Effect of nano-particles on forced convection in sinusoidal-wall channel. Int. Commun. Heat Mass Transf. 2010, 37, 1520–1527. [Google Scholar] [CrossRef]
  41. Sadripour, S. Investigation of Flow Characteristics And Heat Transfer Enhancement in a Nanofluid Flow in a Corrugated Duct. J. Appl. Mech. Tech. Phys. 2018, 59, 1049–1057. [Google Scholar] [CrossRef]
  42. Ho, C.-J.; Wei, L.; Li, Z. An experimental investigation of forced convective cooling performance of a microchannel heat sink with Al2O3/water nanofluid. Appl. Therm. Eng. 2010, 30, 96–103. [Google Scholar] [CrossRef]
  43. Akbari, O.A.; Safaei, M.R.; Goodarzi, M.; Akbar, N.S.; Zarringhalam, M.; Shabani, G.A.S.; Dahari, M. A modified two-phase mixture model of nanofluid flow and heat transfer in a 3-D curved microtube. Adv. Powder Technol. 2016, 27, 2175–2185. [Google Scholar] [CrossRef]
  44. Khoshvaght-Aliabadi, M.; Sahamiyan, M. Performance of nanofluid flow in corrugated minichannels heat sink (CMCHS). Energy Convers. Manag. 2016, 108, 297–308. [Google Scholar] [CrossRef]
  45. Chen, L.; Wu, C.; Sun, F. Finite time thermodynamic optimization or entropy generation minimization of energy systems. J. Non-Equilib. Thermodyn. 1999, 24, 327–359. [Google Scholar] [CrossRef]
  46. Chen, L.; Meng, F.; Sun, F. Thermodynamic analyses and optimization for thermoelectric devices: The state of the arts. Sci. China Technol. Sci. 2016, 59, 442–455. [Google Scholar] [CrossRef]
  47. Sheremet, M.A.; Oztop, H.F.; Pop, I.; Abu-Hamdeh, N. Analysis of entropy generation in natural convection of nanofluid inside a square cavity having hot solid block: Tiwari and Das’ model. Entropy 2016, 18, 9. [Google Scholar] [CrossRef]
  48. Siavashi, M.; Bahrami, H.R.T.; Saffari, H. Numerical investigation of porous rib arrangement on heat transfer and entropy generation of nanofluid flow in an annulus using a two-phase mixture model. Numer. Heat Transf. Part A Appl. 2017, 71, 1251–1273. [Google Scholar] [CrossRef]
  49. Gholamalipour, P.; Siavashi, M.; Doranehgard, M.H. Eccentricity effects of heat source inside a porous annulus on the natural convection heat transfer and entropy generation of Cu-water nanofluid. Int. Commun. Heat Mass Transf. 2019, 109, 104367. [Google Scholar] [CrossRef]
  50. Tayebi, T.; Chamkha, A.J. Entropy generation analysis due to MHD natural convection flow in a cavity occupied with hybrid nanofluid and equipped with a conducting hollow cylinder. J. Therm. Anal. Calorim. 2020, 139, 2165–2179. [Google Scholar] [CrossRef]
  51. Tayebi, T.; Chamkha, A.J. Entropy generation analysis during MHD natural convection flow of hybrid nanofluid in a square cavity containing a corrugated conducting block. Int. J. Numer. Methods Heat Fluid Flow 2019, 30, 1115–1136. [Google Scholar] [CrossRef]
  52. Zhang, L.; Chen, L.; Xia, S.; Wang, C.; Sun, F. Entropy generation minimization for reverse water gas shift (RWGS) reactors. Entropy 2018, 20, 415. [Google Scholar] [CrossRef] [Green Version]
  53. Chen, L.; Yang, A.; Xie, Z.; Sun, F. Constructal entropy generation rate minimization for cylindrical pin-fin heat sinks. Int. J. Therm. Sci. 2017, 111, 168–174. [Google Scholar] [CrossRef]
  54. Chen, L.; Xia, S.; Sun, F. Entropy generation minimization for isothermal crystallization processes with a generalized mass diffusion law. Int. J. Heat Mass Transf. 2018, 116, 1–8. [Google Scholar] [CrossRef]
  55. Feng, H.; Chen, L.; Xie, Z.; Sun, F. “Disc-point” heat and mass transfer constructal optimization for solid–gas reactors based on entropy generation minimization. Energy 2015, 83, 431–437. [Google Scholar] [CrossRef]
  56. Fersadou, I.; Kahalerras, H.; El Ganaoui, M. MHD mixed convection and entropy generation of a nanofluid in a vertical porous channel. Comput. Fluids 2015, 121, 164–179. [Google Scholar] [CrossRef]
  57. Shashikumar, N.; Gireesha, B.; Mahanthesh, B.; Prasannakumara, B.; Chamkha, A.J. Entropy generation analysis of magneto-nanoliquids embedded with aluminium and titanium alloy nanoparticles in microchannel with partial slips and convective conditions. Int. J. Numer. Methods Heat Fluid Flow 2019, 29, 3638–3658. [Google Scholar]
  58. Sheremet, M.; Pop, I.; Öztop, H.F.; Abu-Hamdeh, N. Natural convection of nanofluid inside a wavy cavity with a non-uniform heating. Int. J. Numer. Methods Heat Fluid Flow 2017, 27, 958–980. [Google Scholar] [CrossRef]
  59. Cho, C.-C.; Chen, C.-L. Natural convection heat transfer and entropy generation in wavy-wall enclosure containing water-based nanofluid. Int. J. Heat Mass Transf. 2013, 61, 749–758. [Google Scholar] [CrossRef]
  60. Akbarzadeh, M.; Rashidi, S.; Karimi, N.; Omar, N. First and second laws of thermodynamics analysis of nanofluid flow inside a heat exchanger duct with wavy walls and a porous insert. J Therm. Anal. Calorim. 2019, 135, 177–194. [Google Scholar] [CrossRef] [Green Version]
  61. Esfahani, J.; Akbarzadeh, M.; Rashidi, S.; Rosen, M.; Ellahi, R. Influences of wavy wall and nanoparticles on entropy generation over heat exchanger plat. Int. J. Heat Mass Transf. 2017, 109, 1162–1171. [Google Scholar] [CrossRef]
  62. Dormohammadi, R.; Farzaneh-Gord, M.; Ebrahimi-Moghadam, A.; Ahmadi, M.H. Heat transfer and entropy generation of the nanofluid flow inside sinusoidal wavy channels. J. Mol. Liq. 2018, 269, 229–240. [Google Scholar] [CrossRef]
  63. Siavashi, M.; Bahrami, H.R.T.; Saffari, H. Numerical investigation of flow characteristics, heat transfer and entropy generation of nanofluid flow inside an annular pipe partially or completely filled with porous media using two-phase mixture model. Energy 2015, 93, 2451–2466. [Google Scholar] [CrossRef]
  64. Wang, B.; Wang, L.; Wu, J.; Hong, Y.; Yan, X.; Gong, H.; Xu, Z. Design and characterization of isothermal chambers filled with gradient-porous materials. J. Therm. Sci. Technol. 2017, 12, JTST0008. [Google Scholar] [CrossRef] [Green Version]
  65. Miller, A.; Gidaspow, D. Dense, vertical gas-solid flow in a pipe. AIChE J. 1992, 38, 1801–1815. [Google Scholar] [CrossRef]
  66. Ingham, D.B.; Pop, I. Transport Phenomena in Porous Media III; Elsevier: Amsterdam, The Netherlands, 2005. [Google Scholar]
  67. Wang, B.; Hong, Y.; Wang, L.; Fang, X.; Wang, P.; Xu, Z. Development and numerical investigation of novel gradient-porous heat sinks. Energy Convers. Manag. 2015, 106, 1370–1378. [Google Scholar] [CrossRef]
  68. Maerefat, M.; Mahmoudi, S.Y.; Mazaheri, K. Numerical simulation of forced convection enhancement in a pipe by porous inserts. Heat Transf. Eng. 2011, 32, 45–56. [Google Scholar] [CrossRef]
Figure 1. Schematic geometry of the problem.
Figure 1. Schematic geometry of the problem.
Entropy 22 01008 g001
Figure 2. The generated grid for the current work.
Figure 2. The generated grid for the current work.
Entropy 22 01008 g002
Figure 3. Variations of pressure drop and Nusselt number with diverse cell numbers.
Figure 3. Variations of pressure drop and Nusselt number with diverse cell numbers.
Entropy 22 01008 g003
Figure 4. Comparison between the present numerical data and the numerical data of Maerefat et al. [68].
Figure 4. Comparison between the present numerical data and the numerical data of Maerefat et al. [68].
Entropy 22 01008 g004
Figure 5. Comparison of average Nusselt number for diverse Darcy numbers in various wave amplitudes ( λ = 12   mm ,   φ = 0 ° ) .
Figure 5. Comparison of average Nusselt number for diverse Darcy numbers in various wave amplitudes ( λ = 12   mm ,   φ = 0 ° ) .
Entropy 22 01008 g005
Figure 6. Comparison of average Nusselt number for diverse Darcy numbers with various wavelengths ( a = 0.5   mm ,   φ = 0 ° ) .
Figure 6. Comparison of average Nusselt number for diverse Darcy numbers with various wavelengths ( a = 0.5   mm ,   φ = 0 ° ) .
Entropy 22 01008 g006
Figure 7. Comparison of average Nusselt number for diverse Darcy numbers with various wavelengths and decreasing wave amplitude ( φ = 0 ° ) .
Figure 7. Comparison of average Nusselt number for diverse Darcy numbers with various wavelengths and decreasing wave amplitude ( φ = 0 ° ) .
Entropy 22 01008 g007
Figure 8. Comparison of average Nusselt number for diverse Darcy numbers with various phase shifts and wavelength of ( a = 0.5   mm ,   λ = 10   mm ) .
Figure 8. Comparison of average Nusselt number for diverse Darcy numbers with various phase shifts and wavelength of ( a = 0.5   mm ,   λ = 10   mm ) .
Entropy 22 01008 g008
Figure 9. Comparison of average Nusselt number for diverse Darcy numbers with various phase shifts and increasing wavelength.
Figure 9. Comparison of average Nusselt number for diverse Darcy numbers with various phase shifts and increasing wavelength.
Entropy 22 01008 g009
Figure 10. Comparison of friction factor for diverse Darcy numbers in various wave-amplitudes, ( λ = 12   mm ,   φ = 0 ° ) .
Figure 10. Comparison of friction factor for diverse Darcy numbers in various wave-amplitudes, ( λ = 12   mm ,   φ = 0 ° ) .
Entropy 22 01008 g010
Figure 11. Comparison of friction factor for diverse Darcy numbers with various wavelengths ( a = 0.5   mm ,   φ = 0 ° ) .
Figure 11. Comparison of friction factor for diverse Darcy numbers with various wavelengths ( a = 0.5   mm ,   φ = 0 ° ) .
Entropy 22 01008 g011
Figure 12. Comparison of friction factor for diverse Darcy numbers with various wavelengths and decreasing wave amplitude ( 0.5   mm a 1.5   mm ) .
Figure 12. Comparison of friction factor for diverse Darcy numbers with various wavelengths and decreasing wave amplitude ( 0.5   mm a 1.5   mm ) .
Entropy 22 01008 g012
Figure 13. Comparison of friction factor for diverse Darcy numbers with various phase shifts ( a = 0.5 ,   λ = 10   mm ) .
Figure 13. Comparison of friction factor for diverse Darcy numbers with various phase shifts ( a = 0.5 ,   λ = 10   mm ) .
Entropy 22 01008 g013
Figure 14. Comparison of friction factor for different Darcy numbers with various phase shifts and increasing wavelength.
Figure 14. Comparison of friction factor for different Darcy numbers with various phase shifts and increasing wavelength.
Entropy 22 01008 g014
Figure 15. The variation of performance evaluation criteria (PEC) index against different Darcy numbers for various wave-amplitudes, ( λ = 12   mm ,   φ = 0 ° ) .
Figure 15. The variation of performance evaluation criteria (PEC) index against different Darcy numbers for various wave-amplitudes, ( λ = 12   mm ,   φ = 0 ° ) .
Entropy 22 01008 g015
Figure 16. The variation of PEC index against different Darcy numbers for various wavelengths, ( a = 0.5   mm ,   φ = 0 ° ) .
Figure 16. The variation of PEC index against different Darcy numbers for various wavelengths, ( a = 0.5   mm ,   φ = 0 ° ) .
Entropy 22 01008 g016
Figure 17. The variation of PEC index against different Darcy numbers for various wavelengths and decreasing wave amplitude.
Figure 17. The variation of PEC index against different Darcy numbers for various wavelengths and decreasing wave amplitude.
Entropy 22 01008 g017
Figure 18. The variation of PEC index against different Darcy numbers for various phase shifts, ( a = 0.5   mm ,   λ = 10   mm ) .
Figure 18. The variation of PEC index against different Darcy numbers for various phase shifts, ( a = 0.5   mm ,   λ = 10   mm ) .
Entropy 22 01008 g018
Figure 19. The variation of PEC index against different Darcy numbers for various phase shifts and increasing wavelength, ( a = 0.5   mm ,   9   mm λ 15   mm ) .
Figure 19. The variation of PEC index against different Darcy numbers for various phase shifts and increasing wavelength, ( a = 0.5   mm ,   9   mm λ 15   mm ) .
Entropy 22 01008 g019
Figure 20. Entropy generation in terms of Darcy numbers at various wave-amplitudes for (a) thermal entropy generation and (b) friction entropy generation ( λ = 12   mm ,   φ = 0 ° ) .
Figure 20. Entropy generation in terms of Darcy numbers at various wave-amplitudes for (a) thermal entropy generation and (b) friction entropy generation ( λ = 12   mm ,   φ = 0 ° ) .
Entropy 22 01008 g020
Figure 21. Entropy generation in terms of Darcy numbers at various wavelengths for (a) thermal entropy generation and (b) friction entropy generation ( a = 0.5   mm ,   φ = 0 ° ) .
Figure 21. Entropy generation in terms of Darcy numbers at various wavelengths for (a) thermal entropy generation and (b) friction entropy generation ( a = 0.5   mm ,   φ = 0 ° ) .
Entropy 22 01008 g021
Figure 22. Entropy generation in terms of Darcy numbers at various wavelengths and decreasing wave-amplitude for (a) thermal entropy generation and (b) friction entropy generation.
Figure 22. Entropy generation in terms of Darcy numbers at various wavelengths and decreasing wave-amplitude for (a) thermal entropy generation and (b) friction entropy generation.
Entropy 22 01008 g022
Figure 23. Entropy generation in terms of Darcy numbers at various phase shifts for (a) thermal entropy generation and (b) friction entropy generation, ( a = 0.5 ,   λ = 10   mm ) .
Figure 23. Entropy generation in terms of Darcy numbers at various phase shifts for (a) thermal entropy generation and (b) friction entropy generation, ( a = 0.5 ,   λ = 10   mm ) .
Entropy 22 01008 g023
Figure 24. Entropy generation in terms of Darcy numbers at various phase shifts and increasing wavelength for (a) thermal entropy generation and (b) friction entropy generations.
Figure 24. Entropy generation in terms of Darcy numbers at various phase shifts and increasing wavelength for (a) thermal entropy generation and (b) friction entropy generations.
Entropy 22 01008 g024
Figure 25. Comparison of average Nusselt number for diverse volume fraction in two case of Al 2 O 3 and CuO nanofluid ( λ = 10   mm , a = 0.5   mm , and φ = 60 ° ) .
Figure 25. Comparison of average Nusselt number for diverse volume fraction in two case of Al 2 O 3 and CuO nanofluid ( λ = 10   mm , a = 0.5   mm , and φ = 60 ° ) .
Entropy 22 01008 g025
Figure 26. Comparison of friction factor for diverse volume fraction in two case of Al 2 O 3 and CuO nanofluid, ( λ = 10   mm , a = 0.5   mm , and φ = 60 ° ) .
Figure 26. Comparison of friction factor for diverse volume fraction in two case of Al 2 O 3 and CuO nanofluid, ( λ = 10   mm , a = 0.5   mm , and φ = 60 ° ) .
Entropy 22 01008 g026
Figure 27. Comparison of PEC index for diverse volume fraction in two case of Al 2 O 3 and CuO nanofluid, ( λ = 10   mm , a = 0.5   mm , and φ = 60 ° ) .
Figure 27. Comparison of PEC index for diverse volume fraction in two case of Al 2 O 3 and CuO nanofluid, ( λ = 10   mm , a = 0.5   mm , and φ = 60 ° ) .
Entropy 22 01008 g027
Table 1. Thermo-physical characteristics of water/nanoparticles (water/Al2O3-CuO) [32,63].
Table 1. Thermo-physical characteristics of water/nanoparticles (water/Al2O3-CuO) [32,63].
PropertiesWaterAl2O3CuOUnit
c p 4182773540 J   kg 1   K 1
k0.64018 W   m 1   K 1
ρ998.238806500 kg   m 3
μ0.001003-- kg   m 1 s 1
D p -4025nm
β 2.1 × 10 4 0.85 × 10 5 0.85 × 10 5 K 1
Table 2. Thermo-physical characteristics of the porous medium [67].
Table 2. Thermo-physical characteristics of the porous medium [67].
PropertiesAluminumUnit
c p 871 J   k g 1   K 1
k202 W   m 1   K 1
ρ2719 kg   m 3
ε0.9-

Share and Cite

MDPI and ACS Style

Aminian, E.; Moghadasi, H.; Saffari, H.; Gheitaghy, A.M. Investigation of Forced Convection Enhancement and Entropy Generation of Nanofluid Flow through a Corrugated Minichannel Filled with a Porous Media. Entropy 2020, 22, 1008. https://0-doi-org.brum.beds.ac.uk/10.3390/e22091008

AMA Style

Aminian E, Moghadasi H, Saffari H, Gheitaghy AM. Investigation of Forced Convection Enhancement and Entropy Generation of Nanofluid Flow through a Corrugated Minichannel Filled with a Porous Media. Entropy. 2020; 22(9):1008. https://0-doi-org.brum.beds.ac.uk/10.3390/e22091008

Chicago/Turabian Style

Aminian, Ehsan, Hesam Moghadasi, Hamid Saffari, and Amir Mirza Gheitaghy. 2020. "Investigation of Forced Convection Enhancement and Entropy Generation of Nanofluid Flow through a Corrugated Minichannel Filled with a Porous Media" Entropy 22, no. 9: 1008. https://0-doi-org.brum.beds.ac.uk/10.3390/e22091008

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