Next Article in Journal
Dynamic Contact Networks in Confined Spaces: Synthesizing Micro-Level Encounter Patterns through Human Mobility Models from Real-World Data
Previous Article in Journal
Singular-Value-Decomposition-Based Matrix Surgery
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Entropy Generation via Non-Similar Numerical Approach for Magnetohydrodynamics Casson Fluid Flow with Joule Heating

1
Mathematics Department, Faculty of Science, Northern Border University, Arar 73213, Saudi Arabia
2
Department of Mechanical Engineering, University of Wah, Wah Cantt 47040, Pakistan
3
Department of Mathematics, University of Wah, Wah Cantt 47040, Pakistan
4
Department of Mathematics, College of Sciences and Arts (Majardah), King Khalid University, Magardah 61937, Saudi Arabia
5
Mathematics Department, College of Humanities and Science in Al Aflaj, Prince Sattam Bin Abdulaziz University, Alkharj 11912, Saudi Arabia
*
Author to whom correspondence should be addressed.
Entropy 2024, 26(8), 702; https://doi.org/10.3390/e26080702
Submission received: 24 June 2024 / Revised: 10 August 2024 / Accepted: 14 August 2024 / Published: 19 August 2024
(This article belongs to the Section Multidisciplinary Applications)

Abstract

:
This analysis emphasizes the significance of radiation and chemical reaction effects on the boundary layer flow (BLF) of Casson liquid over a linearly elongating surface, as well as the properties of momentum, entropy production, species, and thermal dispersion. The mass diffusion coefficient and temperature-dependent models of thermal conductivity and species are used to provide thermal transportation. Nonlinear partial differential equations (NPDEs) that go against the conservation laws of mass, momentum, heat, and species transportation are the form arising problems take on. A set of coupled dimensionless partial differential equations (PDEs) are obtained from a set of convective differential equations by applying the proper non-similar transformations. Local non-similarity approaches provide an analytical approximation of the dimensionless non-similar system up to two degrees of truncations. The built-in Matlab (Version: 7.10.0.499 (R2010a)) solver bvp4c is used to perform numerical simulations of the local non-similar (LNS) truncations.

1. Introduction

One kind of non-Newtonian fluid that falls within the category of shear thinning liquids that induce yield stress is Casson fluid. On the other hand, if shear stress is greater than yield stress, it begins to move. Rouleaux are a type of chainlike structure that can be produced in aqueous base plasma due to the presence of certain materials such as fibrinogen, protein, globulin, and human red blood cells. Several investigations have previously been conducted to identify the Casson fluid across flat or cylindrical stretched sheets in various physical conditions. Jelly, toothpaste, concentrated fruit juices, tomato soup, honey, shampoo, blood, etc., are a few examples of Casson fluids [1]. This model has been investigated and examined by numerous scientists and researchers under various rheological circumstances. Mernone et al. [2], for example, looked into the peristaltic transfer of a Casson fluid. Blood flow is explained by the Casson fluid model, which is quite useful for simulating blood hemodialysis and oxygenators. Raju et al. [3] explained the impact of mass transport and heat in Casson fluid flow on the porous stretched surface. Regarding thermally corrugated porous enclosures with Casson liquid suspension, K. U. Rehman et al. have investigated finite element thermal analysis [4]. Hussain et al. provided an investigation of the Casson rheological fluid containing gold nanoparticles under the effect of magnetic and gravitational forces [5]. Nazir et al. propose a theoretical framework for the mathematical modeling of the Casson rheological fluid [6]. Using the shooting technique, Kempannagari et al. [7] have addressed the MHD stagnation-point flow of micropolar fluid over a curved surface under the effects of radiation, Joule heating, and changing thermal conductivity. Runge–Kutta–Fehlberg solutions have been employed by Gopal et al. [8] to describe the effects of MHD and Joule heating on the Casson fluid flow across the stretching sheet. Muthtamilselvan et al. [9] discussed the stagnation-point flow of the Williamson nanofluid containing gyrotactic microorganisms. The results revealed that the density of microorganisms reduces with the stretching velocity.
Radiation effects have significant applications in physics and engineering processes. It is important for processes involving hot materials or liquids and space technology, where the effects of flow temperature are studied for various flows. In the polymer processing business, thermal radiation effects may be crucial in managing heat transmission. Heat transmission via radiation is a phenomenon known as thermal radiation. Its influence on MHD [10,11,12] flows has significance in engineering and industry, involving high thermal effects in the production of paper plates, petroleum pumps, electronic chips, and metallic item cooling. Kumar et al. [13] utilize the nonlinear thermal radiation influence to study MHD shear-thickening fluid dynamics via stretchy geometry.
A point of stagnation is reached when velocity equals zero. It depicts the flow of fluid in a solid surface’s direct zone. As the fluid gets closer to the solid surface, it splits into two streams. Everywhere on the solid surface, the fluid is stationary because of the no-slip condition. N. Khan and Y.M. Chu [14] addressed modeling and dual solutions for magnetized mixed convective stagnation-point flow of a higher convected Maxwell fluid model with second-order velocity slip. The publications (see Refs. [15,16,17,18,19]) contain a thorough description of the stagnation-point flow across either a stretching or contracting sheet. Mahapatra and Gupta [20] investigated the MHD stagnation-point flow past a stretching sheet.
In the field of engineering, minimizing entropy generation is an emerging thermodynamic technique that is crucial to the design of contemporary heat management systems. Bejan [21] is credited with introducing the idea of entropy creation. He came to the conclusion that the process of minimizing entropy formation improves the system’s thermodynamic efficiency and increases the efficiency of a thermal system. Analyzing entropy in Casson fluid flow is essential for understanding various industrial and biological processes involving non-Newtonian fluids. The Casson fluid model is often applied to fluids with yield stress, such as blood flow in arteries [22], as well as in polymer processing and lubrication systems [23]. Conducting entropy analysis in these contexts is crucial for optimizing lubricant performance and extending the operational life of machinery. Additionally, Casson fluids can be utilized in cooling systems for electronic devices, where entropy generation analysis plays a key role in reducing energy loss and enhancing heat removal efficiency. This, in turn, is vital for improving the performance and durability of electronic components. The physical parameters such as Prandtl number and Schmidt number that significantly influence the analysis of entropy generation fall within the appropriate ranges, as documented in [24,25].
A heightened concentration boundary layer was observed for big activation energies, while the contrary was observed for increasing reaction rates, according to the results of a numerical analysis of the mathematical system using the shooting technique. Additionally, entropy generation for Williamson fluid was higher than for Newtonian fluid. Heat transmission of a mixed convection ferrofluid down a vertical duct with variable wall temperatures has been examined by Gul et al. [26]. They discovered that the primary factors influencing the various velocity and temperature outcomes are viscosity and thermal conductivity. The unstable electrically conducting Casson fluid flows across an oscillating vertical porous plate with constant wall temperature have been studied by Khalid et al. [27]. Gul et al. [28] studied the mixed convection time-dependent MHD flow of an incompressible nanofluid (four distinct forms) under three different flow conditions inside a channel filled with saturated porous medium with non-uniform wall temperatures. Entropy analysis plays a key role in enhancing the efficiency of various real-world devices and systems. For instance, in artificial heart valves, stents, and blood pumps, entropy analysis aids in optimizing design to minimize energy losses, improve blood circulation, and reduce the risk of clot formation. In lubrication systems, such as bearings, gears, and hydraulics, it can lead to the development of better lubricant formulations and more efficient system designs, thereby reducing frictional losses and prolonging the life of mechanical components. Additionally, optimizing the flow of Casson fluids in cooling systems, including heat sinks, liquid cooling systems, and thermal management units, through entropy analysis can improve heat dissipation, reduce energy consumption, and enhance the performance and longevity of electronic components.
The entropy number alone cannot adequately distinguish between different irreversibility mechanisms. To address this limitation, the Bejan number ( B e ) is introduced. The Bejan number is calculated by dividing the heat transfer irreversibility by the total entropy production. Its value ranges between 0 and 1: when B e equals 1, heat transfer irreversibility is dominant, whereas when B e equals 0, the irreversibility due to viscous dissipation prevails.

