Next Article in Journal
Mechanical Fault Diagnosis of HVCBs Based on Multi-Feature Entropy Fusion and Hybrid Classifier
Previous Article in Journal
A Simple Explicit Expression for the Flocculation Dynamics Modeling of Cohesive Sediment Based on Entropy Considerations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

MHD Free Convection and Entropy Generation in a Corrugated Cavity Filled with a Porous Medium Saturated with Nanofluids

by
Ali J. Chamkha
1,2 and
Fatih Selimefendigil
3,*
1
Mechanical Engineering Department, Prince Sultan Endowment for Energy and Environment, Prince Mohammad Bin Fahd University, Al-Khobar 31952, Saudi Arabia
2
RAK Research and Innovation Center, American University of Ras Al Khaimah, Ras al Khaimah 10021, UAE
3
Department of Mechanical Engineering, Celal Bayar University, Manisa 45140, Turkey
*
Author to whom correspondence should be addressed.
Entropy 2018, 20(11), 846; https://doi.org/10.3390/e20110846
Submission received: 11 October 2018 / Revised: 23 October 2018 / Accepted: 31 October 2018 / Published: 5 November 2018

Abstract

:
MHD free convection inside a triangular-wave-shaped corrugated porous cavity with Cu-water nanofluid is numerically studied with the finite element method. The influences of the Grashof number ( 10 4 Gr 10 6 ), Hartmann number ( 0 Ha 50 ), Darcy number ( 10 4 Da 10 1 ) and solid volume fraction of the nanoparticle ( 0 ϕ 0.05 ) on the convective flow features are examined. It is observed that increasing the Grashof number and Darcy number enhances the heat transfer, while the effect is opposite for the Hartmann number. As the corrugation frequency of the triangular wave increases, the local and averaged heat transfer rates decrease, which is more effective for higher values of Grashof and Darcy numbers. Normalized total entropy generation increases as the Darcy number and solid volume fraction of the nanoparticles increase and decreases as the Hartmann number increases both for flat and corrugated wall configurations.

1. Introduction

Natural convection in porous media has applications in some engineering problems such as in cooling of MEMs, solar collectors, geothermal applications, building insulation materials and many others. Industrial application examples are given in [1,2,3,4]. In some industrial applications such as MEMs, coolers of nuclear reactors and many other magnetic field effects may be relevant [5]. An external magnetic field has the potential to control the convective heat transfer features [6,7,8]. Nanofluid technology with the addition of nano-sized particles to the base fluid was applied in many thermal engineering problems [9,10]. Nanofluids enhance thermal transport [9]. Nanofluids can also be used in applications with magnetic fields, and both thermal conductivity and electrical conductivity of the base fluid are enhanced with the inclusion of the nanoparticles [11,12]. In the numerical study of Mahmoudi et al. [13], free convection in a triangular cavity with nanofluids was examined under the effects of a magnetic field, and the effects of various parameters such as the Rayleigh number, Hartmann number and nanoparticle volume fraction on the convective flow features were analyzed. Ghasemi et al. [14] numerically examined the free convection with a magnetic field for water: Al 2 O 3 nanofluid. They noted that depending on the Hartmann and Rayleigh number values, the average Nusselt number may be enhanced or deteriorated. In the literature, there exist several studies of convection in porous media filled with a nanofluid. Sun and Pop [15] examined the natural convection in a porous triangular cavity with nanofluids. They performed the study for different nanofluid types, and they reported that as the solid volume fraction of the nanoparticles increases, the heat transfer in the cavity is enhanced only at low Rayleigh numbers. Bourantas et al. [16] numerically examined the nanofluid natural convection in a porous cavity by using a meshless method, while the Darcy–Brinkman equations were used. The effects of nanofluid in the cooling were examined. Other relevant studies for convective heat transfer in a porous medium saturated by a nanofluid can be found in [17,18,19]. MHD flow in porous enclosures was studied by many researchers [20,21,22,23]. In applications for convection, the geometries can be simplified to some basic shapes such as square, triangular or trapezoidal. Surface corrugations can be used to control the convection in cavities. Hydro-thermal performance in a corrugated channel with nanofluid was numerically examined in [24]. As was shown by Hasan et al. [25], the corrugation amplitude and frequency have a significant impact on the heat transfer rate in natural convection for a deferentially-heated cavity. Hussain et al. [26] numerically examined the natural convection in a tilted sinusoidal corrugated cavity with magnetic field effects.
In convective heat transfer applications, second law analysis can be used to evaluate the system performance [27,28,29,30,31]. A review of the recent developments in entropy generation studies relevant to convective heat transfer is given in [32]. Some application examples of second law analysis with magnetic field and nanofluids can be found in [33,34,35,36,37].
As outlined in the above given literature survey, natural convection in a corrugated enclosure filled with a porous medium saturated with nanofluids under the effects of a magnetic field has many important thermal engineering applications. Some of the mentioned methods (surface corrugation, using nanofluids and an external magnetic field) can also be used as efficient control methods to affect the convective heat transfer features of the thermal engineering design for a porous medium. In this study, second law analysis of the present configuration with surface corrugation was also examined, which adds another novelty to the present study. The results of the present study may be useful in the thermal design of systems, which includes porous media, nanofluid and magnetic field interactions.

