Next Article in Journal
High-Temperature Nano-Indentation Creep of Reduced Activity High Entropy Alloys Based on 4-5-6 Elemental Palette
Next Article in Special Issue
Radiative MHD Nanofluid Flow over a Moving Thin Needle with Entropy Generation in a Porous Medium with Dust Particles and Hall Current
Previous Article in Journal
Non-Monogamy of Spatio-Temporal Correlations and the Black Hole Information Loss Paradox
Previous Article in Special Issue
Fluid Flow and Entropy Generation Analysis of Al2O3–Water Nanofluid in Microchannel Plate Fin Heat Sinks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Simulation on Convection and Thermal Radiation of Casson Fluid in an Enclosure with Entropy Generation

1
Department of Mathematics, King Abdulaziz University, Jeddah 21589, Saudi Arabia
2
Department of Mathematics, Kongunadu Polytechnic College, D.Gudalur, Dindigul, Tamilnadu 624620, India
*
Author to whom correspondence should be addressed.
Submission received: 14 January 2020 / Revised: 6 February 2020 / Accepted: 6 February 2020 / Published: 18 February 2020
(This article belongs to the Special Issue Thermal Radiation and Entropy Analysis)

Abstract

:
The goal of the current numerical simulation is to explore the impact of aspect ratio, thermal radiation, and entropy generation on buoyant induced convection in a rectangular box filled with Casson fluid. The vertical boundaries of the box are maintained with different constant thermal distribution. Thermal insulation is executed on horizontal boundaries. The solution is obtained by a finite volume-based iterative method. The results are explored over a range of radiation parameter, Casson fluid parameter, aspect ratio, and Grashof number. The impact of entropy generation is also examined in detail. Thermal stratification occurs for greater values of Casson liquid parameters in the presence of radiation. The kinetic energy grows on rising the values of Casson liquid and radiation parameters. The thermal energy transport declines on growing the values of radiation parameter and it enhances on rising the Casson fluid parameter.

1. Introduction

Since convective flow with heat energy transfer plays a major role in numerous places, research on this topic has been executed over the last few decades [1,2,3,4,5]. However, the research on convective flow of Casson fluids is not active up to the level of current scientific and technological requirements. Therefore, we investigate the energy transport and flow of Casson fluid under different physical environments due to its applications in electronics manufacturing, printing, textile printing, semi-solid casting, etc. [6,7,8]. The Casson liquid is one kind of non-Newtonian liquid and it is a shear thinning (pseudo-plastic) liquid which behaves like an elastic solid [5]. Pop and Sheremet [9] explored the convective stream of Casson fluid in an enclosed space in the occurrence of radiation and viscous dissipation. They established that the Casson parameter favors energy transferal augmentation and stream intensification. Sheremet and Pop [10] explored the convective flow with thermal radiation of a visco-elastic fluid in a square domain. They attained the energy reduction while increasing the elastic number. Srinivas et al. [11] inspected the influence of thermal radiation and chemical reaction on hydromagnetic pulsating Casson liquid flow in a porous channel. The study on buoyant convective stream in the existence of thermal radiation has been enriched in research due to applications in various fields [12,13,14,15]. Mansour et al. [16] observed the impact of thermal radiation on free convection in a wavy box. They established that averaged Nusselt number rises by raising the radiation values. Mahapatra et al. [17] surveyed the combined radiation and convection in a porous box in the occurrence of heat generation. Miroshnichenko et al. [18] inspected the combined radiation and convective motion in enclosures with a heat-generating component.
The entropy generation is the amount of entropy which is produced in any irreversible processes, such as heat/mass transfer, viscous fluid flow, chemical reactions, and Joule heating [19,20]. Every process in nature provides a positive or zero value of entropy generation rate. It is an essential aspect of the 2nd law of thermodynamics [21,22]. The study of entropy generation in closed domain has been analyzed in several studies [23,24,25]. Ahmed et al. [26] analyzed the entropy production due to doubly diffusive convection of Casson liquids in the occurrence of slip boundaries and chemical reactions. They observed that the Casson parameter diminishes the entropy production. Hajji et al. [27] explored the entropy production due to combined free convection and radiative transfer in an inclined box. They concluded that entropy generation rises by increasing wall emissivity due to surface radiation. Daniel et al. [28] explored mixed magneto-convection of chemically reacting nanoliquid with the combined influence of thermal radiation and viscous dissipation. The entropy generation on doubly diffusive convection of Carreau fluid with Dufour and Soret effects and viscous dissipation in a box was investigated by Kefayati and Tang [29]. They found that the growth of Dufour and Soret parameters enhance the entropy production. In another study, Kefayati [30] explored the entropy production due to doubly diffusive convection of Bingham fluid in a tilting box. Chamkha et al. [31] investigated the entropy propagation on MHD convection of nanoliquid in a porous driven box in the existence of partial slip under various heat sinks and sources.
It is understood from the literature that there are several studies involving the convective flow with thermal radiation of Newtonian fluids. In addition, a few studies dealt with the convective stream of Casson liquid in the absence of radiation. This is motivated to initiate the present study to inspect the interaction of thermal radiation and convection of Casson fluid inside a rectangular box. Hence, the numerical model on convective current of Casson liquid in a box in the occurrence of thermal radiation with entropy generation is explored here and it has not been described in the literature.