2. Problem Formulations

As seen in Figure 1, consider a flat vertical plate that has been adjusted parallel to a stationary fluid. The “ x ” and “ y ” coordinates represent the streamwise and crosswise directions, respectively. The surface expands at a linear rate, which causes the magnetized vertical flow to propagate. The “ x ”-momentum equation incorporates the linear buoyancy with respect to concentration and temperature variations. In addition, the energy equation takes thermal radiation and chemical reactions into account. The concentration equation makes the assumption that the chemical reaction is linear.
u x + v y = 0 ,
u u x + v u y = U e d U e d x + ν 1 + 1 β u y y + g β T T T + g β C C C σ B 0 2 ρ u ,
u T x + v T y = α T y y 1 ρ C p q r y + σ B 0 2 ρ C p u 2
u C x + v C y = D C y y k 1 C C .
The associated conditions are as follows:
u = U w x ,         v = 0 ,         T y = h f k f ( T w T ) ,   C = C w   a t   y = 0 , T T ,             u U e ,         C C   a s     y .
With the help of Rosseland estimation, the expression for q r is given by
q r = 4 σ 3 k T y 4 .
Ignoring T 4 and all higher terms, the expression may be expressed as
T 4 4 T 3 T 3 T 4 .
By using (6) and (7) in (4), we obtain
u T x + v T y = α + 16 σ T 3 3 k ρ C p T y y
u and v are the component of velocity in the x and y directions. U w x = u w e x l is the velocity at the wall, U e = u e e x l is the velocity of the fluid away from the plate at free stream, β is the Casson parameter, B 0 2 is the magnetic field parameter, C p is the specific heat capacity, and q r is the radiative heat flux.

2.1. Non-Similar Transformation

The transformation of local non-similarity was first presented in 1970 by Sparrow et al. [29]. Two independent variables are involved in this transformation: one is non-similar and the other is pseudo-similar. By using this method, problems become dimensionless and partial differential equations remain partial differential equations. The local non-similarity method yields highly accurate results at all streamwise locations, except near the point of separation.
ξ = e x l , η = y U e 2 l υ , ψ x , y = 2 υ l U e f ( ξ , η ) , u = U e f η ,   v = υ U e 2 l f 2 υ U e l ξ f ξ y U e 2 l f η f , U e = u e e x l   and   θ ξ , η = T T T w T , φ ξ , η = C C C w C .
After applying non-similarity transformations on Equations (2), (4), and (8), they are reduced to (10)–(12):
f η 2 + ξ 2 f ξ η f η f 2 2 f η 2 ξ f ξ 2 f η 2 = 1 + 1 β 3 f η 3 + G r l R e 2 θ + G r c R e 2 φ M f η ,
P r f θ η 2 ξ θ ξ f η + 2 ξ f ξ θ η + J f η 2 = 1 + 4 R d 3 2 θ η 2 ,
f φ η 2 ξ φ ξ f η + 2 ξ f ξ φ η = 1 2 S c 2 φ η 2 + γ φ .
Associated conditions are
f ξ , η = ϵ ,         θ ξ , η = 2 R e B i 1 θ ,         f ξ , η = ξ f ξ , φ ξ , η = 1   a t   η = 0 , θ ξ , η 0 ,   f ξ , η 1 ,   φ ξ , η 0   a t   η .
G r l = g β T T f T l 3 v 2 , G r C = g β C C w C l 3 v 2 is the Grashof number for heat and mass transfer, R e = U e l υ is the Reynolds number, B i = h f k f l is the Biot number, and ϵ = U w U e is the ratio of stretching to free stream velocities.