2. Mathematical Formulation

A schematic diagram of the problem for natural convection in a triangular wave-shaped corrugated cavity is shown in Figure 1. The enclosure size is H, while the corrugated wave has a length and height as a and b = 0.2 H . The cavity is deferentially heated with isothermal temperatures of T c and T h for the right and left corrugated vertical walls with T h > T c . Other walls are adiabatic with no-slip boundary conditions. A uniform vertical magnetic field is imposed. The nanofluid is assumed to saturate the solid matrix, and both are in thermodynamic equilibrium. Thermo-physical properties are shown in Table 1 [13]. Various effects such as joule heating, induced magnetic field and displacement currents are neglected. We assume that the magnetic Reynolds number is much smaller than the induced magnetic field so that it can be neglected when compared to the applied magnetic field.
Conservation equations and entropy generation in a two-dimensional system can be stated as follows [38,39]:
u x + v y = 0 ,
u u x + v u y = 1 ρ nf p x + ν nf 2 u x 2 + 2 u y 2 ν nf κ u ,
u v x + v v y = 1 ρ nf p y + ν nf 2 v x 2 + 2 v y 2 ν nf κ v + β nf g ( T T c ) σ nf B 0 2 ρ nf v ,
u T x + v T y = α nf 2 T x 2 + 2 T y 2 ,
S = k n f T 0 2 T x 2 + T y 2 + μ n f T 0 2 u x 2 + v y 2 + u x + v y 2 + μ n f T 0 κ u 2 + v 2 + σ n f B 0 2 T 0 u 2 .
In the above representation, entropy generation due to heat transfer, viscous dissipation and magnetic field is shown.
The effective density and specific heat of the nanofluid are defined as:
ρ n f = ( 1 ϕ ) ρ b f + ϕ ρ p , ( ρ c p ) n f = ( 1 ϕ ) ( ρ c p ) b f + ϕ ( ρ c p ) p .
The thermal expansion coefficient of the nanofluid is given as:
( ρ β ) n f = ( 1 ϕ ) ( ρ β ) b f + ϕ ( ρ β ) p .
The Maxwell–Garnett model and Maxwell model were used for the definition of the thermal conductivity and electrical conductivity of the nanofluid as:
k n f = k f ( k p + 2 k f ) 2 ϕ ( k f k p ) ( k p + 2 k f ) + ϕ ( k f k p ) , σ n f = σ f 1 + 3 σ p σ f 1 ϕ σ p σ f + 2 σ p σ f 1 ϕ .
The dynamic viscosity of the nanofluid is defined as [40]:
μ n f = μ b f ( 1 ϕ ) 0.25 .
Neither the viscosity nor the thermal conductivity of the nanofluid as described above include the effects of Brownian motion. Thermophoresis and Brownian motion in nanofluid flow were considered in many studies [41,42,43].
The following dimensionless parameters are used to convert the nondimensional form of equations:
X = x H , Y = y H , U = u H α , V = v H α , P = p ρ n f α 2 , θ = T T c T h T c , Gr = g β f ( T h T c ) H 3 ν f 2 , Pr = ν f α f , Ha = B 0 H σ n f ρ n f ν f , Da = κ H 2 .
Non-dimensional equations are written as follows:
U X + V Y = 0 ,
U U X + V U Y = P X + 1 Pr ν n f ν f 2 U X 2 + 2 U Y 2 Pr Da ν n f ν f U ,
U V X + V V Y = P Y + 1 Pr ν n f ν f 2 V X 2 + 2 V Y 2 , Pr Da ν n f ν f V + β n f β f Pr 2 Gr θ Ha 2 Pr V ,
U θ X + V θ Y = α n f α f 2 θ X 2 + 2 θ Y 2 .
Non-dimensional boundary conditions as:
  • At the bottom and top walls:
    U = V = 0 , θ Y = 0 ,
  • At the right vertical wall:
    U = V = 0 , θ = 0 ,
  • At the corrugated left wall:
    U = V = 0 , θ = 1 .