2. Mathematical Model

Buoyant convective laminar incompressible flow in a 2-dimensional box of height H and width L filled with Casson liquid is considered in the current study, as given in Figure 1. The left wall has the constant higher temperature ( T ( y ) = T h ) and right wall has constant lower temperature ( T ( y ) = T c ) . The thermal insulation is imposed on horizontal walls. The gravity acts in the negative y-direction. The thermal radiation is taken into account. The medium is optically thick. The thermal properties of the liquid are kept invariable except the density in the buoyancy term (Boussinesq approximation). The density is defined as ρ = ρ 0 [ 1 β T ( T T c ) ] , where βT is thermal expansion coefficient. No viscous dissipation or Joule heating are counted.
The rheological model of Casson fluid is derived as below [9,10]:
τ i j = { 2 ( μ B + P y 2 π ) e i j , π > π C 2 ( μ B + P y 2 π C ) e i j , π < π C
where π = e i j . e i j and e i j are (i, j)th component of the deformation rate, π c is a critical rate of the product of deformation rate, μ B and P y is dynamic viscosity, yield stress of liquid. The Casson liquid satisfies the Boussinesq approximation. According to aforementioned assumptions, the principal equations are
u   x + v y = 0
u t + [ u u x + v u y ] = 1 ρ 0 p x + ν ( 1 + 1 β ) [ 2 u x 2 + 2 u y 2 ]
v t + [ u v x + v v y ] = 1 ρ 0 p y + ν ( 1 + 1 β ) [ 2 v x 2 + 2 v y 2 ]   + g β ( T T c )
T t + u T x + v T y = k ρ 0 c p [ 2 T x 2 + 2 T y 2 ] 1 ρ 0 c p [ q r x + q r y ]
where β is the Casson liquid parameter. The heat flux due to radiation along x and y directions are given by
q r x = 4 σ * 3 K T 4 x   and   q r y = 4 σ * 3 K T 4 y
Using Rosseland approximation for radiation, the temperature differences within the stream are too small. Expanding T 4 about T c using Taylor series and approximating it. Finally, we get
T 4 4 T c 3 T 3 T c 4
Therefore, heat flux due to radiation along x and y directions (Equation (6)) reduces to
q r x = 16 σ * T c 3 3 K T x   and     q r y = 16 σ * T c 3 3 K T y
Substituting Equation (7) into Equation (5), we get
T t + u T x + v T y = k ρ 0 c p [ 2 T x 2 + 2 T y 2 ] + 1 ρ 0 c p 16 σ * T c 3 3 K [ 2 T x 2 + 2 T y 2 ]
The relevant initial and boundary settings are
t = 0 :   u = v = 0 ,   T = T 0 0 x L ,   0 y H t > 0 :   u = v = 0 , x = 0 & L , y = 0 & H T ( y ) = T h x = 0 T ( y ) = T c x = L T y = 0 y = 0 & L
The dimensionless form of governing model is done by the following variables:
X = x L ,   Y = y L ,   U = u L ν ,   V = v L ν ,   θ = T T c ( T h T c ) ,   τ = t ν L 2 ,   and   P = p L 2 ρ 0 ν 2
The non-dimensional equations are
U X + V Y = 0
U τ + [ U U X + V U Y ]   =   P X +   ( 1 + 1 β ) 2 U
V τ + [ U V X + V V Y ]   =   P Y +   ( 1 + 1 β ) 2 V + G r   θ
θ τ + U θ X + V θ Y = 1 P r ( 1 + 4 3 R d ) [ 2 θ X 2 + 2 θ Y 2 ]
The non-dimensional parameters in the equations are A r = H L aspect ratio, G r = g β ( T h T c ) L 3 ν 2 the Grashof number, P r = ν α   ( = 10 ) , the Prandtl number, and the radiation parameter, R d = 4 σ * T 0 3 k K .
The initial and border settings (in dimensionless form) are
τ = 0 :   V = U = 0   θ = 0 0 X 1 ,   0 Y A r , τ > 0 :   V = U = 0   θ Y = 0 Y = 0   &   A r , V = U = 0     θ = 1 X = 0 , V = U = 0   θ = 0 X = 1 .
The thermal energy transferal rate across the box is an essential factor in the heat removal applications. The Nusselt number is the ratio of convection to conduction thermal transport across the box. The local Nusselt number along the left-wall of the box is derived as
N u = q w L k ( Δ T )   where   q w = ( k T x + 16 σ * T c 3 3 K T x ) .
The non-dimensional form of local Nusselt number is derived as follows:
N u = ( 1 + 4 3 R d ) ( θ X ) X = 0 .
The total heat energy transfer across the box is computed by the averaged Nusselt number, which is derived as follows:
N u ¯ = 1 A r 0 A r N u d Y .
The drag force on the surface is explored by the skin friction and it is defined as C f l o c = u x . The averaged skin friction factor at left-wall is estimated by integrating C f l o c over the length of the left-wall as follows:
C f y = 1 A r 0 A r C f l o c d Y .

3. Entropy Generation Analysis

The buoyant convection in a closed box presents significant possibilities for thermal industrial applications. Though, the usage of entropy generation analysis assists to detect the optimal conditions for real-world applications. Since the entropy production is due to the irreversible process of heat transfer and viscosity effects, entropy production can be calculated from the known velocity and thermal fields.
The local entropy generation per unit area can be expressed by two quantities. That is, the local entropy generation subject to heat transference and due to fluid friction [24,25,26,27,28].
S h e a t = k T c 2 [ ( T x ) 2 + ( T y ) 2 ]
S f l u i d = ( μ T c ) ( 1 + 1 β ) { 2 [ ( u x ) 2 + ( v y ) 2 ] + ( u y + v x ) 2 }
That is, the total local entropy production is sum of the above two quantities.
S G e n = k T c 2 [ ( T x ) 2 + ( T y ) 2 ] + ( μ T c ) ( 1 + 1 β ) { 2 [ ( u x ) 2 + ( v y ) 2 ] + ( u y + v x ) 2 }
The dimensionless form for entropy generation is obtained in a usual way by using Equation (10)
S t o t a l = S h e a t * + S f l u i d *
S h e a t * = ( θ X ) 2 + ( θ Y ) 2
S f l u i d * = ϕ 2 ( 1 + 1 β ) { 2 [ ( U X ) 2 + ( V Y ) 2 ] + ( U Y + V X ) 2 }
The overall entropy generation is attained by integrating the local entropy production over the box.
S G t o t a l = V S t o t a l ( X , Y ) d A
The local Bejan number specifies the strength of the entropy production subject to heat transference irreversibility. It is defined as
B e l o c = S h e a t * S t o t a l
For any point inside the box, when B e l o c > 1 2 , the heat transfer irreversibility is dominating. When B e l o c < 1 2 , the liquid friction irreversibility dominates. If B e l o c = 1 2 , the viscous and thermal irreversibilities are equal. The averaged Bejan number is used to define the relative significance of the heat energy transfer irreversibility for the entire cavity.
B e = A B e l o c ( X , Y ) d A A d A
The computer program built-up for the flow and thermal fields is extended to calculate the entropy generation inside the box.

4. Cup Mixing Temperature and RMSD

To explore the thermal mixing in the enclosed box, cup mixing temperature is defined. The velocity-weighted averaged temperature is more suitable for convective flow than spatial average temperature. The cup mixing temperature (Tcup) and area averaged temperature (Tavg) are given as [32]
T C u p = V ^ ( X , Y )   θ ( X , Y ) d X d Y V ^ ( X , Y ) d X d Y
where V ^ ( X , Y ) = U 2 + V 2 and
T a v g = θ ( X , Y ) d X d Y d X d Y
Also, root mean square deviation (RMSD) is derived to measure the degree of temperature regularity in all cases. The RMSDs are derived based on cup mixing temperature and area averaged temperature as follows:
R M S D T c u p = i = 1 N ( θ i T C u p ) 2 N
R M S D T a v g = i = 1 N ( θ i T a v g ) 2 N
The higher values of RMSD indicate lower temperature uniformity inside the box and vice-versa. Also, RMSD should not exceed 1, because the dimensionless temperature varies between 0 and 1. The above parameters are calculated using the obtained values of thermal and flow fields in the same computer program.

5. Numerical Technique and Validation

The non-dimensional model Equations (11)–(14) subject to the boundary settings (Equation (15)), are solved by the control volume method [2] using a 122 × 122 non-uniform grid. The discretization of diffusive and convective terms are taken by central difference and QUICK schemes. The iterative process is executed by using a tri-diagonal matrix. The convergence condition of the iteration is kept as 10 −6. The justification of current computer program is tested against the existing results for buoyant convection in a box [33,34], see Table 1. An agreement between the results provides assurance in the accurateness of the current code to explore the problem.

6. Results and Discussion

The numerical investigation is done to conclude the influence of the thermal radiation on natural convection flow in a Casson fluid saturated porous square box. The controlling parameters for this investigation are the Grashof number ( G r = 10 4 & 10 6 ), the aspect ratio ( 0.25 A r 5 ), the Casson fluid parameter ( 0.01 β 1 ), and the radiation parameter ( 0 R d 10 ). The calculations have been computed for diverse combinations of above said parameters to analyze various effects.
The stream pattern for various values of Casson liquid parameter and thermal radiation parameter with G r = 10 6 , A r = 1 is displayed in Figure 2. A clockwise circling eddy occupies the whole box for all values of the considered parameters. The core region of the eddy is elongated horizontally and lies in the middle of the box in the nonappearance of thermal radiation for β = 0.01. On raising the thermal radiation intensity, the shape of the core region turns from elliptical shape to circular shape. However, this pattern is not found on raising the values of Casson fluid parameter to β = 0.1 to β = 1. That is, the viscous diffusion dominates and it slowed down the flow. The core region is divided into two parts for all values of thermal radiation parameter when β = 1 case. However, the similar trend exists for weak thermal radiation case when β = 0.1 and it disappears on strong thermal radiation case. That is, the stream speed is reduced on growing the values of β. Figure 3 displays the thermal distributions inside the box for several values of β and Rd with G r = 10 6 ,     A r = 1 . The thermal layers at boundary are molded along the isothermal walls for all given values of Casson liquid parameter in the nonappearance of radiation parameter. The thermal boundary layers are feeble in the presence of radiation parameter and they disappear for strong values of radiation parameter (Rd = 10). The vertical temperature stratification is found for all Rd values when β = 1. The thermal layers near boundary are stronger on rising the values of β.
Figure 4a–d demonstrates the impact of skin friction ( C f y ) of the liquid for various values of thermal radiation, aspect ratio, and Casson fluid parameters. It is visibly understood from Figure 4a,b that skin friction rises on growing the values of Casson fluid parameter. Since the liquid friction inside the liquid rises on growing the Casson liquid parameter, it produces more friction between the liquid molecules and it results in resistant force being exerted on moving liquid due to drag caused by the viscosity of liquids. Comparing Figure 4a and Figure 4b, the skin friction rises in large amount by increasing Gr. The influence of aspect ratio on skin friction is clearly demonstrated in Figure 4c,d. The drag force rises while the aspect ratio of the box rises. The friction also grows when rising the values of radiation parameter. Figure 5a–d demonstrates the impact of kinematic energy (KE) of the liquid for various values of thermal radiation, aspect ratio, and Casson fluid parameters. It is visibly understood from Figure 5a,b that the kinematic energy rises on growing the values of Casson liquid parameter. However, the increment of KE augments on rising the values of the Gr number. When rising the radiation parameter to Rd = 1, 5, and 10, the KE of the liquid particle increases.
Figure 6 indicates the local Nusselt number for various radiation parameter with β = 0.1 and 1, G r = 10 6 , and A r = 1 . The shapes of local Nusselt number curve visibly establish the impact of rate of the local thermal transport from the wall. The maximum local energy transfer is attained near the upper portion of the left wall. It is spotted that the local energy transference enhances on rising the values of radiation parameter (Rd). The local energy transport raises first and then declines along the wall for a given Rd value. Figure 7 depicts the impact of local thermal energy transport for several values of Casson liquid parameter (β) with Rd=0 & 5, G r = 10 6 , and Ar = 1. The local energy transfer increases first along the height of warm wall up to y=0.05 and then it drops along y direction. It is also achieved the highest local energy transfer at the topmost area of the hot wall. The local energy transfer increases with β.
To determine the influence of thermal energy transference across the domain, the mean Nusselt number is graphed in Figure 8a–d against Rd and Ar for several values of radiation and Casson liquid parameters. It is clearly seen from these figures that the thermal energy transmission enhances on rising the values of Rd. The averaged Nusselt number rises rapidly for lower values of Ar and then falls slightly for higher values of Ar for a given Rd and β. It is also acquired that the mean energy transfer enhances on rising the values of Casson liquid parameter for all Rd values. The mean energy transport enhances on growing the values of radiation parameter. At Rd = 0, the averaged Nusselt number is almost constant when changing the values of Ar. Figure 9a–d depicted the influence of Bejan number for diverse values of Casson fluid and radiation parameters and aspect ratio. The Bejan number behaves nonlinearly with Rd for lower values of Gr = 104. However, the Bejan number declines when raising the values of Rd for Gr = 106. There is no difference on Bejan number when changing the values of Casson fluid parameter at high aspect ratios (Ar ≥ 3).
Figure 10a–d depicted the impact of cup mixing temperature inside the box with various values of β, Rd, and Ar. The cup mixing temperature enhances while reducing the Casson fluid parameter for Gr = 106. However, the opposite trend is found for Gr = 104 in Figure 10a. The cup mixing temperature augments when increasing the thermal radiation for Gr = 106, however, it declines when increasing the radiation for low values of Gr and low values of β (≤0.1). It is seen that the cup mixing temperature behaves nonlinearly with aspect ratio. That is, it is increasing rapidly initially and then declines gradually on raising the size of the box (aspect ratio). Figure 11a–d portrayed the impression of average temperature inside the box with various values of β, Rd, and Ar. The average temperature enhances while reducing the Casson fluid parameter. The average temperature augments when increasing the thermal radiation for Gr = 104 for all values of Casson fluid parameter, however, it declines when increasing the radiation for higher values of Gr except β = 0.01. It is seen that the average temperature behaves nonlinearly with aspect ratio. That is, it is declining rapidly first (until Ar = 1) and then rises gradually on raising the aspect ratio of the box. The tall and slender cavities behave in different manner.
Figure 12a–d depicted the influence of RMSDTcup inside the box with various values of β, Rd, and Ar. The RMSDTcup enhances while reducing the Casson fluid parameter. The RMSDTcup augments by increasing the thermal radiation for all values of Casson fluid parameter and Gr. It is perceived that the RMSDTcup temperature performs nonlinearly with aspect ratio. RMSDTcup declines rapidly first (until Ar = 1/2) for some cases and then rises gradually on raising the aspect ratio of the box. Figure 13a–d demonstrated the impact of RMSDTavg inside the box with various values of β, Rd, and Ar. The RMSDTavg augments while the Casson fluid parameter falls. The RMSDTavg enhances by rising the thermal radiation for all values of Casson fluid parameter and Gr. It is perceived that the RMSDTavg achieves nonlinear fashion with aspect ratio. RMSDTavg declines rapidly first (until Ar = 0.5) for some cases and then rises gradually on raising the aspect ratio of the box.

7. Conclusions

The impact of aspect ratio and entropy generation analysis on buoyant convective flow of Casson liquid in a rectangular box has been performed numerically. The following conclusions are arrived from this study:
Strong thermal layers at boundary are formed along the thermal walls.
Thermal stratification found for higher values of β (= 1) for all values of radiation parameter.
Skin friction develops with aspect ratio, thermal radiation, and Casson fluid parameter.
The kinetic energy enhances with aspect ratio, thermal radiation, and Casson fluid parameter.
Averaged heat transfer enhances with thermal radiation and Casson fluid parameter. However, it increases first and then declines when growing the aspect ratio of the box.
The Bejan number enhances with Casson fluid parameter and declines with Ar.
The cup mixing and average temperature behaves in a nonlinear fashion with aspect ratio of the box.
The RMSDTavg augments while the Casson fluid parameter falls and it enhances by rising the thermal radiation.
The RMSDTcup enhances while reducing the Casson fluid parameter and it augments by growing the thermal radiation.

Author Contributions

Conceptualization, S.S.; methodology, S.S.; software, S.S.; validation, S.S. and M.B.; investigation, S.S. and M.B.; writing—original draft preparation, S.S. and M.B.; writing—review and editing, S.S. and M.B.; visualization, S.S. and M.B.; project administration, S.S. and A.K.A.; All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Deanship of Scientific Research (DSR), King Abdulaziz University, Jeddah, Saudi Arabia, under Grant No. (G-397-130-1439).

Acknowledgments

This project was funded by the Deanship of Scientific Research (DSR), King Abdulaziz University, Jeddah, Saudi Arabia, under Grant No. (G-397-130-1439) The authors, therefore, acknowledge with thanks DSR technical and financial support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Revnic, C.; Ghalambaz, M.; Grosan, T.; Sheremet, M.A.; Pop, I. Impacts of non-uniform border temperature variations on time-dependent nanofluid free convection within a trapezium: Buongiorno’s Nanofluid Model. Energies 2019, 12, 1461. [Google Scholar] [CrossRef] [Green Version]
  2. Cheong, H.T.; Sivasankaran, S.; Bhuvaneswari, M. Natural convection in a wavy porous cavity with sinusoidal heating and internal heat generation. Int. J. Numer. Methods Heat Fluid Flow 2017, 27, 287–309. [Google Scholar] [CrossRef]
  3. Janagi, K.; Sivasankaran, S.; Bhuvaneswari, M.; Eswaramurthi, M. Numerical study on free convection of cold water in a square porous cavity with sinusoidal wall temperature. Int. J. Numer. Methods Heat Fluid Flow 2017, 27, 1000–1014. [Google Scholar] [CrossRef]
  4. Sivasankaran, S.; Cheong, H.T.; Bhuvaneswari, M.; Ganesan, P. Effect of moving wall direction on mixed convection in an inclined lid-driven square cavity with sinusoidal heating. Numer. Heat Transf. A 2016, 69, 630–642. [Google Scholar] [CrossRef]
  5. Eswaramoorthi, S.; Bhuvaneswari, M.; Sivasankaran, S.; Rajan, S.; Alshomrani, A.S. Effect of viscous dissipation and Newtonian heating on convective flow of a second grade fluid over a stretching surface: Analytical and numerical study. Scientia Iranica-Trans. B Mech. Eng. 2019, 26, 1350–1357. [Google Scholar]
  6. Abbas, Z.; Sheikh, M.; Motsa, S.S. Numerical solution of binary chemical reaction on stagnation-point flow of Casson fluid over a stretching/shrinking sheet with thermal radiation. Energy 2016, 95, 12–20. [Google Scholar] [CrossRef]
  7. Rehman, K.U.; Malik, A.A.; Malik, M.Y.; Sandeep, N.; Saba, N.U. Numerical study of double stratification in Casson fluid flow in the presence of mixed convection and chemical reaction. Results Phys. 2017, 7, 2997–3006. [Google Scholar] [CrossRef]
  8. Malik, M.Y.; Khan, M.; Salahuddin, T.; Khan, I. Variable viscosity and MHD flow in Casson fluid with Cattaneo–Christov heat flux model: Using Keller box method. Eng. Sci. Technol. Int. J. 2016, 19, 1985–1992. [Google Scholar] [CrossRef] [Green Version]
  9. Pop, I.; Sheremet, M.A. Free convection in a square cavity filled with a Casson fluid under the effects of thermal radiation and viscous dissipation. Int. J. Numer. Methods Heat Fluid Flow 2017, 27, 2318–2332. [Google Scholar] [CrossRef]
  10. Sheremet, M.A.; Pop, I. Natural convection combined with thermal radiation in a square cavity filled with a viscoelastic fluid. Int. J. Numer. Methods Heat Fluid Flow 2018, 28, 624–640. [Google Scholar] [CrossRef]
  11. Srinivas, S.; Kumar, C.K.; Subramanyam Reddy, A. Pulsating flow of Casson fluid in a porous channel, with thermal radiation, chemical reaction and applied magnetic field. Nonlinear Anal. Model. Control 2018, 23, 213–233. [Google Scholar] [CrossRef]
  12. Niranjan, H.; Sivasankaran, S.; Bhuvaneswari, M. Chemical reaction, soret and dufour effects on MHD mixed convection stagnation point flow with radiation and slip condition. Sci. Iranica-Trans. B Mech. Eng. 2017, 24, 698–706. [Google Scholar] [CrossRef] [Green Version]
  13. Sivasankaran, S.; Niranjan, H.; Bhuvaneswari, M. Chemical reaction, radiation and slip effects on MHD mixed convection stagnation-point flow in a porous medium with convective boundary condition. Int. J. Numer. Methods Heat Fluid Flow 2017, 27, 454–470. [Google Scholar] [CrossRef]
  14. Zheng, S.; Liang, W.; Chu, H.; Zhou, H. Effects of radiation reabsorption of C1-C6 hydrocarbon flames at normal and elevated pressures. Fuel 2020, 266, 117061. [Google Scholar] [CrossRef]
  15. Sun, Y.; Zheng, S.; Jiang, B.; Tang, J.; Liu, F. One-dimensional P1 method for gas radiation heat transfer in spherical geometry. Int. J. Heat Mass Transf. 2019, 145, 118777. [Google Scholar] [CrossRef]
  16. Mansour, M.A.; Abd El-Aziz, M.M.; Mohamed, R.A.; Ahmed, S.E. Numerical simulation of natural convection in wavy porous cavities under the influence of thermal radiation using a thermal non equilibrium model. Transp. Porous Media 2011, 86, 585–600. [Google Scholar] [CrossRef]
  17. Mahapatra, T.R.; Pulal, D.; Mondal, S. Combined effects of thermal radiation and heat generation on natural convection in a square cavity filled with Darcy-Forchheimer porous medium. Int. J. Appl. Math. Comput. 2012, 4, 359–368. [Google Scholar]
  18. Miroshnichenko, I.V.; Sheremet, M.A.; Mohamad, A.A. The influence of surface radiation on the passive cooling of a heat-generating element. Energies 2019, 12, 980. [Google Scholar] [CrossRef] [Green Version]
  19. 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]
  20. Sheremet, M.A.; Grosan, T.; Pop, I. Natural convection and entropy generation in a square cavity with variable temperature side walls filled with a nanofluid: Buongiorno’s Mathematical Model. Entropy 2017, 19, 337. [Google Scholar] [CrossRef]
  21. 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. [Google Scholar] [CrossRef] [Green Version]
  22. Chamkha, A.J.; Selimefendigil, F.; Oztop, H.F. MHD mixed convection and entropy generation in a lid-driven triangular cavity for various electrical conductivity models. Entropy 2018, 20, 903. [Google Scholar] [CrossRef] [Green Version]
  23. Alsabery, A.I.; Ismael, M.A.; Chamkha, A.J.; Hashim, I. Numerical investigation of mixed convection and entropy generation in a wavy-walled cavity filled with nanofluid and involving a rotating cylinder. Entropy 2018, 20, 664. [Google Scholar] [CrossRef] [Green Version]
  24. Oztop, H.F.; A Almeshaal, M.; Kolsi, L.; Rashidi, M.M.; E Ali, M. Natural Convection and irreversibility evaluation in a cubic cavity with partial opening in both top and bottom sides. Entropy 2019, 21, 116. [Google Scholar] [CrossRef] [Green Version]
  25. Khan, N.S.; Shah, Z.; Islam, S.; Khan, I.; Alkanhal, T.A.; Tlili, I. Entropy generation in MHD mixed convection non-Newtonian second-grade nanoliquid thin film flow through a porous medium with chemical reaction and stratification. Entropy 2019, 21, 139. [Google Scholar] [CrossRef] [Green Version]
  26. Ahmed, S.E.; Mansour, M.A.; Mahdy, A.; Mohamed, S.S. Entropy generation due to double diffusive convective flow of Casson fluids over nonlinearity stretching sheets with slip conditions. Eng. Sci. Technol. Int. J. 2017, 20, 1553–1562. [Google Scholar] [CrossRef]
  27. Hajji, F.; Mazgar, A.; Sakly, A.; Nejma, F.B. Entropy generation due to combined natural convection and thermal radiation within a rectangular enclosure. Heat Transf. Eng. 2018, 39, 1698–1714. [Google Scholar] [CrossRef]
  28. Daniel, Y.S.; Aziz, Z.A.; Ismail, Z.; Salah, F. Entropy analysis in electrical magnetohydro-dynamic (MHD) flow of nanofluid with effects of thermal radiation, viscous dissipation, and chemical reaction. Theor. Appl. Mech. Lett. 2017, 7, 235–242. [Google Scholar] [CrossRef]
  29. Kefayati, G.H.R.; Tang, H. Double-diffusive natural convection and entropy generation of Carreau fluid in a heated enclosure with an inner circular cold cylinder (Part I: Entropy generation). Int. J. Heat Mass Transf. 2018, 120, 683–713. [Google Scholar] [CrossRef]
  30. Kefayati, G.H.R. Double-diffusive natural convection and entropy generation of Bingham fluid in an inclined cavity. Int. J. Heat Mass Transf. 2018, 116, 762–812. [Google Scholar] [CrossRef]
  31. Chamkha, A.J.; Rashad, A.M.; Mansour, M.A.; Armaghani, T.; Ghalambaz, M. Effects of heat sink and source and entropy generation on MHD mixed convection of a Cu-water nanofluid in a lid-driven square porous enclosure with partial slip. Phys. Fluids 2017, 29, 052001. [Google Scholar] [CrossRef]
  32. Basak, T.; Kaluri, R.S.; Balakrishnan, A.R. Effects of thermal boundary conditions on entropy generation during natural convection. Numer. Heat Transf. Part A 2011, 59, 372–402. [Google Scholar] [CrossRef]
  33. Ho, C.J.; Chen, M.W.; Li, Z.W. Numerical simulation of natural convection of nanofluid in a square enclosure: Effects due to uncertainties of viscosity and thermal conductivity. Int. J. Heat Mass Transf. 2008, 51, 4506–4516. [Google Scholar] [CrossRef]
  34. Fusegi, T.; Hyun, J.M.; Kuwahara, K.; Farouk, B. A numerical study of three-dimensional natural convection in a differentially heated cubical enclosure. Int. J. Heat Mass Transf. 1991, 34, 1543–1557. [Google Scholar] [CrossRef]