2.2. Non-Similarity Technique

For the first and second level of truncation, the following equations are used.

2.2.1. Truncation Level 1

We consider ξ ( . ) ξ = 0 at this level to obtain the required ordinary differential equations:
f = 1 1 + 1 β f 2 f 2 f G r l R e 2 θ G r c R e 2 φ + M f ,
θ = P r 1 + 4 R d 3 f θ + J ( f ) 2
φ = 2 S c γ φ f φ
Associated conditions are
f ξ , η = ϵ ,       θ ξ , η = 2 R e B i 1 θ ,                 f ξ , η = 0 ,   φ ξ , η = 1   a t   η = 0 θ ξ , η 0 ,   f ξ , η 1 ,   φ ξ , η 0   a t   η .  

2.2.2. Truncation Level 2

By introducing the new variables, the derivative of the following functions will be written as
Let   f ξ = p ,   2 f ξ η = f ξ = p , θ ξ = q , 2 θ ξ η = θ ξ = q , φ ξ = g , 2 φ ξ η = φ ξ = g f = 1 1 + 1 β f 2 + ξ p f f 2 f ξ p f G r l R e 2 θ G r c R e 2 φ + M f ,
θ = P r 1 + 4 R d 3 2 ξ q f f θ 2 ξ p θ f θ J ( f ) 2 ,
φ = 2 S c 2 ξ   g f f φ 2 ξ p φ + γ φ .
With conditions,
f ξ , η = ϵ ,             θ ξ , η = 2 R e B i 1 θ ,     f ξ , η = 0 ,           φ ξ , η = 1   a t   η = 0 ,           θ ξ , η 0 ,           f ξ , η 1 ,           φ ξ , η 0   a t   η .
Taking the derivative with respect to ξ for (18)–(20), we obtain
p = 1 1 + 1 β 3 p f + ξ p 2 f 3 2 p f f 2 p ξ p p G r l R e 2 q G r c R e 2 g + M p ,
q = P r 1 + 4 R d 3 2 ξ q p 3 p θ f q + 2 q f 2 ξ p q 2 J f p ,
g = 2 S c γ g 3 p φ f g + 2 g f + 2 ξ g p 2 ξ p g ξ γ φ .
With conditions,
p ξ , η = 0 ,         q ξ , η = 2 R e B i . q ,         p ξ , η = 0 ,         g ξ , η = 0   a t   η = 0 , q ξ , η 0 ,         p ξ , η 0 ,         g ξ , η 0   a t   η .
The mathematical expression of physical interest is specified by
C f = τ w ρ U e 2 , N u x = x q w k T w T , S h x = x q m D B C w C
τ w = μ u y y = 0 ,   q w = k + 16 σ * T 3 3 k * T y y = 0 , and q m = D C y y = 0 are the shear stress, the wall heat flux, and wall mass flux, respectively.
The significant transport quantities [30] are given as
2 1 2 R e x 1 2 C f x = 1 + 1 β 1 ϵ 2 2 f η 2 η = 0
2 R e x 1 2 N u x = 1 + 4 3 R d θ η η = 0
2 R e x 1 2 S h x = φ η η = 0

2.3. Entropy Analysis

The loss of productive work in an engineering system is directly correlated with the generation of entropy. According to But et al. [31], the volumetric entropy expression for a viscous incompressible Casson fluid in the Cartesian coordinate system is as follows:
S G = κ T 2   1 + 16 σ 1 T 3 3 k k T y 2 + μ T     1 + 1 β u y 2 + D C C y 2 + D T T y C y + σ B 0 2 T u
S G = κ T 2   1 + 16 σ 1 T 3 3 k k T w T 2 U e 2 l υ θ 2 + μ T     1 + 1 β U e 3 2 l v f 2 + D B C C w C 2 U e 2 l υ φ 2 + D B T T w T C w C U e 2 l υ θ φ + σ T B 0 2 U 0 2 f 2
The irreversibility of heat transfer causes the first part of the statement above to arise; viscous dissipation causes the second part to appear; and the concentration of nanoparticles causes the third component to generate local entropy. Using Equation (9) in the above equations yields the entropy number expression in a dimensionless form:
N s = S G S G 0 = R e l 1 2 1 + 4 3 R d θ 2 + B r   Ω T R e l 2 1 + 1 β f 2 + ϵ Ω C   Ω T θ φ + ϵ R e l 2 Ω C   Ω T 2 φ 2 + M B r   Ω T R e l f 2
where R d = 4 σ 1 T 3 k * k ,   Ω T = T T ,   Ω C = C C ,   B r = μ U e 2 / κ T ,   M = σ B 0 2 l ρ U e and ϵ = R D B C / κ are the characteristic entropy, the dimensionless temperature, concentration differences, the Brinkman number, and the nanoparticle mass transfer parameter, respectively. Classifying the dominance of the irreversibility mechanism is not possible with the entropy number. The Bejan number is created to examine the irreversibility method in order to solve this problem. The formation of the Bejan number B e involves dividing the irreversibility of heat transfer by the overall creation of entropy. The Bejan number ( B e ) for the current analysis has the following form (see, for instance, [32]):
B e = R e l 1 2 1 + 4 3 R d θ 2 + ϵ Ω C   Ω T θ φ + ϵ R e l 2 Ω C   Ω T 2 φ 2 R e l 1 2 1 + 4 3 R d θ 2 + B r   Ω T R e l 2 1 + 1 β f 2 + ϵ Ω C   Ω T θ φ + ϵ R e l 2 Ω C   Ω T 2 φ 2 + M B r   Ω T R e l f 2 .
The aforementioned equation unequivocally demonstrates that the range of B e is between 0 and 1 : at B e = 1 , the irreversibility mechanism resulting from viscous dissipation and heat mass diffusion dominates the entropy generation; at B e = 1 / 2 , the involvement in entropy production due to heat transfer and the sum of viscous dissipation and heat and mass diffusion are equal; and at B e = 0 , the entropy production is dominated by heat transfer.