Local and average Nusselt numbers are calculated as [39]:
Nu s = k n f k f θ n n = 0 , Nu m = 0 1 Nu s d s .

3. Solver and Code Validation

The finite element method was used to obtain the solution of governing equations with boundary conditions as described above. The weak form of the governing equations is obtained, and when the approximated flow variables are substituted in the governing equations, residual R results, and the weighted average of R will be forced to be zero:
Ω W R d v = 0 ,
where W is the weight function, which is chosen from the same set of functions as of the trial functions. A system of nonlinear ordinary differential equations is achieved, which is solved with the Newton–Raphson method. Grid independence tests were performed, and the grid was refined near the walls to resolve the significant gradients in those locations. Table 2 shows the average Nusselt number for various grid sizes ( Gr = 10 6 , Ha = 50 , ϕ = 0.05 , Da = 10 2 ). Figure 2 demonstrates the grid distribution. The present solver is validated by using the results of Oztop and Abu-Nada [9] for cavity flow at Ra = 10 5 and the numerical results of Ahmed et al. [24] for corrugated channel flow at Re = 500 . The results in Figure 3 show sufficient accuracy of the current solver.

4. Results and Discussion

The effects of varying the Grashof number on the variation of flow and thermal patterns are shown in Figure 4 and Figure 5 ( Ha = 20 , Da = 10 2 , ϕ = 0.02 ). A single main recirculation cell is observed for all values of Grashof numbers. When the Grashof number value is augmented, the flow strength increases. At a low Grashof number, the isotherms are parallel to the vertical walls, indicating the characteristics of the conduction-dominated regime. By increasing the Grashof number, the flow strength increases, and temperature gradients near the hot and cold wall become significant. The effects of the corrugated wall on the streamlines and isotherms are also shown in Figure 4 and Figure 5. The global features of the flow and thermal patterns are not disturbed, except in the vicinity of the left wall. The streamlines and isotherms mimic the wall’s profile. The value of the maximum stream function increases with surface corrugation. The variation of the local and averaged Nusselt numbers for various Grashof numbers is shown in Figure 6 for flat and corrugated walls. As the Grashof number increases, heat transfer is locally enhanced along the heated wall for all configurations. The effect of the frequency of the triangular wave of the corrugated surface is also demonstrated in Figure 6. As is seen from the figure, the imposed frequency of the triangular wave affects the variation of the local heat transfer for the corrugated wall. As the frequency of the triangular wave increases, averaged heat transfer increases, and this effect is significant at higher Grashof numbers. This is due to the variation of the thermal and flow patterns adjacent to the boundary with respect to the change in the frequency. The increase in frequency reduces the space for the circulation adjacent to the boundary of the heated wall. This result is supported by Hasan et al. [25].
The influences of varying the Hartmann number on streamlines and isotherms ( Gr = 5 × 10 4 , Da = 10 2 , ϕ = 0.02 ) are demonstrated in Figure 7 and Figure 8. The strength of convection decreases as the Hartmann number increases for flat and corrugated wall configurations. When the value of the Hartmann number rises, a less dense arrangement of isotherms is obtained along the corrugated wall. The influence of the Hartmann number on the variation of the local and average Nusselt number is shown in Figure 9. Nusselt numbers decrease as the value of the Hartman number rises, which is due to the convective heat transfer suppression with a magnetic field, both for a flat wall and a corrugated wall configuration. These results are supported by [13,14].
Figure 10 and Figure 11 demonstrate the effects of varying the Darcy number on streamlines and isotherms for the flat wall and corrugated wall cases ( Gr = 10 4 , Ha = 10 , ϕ = 0.02 ). For low values of the Darcy number, fluid motion is not signification, and this can contribute to the high resistance of the medium. Convection becomes significant for higher values of the Darcy number with flat and corrugated wall configurations of the cavity. Conduction is dominant for low values of the Darcy number, which is demonstrated in Figure 11, while the isotherms are distributed approximately parallel to the vertical walls. The main cell is elongated horizontally as the Da number increases, as shown in Figure 10. An increment in Da results in increased temperature gradients near the heated wall, which can contribute to the enhanced permeability of the medium. In this case, the resistance from the boundary friction has been reduced with increasing Da number. Local and averaged heat transfer distributions for various Da numbers are shown in Figure 12. Local heat transfer is augmented for higher values of the Darcy number when both flat and corrugated wall configurations are considered. The effect of Da number on the average Nusselt number is less for the flat wall and corrugated wall with a = H / 4 and a = H / 8 for a low Da number due to the dominance of conduction. As the Da number increases, the effect of corrugation on the averaged heat transfer is more pronounced. The effect of nanofluid on the averaged heat transfer is demonstrated in Figure 13 ( Gr = 5 × 10 4 , Ha = 20 , Da = 10 2 ). Heat transfer enhancements of 34.4%, 30.7%, 30.4% and 33.5% are achieved for the nanofluid with a volume fraction of ϕ = 0.05 when it is compared to the base fluid ( ϕ = 0 ) for the flat wall and corrugated walls with a = H / 4 , a = H / 8 and a = H / 16 , respectively. The distribution of velocity and temperature along a horizontal plane located at y = 0.5 H are shown in Figure 14 ( Gr = 10 5 , Ha = 20 , Da = 10 2 , ϕ = 0.02 ) for flat and corrugated wall configurations. These profiles are affected from the corrugation frequency in the vicinity of the heated wall. Temperature increases with wall corrugation and its frequency. A shift and reduction in the peak values of the velocity components are seen with wall corrugation.
Normalized total entropy generation (normalized by the total entropy generation value for the first parameter) is shown in Figure 15 for various Darcy numbers, Hartmann numbers and nanoparticle volume fractions for both a flat wall and a corrugated wall with the a = H / 4 case. The total entropy generation increases as the value of the Darcy number increases due to the fluid friction irreversibility. When the value of the Hartmann number rises, the total entropy generation deteriorates, which is due to the suppression of the convection with the magnetic field strength. An increment in the solid volume fraction of nanoparticles resulted in the enhancement of the total entropy generation rate, which is due to the rise of heat transfer and fluid friction irreversibility. The discrepancy between the normalized entropy generation of the flat and corrugated wall cases increases as the values of the Darcy and Hartmann numbers rise.