Figure 1. Physical domain.
Figure 1. Physical domain.
Entropy 22 00229 g001
Figure 2. Streamlines for diverse values of radiation and Casson fluid parameters with Gr = 106, Ar = 1.
Figure 2. Streamlines for diverse values of radiation and Casson fluid parameters with Gr = 106, Ar = 1.
Entropy 22 00229 g002
Figure 3. Isotherms for diverse values of radiation and Casson fluid parameters with Gr = 106, Ar = 1.
Figure 3. Isotherms for diverse values of radiation and Casson fluid parameters with Gr = 106, Ar = 1.
Entropy 22 00229 g003
Figure 4. Skin friction vs. radiation (a,b) and Ar (c,d) for different values of β and Rd.
Figure 4. Skin friction vs. radiation (a,b) and Ar (c,d) for different values of β and Rd.
Entropy 22 00229 g004
Figure 5. Kinetic energy vs. radiation (a,b) and Ar (c,d) for different values of β and Rd.
Figure 5. Kinetic energy vs. radiation (a,b) and Ar (c,d) for different values of β and Rd.
Entropy 22 00229 g005aEntropy 22 00229 g005b
Figure 6. Local Nusselt number for various values of Rd with β = 0.01, and 1, Gr = 106, Ar = 1.
Figure 6. Local Nusselt number for various values of Rd with β = 0.01, and 1, Gr = 106, Ar = 1.
Entropy 22 00229 g006
Figure 7. Local Nusselt number for different values of Casson liquid parameter with Rd = 0 and 5 Gr = 106, Ar = 1.
Figure 7. Local Nusselt number for different values of Casson liquid parameter with Rd = 0 and 5 Gr = 106, Ar = 1.
Entropy 22 00229 g007
Figure 8. Averaged Nusselt number vs. radiation (a,b) and Ar (c,d) for various values of β and Rd.
Figure 8. Averaged Nusselt number vs. radiation (a,b) and Ar (c,d) for various values of β and Rd.
Entropy 22 00229 g008aEntropy 22 00229 g008b
Figure 9. Bejan number vs. radiation (a,b) and Ar (c,d) for different values of β and Rd.
Figure 9. Bejan number vs. radiation (a,b) and Ar (c,d) for different values of β and Rd.
Entropy 22 00229 g009
Figure 10. Tcup vs. radiation (a,b) and Ar (c,d) for different values of β and Rd.
Figure 10. Tcup vs. radiation (a,b) and Ar (c,d) for different values of β and Rd.
Entropy 22 00229 g010
Figure 11. TAvg vs radiation (a,b) and Ar (c,d) for different values of β and Rd.
Figure 11. TAvg vs radiation (a,b) and Ar (c,d) for different values of β and Rd.
Entropy 22 00229 g011aEntropy 22 00229 g011b
Figure 12. RMSDTcup vs. radiation (a,b) and Ar (c,d) for different values of β and Rd.
Figure 12. RMSDTcup vs. radiation (a,b) and Ar (c,d) for different values of β and Rd.
Entropy 22 00229 g012
Figure 13. RMSDTavg vs radiation (a,b) and Ar (c,d) for different values of β and Rd.
Figure 13. RMSDTavg vs radiation (a,b) and Ar (c,d) for different values of β and Rd.
Entropy 22 00229 g013
Table 1. Comparison of N u ¯ for a square box with Pr = 0.71, Rd = 0, β = .
Table 1. Comparison of N u ¯ for a square box with Pr = 0.71, Rd = 0, β = .
N u ¯
RaHo et al. [33]Fusegi et al. [34]Present
1031.1181.1061.103
1042.2462.3022.292
1054.5224.6464.628
1068.8259.0128.935

Share and Cite

MDPI and ACS Style

Alzahrani, A.K.; Sivasankaran, S.; Bhuvaneswari, M. Numerical Simulation on Convection and Thermal Radiation of Casson Fluid in an Enclosure with Entropy Generation. Entropy 2020, 22, 229. https://0-doi-org.brum.beds.ac.uk/10.3390/e22020229

AMA Style

Alzahrani AK, Sivasankaran S, Bhuvaneswari M. Numerical Simulation on Convection and Thermal Radiation of Casson Fluid in an Enclosure with Entropy Generation. Entropy. 2020; 22(2):229. https://0-doi-org.brum.beds.ac.uk/10.3390/e22020229

Chicago/Turabian Style

Alzahrani, A. K., S. Sivasankaran, and M. Bhuvaneswari. 2020. "Numerical Simulation on Convection and Thermal Radiation of Casson Fluid in an Enclosure with Entropy Generation" Entropy 22, no. 2: 229. https://0-doi-org.brum.beds.ac.uk/10.3390/e22020229

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