2.4. Solution Procedure

A numerical scheme called bvp4c is taken into consideration for solving ordinary differential equations that are derived by applying transformations on the governing equations of the specified flow issue. The formula known as three-stage Lobatto III, which yields fourth-order accuracy when solving a system of ordinary differential equations (ODEs), serves as the foundation for the solver. First-order ordinary differential equations (ODEs), boundary data, and an initial guess are required by the solver. We receive a system of initial value issues for the solver to solve if the differential equations have more than one order.
f = f 2 ,   f = f 3
f 3 = f = 1 / 1 + 1 / β f 2 2 + ξ f 9 f 2 f 1 2 f 3 ξ f 8 f 3 G r l R e 2 f 4 G r c R e 2 f 6 + M f 2 ,
θ = f 5 ,
                                                                                                f 5 = θ = P r / ( 1 + 4 R d / 3 ) 2 ξ   f 11 f 2 f 1 f 5 2 ξ f 8 f 7 J f 2 2 ,
φ = f 7 ,
                                                                                                                  f 7 = φ = 2 S c 2 ξ   f 13 f 2 f 1 f 7 2 ξ f 7 f 8 + γ f 6 ,
p = f 9 ,   p = f 10 ,
f 11 = p = 1 / 1 + 1 / β 3 f 9 f 2 + ξ f 9 2 f 2 3 2 f 8 f 3 f 1 2 f 10 ξ f 8 f 10 G r l R e 2 f 11 G r c R e 2 f 13 + M f 9 ,
q = f 12 ,
f 13 = q = P r / ( 1 + 4 R d / 3 ) 2 ξ f 11 f 9 3 f 8 f 5 f 1 f 12 + 2 f 11 f 2 2 ξ f 8 f 12 2 J f 2 f 8 ,
                                                                                                                                                                      g = f 14 ,
f 15 = g = 2 S c γ f 13 3 f 8 f 7 f 1 f 14 + 2 f 13 f 2 + 2 ξ f 13 f 9 2 ξ f 8 f 14 ξ γ f 6 .
where f = f 1 , θ = f 4 , φ = f 6 , p = f 8 , q = f 11 and g = f 13 and the related boundary data read as
f 1 0 + 2 ξ f 8 , f 2 0 ϵ , f 2 i n f 1 , f 5 0 + 2 R e B i 1 f 4 , f 4 i n f ,   f 6 0 1 ,   f 6 i n f ,   f 8 0 + 2 f 8 ,   f 9 0 ,   f 9 i n f ,   f 12 0 2 R e B i f 11 ,   f 11 i n f ,   f 13 0 ,   f 13 i n f t
The use of a non-similar numerical approach allows for the analysis of complex flow phenomena that do not conform to traditional similarity solutions. This method provides more flexibility and accuracy in solving partial differential equations governing the fluid flow. The non-similar numerical approach fills the gap left by traditional methods in terms of providing sufficient predictive capabilities for real-world problems that are complex. It makes it possible to accurately simulate intricate fluid behaviors in a variety of situations.

3. Results and Discussions