5. Conclusions

A numerical study of MHD free convection in a triangular wave-shaped corrugated cavity filled with a porous medium saturated with nanofluids was analyzed. The following conclusions can be stated as:
  • Increasing the Grashof number and Darcy number results in heat transfer enhancement while keeping other parameters fixed.
  • The imposed frequency of the triangular wave affects the variation of the local Nusselt number for the corrugated wall. The velocity and temperature distributions are very sensitive to the corrugation frequency near the heated wall.
  • Increasing the corrugation frequency reduces the local and average Nusselt number.
  • The effect of the frequency of the triangular wave on the average Nusselt number is significant at higher values of the Grashof and Darcy numbers.
  • The magnetic field reduced the convective heat transfer.
  • Thirty-four-point-four percent of the average Nusselt number enhancement is obtained for the nanofluid with a nanoparticle volume fraction of ϕ = 0.05 when it is compared to the base fluid ( ϕ = 0 ) for the flat wall.
  • The normalized total entropy generation increases as the values of the Da number (due to the fluid friction irreversibility) and ϕ increase and the Ha number decreases (suppression of convection). The discrepancy between the normalized entropy generation for a flat and a corrugated wall is more pronounced for the highest values of the Darcy and Hartmann numbers.
An extension of this study may be considered to include the influence of corrugation amplitude, various corrugation types (rectangular, trapezoidal, sinusoidal, etc.), different nanofluid types and transient flow effects.

Author Contributions

F.S. performed the numerical simulations and wrote some sections of the manuscript. A.J.C. prepared some other sections of the paper and analyzed the results. All of the authors contributed equally to reviewing and revising the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations/Nomenclature