Equations (18)–(20) and (22)–(24) with boundary conditions (21) and (25) were solved numerically using the built-in MATLAB package, bvp4c, to obtain a numerical solution. Table 1 presents the numerical values for the skin friction coefficient, local Nusselt number, and local Sherwood number. The data indicate that an increase in the Prandtl number and radiation parameter leads to an enhancement in the local Nusselt number. This is because the Prandtl number is inversely related to thermal diffusivity, which results in a decrease in thermal diffusivity. As a result, the temperature slows down, causing the thermal boundary layer to contract, and thereby increasing the local Nusselt number. The skin friction coefficient rises with an increase in the parameter β , which is due to the more rapid dispersion of momentum in non-Newtonian fluids compared to Newtonian fluids. This increased dispersion occurs because non-Newtonian fluids sustain a greater applied force than Newtonian fluids. Additionally, the local Sherwood number decreases with an increase in the Schmidt number, and the skin friction coefficient decreases as the Grashof number in heat transfer declines. Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11 explore the effects of the parameters involved on the entropy generation rate and Bejan number. The graphical results indicate that the stretching plate surface is a significant contributor towards the production of entropy. Figure 2 examines the effects of the Brinkman number ( B r ) on local entropy generation. The figure shows that entropy generation rises as B r increases. This increase is clearly attributed to fluid friction, leading to an elevation in entropy production. In physical terms, a higher Brinkman number indicates that a greater portion of the mechanical energy is being converted into heat due to internal friction within the fluid and consequently decreases the performance of the thermal system. Figure 3 illustrates the impact of the Brinkman number on the local Bejan number B e . As the Brinkman number increases, the local Bejan number decreases. Elevated B r values enhance the viscous effects, ultimately resulting in a dominance of fluid friction irreversibility. However, for a fixed value of B r , this figure illustrates that the major contributors to entropy production are diffusivity and heat transfer at the plate, and in the distant region (far from the boundary), fluid friction and magnetic field entropy become stronger. Table 2 presents the quantitative results for the average entropy generation number ( N s ¯ ) and average Bejan number B e ¯ as the Brinkman number ( B r ) varies. The average entropy generation demonstrates an upward trend with increasing B r , while the Bejan number shows the opposite trend. Increasing B r from 0.1 to 1.1, it is observed that N s ¯ is increased up to 10.43 times and B e ¯ decreases up to 3.60 times. Figure 4 explores the impact of the Casson parameter β on the entropy profile. The decreasing trend has been observed in the entropy profile due to the increasing value of β , but this depreciation is not so prominent, which can also be observed from Table 3 in the case of average entropy ( N s ¯ ) . A similar trend can be seen in Figure 5 for the case of Bejan numbers for different β values. Increasing the value of the mass transfer parameter ε results in decreasing behavior in the entropy generation profile, as shown in Figure 6. On the other hand, in Figure 7, increasing the value of ε also increases the Bejan number, which indicates that thermal irreversibility and diffusion entropy profile are more dominant. The impact of the magnetic parameter M on the entropy profile is presented in Figure 8, which shows an increment in N s . When M changes from 0.1 to 1.1, N s ¯ augments 4.81 times (see in Table 4). The Lorentz force acts as a damping force on fluid motion, converting kinetic energy into thermal energy. This conversion increases the process’s irreversibility. Due to the strengthened magnetic field parameter M , the Bejan number profile decreases, which ultimately depicts the dominance of fluid friction and magnetic field entropy. Figure 10 reveals the influence of the radiation parameter R d on the entropy profile. It is seen that entropy increases as R d increases. Enhanced radiative heat transfer results in increased energy dissipation rates. Because radiation tends to be directional and can cause uneven temperature distributions, it amplifies the system irreversibly, thus boosting entropy generation. Table 5 also shows the same impact upon average entropy. The local Bejan number is not influenced by the radiation parameter R d in the vicinity of the plate, but in the main flow region, the Bejan profile is enhanced, as one can see in Figure 11. Figure 12, Figure 13 and Figure 14 were drawn using software that uses the finite element method. The two-dimensional region is used as the geometry of the problem. The left and right walls are considered as outlets, the lower wall is considered to be a moving wall, part of the top wall is considered with no-slip conditions, and the rest of the wall is considered to be an inlet with velocity along the y-axis. The mixed convection and MHD effects are incorporated with different velocities of the plate. The solution is displayed in different graphs shown in Figure 12, Figure 13 and Figure 14. The impact of the physical parameters Brinkman number, Casson parameter, and magnetic parameter on entropy generation and Bejan number follows a trend similar to that reported by Butt et al. [33].

4. Conclusions

This paper mathematically models the heat and mass transfer of a thermally radiative, chemically reactive mixed convection flow of Casson fluid over a stretching sheet. By applying the local non-similarity method, the governing equations of the flow phenomena are transformed into their dimensionless form. The resulting partial differential equations are transformed to an ordinary differential equation, and then solved numerically using the bvp4c built-in solver in MATLAB. The results obtained with different physical parameters are plotted and discussed. The outcomes of this study are as follows:
  • The local entropy generation number N s is enhanced by increases in the Brinkman number, magnetic field, and radiation parameter.
  • An increase in the value of the Casson parameter and mass transfer parameter causes the local entropy generation number to decrease in the vicinity of the convective stretching surface.
  • The average entropy generation number ( N s ¯ ) increases up to 10.43 times as B r changes from 0.1 to 1.1, while the average Bejan number B e ¯ decreases up to 3.60 times.
  • The Bejan number ( B e ) graphs indicate that near the stretching surface, the effects of heat transfer and diffusion entropy are predominant, while in regions farther away, the entropy effects due to fluid friction and the magnetic field become more significant.
  • At the surface of the plate, the average Bejan number decreased up to 3.21 times as the magnetic parameter varied from 0.1 to 1.1.
The current approach towards the entropy analysis can further be extended to the problem related to Casson nanofluids in the presence of a magnetic field.

Author Contributions

Methodology, H.L.; Software, S.K.; Validation, M.M.; Formal analysis, S.O.H.; Investigation, A.G. All authors have read and agreed to the published version of the manuscript.

Funding

The authors extend their appreciation to the Deanship of Research and Graduate Studies at King Khalid University for funding this work through Large Research Project under grant number RGP2/115/45. The authors extend their appreciation to the Deanship of Scientific Research at Northern Border University, Arar, KSA for funding this research work through the project number NBU-FFR-2024-2920-03 and this study is supported via funding from Prince Sattam bin Abdulaziz University project number (PSAU/2024/R/1445).

Institutional Review Board Statement

Not applicable.

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Vijayaragavan, R. Magnetohydrodynamic radiative casson fluid flow over a stretching sheet with heat source/sink. Adv. Phys. Theo. Appl. 2016, 55, 13–23. [Google Scholar]
  2. Mernone, A.V.; Mazumdar, J.N. A mathematical study of peristaltic transport of a casson fluid. Math. Comput. Model. 2002, 35, 895–912. [Google Scholar] [CrossRef]
  3. Raju, C.; Sandeep, N.; Sugunamma, V.; Babu, M.J.; Reddy, J.R. Heat and mass transfer in magnetohydrodynamic Casson fluid over an exponentially permeable stretching surface. Eng. Sci. Technol. Int. J. 2016, 19, 45–52. [Google Scholar] [CrossRef]
  4. Rehman, K.U.; Algehyne, E.A.; Shahzad, F.; Sherif, E.-S.M.; Chu, Y.-M. On thermally corrugated porous enclosure (TCPE) equipped with casson liquid suspension: Finite element thermal analysis. Case Stud. Therm. Eng. 2021, 25, 100873. [Google Scholar] [CrossRef]
  5. Hussain, F.; Nazeer, M.; Altanji, M.; Saleem, A.; Ghafar, M. Thermal analysis of Casson rheological fluid with gold nanoparticles under the impact of gravitational and magnetic forces. Case Stud. Therm. Eng. 2021, 28, 101433. [Google Scholar] [CrossRef]
  6. Nazeer, M.; Hussain, F.; Hameed, M.K.; Khan, M.I.; Ahmad, F.; Malik, M.Y.; Shi, Q.H. Development of mathematical modeling of multi-phase flow of Casson rheological fluid: Theoretical approach. Chaos Solitons Fractals 2021, 150, 111198. [Google Scholar] [CrossRef]
  7. Kempannagari, A.K.; Buruju, R.R.; Naramgari, S.; Vangala, S. Effect of Joule heating on MHD non-Newtonian fluid flow past an exponentially stretching curved surface. Heat Transf. 2020, 49, 3575–3592. [Google Scholar] [CrossRef]
  8. Gopal, D.; Kishan, N.; Raju, C.S.K. Viscous and Joule’s dissipation on Casson fluid over a chemically reacting stretching sheet with inclined magnetic field and multiple slips. Inform. Med. Unlocked 2017, 9, 154–160. [Google Scholar] [CrossRef]
  9. Muthtamilselvan, M.; Suganya, S.; Al-Mdallal, Q.M. Stagnation-Point Flow of the Williamson Nanofluid Containing Gyrotactic Micro-organisms. Proc. Natl. Acad. Sci. India Sect. A Phys. Sci. 2021, 91, 633–648. [Google Scholar] [CrossRef]
  10. Al-Zubaidi, A.; Nazeer, M.; Subia, G.S.; Hussain, F.; Saleem, S.; Ghafar, M. Mathematical modeling and simulation of MHD electro-osmotic flow of Jeffrey fluid in convergent geometry. Waves Random Complex Media 2021, 1–17. [Google Scholar] [CrossRef]
  11. Nazeer, M.; Hussain, F.; Shabbir, L.; Saleem, A.; Khan, M.I.; Malik, M.Y.; Sun, T.-C.; Hussain, A. A comparative study of MHD fluid-particle suspension induced by metachronal wave under the effects of lubricated walls. Int. J. Mod. Phys. B 2021, 35, 2150204. [Google Scholar] [CrossRef]
  12. Nazeer, M.; Saleem, S.; Hussain, F.; Zia, Z.; Khalid, K.; Feroz, N. Heat transmission in a magnetohydrodynamic multiphase flow induced by metachronal propulsion through porous media with thermal radiation. Proc. Inst. Mech. Eng. Part E J. Process Mech. Eng. 2022, 09544089221075299. [Google Scholar] [CrossRef]
  13. Kumar, A.; Sugunamma, V.; Sandeep, N. Impact of non-linear radiation on MHD non-aligned stagnation point flow of micropolar fluid over a convective surface. J. Non-Equilib. Thermodyn. 2018, 43, 327–345. [Google Scholar]
  14. Khan, N.M.; Chu, Y.M.; Khan, M.I.; Kadry, S.; Qayyum, S. Modeling and dual solutions for magnetized mixed convective stagnation point flow of upper convected Maxwell fluid model with second-order velocity slip. Math. Methods Appl. Sci. 2020. [Google Scholar] [CrossRef]
  15. Bhattacharyya, K. Dual solutions in boundary layer stagnation-point flow and mass transfer with chemical reaction past a stretching/shrinking sheet. Int. Commun. Heat Mass Transf. 2011, 38, 917–922. [Google Scholar] [CrossRef]
  16. Bhattacharyya, K. Dual solutions in unsteady stagnation-point flow over a shrinking sheet. Chin. Phys. Lett. 2011, 28, 084702. [Google Scholar] [CrossRef]
  17. Bhattacharyya, K. MHD stagnation-point flow of Casson fluid and heat transfer over a stretching sheet with thermal radiation. J. Thermodyn. 2013, 2013, 169674. [Google Scholar] [CrossRef]
  18. Bhattacharyya, K. Heat transfer analysis in unsteady boundary layer stagnation-point flow towards a shrinking/stretching sheet. Ain Shams Eng. J. 2013, 4, 259–264. [Google Scholar] [CrossRef]
  19. Bhattacharyya, K.; Mukhopadhyay, S.; Layek, G. Slip effects on an unsteady boundary layer stagnation-point flow and heat transfer towards a stretching sheet. Chin. Phys. Lett. 2011, 28, 094702. [Google Scholar] [CrossRef]
  20. Mahapatra, T.R.; Gupta, A. Magnetohydrodynamic stagnation-point flow towards a stretching sheet. Acta Mech. 2001, 152, 191–196. [Google Scholar] [CrossRef]
  21. Bejan, A. Entropy Generation Through Heat and Fluid Flow; Wiley and Sons: Hoboken, NJ, USA, 1994. [Google Scholar]
  22. Usha, R.; Uma, R. Entropy generation in a Casson fluid flow in a stenosed artery. J. Mech. Med. Biol. 2013, 13, 1350045. [Google Scholar]
  23. Rizvi, S.A.A.; Ansari, A.K. Heat transfer and entropy generation in Casson fluid flow over a stretching sheet with thermal radiation and viscous dissipation. J. Mol. Liq. 2016, 223, 1202–1210. [Google Scholar]
  24. Shevchuk, I.V. Convective Heat and Mass Transfer in Rotating Disk Systems; Springer: Berlin/Heidelberg, Germany, 2009. [Google Scholar]
  25. Shevchuk, I.V. A Self-Similar Solution of Navier–Stokes and Energy Equations for Rotating Flows Between a Cone and a Disk. High Temp. 2004, 42, 104–110. [Google Scholar] [CrossRef]
  26. Gul, A.; Khan, I.; Shafie, S.; Khalid, A.; Khan, A. Heat transfer in MHD mixed convection flow of a ferrofluid along a vertical channel. PLoS ONE 2015, 10, e0141213. [Google Scholar] [CrossRef]
  27. Khalid, A.; Khan, I.; Khan, A.; Shafie, S. Unsteady MHD free convection flow of Casson fluid past over an oscillating vertical plate embedded in a porous medium. Eng. Sci. Technol. Int. J. 2015, 18, 309–317. [Google Scholar] [CrossRef]
  28. Gul, A.; Khan, I.; Shafie, S. Energy transfer in mixed convection MHD flow of nanofluid containing different shapes of nanoparticles in a channel filled with saturated porous medium. Nanoscale Res. Lett. 2015, 10, 490. [Google Scholar]
  29. Sparrow, E.; Quack, H.; Boerner, C. Local nonsimilarity boundary-layer solutions. AIAA J. 1970, 8, 1936–1942. [Google Scholar] [CrossRef]
  30. Shevchuk, I.V.; Saniei, N.; Yan, X.T. Impingement Heat Transfer over a Rotating Disk: Integral Method. J. Thermophys. Heat Transf. 2003, 17, 291–293. [Google Scholar] [CrossRef]
  31. Butt, A.S.; Ali, A.; Mehmood, A. Entropy analysis in MHD nanofluid flow near a convectively heated stretching surface. Int. J. Exergy 2016, 20, 318–342. [Google Scholar] [CrossRef]
  32. Kefayati, G.H.R. Simulation of natural convection and entropy generation of non-Newtonian nanofluid in a porous cavity using Buongiorno’s mathematical model. Int. J. Heat Mass Transf. 2017, 112, 709–744. [Google Scholar] [CrossRef]
  33. Butt, A.S.; Maqbool, K.; Imran, S.M.; Ahmad, B. Entropy generation effects in MHD Casson nanofluid past a permeable stretching surface. Int. J. Exergy 2020, 31, 150–171. [Google Scholar] [CrossRef]