alength of triangular wave (m)
bheight of triangular wave (m)
B 0 magnetic field strength
DaDarcy number, κ H 2
GrGrashof number, g β f ( T h T c H 3 ν f 2
hlocal heat transfer coefficient (W/m 2 K)
HaHartmann number, B 0 H σ n f ρ n f ν f
kthermal conductivity (W/m·K)
Hlength of the enclosure (m)
nunit normal vector
Nulocal Nusselt number
ppressure (Pa)
Pnon-dimensional pressure
PrPrandtl number, ν f α f
Ttemperature (K)
u, vx-y velocity components (m/s)
U, Vdimensionless velocity components
x, yCartesian coordinates (m)
X, Ydimensionless coordinates
 Greek characters:
α thermal diffusivity (m 2 /s)
β expansion coefficient (1/K)
ϕ nanoparticle volume fraction
θ non-dimensional temperature,
T T c T h T c
ν kinematic viscosity (m 2 /s)
ρ density of the fluid (kg/m 3 )
σ electrical conductivity (S/m)
 Subscripts:
ccold wall
maverage
hhot wall

References

  1. Pop, I.; Ingham, D. Convective Heat Transfer: Mathematical and Computational Modeling of Viscous Fluids and Porous Media; Pergamon: Oxford, UK, 2001. [Google Scholar]
  2. Bejan, A.; Dincer, I.; Lorente, S.; Miguel, A.; Reis, A. Porous and Complex Flow Structures in Modern Technologies; Springer: New York, NY, USA, 2004. [Google Scholar]
  3. Vafai, K. Porous Media: Applications in Biological Systems and Biotechnology; CRC Press: New York, NY, USA, 2010. [Google Scholar]
  4. Sankar, M.; Kiran, S.; Ramesh, G.; Makinde, O. Natural convection in a non-uniformly heated vertical annular cavity. Defect Diffus. Forum 2017, 377, 189–199. [Google Scholar] [CrossRef]
  5. Nkurikiyimfura, I.; Wang, Y.; Pan, Z. Heat transfer enhancement by magnetic nanofluids—A review. Renew. Sustain. Energy Rev. 2013, 21, 548–561. [Google Scholar] [CrossRef]
  6. Hasanuzzaman, M.; Oztop, H.F.; Rahman, M.; Rahim, N.; Saidur, R.; Varol, Y. Magnetohydrodynamic natural convection in trapezoidal cavities. Int. Commun. Heat Mass Transf. 2012, 39, 1384–1394. [Google Scholar] [CrossRef]
  7. Selimefendigil, F.; Oztop, H.F. Effect of a rotating cylinder in forced convection of ferrofluid over a backward facing step. Int. J. Heat Mass Transf. 2014, 71, 142–148. [Google Scholar] [CrossRef]
  8. Yu, P.; Qiu, J.; Qin, Q.; Tian, Z.F. Numerical investigation of natural convection in a rectangular cavity under different directions of uniform magnetic field. Int. J. Heat Mass Transf. 2013, 67, 1131–1144. [Google Scholar] [CrossRef]
  9. Oztop, H.F.; Abu-Nada, E. Numerical study of natural convection in partially heated rectangular enclosures filled with nanofluids. Int. J. Heat Fluid Flow 2008, 29, 1326–1336. [Google Scholar] [CrossRef]
  10. Selimefendigil, F.; Oztop, H.F. Pulsating nanofluids jet impingement cooling of a heated horizontal surface. Int. J. Heat Mass Transf. 2014, 69, 54–65. [Google Scholar] [CrossRef]
  11. Sheikholeslami, M.; Bandpy, M.G.; Ganji, D. Numerical investigation of MHD effects on Al2O3-water nanofluid flow and heat transfer in a semi-annulus enclosure using LBM. Energy 2013, 60, 501–510. [Google Scholar] [CrossRef]
  12. Hamad, M.; Ismail, I.P.A. Magnetic field effects on free convection flow of a nanofluid past a vertical semi-infinite flat plate. Nonlinear Anal. Real World Appl. 2010, 12, 1338–1346. [Google Scholar] [CrossRef]
  13. Mahmoudi, A.H.; Pop, I.; Shahi, M. Effect of magnetic field on natural convection in a triangular enclosure filled with nanofluid. In. J. Therm. Sci. 2012, 59, 126–140. [Google Scholar] [CrossRef]
  14. Ghasemi, B.; Aminossadati, S.; Raisi, A. Magnetic field effect on natural convection in a nanofluid-filled square enclosure. Int. J. Therm. Sci. 2011, 50, 1748–1756. [Google Scholar] [CrossRef]
  15. Sun, Q.; Pop, I. Free convection in a triangle cavity filled with a porous medium saturated with nanofluids with flush mounted heater on the wall. Int. J. Therm. Sci. 2011, 50, 2141–2153. [Google Scholar] [CrossRef]
  16. Bourantas, G.; Skouras, E.; Loukopoulos, V.; Burganos, V. Heat transfer and natural convection of nanofluids in porous media. Eur. J. Mech. B-Fluids 2014, 43, 45–56. [Google Scholar] [CrossRef]
  17. Nielda, D.; Kuznetsov, A. Forced convection in a parallel-plate channel occupied by a nanofluid or a porous medium saturated by a nanofluid. Int. J. Heat Mass Transf. 2014, 70, 430–433. [Google Scholar] [CrossRef]
  18. Cheng, C.Y. Natural convection boundary layer flow over a truncated cone in a porous medium saturated by a nanofluid. Int. Commun. Heat Mass Transf. 2012, 39, 231–235. [Google Scholar] [CrossRef]
  19. Cimpean, D.S.; Pop, I. Fully developed mixed convection flow of a nanofluid through an inclined channel filled with a porous medium. Int. J. Heat Mass Transf. 2012, 55, 907–914. [Google Scholar] [CrossRef]
  20. Pekmen, B.; Sezgin, M.T. MHD flow and heat transfer in a lid-driven porous enclosure. Comput. Fluids 2014, 89, 191–199. [Google Scholar] [CrossRef]
  21. Guedda, M.; Ouahsine, A. Similarity solutions of MHD flows in a saturated porous medium. Eur. J. Mech. B-Fluids 2012, 33, 87–94. [Google Scholar] [CrossRef]
  22. Khan, M.; Fetecau, C.; Hayat, T. MHD transient flows in a channel of rectangular cross-section with porous medium. Phys. Lett. A 2007, 369, 44–54. [Google Scholar] [CrossRef]
  23. Makinde, O.; Aziz, A. MHD mixed convection from a vertical plate embedded in a porous medium with a convective boundary condition. Int. J. Therm. Sci. 2010, 49, 1813–1820. [Google Scholar] [CrossRef]
  24. Ahmed, M.; Shuaib, N.; Yusoff, M.; Al-Falahi, A. Numerical investigations of flow and heat transfer enhancement in a corrugated channel using nanofluid. Int. Commun. Heat Mass Transf. 2011, 38, 1368–1375. [Google Scholar] [CrossRef]
  25. Hasan, M.N.; Saha, S.C.; Gu, Y. Unsteady natural convection within a differentially heated enclosure of sinusoidal corrugated side walls. Int. J. Heat Mass Transf. 2012, 55, 5696–5708. [Google Scholar] [CrossRef] [Green Version]
  26. Hussain, S.H.; Hussein, A.K.; Mohammed, R.N. Studying the effects of a longitudinal magnetic field and discrete isoflux heat source size on natural convection inside a tilted sinusoidal corrugated enclosure. Comput. Math. Appl. 2012, 64, 476–488. [Google Scholar] [CrossRef]
  27. Rashidi, S.; Akar, S.; Bovand, M.; Ellahi, R. Volume of fluid model to simulate the nanofluid flow and entropy generation in a single slope solar still. Renew. Energy 2018, 115, 400–410. [Google Scholar] [CrossRef]
  28. Basak, T.; Gunda, P.; Anandalakshmi, R. Analysis of entropy generation during natural convection in porous right-angled triangular cavities with various thermal boundary conditions. Int. J. Heat Mass Transf. 2012, 55, 4521–4535. [Google Scholar] [CrossRef]
  29. Eegunjobi, A.S.; Makinde, O.D. Inherent irreversibility in a variable viscosity Hartmann flow through a rotating permeable channel with Hall effects. Defect Diffus. Forum 2017, 377, 180–188. [Google Scholar] [CrossRef]
  30. Qasim, M.; Khan, Z.H.; Khan, I.; Al-Mdallal, Q.M. Analysis of Entropy Generation in Flow of Methanol-Based Nanofluid in a Sinusoidal Wavy Channel. Entropy 2017, 19, 490. [Google Scholar] [CrossRef]
  31. Abbas, T.; Ayub, M.; Bhatti, M.M.; Rashidi, M.M.; Ali, M.E.S. Entropy Generation on Nanofluid Flow through a Horizontal Riga Plate. Entropy 2016, 18, 223. [Google Scholar] [CrossRef]
  32. Oztop, H.F.; Al-Salem, K. A review on entropy generation in natural and mixed convection heat transfer for energy systems. Renew. Sustain. Energy Rev. 2012, 16, 911–920. [Google Scholar] [CrossRef]
  33. Bhatti, M.M.; Abbas, T.; Rashidi, M.M.; Ali, M.E.S.; Yang, Z. Entropy Generation on MHD Eyring-Powell Nanofluid through a Permeable Stretching Surface. Entropy 2016, 18, 224. [Google Scholar] [CrossRef]
  34. Ellahi, R.; Alamri, S.Z.; Basit, A.; Majeed, A. Effects of MHD and slip on heat transfer boundary layer flow over a moving plate based on specific entropy generation. J. Taibah Univ. Sci. 2018, 12, 476–482. [Google Scholar] [CrossRef]
  35. Mahian, O.; Oztop, H.; Pop, I.; Mahmud, S.; Wongwises, S. Entropy generation between two vertical cylinders in the presence of MHD flow subjected to constant wall temperature. Int. Commun. Heat Mass Transf. 2013, 44, 87–92. [Google Scholar] [CrossRef]
  36. Mahian, O.; Kianifar, A.; Kleinstreuer, C.; Al-Nimr, M.A.; Pop, I.; Sahin, A.Z.; Wongwises, S. A review of entropy generation in nanofluid flow. Int. J. Heat Mass Transf. 2013, 65, 514–532. [Google Scholar] [CrossRef]
  37. Bhatti, M.M.; Abbas, T.; Rashidi, M.M.; Ali, M.E.S. Numerical Simulation of Entropy Generation with Thermal Radiation on MHD Carreau Nanofluid towards a Shrinking Sheet. Entropy 2016, 18, 200. [Google Scholar] [CrossRef]
  38. Rashad, A.M.; Armaghani, T.; Chamkha, A.J.; Mansoure, M.A. Entropy generation and MHD natural convection of a nanofluid in an inclined square porous cavity: Effects of a heat sink and source size and location. Chin. J. Phys. 2018, 56, 193–211. [Google Scholar] [CrossRef]
  39. Sheikholeslami, M.; Li, Z.; Shamlooei, M. Nanofluid MHD natural convection through a porous complex shaped cavity considering thermal radiation. Phys. Lett. A 2018, 382, 1615–1632. [Google Scholar] [CrossRef]
  40. Brinkman, H. The viscosity of concentrated suspensions and solutions. J. Chem. Phys. 1952, 20, 571–581. [Google Scholar] [CrossRef]
  41. Sheikholeslami, M.; Gorji-Bandpy, M.; Ganji, D.D.; Rana, P.; Soleimani, S. Magnetohydrodynamic free convection of Al2O3-water nanofluid considering Thermophoresis and Brownian motion effects. Comput. Fluids 2014, 94, 147–160. [Google Scholar] [CrossRef]
  42. Makinde, O.D.; Animasaun, I.L. Thermophoresis and Brownian motion effects on MHD bioconvection of nanofluid with nonlinear thermal radiation and quartic chemical reaction past an upper horizontal surface of a paraboloid of revolution. J. Mol. Liquids 2016, 221, 733–743. [Google Scholar] [CrossRef]
  43. Sheikholeslami, M.; Ganji, D.D. Nanofluid flow and heat transfer between parallel plates considering Brownian motion using DTM. Comput. Methods Appl. Mech. Eng. 2015, 283, 651–663. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of the physical model with boundary conditions.
Figure 1. Schematic diagram of the physical model with boundary conditions.
Entropy 20 00846 g001
Figure 2. Grid distribution.
Figure 2. Grid distribution.
Entropy 20 00846 g002
Figure 3. Code verification with the results of Oztop and Abu-Nada [9] (a,b) and Ahmed et al. [24] (c).
Figure 3. Code verification with the results of Oztop and Abu-Nada [9] (a,b) and Ahmed et al. [24] (c).
Entropy 20 00846 g003
Figure 4. The effect of varying the Grashof number on the streamlines for flat and corrugated walls ( Ha = 20 , Da = 10 2 , ϕ = 0.02 ).
Figure 4. The effect of varying the Grashof number on the streamlines for flat and corrugated walls ( Ha = 20 , Da = 10 2 , ϕ = 0.02 ).
Entropy 20 00846 g004
Figure 5. The effect of varying the Grashof number on the isotherms for flat and corrugated walls ( Ha = 20 , Da = 10 2 , ϕ = 0.02 ).
Figure 5. The effect of varying the Grashof number on the isotherms for flat and corrugated walls ( Ha = 20 , Da = 10 2 , ϕ = 0.02 ).
Entropy 20 00846 g005aEntropy 20 00846 g005b
Figure 6. Local and averaged Nusselt number along the flat and corrugated walls for various Grashof numbers ( Ha = 20 , Da = 10 2 , ϕ = 0.02 ).
Figure 6. Local and averaged Nusselt number along the flat and corrugated walls for various Grashof numbers ( Ha = 20 , Da = 10 2 , ϕ = 0.02 ).
Entropy 20 00846 g006
Figure 7. The effect of varying the Hartmann number on the streamlines for flat and corrugated walls ( Gr = 5 × 10 4 , Da = 10 2 , ϕ = 0.02 ).
Figure 7. The effect of varying the Hartmann number on the streamlines for flat and corrugated walls ( Gr = 5 × 10 4 , Da = 10 2 , ϕ = 0.02 ).
Entropy 20 00846 g007
Figure 8. The effect of varying the Hartmann number on the isotherms for flat and corrugated walls ( Gr = 5 × 10 4 , Da = 10 2 , ϕ = 0.02 ).
Figure 8. The effect of varying the Hartmann number on the isotherms for flat and corrugated walls ( Gr = 5 × 10 4 , Da = 10 2 , ϕ = 0.02 ).
Entropy 20 00846 g008aEntropy 20 00846 g008b
Figure 9. Local and averaged Nusselt number along the flat and corrugated wall for various Hartmann numbers ( Gr = 5 × 10 4 , Da = 10 2 , ϕ = 0.02 ).
Figure 9. Local and averaged Nusselt number along the flat and corrugated wall for various Hartmann numbers ( Gr = 5 × 10 4 , Da = 10 2 , ϕ = 0.02 ).
Entropy 20 00846 g009
Figure 10. The effect of varying the Darcy number on the streamlines for flat and corrugated walls ( Gr = 10 4 , Ha = 10 , ϕ = 0.02 ).
Figure 10. The effect of varying the Darcy number on the streamlines for flat and corrugated walls ( Gr = 10 4 , Ha = 10 , ϕ = 0.02 ).
Entropy 20 00846 g010aEntropy 20 00846 g010b
Figure 11. The effect of varying the Darcy number on the isotherms for flat and corrugated walls ( Gr = 10 4 , Ha = 10 , ϕ = 0.02 ).
Figure 11. The effect of varying the Darcy number on the isotherms for flat and corrugated walls ( Gr = 10 4 , Ha = 10 , ϕ = 0.02 ).
Entropy 20 00846 g011
Figure 12. Local and averaged Nusselt number along the flat and corrugated walls for various Darcy numbers ( Gr = 10 4 , Ha = 10 , ϕ = 0.02 ).
Figure 12. Local and averaged Nusselt number along the flat and corrugated walls for various Darcy numbers ( Gr = 10 4 , Ha = 10 , ϕ = 0.02 ).
Entropy 20 00846 g012
Figure 13. Averaged Nusselt number along the flat and corrugated walls for various nanoparticle volume fractions ( Gr = 5 × 10 4 , Ha = 20 , Da = 10 2 ).
Figure 13. Averaged Nusselt number along the flat and corrugated walls for various nanoparticle volume fractions ( Gr = 5 × 10 4 , Ha = 20 , Da = 10 2 ).
Entropy 20 00846 g013
Figure 14. Velocity and temperature profiles along the horizontal plane at y = 0.5 H for flat and corrugated walls ( Gr = 10 5 , Ha = 20 , Da = 10 2 , ϕ = 0.02 ).
Figure 14. Velocity and temperature profiles along the horizontal plane at y = 0.5 H for flat and corrugated walls ( Gr = 10 5 , Ha = 20 , Da = 10 2 , ϕ = 0.02 ).
Entropy 20 00846 g014
Figure 15. The effects of varying the Darcy number (a) Hartmann number (b) nanoparticle volume fraction (c) on the normalized total entropy generation within the cavity.
Figure 15. The effects of varying the Darcy number (a) Hartmann number (b) nanoparticle volume fraction (c) on the normalized total entropy generation within the cavity.
Entropy 20 00846 g015
Table 1. Thermophysical properties.
Table 1. Thermophysical properties.
PropertyWaterCu
ρ 997.18954
c p 4179383
k0.6400
β 2 . 1 × 10 4 1 . 67 × 10 5
σ 0.05 5 . 97 × 10 7
Table 2. Grid independence test.
Table 2. Grid independence test.
Grid NameNumber of ElementsAveraged Nusselt Number
G111369.232
G227969.983
G3454410.483
G418,17611.087
G552,62411.184

Share and Cite

MDPI and ACS Style

Chamkha, A.J.; Selimefendigil, F. MHD Free Convection and Entropy Generation in a Corrugated Cavity Filled with a Porous Medium Saturated with Nanofluids. Entropy 2018, 20, 846. https://doi.org/10.3390/e20110846

AMA Style

Chamkha AJ, Selimefendigil F. MHD Free Convection and Entropy Generation in a Corrugated Cavity Filled with a Porous Medium Saturated with Nanofluids. Entropy. 2018; 20(11):846. https://doi.org/10.3390/e20110846

Chicago/Turabian Style

Chamkha, Ali J., and Fatih Selimefendigil. 2018. "MHD Free Convection and Entropy Generation in a Corrugated Cavity Filled with a Porous Medium Saturated with Nanofluids" Entropy 20, no. 11: 846. https://doi.org/10.3390/e20110846

APA Style

Chamkha, A. J., & Selimefendigil, F. (2018). MHD Free Convection and Entropy Generation in a Corrugated Cavity Filled with a Porous Medium Saturated with Nanofluids. Entropy, 20(11), 846. https://doi.org/10.3390/e20110846

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