Figure 1. Geometry of the problem.
Figure 1. Geometry of the problem.
Entropy 26 00702 g001
Figure 2. Local entropy generation profile for different Br.
Figure 2. Local entropy generation profile for different Br.
Entropy 26 00702 g002
Figure 3. B e plotted against η for different Brinkman numbers.
Figure 3. B e plotted against η for different Brinkman numbers.
Entropy 26 00702 g003
Figure 4. Dependence of local entropy generation number on β .
Figure 4. Dependence of local entropy generation number on β .
Entropy 26 00702 g004
Figure 5. Impact of β on local Bejan number.
Figure 5. Impact of β on local Bejan number.
Entropy 26 00702 g005
Figure 6. Variations in local entropy generation due to ε .
Figure 6. Variations in local entropy generation due to ε .
Entropy 26 00702 g006
Figure 7. Local Bejan number plotted for different values of ε .
Figure 7. Local Bejan number plotted for different values of ε .
Entropy 26 00702 g007
Figure 8. Effects of M on local entropy generation.
Figure 8. Effects of M on local entropy generation.
Entropy 26 00702 g008
Figure 9. Local Bejan number profile for various values of M.
Figure 9. Local Bejan number profile for various values of M.
Entropy 26 00702 g009
Figure 10. Local entropy generation profile for various values of R d .
Figure 10. Local entropy generation profile for various values of R d .
Entropy 26 00702 g010
Figure 11. Influence of R d on local Bejan number profile.
Figure 11. Influence of R d on local Bejan number profile.
Entropy 26 00702 g011
Figure 12. Surface and streamlines for U e = 0.1 . (a) Surface plot of the velocity profile, (b) streamlines, (c) surface plot for the temperature profile (d) Isothermal contours.
Figure 12. Surface and streamlines for U e = 0.1 . (a) Surface plot of the velocity profile, (b) streamlines, (c) surface plot for the temperature profile (d) Isothermal contours.
Entropy 26 00702 g012
Figure 13. Surface and streamlines for U e = 0.3 . (a) Surface plot of the velocity profile, (b) streamlines, (c) surface plot for the temperature profile (d) Isothermal contours.
Figure 13. Surface and streamlines for U e = 0.3 . (a) Surface plot of the velocity profile, (b) streamlines, (c) surface plot for the temperature profile (d) Isothermal contours.
Entropy 26 00702 g013
Figure 14. Surface and streamlines for U e = 0.5 . (a) Surface plot of the velocity profile, (b) streamlines, (c) surface plot for the temperature profile (d) Isothermal contours.
Figure 14. Surface and streamlines for U e = 0.5 . (a) Surface plot of the velocity profile, (b) streamlines, (c) surface plot for the temperature profile (d) Isothermal contours.
Entropy 26 00702 g014
Table 1. List of numerical values for skin friction coefficient and local Nusselt and Sherwood numbers using γ = 0.9 , B i = 0.9 and R e = 0.7 .
Table 1. List of numerical values for skin friction coefficient and local Nusselt and Sherwood numbers using γ = 0.9 , B i = 0.9 and R e = 0.7 .
P r G r L G r C β S c R d ϵ 2 1 2 R e x 1 2 C f x 2 R e x 1 2 N u x 2 R e x 1 2 S h x
0.7 1.5126 0.5854 1.2851
1.0 1.4793 0.6245 1.2841
0.4 1.3716 0.6225 1.2819
0.6 1.5862 0.6264 1.2863
0.5 1.3819 0.6227 1.2824
0.6 1.4793 0.6245 1.2841
0.6 1.5691 0.6188 1.2798
0.7 1.5307 0.6207 1.2814
0.4 1.5759 0.6222 1.1487
0.6 1.4947 0.6196 1.4012
0.6 1.5101 0.6509 1.4016
0.7 1.5245 0.6814 1.4020
0.9 1.4947 0.6196 1.4012
1.1 0.5685 0.6388 1.4402
Table 2. Average entropy number ( N s ¯ ) and average Bejan number ( B e ¯ ) with variation of Br at fixed γ = 0.3 , P r = 0.7 , S c = 0.9 , R d = 0.3 , G r c = 0.8 , G r l = 0.5 , β = 1.0 , B i = 0.1 , ε = 0.1 , R e = 0.7 , Ω c = Ω t = 0.1 , M = 0.1 .
Table 2. Average entropy number ( N s ¯ ) and average Bejan number ( B e ¯ ) with variation of Br at fixed γ = 0.3 , P r = 0.7 , S c = 0.9 , R d = 0.3 , G r c = 0.8 , G r l = 0.5 , β = 1.0 , B i = 0.1 , ε = 0.1 , R e = 0.7 , Ω c = Ω t = 0.1 , M = 0.1 .
B r 0.10.30.50.70.91.1
N s ¯ 0.07530.21740.35940.50150.64350.7856
B e ¯ 0.08910.05630.04230.03410.02870.0247
Table 3. Average entropy number ( N s ¯ ) and average Bejan number ( B e ¯ ) with variation of β at fixed γ = 0.3 , P r = 0.7 , S c = 0.9 , R d = 0.3 , G r c = 0.8 , G r l = 0.5 , B r = 0.1 , B i = 0.1 , ε = 0.1 , R e = 0.7 , Ω c = Ω t = 0.1 , M = 0.1 .
Table 3. Average entropy number ( N s ¯ ) and average Bejan number ( B e ¯ ) with variation of β at fixed γ = 0.3 , P r = 0.7 , S c = 0.9 , R d = 0.3 , G r c = 0.8 , G r l = 0.5 , B r = 0.1 , B i = 0.1 , ε = 0.1 , R e = 0.7 , Ω c = Ω t = 0.1 , M = 0.1 .
β 123456
N s ¯ 0.07530.07320.07240.07190.07170.0715
B e ¯ 0.08910.08810.08770.08750.08730.0872
Table 4. Average entropy number ( N s ¯ ) and average Bejan number ( B e ¯ ) with variation of M at fixed γ = 0.3 , P r = 0.7 , S c = 0.9 , R d = 0.3 , G r c = 0.8 , G r l = 0.5 , B r = 0.1 , B i = 0.1 , ε = 0.1 , R e = 0.7 , Ω c = Ω t = 0.1 , B r = 0.1 .
Table 4. Average entropy number ( N s ¯ ) and average Bejan number ( B e ¯ ) with variation of M at fixed γ = 0.3 , P r = 0.7 , S c = 0.9 , R d = 0.3 , G r c = 0.8 , G r l = 0.5 , B r = 0.1 , B i = 0.1 , ε = 0.1 , R e = 0.7 , Ω c = Ω t = 0.1 , B r = 0.1 .
M 0.10.30.50.70.91.1
N s ¯ 0.07530.13270.19010.24750.30490.3623
B e ¯ 0.08910.05840.04520.03720.03170.0277
Table 5. Average entropy number ( N s ¯ ) and average Bejan number ( B e ¯ ) with variation of R d at fixed γ = 0.3 , P r = 0.7 , S c = 0.9 , B r = 0.1 , G r c = 0.8 , G r l = 0.5 , B r = 0.1 , B i = 0.1 , ε = 0.1 , R e = 0.7 , Ω c = Ω t = 0.1 , M = 0.1 .
Table 5. Average entropy number ( N s ¯ ) and average Bejan number ( B e ¯ ) with variation of R d at fixed γ = 0.3 , P r = 0.7 , S c = 0.9 , B r = 0.1 , G r c = 0.8 , G r l = 0.5 , B r = 0.1 , B i = 0.1 , ε = 0.1 , R e = 0.7 , Ω c = Ω t = 0.1 , M = 0.1 .
R d 0.52.54.56.58.510.5
N s ¯ 0.07530.08260.08780.09250.09680.1010
B e ¯ 0.08910.10790.12960.15000.16820.1841
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Louati, H.; Khan, S.; Mansoor, M.; Hilali, S.O.; Gargouri, A. Analysis of Entropy Generation via Non-Similar Numerical Approach for Magnetohydrodynamics Casson Fluid Flow with Joule Heating. Entropy 2024, 26, 702. https://doi.org/10.3390/e26080702

AMA Style

Louati H, Khan S, Mansoor M, Hilali SO, Gargouri A. Analysis of Entropy Generation via Non-Similar Numerical Approach for Magnetohydrodynamics Casson Fluid Flow with Joule Heating. Entropy. 2024; 26(8):702. https://doi.org/10.3390/e26080702

Chicago/Turabian Style

Louati, Hanen, Sajid Khan, Muavia Mansoor, Shreefa O. Hilali, and Ameni Gargouri. 2024. "Analysis of Entropy Generation via Non-Similar Numerical Approach for Magnetohydrodynamics Casson Fluid Flow with Joule Heating" Entropy 26, no. 8: 702. https://doi.org/10.3390/e26080702

APA Style

Louati, H., Khan, S., Mansoor, M., Hilali, S. O., & Gargouri, A. (2024). Analysis of Entropy Generation via Non-Similar Numerical Approach for Magnetohydrodynamics Casson Fluid Flow with Joule Heating. Entropy, 26(8), 702. https://doi.org/10.3390/e26080702

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