Next Article in Journal / Special Issue
Dissolution-Driven Convection in a Porous Medium Due to Vertical Axis of Rotation and Magnetic Field
Previous Article in Journal
A Note on Gerber–Shiu Function with Delayed Claim Reporting under Constant Force of Interest
Previous Article in Special Issue
Magneto Mixed Convection of Williamson Nanofluid Flow through a Double Stratified Porous Medium in Attendance of Activation Energy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exploration of Darcy–Forchheimer Flows of Non-Newtonian Casson and Williamson Conveying Tiny Particles Experiencing Binary Chemical Reaction and Thermal Radiation: Comparative Analysis

by
Sheniyappan Eswaramoorthi
1,
S. Thamaraiselvi
1 and
Karuppusamy Loganathan
2,3,*
1
Department of Mathematics, Dr. N.G.P. Arts and Science College, Coimbatore 641035, India
2
Department of Mathematics and Statistics, Manipal University Jaipur, Jaipur 303007, India
3
Research and Development Wing, Live4Research, Tiruppur 638106, India
*
Author to whom correspondence should be addressed.
Math. Comput. Appl. 2022, 27(3), 52; https://doi.org/10.3390/mca27030052
Submission received: 23 May 2022 / Revised: 15 June 2022 / Accepted: 17 June 2022 / Published: 20 June 2022

Abstract

:
This discussion intends to scrutinize the Darcy–Forchheimer flow of Casson–Williamson nanofluid in a stretching surface with non-linear thermal radiation, suction and heat consumption. In addition, this investigation assimilates the influence of the Brownian motion, thermophoresis, activation energy and binary chemical reaction effects. Cattaneo–Christov heat-mass flux theory is used to frame the energy and nanoparticle concentration equations. The suitable transformation is used to remodel the governing PDE model into an ODE model. The remodeled flow problems are numerically solved via the BVP4C scheme. The effects of various material characteristics on nanofluid velocity, nanofluid temperature and nanofluid concentration, as well as connected engineering aspects such as drag force, heat, and mass transfer gradients, are also calculated and displayed through tables, charts and figures. It is noticed that the nanofluid velocity upsurges when improving the quantity of Richardson number, and it downfalls for larger magnitudes of magnetic field and porosity parameters. The nanofluid temperature grows when enhancing the radiation parameter and Eckert number. The nanoparticle concentration upgrades for larger values of activation energy parameter while it slumps against the reaction rate parameter. The surface shear stress for the Williamson nanofluid is greater than the Casson nanofluid. There are more heat transfer gradient losses the greater the heat generation/absorption parameter and Eckert number. In addition, the local Sherwood number grows when strengthening the Forchheimer number and fitted rate parameter.

1. Introduction

Nowadays, heat transfer enrichment is a fascinating topic because of its numerous applications in engineering and industry. In many industrial processes, regular fluids (water, oil and ethylene glycol) are often employed. However, these fluids have a low heat transfer phenomenon because of their low thermal conductivity. To address this shortcoming, the nanometer-sized particle was mixed with regular fluids and enriches the regular fluid thermal conductivity; see [1,2,3]. This is the way of preparing the nanofluid, and this fluid has played an essential role in many fields such as solar water heating, heat exchangers, transformer cooling, cancer therapy, etc. Choi [4] was the first to publish the characteristics of nanoparticles, which were coupled with experimental evidence data. The nanofluid flow over a cylinder with suction was explored by Sheikholeslami [5]. It was uncovered that the local Nusselt number elevates when mounting the nanoparticle volume fraction values. Ramana Reddy et al. [6] addressed the time-dependent MHD flow of nanofluid past a slendering surface. It was detected that the fluid temperature progressed when enhancing the Brownian motion parameter. Makinde et al. [7] employed the impact of Brownian motion and thermophoresis effects of MHD flow of nanofluid past a heated surface. It was noted that the heat transfer rate decays when upturning the quantity of the Brownian motion parameter. The consequences of Brownian motion and thermophoresis of stagnation point flow of nanofluid past a non-uniform cylinder were presented by Shafey et al. [8]. It was noticed that the heat transfer gradient slumps when rising the thermophoresis quantity. Rasheed et al. [9] addressed the MHD flow of water-based nanofluids with convective heating conditions. It was seen that the thermophoretic parameter improves the thickness of the thermal boundary layer.
The non-Newtonian fluid has stimulated various scientists to investigate the events of heat-mass transport because of its necessary part in industrial and engineering processes, such as drilling muds, polymer extrusion, optical fibers, polymer production, etc. The non-Newtonian fluid defies the Newton’s viscosity law. To deal with the huge nature of the rheological behavior of such fluids, several non-Newtonian models have been devised. Casson fluid is one of the non-Newtonian type models, and at the infinite non-linear shear rate, the fluid material’s yield stress does not push flow, and it has zero viscosity. The MHD flow of Casson nanofluid past a heated surface with viscous dissipation was analyzed by Alotaibi et al. [10]. It was found that the drag force coefficient decays when strengthening the Casson parameter. Nayak et al. [11] provide the impact of the triple diffusive bioconvective flow of Casson nanofluid past a sheet. It was observed that the wall motile micro-organism decimates when developing the Casson parameter. Entropy optimization of MHD flow of Casson nanofluid over a stretching surface with convective heating and mass conditions was illustrated by Butt et al. [12]. It was proved that the Casson fluid parameter leads to a slow down of the entropy production. Ibrahim et al. [13] discovered the chemically reactive MHD flow of Casson nanofluid past a stretching surface with viscous dissipation. It was noticed that the nanoparticle concentration profile decreases when raising the Casson parameter. The multiple slip effects of a Casson nanofluid on a stretching surface were numerically performed by Afify [14], and he proved that the mass transfer gradient enriches when enhancing the Casson parameter. Varun Kumar et al. [15] scrutinized the MHD chemically reactive flow of Casson nanoliquid past a curved stretching sheet. It was noted that liquid velocity depresses when enhancing the Casson parameter. The 2D flow of Casson nanofluid on a thin moving needle was examined by Naveen Kumar et al. [16], and they proved that the thermophoresis parameter improves the mass transfer rate. Gohar et al. [17] studied the Darcy–Forchheimer flow of Casson hybrid nanofluid on a curved surface. It was detected that the Casson parameter suppresses the hybrid nanofluid motion.
Williamson fluid is also the non-Newtonian division model, which exhibits the shear thinning property; that is, the fluid viscosity decays when rising the shear stress rate. Waqas et al. [18] examined the MHD flow of Williamson nanofluid past a heated wedge. It was revealed that the wall shear stress downfalls when mounting the Weissenberg number. The MHD flow of Williamson nanofluid past a porous stretching surface with suction was presented by Li et al. [19]. It was noted that the Williamson parameter leads to depromoting the friction drag. Ahmed et al. [20] presented the consequences of MHD Williamson nanofluid flow on an exponentially porous stretching surface. It was uncovered that the fluid speed depresses when escalating the Williamson parameter. The 2D flow of Williamson fluid over a cylinder was addressed by Iqbal et al. [21], and it was acknowledged that the skin friction coefficient decreases as the Weissenberg number increases. Gorla and Gireesha [22] demonstrated the convective heat transport analysis of a Williamson nanofluid past the stretching surface. It was noticed that the nanofluid volume fraction intensifies when heightening the Williamson parameter. The MHD flow of Williamson nanofluid past a heated stretching surface was examined by Srinivasulu and Goud [23]. It was concluded that the heat transfer gradient downfalls when promoting the quantity of the Williamson parameter.
In the last few decades, many researchers have focused on studying the thermal radiation effect because the consequences of thermal radiation in flow structures are helpful in atomic reactors, spacecraft, ship compressors, and solar radiation. Most of the investigation is based on linearized Rosseland approximation; however, this concept is applicative when the temperature distinction between ambient and fluid is small. However, on many industrial occasions, this difference is enormous. So, a non-linearized Rosseland approximation is introduced to overcome this restriction. MHD Casson nanofluid in a bi-directional heated stretching surface with non-linear radiation was deliberated by Mahanta et al. [24]. It was detected that the temperature ratio parameter leads to enriching the entropy generation profile. Humane et al. [25] scrutinized the thermally radiative MHD Casson–Williamson nanofluid flow on a porous stretching surface with a chemical reaction. MHD heat-generating Casson nanofluid through a thin needle with non-linear thermal radiation was examined by Akinshilo et al. [26]. Ghasemi et al. [27] numerically solved the non-linear thermal radiative flow of nanofluid with a magnetic field via the spectral relaxation method. It was noted that the nanofluid concentration upsurges when enhancing the thermal radiation parameter. The bio-convective flow of Carreau nanofluid with non-linear thermal radiation with a magnetic dipole was presented by Imran et al. [28]. It was shown that the thermal boundary layer thickens when the temperature ratio parameter is large. Bhatti et al. [29] demonstrated the impact of MHD flow of Williamson nanofluid through a shrinking porous sheet. The problem of non-linear radiative flow of nanofluid with the inclined magnetic field was numerically solved via the finite difference method by Mahanthesh and Thriveni [30]. Their results clearly explain that the fluid temperature ascends when upgrading the quantity of the thermal radiation parameter. Cao et al. [31] investigated the non-linear thermal radiative flow of a ternary-hybrid nanofluid with partial slip. The 3D radiative flow of C u / A g -water-based nanofluid with entropy optimization was illustrated by Eswaramoorthi et al. [32], and they detected that the Bejan number rises as the radiation parameter enhances.
The smallest amount of energy necessary to start a chemical reaction is known as activation energy. This conception was initiated by Arrhenius in 1889, and this incident has plentiful appliances in geothermal engineering, water emulsions, oil emulsion and food processing. Shah et al. [33] addressed the chemically reactive flow of Casson nanofluid with activation energy and radiation, and they found that activation energy leads to magnifying the nanoparticle concentration. The 3D time-dependent flow of Williamson nanofluid with heat generation and the activation energy was inspected by Aziz et al. [34]. Their findings show that the higher chemical reaction parameter suppresses the nanofluid concentration. Kalaivanan et al. [35] discussed the Arrhenius activation energy and non-linear thermal impacts of second-grade nanofluid past a stretching surface. It was exposed that the heat transfer gradient weakens when strengthening the exponential fitted rate. The MHD flow of Casson nanofluid over a stretching cylinder with Arrhenius activation energy was examined by Zeeshan et al. [36]. It was seen that the nanoparticle concentration enhances for strengthening the activation energy parameter. Tayyab et al. [37] securitized the consequences of Darcy–Forchheimer flow of 3D nanofluid on a sheet with activation energy. The 3D Darcy–Forchheimer flow past a porous space with the presence of Arrhenius activation energy was presented by Rashid et al. [38]. It was revealed that the reaction rate leads to a decline in the nanoparticle concentration profile. Alsaadi et al. [39] elucidated the flow of MHD WNF with the influence of Arrhenius activation energy. The impact of activation energy of a second-grade nanofluid on a surface with heat source/sink was analyzed by Punith Gowda et al. [40]. Varun Kumar et al. [41] studied the impact of Arrhenius activation energy on a hybrid nanofluid past a curved surface. It was proved that the nanofluid concentration improves when escalating the activation energy parameter. The MHD flow of Williamson nanofluid with activation energy was investigated by Tamilzharasan et al. [42], and they found that the activation energy parameter improves the heat transfer rate.
In light of the above analysis, no research articles provide the impact of the non-linear thermal radiative flow of Casson–Williamson nanofluid over a heated stretchy plate with activation energy. In addition, zero nanoparticle mass flux and Cattaneo–Christov heat-mass flux conditions are included in our study. This research has implications for thermal sciences, food processing, chemical engineering, polymer extrusion, and many other fields in which heat conduction and convection are improved. In the limiting scenarios, the calculated values derived from specified parameters are consistent with existing findings in the literature, while tables and graphs have been built and explained to spread the responses of dimensionless quantities. This type of flow model is used in food processing, heavy mechanical apparatus, enzymes, ceramic processing, heating/cooling processes, etc.

2. Mathematical Formulation

Let us consider the chemically reactive 2D Darcy–Forchheimer flow of Casson–Williamson nanofluid past a stretchy plate. The x-axis is in the stretching direction and the y-axis is perpendicular to it. The uniform magnetic field of strength B 0 is applied in the y-direction, and the induced magnetic effect is neglected because of the small quantity of Reynolds number. Moreover, the flow is disclosed under the consequences of Arrhenius activation energy; suction/injection, viscous dissipation and non-linear thermal radiation are all taken into our account. The zero nanoparticle mass flux condition is assumed on the surface of the sheet. Let T w and C w represent the temperature and nanoparticle concentration, which are higher than the free-stream temperature ( T ) and nanoparticle concentration ( C ) , see Figure 1.
The Cauchy stress tensor of Williamson fluid is expressed as S = p I + τ 1 , where τ 1 = μ + μ 0 μ 1 Γ 1 γ 1 * A 1 ; here, τ 1 is the extra stress tensor, μ 0 is the limiting viscosity at zero shear rate, μ is the limiting viscosity at infinity shear rate, Γ 1 > 0 is the time constant and A 1 is the Rivlin–Ericson tensor. The simplified form of the extra stress tensor is τ 1 = + μ 0 1 Γ 1 γ 1 * A 1 ; see [43].
Similarly, Casson fluid flow is
τ i j = 2 μ n f + Q y 2 π k i j , π > π c 2 μ n f + Q y 2 π c k i j , π < π c
here, Q y is the yield stress of fluid, k i j is the ( i , j ) th laceration direction component rate, π = k i j k i j is the product of the component of rate of deformation with itself and π c is the critical value of the product of the component of the strain tensor rate with itself; see [44].
The flow model may be described as follows using the given assumptions; see Mustafa et al. [45].
u ¯ x ¯ + v ¯ y ¯ = 0
u u ¯ x ¯ + v u ¯ y ¯ = ν 1 + 1 β u ¯ y y ¯ + 2 Γ ν u ¯ y ¯ u ¯ y y ¯ ν k 1 u ¯ C b x ¯ k 1 u ¯ 2 σ B 0 2 u ¯ ρ f + 1 ρ f ( 1 C ¯ ) ρ f β ( T ¯ T ¯ ) ( ρ p ρ f ) ( C ¯ C ¯ ) g ,
u ¯ T ¯ x ¯ + v ¯ T ¯ y ¯ + λ T Ω T = α T ¯ y y ¯ + 1 ρ C p 16 σ * 3 k * y ¯ T ¯ 3 T ¯ y ¯ + Q ρ f C p ( T ¯ T ¯ ) + τ D B T ¯ y ¯ C ¯ y ¯ + D T T ¯ T ¯ y ¯ 2 + μ ρ C p 1 + 1 β u ¯ y ¯ 2
u ¯ C ¯ x ¯ + v ¯ C ¯ y ¯ + λ C Ω C = D B C ¯ y y ¯ + D T T ¯ C ¯ y y ¯ k r 2 ( C ¯ C ¯ ) T ¯ T ¯ n e x p E a κ T ¯
where
Ω T = u u ¯ x ¯ T ¯ x ¯ + v v ¯ y ¯ T ¯ y ¯ + u ¯ 2 T ¯ x x ¯ + v ¯ 2 T ¯ y y ¯ + 2 u v ¯ T ¯ x y ¯ + u v ¯ x ¯ T ¯ y ¯ + v u ¯ y ¯ T ¯ x ¯ Ω C = u u ¯ x ¯ C ¯ x ¯ + v v ¯ y ¯ C ¯ y ¯ + u ¯ 2 C ¯ x x ¯ + v ¯ 2 C ¯ y y ¯ + 2 u v ¯ C ¯ x y ¯ + u v ¯ x ¯ C ¯ y ¯ + v u ¯ y ¯ C ¯ x ¯
The boundary conditions are
u ¯ = U ¯ ω + L 1 + 1 β + Γ u ¯ y ¯ u ¯ y ¯ ; v ¯ = V ¯ ω T ¯ = T ¯ ω , D B C ¯ y ¯ + D T T ¯ T ¯ y ¯ = 0 a s y ¯ = 0 u ¯ 0 , T ¯ T ¯ , C ¯ C ¯ a s y ¯
From the above Equations (Equations (2) and (5)), β & Γ 0 is treated as a Williamson fluid model and β & Γ = 0 is treated as a Casson fluid model.
The dimensionless parameters are
ϖ = a ν y ¯ ; u ¯ = a x ¯ υ 1 ; v ¯ = a ν υ 1 ( ϖ ) υ 2 ( ϖ ) = T ¯ T ¯ T ¯ ω T ¯ ; υ 3 ( ϖ ) = C ¯ C ¯ C ¯
Substituting Equation (6) into Equations (2)–(4), we obtain
1 + 1 β υ 1 ( ϖ ) υ 1 2 ( ϖ ) + υ 1 ( ϖ ) υ 1 ( ϖ ) + W e υ 1 ( ϖ ) υ 1 ( ϖ ) λ υ 1 ( ϖ ) F r υ 1 2 ( ϖ ) M υ 1 ( ϖ ) + R i ( υ 2 ( ϖ ) N r υ 3 ( ϖ ) ) = 0
1 P r υ 2 ( ϖ ) + 1 P r 4 3 R [ θ w 1 3 { υ 2 ( ϖ ) υ 2 3 ( ϖ ) + 3 υ 2 2 ( ϖ ) υ 2 2 ( ϖ ) } + 3 θ w 1 2 { υ 2 ( ϖ ) υ 2 2 ( ϖ ) + 2 υ 2 ( ϖ ) υ 2 ( ϖ ) 2 } + 3 θ w 1 { υ 2 ( ϖ ) υ 2 ( ϖ ) + υ 2 2 ( ϖ ) } + υ 2 ( ϖ ) ] + υ 1 ( ϖ ) υ 2 ( ϖ ) Γ T { υ 1 ( ϖ ) υ 1 ( ϖ ) υ 2 ( ϖ ) + υ 1 2 ( ϖ ) υ 2 ( ϖ ) } + H g υ 2 + 1 + 1 β E c υ 2 2 + N b υ 2 υ 3 + N t υ 2 2 = 0
1 S c υ 3 ( ϖ ) + υ 1 ( ϖ ) υ 3 ( ϖ ) Γ C { υ 1 ( ϖ ) υ 1 ( ϖ ) υ 3 ( ϖ ) + υ 1 2 ( ϖ ) υ 3 ( ϖ ) } + 1 S c N t N b υ 2 ( ϖ )
σ * * ( 1 + δ υ 2 ( ϖ ) ) n υ 3 ( ϖ ) e x p E 1 + δ υ 2 ( ϖ ) = 0
The covered boundary conditions are
υ 1 ( 0 ) = f w , υ 1 ( 0 ) = 1 + K 1 + 1 β + W e 2 υ 1 ( 0 ) υ 1 ( 0 ) , υ 2 ( 0 ) = 1 , N b υ 3 ( 0 ) + N t υ 2 ( 0 ) = 0 υ 1 ( ) 0 , υ 2 ( ) 0 , υ 3 ( ) 0
The non-dimensional form of skin friction coefficient, local Nusselt number and local Sherwood number are expressed as
1 2 C f R e = 1 + 1 β υ 1 ( 0 ) + W e 2 υ 1 2 ( 0 ) ; N u R e = 1 + 4 3 R 1 + θ w 1 υ 2 ( 0 ) 3 υ 2 ( 0 ) S h R e = N b N t υ 2 ( 0 )

3. Numerical Solution

The derived ODE models (7)–(9) along with the conditions (10) are numerically solved by applying the MATLAB bvp4c scheme. Initially, the higher-order terms are converted into first-order terms, see [46,47].
Let υ 1 = y 1 , υ 1 = y 2 , υ 1 = y 3 , υ 1 = y 3 , υ 2 = y 4 , υ 2 = y 5 , υ 2 = y 5 , υ 3 = y 6 , υ 3 = y 7 , υ 3 = y 7 .
y 1 = y 2 y 2 = y 3 y 3 = y 2 2 y 1 y 3 + λ y 2 + F r y 2 2 + M y 2 R i [ y 4 N r y 6 ] ( 1 + 1 β ) + W e y 3 y 4 = y 5 A = y 1 y 5 + Γ T y 1 y 2 y 5 1 P r 4 3 R [ ( 3 θ n 1 ) 3 y 4 2 y 5 2 + 6 ( θ n 1 ) 2 y 4 y 5 2 + 3 ( θ n 1 ) y 5 2 ] B = H g y 4 N b y 5 y 7 N t y 5 2 1 + 1 β E c y 3 2 y 5 = A + B 1 P r [ 1 + 4 3 R [ 1 + ( θ n 1 ) 3 y 4 3 + 3 ( θ n 1 ) 2 y 4 2 + 3 ( θ n 1 ) y 4 ] ] Γ T y 1 2 y 6 = y 7 C = y 1 y 5 + Γ T y 1 y 2 y 5 1 P r 4 3 R [ ( 3 θ n 1 ) 3 y 4 2 y 5 2 + 6 ( θ n 1 ) 2 y 4 y 5 2 + 3 ( θ n 1 ) y 5 2 ] D = H g y 4 N b y 5 y 7 N t y 5 2 1 + 1 β E c y 3 2 E = 1 P r [ 1 + 4 3 R [ 1 + ( θ n 1 ) 3 y 4 3 + 3 ( θ n 1 ) 2 y 4 2 + 3 ( θ n 1 ) y 4 ] ] Γ T y 1 2 y 7 = y 1 y 7 + Γ C y 1 y 2 y 7 1 S c N t N b C + D E + σ * * ( 1 + δ y 4 ) n y 6 e x p E 1 + δ y 4 1 S c Γ C y 1 2
With the conditions
y 1 ( 0 ) = f w , y 2 ( 0 ) = 1 + K 1 + 1 β + W e 2 y 3 ( 0 ) y 3 ( 0 ) , y 2 ( ) = 0 , y 4 ( 0 ) = 1 , N b y 7 ( 0 ) + N t y 5 ( 0 ) = 0 , y 4 ( ) = 0 , y 6 ( ) = 0
We implemented the MATLAB bvp4c scheme to find the numerical solution for the above problem with maximum error is 10 5 and step size is 0.05.

4. Results and Discussion

This segment provides the details about the changes of velocity, temperature, naofluid concentration, skin friction coefficient (SFC), local Nusselt number (LNN) and local Sherwood number (LSN) for different flow parameters through graphs and tables. The consequences of f w , λ , F r , M, R i and N r (Table 1), R, H g , E c , Γ T , N b , N t and θ n (Table 2) and Γ c , σ * * , δ , n and E (Table 3) on SFC, LNN and LSN for Casson–Williamson nanofluid are deliberated in Table 1, Table 2 and Table 3. Table 4 provides the comparison of our numerical results to Mustafa et al. [45] and found excellent agreement. It is detected that there is SFC shrinkage when boosting the f w , λ , F r , M, N r , Γ T , N t Γ C and E, and it upturns when enriching the R i R, H g , E c , N b , θ n , σ * * , δ and n. The heat transfer gradient (HTG) loses when strengthening the quantity of λ , F r , M, N r , H g , E c , N t , Γ C , δ , n and E, and it upturns when enhancing the amount of f w , R i , R, Γ T , N b , θ n and σ * * . The LSN proliferate when mounting the quantity of λ , F r , M, N r , R, H g , N b , θ n , Γ c and δ . The quite opposite trend is obtained when changing the presence of f w , R i , E c , Γ T , N t , σ * * , n and E.
Figure 2a–d display the variances of fluid velocity versus R i (a), f w (b), λ (c) and M (d). It is clearly shown that the fluid speed enhances when heightening the quantity of R i and it depresses when mounting the quantity of f w , λ and M for both fluids. In addition, the velocity of the Casson nanofluid is low near the plate and high away from the plate compared to Williamson nanofluid. Physically, a larger quantity of M generates a drag force named the Lorentz force. This force leads to suppressing the fluid movement on a plate surface, and this causes a decline in the fluid speed and thinner momentum boundary layer. The fluid temperature variations on R(a), H g (b), Γ T (c) and E c (d) for both fluids are illustrated in Figure 3a–d. It is acknowledged that the fluid warmness escalates when enhancing the R and H g values, and it suppresses when the Γ T and E c values are rising. Physically, the presence of a radiation parameter has enriched the fluid thermal state, thereby strengthening the fluid warmness and thicker thermal boundary layer thickness. In addition, the greater availability of Eckert number creates a more robust viscous dissipation effect, which enriches the fluid warmness. Figure 4a–d show the consequences of σ * * (a), f w (b), E(c) and N t (d) on nanoparticle concentration profile. It is seen that the nanoparticle concentration reduces when raising the values σ * * . A opposite behavior occurs for varying the values of f w , E and N t . The skin friction coefficient for a distinct combination of M, λ and f w is presented in Figure 5a–d. It is found that the surface shear stress decays when enhancing the magnetic field and porosity parameter for both f w values. In addition, the Williamson nanofluid has a greater skin friction coefficient value than the Casson nanofluid. Figure 6a–d portrayed the changes of local Nusselt number for a distinct combination of values of M, λ and f w . It is concluded that the heat transfer gradient slowly depresses when increasing the magnetic field and porosity parameters for both f w values. The local Sherwood number for various combination of values of M, λ and f w is shown in Figure 7a–d. It is seen from these figures that the local Sherwood number slowly depresses when increasing the magnetic field and porosity parameters for both f w values.

5. Conclusions

Here, Brownian motion and the thermophoresis impact of the non-linear radiative flow of C-WNF in a Darcy–Forchheimer porous space with suction and heat consumption is investigated. The present investigation includes the consequences of activation energy and binary chemical reaction. The governing mathematical models are numerically solved by the bvp4c algorithm with MATLAB. The main outcomes of our discussion are as follows:
  • The fluid speed enhances for Richardson number but it slows against porosity, suction/injection and magnetic field parameters.
  • The fluid becomes more warmed as the radiation, heat generation parameters and Eckert number increase.
  • The nanoparticle concentration enhances upon strengthening the suction/injection and thermophoresis parameters and it downfalls upon escalating the reaction rate.
  • The skin friction reduces after enriching the Forchheimer number, porosity and magnetic field parameters.
  • The heat transfer gradient increases when escalating the values of radiation parameter and it downturns against radiation and heat generation parameters.
  • The mass transfer gradient enhances upon heightening the Brownian motion parameter and it weakens against the thermophoresis parameter.
  • In the future, we extend this flow model through the Riga plate with the convective heating condition.

Author Contributions

Data curation, S.E.; Formal analysis, K.L.; Investigation, S.T.; Methodology, S.E. and S.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:

Nomenclature

Symbols Description
a , b Positive   constants
B 0 Magnetic   field   strength ( T )
C Fluid   concentration ( m o l L 1 )
C b Drag   coefficient
C p Specific   heat ( J k g 1 K 1 )
C Ambient   fluid   concentration
D B Brownian   diffusion   coefficient ( m 2 s 1 )
D T Thermophoretic   diffusion   coefficient
E ( = E a / k T ) Activation   energy   parameter
E a Activation   energy
F r ( = C b / k 1 ) Forchheimer   number
f Dimensionless   velocity
g Acceleration   due   to   gravity ( m s 2 )
G r x ( = ( g β ( 1 C ) ( T ω T ) x 3 / ν 2 ) Local   Grashof   number
H g ( = Q ρ f C p a ) Heat   generation / absorption   parameter
k Thermal   conductive ( W m 1 K 1 )
k 1 Permeability   of   porous   medium ( m 2 )
k * Mean   absorption   coefficient
k r Reaction   rate
M ( = σ B 0 2 ρ f a ) Magnetic   parameter
n Fitted   rate   or   stretching   sheet   index   parameter
E c Eckert   number
N b ( = τ D B ( C w C ) ν ) Brownian   diffusion   parameter
N t ( = τ D T ( T ω T ) T ν ) Thermophoresis   parameter
P r ( = ν α = m 2 s 1 m 2 s 1 = 1 ) Prandtl   number
QHeat generation/absorption coefficient
( J M 1 m 3 s 1 )
R ( = 4 σ * T 3 k k * ) Thermal radiation
R e x ( = U ω x ν ) Local Reynolds number
R i ( = G r x R e x 2 = g β ( 1 C ) ( T ω T ) a 2 x ) Richardson number
S c ( = ν D B ) Schmith number
TFluid temperature ( K )
T Ambient temperature ( K )
u , v Velocity components ( m s 1 )
U ω Stretching surface velocity ( m s 1 )
W e ( = Γ x 2 a 3 / ν ) = s m 2 s 3 m 2 s 1 = c o n s t a n t ) Weissenberg number
x , y Direction coordinates ( m )
Greek Symbols Description
α Thermal diffusivity ( m 2 s 1 )
β Casson parameter
δ ( = T ω T T ) Temperature difference parameter
Γ Williamson parameter or time constant
Γ T ( = a λ T ) Thermal relaxation parameter
Γ C ( = a λ C ) Solute relaxation parameter
λ Local porosity parameter
λ C Relaxation time of mass flux
λ T Relaxation time of heat flux
ν Kinetic viscosity ( m 2 s 1 )
ϕ Non-dimensional nanofluid concentration
ρ f Fluid density ( k g m 2 )
σ Electrical conductivity ( S m 1 )
σ * Stefan-Boltzmann constant ( W m 2 K 4 )
σ * * ( = k r 2 a ) Dimensionless reaction rate
τ Heat capacity ratio
θ Non-dimensional temperature
θ n Temperature ratio parameter
η Similarity variable
μ Dynamic viscosity ( K g m 1 s 1 )
pDust phase
Fluid properties at ambient condition

References

  1. Elnaqeeb, T.; Animasaun, I.L.; Shah, N.A. Ternary-hybrid nanofluids: Significance of suction and dual-stretching on three-dimensional flow of water conveying nanoparticles with various shapes and densities. Z. Naturforsch. A 2021, 76, 231–243. [Google Scholar] [CrossRef]
  2. Islam, T.; Yavuz, M.; Parveen, N.; Fayz-Al-Asad, M. Impact of non-uniform periodic magnetic field on unsteady natural convection flow of nanofluids in square enclosure. Fractal Fract 2022, 6, 101. [Google Scholar] [CrossRef]
  3. Saleem, S.; Animasaun, I.L.; Yook, S.J.; Al-Mdallal, Q.M.; Shah, N.A.; Faisal, M. Insight into the motion of water conveying three kinds of nanoparticles shapes on a horizontal surface: Significance of thermo-migration and Brownian motion. Surf. Interfaces 2022, 30, 101854. [Google Scholar] [CrossRef]
  4. Choi, S.U.S. Enhancing thermal conductivity of fluids with nanoparticles. ASME-Publications-Fed 1995, 231, 99–106. [Google Scholar]
  5. Sheikholeslami, M. Effect of uniform suction on nanofluid flow and heat transfer over a cylinder. J. Braz. Soc. Mech. Sci. Eng. 2015, 37, 1623–1633. [Google Scholar] [CrossRef]
  6. Ramana Reddy, J.V.; Sugunamma, V.; Sandeep, N. Thermophoresis and Brownian motion effects on unsteady MHD nanofluid flow over a slendering stretching surface with slip effects. Alex. Eng. J. 2018, 57, 2465–2473. [Google Scholar] [CrossRef]
  7. Makinde, O.D.; Mabood, F.; Khan, W.A.; Tshehla, M.S. MHD flow of a variable viscosity nanofluid over a radially stretching convective surface with radiative heat. J. Mol. Liq. 2016, 219, 624–630. [Google Scholar] [CrossRef]
  8. Shafey, A.M.E.; Alharbi, F.M.; Javed, A.; Abbas, N.; ALrafai, H.A.; Nadeem, S.; Issakhov, A. Theoretical analysis of Brownian and thermophoresis motion effects for Newtonian fluid flow over nonlinear stretching cylinder. Case Stud. Therm. Eng. 2021, 28, 101369. [Google Scholar] [CrossRef]
  9. Rasheed, H.U.; Islam, S.; Khan, Z.; Khan, J.; Mashwani, W.K.; Abbas, T.; Shah, Q. Computational analysis of hydromagnetic boundary layer stagnation point flow of nano liquid by a stretched heated surface with convective conditions and radiation effect. Adv. Mech. Eng. 2021, 13, 16878140211053142. [Google Scholar] [CrossRef]
  10. Alotaibi, H.; Althubiti, S.; Eid, M.R.; Mahny, K.L. Numerical treatment of MHD flow of Casson nanofluid via convectively heated non-linear extending surface with viscous dissipation and suction/injection effects. Comput. Mater. Continua 2021, 66, 229–245. [Google Scholar] [CrossRef]
  11. Nayak, M.K.; Prakash, J.; Tripathi, D.; Pandey, V.S.; Shaw, S. 3D Bioconvective multiple slip flow of chemically reactive Casson nanofluid with gyrotactic microorganisms. Heat Transf. Res. 2020, 49, 135–153. [Google Scholar] [CrossRef]
  12. 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]
  13. Ibrahim, S.M.; Lorenzini, G.; Vijaya Kumar, P.; Raju, C.S.K. Influence of chemical reaction and heat source on dissipative MHD mixed convection flow of a Casson nanofluid over a nonlinear permeable stretching sheet. Int. J. Heat Mass Transf. 2017, 111, 346–355. [Google Scholar] [CrossRef]
  14. Afify, A.A. The influence of slip boundary condition on Casson nanofluid flow over a stretching sheet in the presence of viscous dissipation and chemical reaction. Math. Probl. Eng. 2017, 2017, 3804751. [Google Scholar] [CrossRef] [Green Version]
  15. Varun Kumar, R.S.; Dhananjaya, P.G.; Naveen Kumar, R.; Punith Gowda, R.J.; Prasannakumara, B.C. Modeling and theoretical investigation on Casson nanofluid flow over a curved stretching surface with the influence of magnetic field and chemical reaction. Int. J. Comput. Methods Eng. Sci. Mech. 2022, 23, 12–19. [Google Scholar] [CrossRef]
  16. Naveen Kumar, R.; Punith Gowda, R.J.; Madhukesh, J.K.; Prasannakumara, B.C.; Ramesh, G.K. Impact of thermophoretic particle deposition on heat and mass transfer across the dynamics of Casson fluid flow over a moving thin needle. Phys. Scr. 2021, 96, 075210. [Google Scholar] [CrossRef]
  17. Khan, T.S.; Sene, N.; Mouldi, A.; Brahmia, A. Heat and mass transfer of the Darcy–Forchheimer Casson hybrid nanofluid flow due to an extending curved surface. J. Nanomat. 2022, 2022, 3979168. [Google Scholar]
  18. Waqas, S.H.; Khan, S.U.; Imran, M.; Bhatti, M.M. Thermally developed Falkner–Skan bioconvection flow of a magnetized nanofluid in the presence of a motile gyrotactic microorganism: Buongiorno’s nanofluid model. Phys. Scr. 2019, 94, 115304. [Google Scholar] [CrossRef]
  19. Li, Y.X.; Alshbool, M.H.; Lv, Y.P.; Khan, I.; Khan, M.R.; Issakhov, A. Heat and mass transfer in MHD Williamson nanofluid flow over an exponentially porous stretching surface. Case Stud. Therm. Eng. 2021, 26, 100975. [Google Scholar] [CrossRef]
  20. Ahmed, K.; Akbar, T. Numerical investigation of magnetohydrodynamics Williamson nanofluid flow over an exponentially stretching surface. Adv. Mech. Eng. 2021, 13, 168781402110198. [Google Scholar] [CrossRef]
  21. Iqbal, W.; Naeem, M.N.; Jalil, M. Numerical analysis of Williamson fluid flow along an exponentially stretching cylinder. AIP Adv. 2019, 9, 055118. [Google Scholar] [CrossRef] [Green Version]
  22. Gorla, R.S.R.; Gireesha, B.J. Dual solutions for stagnation-point flow and convective heat transfer of a Williamson nanofluid past a stretching/shrinking sheet. Heat Mass Transf. 2016, 52, 1153–1162. [Google Scholar] [CrossRef]
  23. Srinivasulu, T.; Goud, B.S. Effect of inclined magnetic field on flow, heat and mass transfer of Williamson nanofluid over a stretching sheet. Case Stud. Therm. Eng. 2021, 23, 100819. [Google Scholar] [CrossRef]
  24. Mahanta, G.; Das, M.; Nayak, M.K.; Shaw, S. Irreversibility Analysis of 3D Magnetohydrodynamic Casson Nanofluid Flow Past Through Two Bi-Directional Stretching Surfaces with Nonlinear Radiation. J. Nanofluids 2021, 10, 316–326. [Google Scholar] [CrossRef]
  25. Humane, P.P.; Patil, V.S.; Patil, A.B. Chemical reaction and thermal radiation effects on magnetohydrodynamics flow of Casson–Williamson nanofluid over a porous stretching surface. Proc. Inst. Mech. Eng. E: J. Process Mech. Eng. 2021, 235, 2008–2018. [Google Scholar] [CrossRef]
  26. Akinshilo, A.T.; Mabood, F.; Ilegbusi, A.O. Heat generation and nonlinear radiation effects on MHD Casson nanofluids over a thin needle embedded in porous medium. Int. Commun. Heat Mass Transf. 2021, 127, 105547. [Google Scholar] [CrossRef]
  27. Ghasemi, S.E.; Mohsenian, S.; Gouran, S.; Zolfagharian, A. A novel spectral relaxation approach for nanofluid flow past a stretching surface in presence of magnetic field and nonlinear radiation. Results Phys. 2022, 32, 105141. [Google Scholar] [CrossRef]
  28. Imran, M.; Farooq, U.; Muhammad, T.; Khan, S.U.; Waqas, H. Bioconvection transport of Carreau nanofluid with magnetic dipole and nonlinear thermal radiation. Case Stud. Therm. Eng. 2021, 26, 101129. [Google Scholar] [CrossRef]
  29. Bhatti, M.M.; Abbas, T.; Rashidi, M.M. Numerical study of entropy generation with nonlinear thermal radiation on magnetohydrodynamics non-Newtonian nanofluid through a porous shrinking sheet. J. Magn. 2016, 21, 468–475. [Google Scholar] [CrossRef] [Green Version]
  30. Mahanthesh, B.; Thriveni, K. Significance of inclined magnetic field on nano-bioconvection with nonlinear thermal radiation and exponential space based heat source: A sensitivity analysis. Eur. Phys. J. Spec. Top. 2021, 230, 1487–1501. [Google Scholar] [CrossRef]
  31. Cao, W.; Animasaun, I.L.; Yook, S.J.; Oladipupo, V.A.; Ji, X. Simulation of the dynamics of colloidal mixture of water with various nanoparticles at different levels of partial slip: Ternary-hybrid nanofluid. Int. Commun. Heat Mass Transf 2022, 135, 106069. [Google Scholar] [CrossRef]
  32. Eswaramoorthi, S.; Divya, S.; Faisal, M.; Namgyel, N. Entropy and heat transfer analysis for MHD flow of-water-based nanofluid on a heated 3D plate with nonlinear radiation. Math. Probl. Eng. 2022, 2022, 7319988. [Google Scholar] [CrossRef]
  33. Shah, Z.; Kumam, P.; Deebani, W. Radiative MHD Casson Nanofluid Flow with Activation energy and chemical reaction over past nonlinearly stretching surface through Entropy generation. Sci. Rep. 2020, 10, 4402. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Aziz, S.; Ahmad, I.; Khan, S.U.; Ali, N. A three-dimensional bioconvection Williamson nanofluid flow over bidirectional accelerated surface with activation energy and heat generation. Int. J. Mod. Phys. B 2021, 35, 2150132. [Google Scholar] [CrossRef]
  35. Kalaivanan, R.; Vishnu Ganesh, N.; Al-Mdallal, Q.M. An investigation on Arrhenius activation energy of second grade nanofluid flow with active and passive control of nanomaterials. Case Stud. Therm. Eng. 2020, 22, 100774. [Google Scholar] [CrossRef]
  36. Zeeshan, A.; Mehmood, O.U.; Mabood, F.; Alzahrani, F. Numerical analysis of hydromagnetic transport of Casson nanofluid over permeable linearly stretched cylinder with Arrhenius activation energy. Int. Commun. Heat Mass Transf. 2022, 130, 105736. [Google Scholar] [CrossRef]
  37. Tayyab, M.; Siddique, I.; Jarad, F.; Ashraf, M.K.; Ali, B. Numerical solution of 3D rotating nanofluid flow subject to Darcy–Forchheimer law, bio-convection and activation energy. S. Afr. J. Chem. Eng. 2022, 40, 48–56. [Google Scholar] [CrossRef]
  38. Rashid, S.; Hayat, T.; Qayyum, S.; Ayub, M.; Alsaedi, A. Three dimensional rotating Darcy–Forchheimer flow with activation energy. Int. J. Numer. Methods Heat Fluid Flow 2018, 29, 935–948. [Google Scholar] [CrossRef]
  39. Alsaadi, F.E.; Hayat, T.; Khan, M.I.; Alsaadi, F.E. Heat transport and entropy optimization in flow of magneto-Williamson nanomaterial with Arrhenius activation energy. Comput. Methods Progr. Biomed. 2020, 183, 105051. [Google Scholar] [CrossRef]
  40. Punith Gowda, R.J.; Naveen Kumar, R.; Jyothi, A.M.; Prasannakumara, B.C.; Sarris, I.E. Impact of binary chemical reaction and activation energy on heat and mass transfer of marangoni driven boundary layer flow of a non-Newtonian nanofluid. Processes 2021, 9, 702. [Google Scholar] [CrossRef]
  41. Varun Kumar, R.S.; Alhadhrami, A.; Punith Gowda, R.J.; Naveen Kumar, R.; Prasannakumara, B.C. Exploration of Arrhenius activation energy on hybrid nanofluid flow over a curved stretchable surface. ZAMM-J. Appl. Math. Mech./Z. Angew. Math. Mech. 2021, 101, e202100035. [Google Scholar] [CrossRef]
  42. Tamilzharasan, B.M.; Karthikeyan, S.; Kaabar, M.K.; Yavuz, M.; Özköse, F. Magneto Mixed Convection of Williamson Nanofluid Flow through a Double Stratified Porous Medium in Attendance of Activation Energy. Math. Comput. Appl. 2022, 27, 46. [Google Scholar] [CrossRef]
  43. Nadeem, S.; Hussain, S.T. Flow and heat transfer analysis of Williamson nanofluid. Appl. Nanosci. 2014, 4, 1005–1012. [Google Scholar] [CrossRef] [Green Version]
  44. Raju, C.S.K.; Sandeep, N.; Ali, M.E.; Nuhait, A.O. Heat and mass transfer in 3-D MHD Williamson–Casson fluids flow over a stretching surface with non-uniform heat source/sink. Therm. Sci. 2019, 23, 281–293. [Google Scholar] [CrossRef]
  45. Mustafa, M.; Khan, J.A.; Hayat, T.; Alsaedi, A. Buoyancy effects on the MHD nanofluid flow past a vertical surface with chemical reaction and activation energy. Int. J. Heat Mass Transf. 2017, 108, 1340–1346. [Google Scholar] [CrossRef]
  46. Ali, F.; Loganathan, K.; Eswaramoorthi, S.; Prabu, K.; Zaib, A.; Chaudhary, D.K. Heat transfer analysis on Carboxymethyl cellulose water-based cross hybrid nanofluid flow with entropy generation. J. Nanomater. 2022, 2022, 5252918. [Google Scholar] [CrossRef]
  47. Eswaramoorthi, S.; Loganathan, K.; Jain, R.; Gyeltshen, S. Darcy–Forchheimer 3D flow of glycerin-based carbon nanotubes on a Riga plate with nonlinear thermal radiation and Cattaneo–Christov heat flux. J. Nanomater. 2022, 2022, 5286921. [Google Scholar] [CrossRef]
Figure 1. Physical model of flow.
Figure 1. Physical model of flow.
Mca 27 00052 g001
Figure 2. The nanofluid velocity for various values R i (a), f w (b), λ (c) and M (d).
Figure 2. The nanofluid velocity for various values R i (a), f w (b), λ (c) and M (d).
Mca 27 00052 g002
Figure 3. The nanofluid temperature profile for various vales of R (a), H g (b), Γ T (c) and E c (d).
Figure 3. The nanofluid temperature profile for various vales of R (a), H g (b), Γ T (c) and E c (d).
Mca 27 00052 g003
Figure 4. The nanoparticle concentration for various values of σ * * (a), f w (b), E (c) and N t (d).
Figure 4. The nanoparticle concentration for various values of σ * * (a), f w (b), E (c) and N t (d).
Mca 27 00052 g004
Figure 5. The skin friction coefficient (SFC) for different values of M with f w = 0.6 (a), M with f w = 0.6 (b), λ with f w = 0.6 (c) and λ with f w = 0.6 (d) for Casson nanofluid (CF) and Williamson nanofluid (WF).
Figure 5. The skin friction coefficient (SFC) for different values of M with f w = 0.6 (a), M with f w = 0.6 (b), λ with f w = 0.6 (c) and λ with f w = 0.6 (d) for Casson nanofluid (CF) and Williamson nanofluid (WF).
Mca 27 00052 g005
Figure 6. The local Nusselt number (LNN) for different values of M with f w = 0.6 (a), M with f w = 0.6 (b), λ with f w = 0.6 (c) and λ with f w = 0.6 (d) for Casson nanofluid (CF) and Williamson nanofluid (WF).
Figure 6. The local Nusselt number (LNN) for different values of M with f w = 0.6 (a), M with f w = 0.6 (b), λ with f w = 0.6 (c) and λ with f w = 0.6 (d) for Casson nanofluid (CF) and Williamson nanofluid (WF).
Mca 27 00052 g006
Figure 7. The local Sherwood number (LSN) for different values of M with f w = 0.6 (a), M with f w = 0.6 (b), λ with f w = 0.6 (c) and λ with f w = 0.6 (d) for Casson nanofluid (CF) and Williamson nanofluid (WF).
Figure 7. The local Sherwood number (LSN) for different values of M with f w = 0.6 (a), M with f w = 0.6 (b), λ with f w = 0.6 (c) and λ with f w = 0.6 (d) for Casson nanofluid (CF) and Williamson nanofluid (WF).
Mca 27 00052 g007
Table 1. The skin friction coefficient, local Nusselt number and local Sherwood number for f w , λ , F r , M, R i and N r for both fluids.
Table 1. The skin friction coefficient, local Nusselt number and local Sherwood number for f w , λ , F r , M, R i and N r for both fluids.
fw λ FrMRiNr 1 2 C f R e N u / R e S h / R e
CassonWilliamsonCassonWilliamsonCassonWilliamson
−0.60.20.40.50.50.5−0.509450−0.3917090.5793210.583070−0.325537−0.327817
−0.3 −0.540008−0.4279640.6418090.642974−0.363858−0.364580
0.0 −0.571635−0.4682550.7086730.706679−0.405644−0.404386
0.3 −0.603205−0.5116660.7791390.773945−0.450585−0.447240
0.6 −0.633606−0.5566340.8522620.844421−0.498239−0.493078
0.40.00.40.50.50.5−0.594642−0.5080550.8064990.800459−0.468292−0.464370
0.4 −0.629919−0.5431880.8005620.794198−0.464437−0.460312
0.8 −0.657162−0.5718830.7962830.789413−0.461663−0.457217
1.2 −0.679033−0.5959490.7930560.785643−0.459573−0.454781
0.40.20.00.50.50.5−0.603953−0.5146060.8042720.798227−0.466845−0.462923
0.6 −0.617898−0.5320820.8028080.796584−0.465894−0.461858
1.2 −0.629722−0.5470220.8015600.795165−0.465084−0.460939
1.8 −0.639957−0.5600480.8004800.737920−0.464384−0.460133
0.40.20.40.00.50.5−0.560068−0.4757200.8127210.806599−0.472339−0.468357
0.5 −0.613519−0.5265770.8032680.797104−0.466193−0.462195
1.0 −0.650934−0.5651990.7972370.790500−0.462281−0.457920
1.5 −0.679033−0.5959490.7937860.785643−0.459573−0.454781
0.40.20.40.50.00.5−0.660336−0.5803490.7984240.788312−0.463051−0.456505
0.4 −0.622571−0.5368480.8024120.795557−0.465637−0.461193
0.7 −0.595805−0.5066270.8048360.799939−0.467212−0.464033
1.0 −0.570084−0.4779800.8068650.803636−0.468530−0.466432
0.40.20.40.50.5−1.0−0.606275−0.5227590.8083620.804310−0.469503−0.466870
−0.5−0.608591−0.5238960.8067520.802074−0.468456−0.465418
0.0−0.611002−0.5251650.8050550.799680−0.467353−0.463865
0.5−0.613519−0.5265770.8032680.797104−0.466193−0.462195
1.0−0.616151−0.5281550.8013810.794318−0.464969−0.460390
Table 2. The skin friction coefficient, local Nusselt number and local Sherwood number for R, H g , E c , Γ T , N b , N t and θ n for both fluids.
Table 2. The skin friction coefficient, local Nusselt number and local Sherwood number for R, H g , E c , Γ T , N b , N t and θ n for both fluids.
RHgEc Γ T NbNt θ n 1 2 C f R e N u / R e S h / R e
CassonWilliamsonCassonWilliamsonCassonWilliamson
0.0−0.50.40.10.50.51.2−0.629685−0.54437590.5523860.545982−0.552386−0.545982
0.4 −0.613519−0.5265770.8032680.797104−0.466193−0.462195
0.8 −0.600650−0.5135801.0233771.017386−0.411036−0.408267
1.2 −0.590235−0.5034701.2210391.215175−0.372078−0.370007
0.4−0.40.40.10.50.51.2−0.609973−0.5228670.7765460.770393−0.448914−0.444955
−0.2 −0.600663−0.5132120.7089010.702637−0.405788−0.401838
0.0 −0.586188−0.4983330.6083280.600989−0.343240−0.338747
0.2 −0.560224−0.4716860.4311090.417147−0.237329−0.229206
0.4 −0.548398−0.4436440.2487460.158094−0.133635−0.083938
0.4−0.50.00.10.50.51.2−0.616797−0.5301650.8289580.830200−0.482937−0.483750
0.5 −0.612698−0.5256860.7968340.788861−0.462020−0.456860
1.0 −0.608584−0.5212710.7645980.747832−0.441233−0.430502
1.5 −0.604457−0.5169180.7322660.707121−0.420586−0.404665
0.4−0.50.40.00.50.51.2−0.612511−0.5256200.7992750.793245−0.463602−0.459696
0.2 −0.614542−0.5275550.8073560.801051−0.468849−0.464754
0.4 −0.616637−0.5295760.8158340.809232−0.474367−0.470068
0.6 −0.618798−0.5316900.8247580.817834−0.480190−0.475670
0.4−0.50.40.10.10.51.2−0.624871−0.5342660.7950020.784238−2.304167−2.269369
0.5 −0.613519−0.5265770.8032680.797104−0.466193−0.462195
1.0 −0.612246−0.5258510.8041740.798416−0.233391−0.231523
1.5 −0.611828−0.5256180.8044700.798843−0.155658−0.154441
0.4−0.50.40.10.50.51.2−0.613519−0.5265770.8032680.797104−0.466193−0.46295
1.0 −0.614246−0.5261340.7903450.783132−0.915639−0.906320
1.5 −0.614786−0.5255090.7771960.768887−1.347998−1.331962
2.0 −0.615136−0.5247130.7638700.754421−1.763066−1.738850
0.4−0.50.40.10.50.51.0−0.617049−0.5303980.7480780.741525−0.487877−0.483603
1.2−0.613519−0.5265770.8032680.797104−0.466193−0.462195
1.4−0.609177−0.5219230.8716680.865948−0.440864−0.437192
1.6−0.603889−0.5163280.9560340.950777−0.412023−0.408740
Table 3. The skin friction coefficient, local Nusselt number and local Sherwood number for Γ c , σ * * , δ , n and E for both fluids.
Table 3. The skin friction coefficient, local Nusselt number and local Sherwood number for Γ c , σ * * , δ , n and E for both fluids.
Γ C σ * * δ nE 1 2 C f R e N u / R e S h / R e
CassonWilliamsonCassonWilliamsonCassonWilliamson
0.01.01.00.51.0−0.613044−0.5260390.8035080.797397−0.466349−0.462385
0.1 −0.613519−0.5265770.8032680.797104−0.466193−0.462195
0.2 −0.613985−0.5271120.8030240.796807−0.466035−0.462003
0.3 −0.614443−0.5281720.8027750.796199−0.465873−0.461608
0.10.01.00.51.0−0.621718−0.5373160.8020430.793644−0.465398−0.459954
0.4 −0.616422−0.5298960.8033440.796616−0.466243−0.461879
0.8 −0.614204−0.5273240.8033930.797110−0.466274−0.462199
1.2 −0.612990−0.5260160.8030970.797019−0.466082−0.462140
0.11.00.00.51.0−0.613684−0.5267840.8037130.797484−0.466482−0.462441
0.3 −0.613634−0.5267190.8035720.797365−0.466391−0.462364
0.6 −0.613585−0.5266570.8034380.797250−0.466303−0.462290
1.0 −0.613519−0.5265770.8032680.797104−0.466193−0.462195
0.11.01.00.01.0−0.613558−0.5266290.8034680.797283−0.466323−0.462311
0.5 −0.613519−0.5265770.8032680.797104−0.466193−0.462195
1.0 −0.613492−0.5265400.8030290.796883−0.466038−0.462052
1.5 −0.613478−0.5265200.8027450.796615−0.465854−0.461878
0.11.01.00.50.0−0.611041−0.5240710.8022640.796452−0.465541−0.461773
1.0−0.613519−0.5265770.8032680.797104−0.466193−0.462195
2.0−0.616499−0.5299620.8032400.796521−0.466175−0.461817
3.0−0.618951−0.5331180.8028020.795459−0.465891−0.461129
Table 4. Comparison of local Nusselt number when W e = λ = F r = R = Γ T = H g = Γ C = 0 and M = N r = 0 : 5 , S c = 5 a n d δ = 1 with Mustafa et al. [45].
Table 4. Comparison of local Nusselt number when W e = λ = F r = R = Γ T = H g = Γ C = 0 and M = N r = 0 : 5 , S c = 5 a n d δ = 1 with Mustafa et al. [45].
P r N t E σ * * n R i N u / R e
Ref. [45]Present
2.00.51.01.00.50.50.7066050.706604
4.0 0.9359520.935955
7.0 1.1327871.132788
10.0 1.2574761.257482
5.00.11.01.00.50.51.4262671.426269
0.5 1.0139391.013938
0.7 0.8469430.846928
1.0 0.6499400.649939
5.00.50.01.00.50.50.9412010.941209
1.0 1.0139391.013943
2.0 1.0645511.064563
4.0 1.1145491.114191
5.00.51.00.00.50.51.1453041.145301
1.0 1.0139391.013938
2.0 0.9262820.926281
5.0 0.7986710.798669
5.00.51.02.0−1.00.51.0308051.030804
−0.5 0.9994700.999468
0.0 0.9642860.964285
1.0 0.8868300.886830
10.00.51.02.00.50.01.0322811.032280
0.51.0567041.056706
3.01.1545391.154538
5.01.2159371.215938
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Eswaramoorthi, S.; Thamaraiselvi, S.; Loganathan, K. Exploration of Darcy–Forchheimer Flows of Non-Newtonian Casson and Williamson Conveying Tiny Particles Experiencing Binary Chemical Reaction and Thermal Radiation: Comparative Analysis. Math. Comput. Appl. 2022, 27, 52. https://doi.org/10.3390/mca27030052

AMA Style

Eswaramoorthi S, Thamaraiselvi S, Loganathan K. Exploration of Darcy–Forchheimer Flows of Non-Newtonian Casson and Williamson Conveying Tiny Particles Experiencing Binary Chemical Reaction and Thermal Radiation: Comparative Analysis. Mathematical and Computational Applications. 2022; 27(3):52. https://doi.org/10.3390/mca27030052

Chicago/Turabian Style

Eswaramoorthi, Sheniyappan, S. Thamaraiselvi, and Karuppusamy Loganathan. 2022. "Exploration of Darcy–Forchheimer Flows of Non-Newtonian Casson and Williamson Conveying Tiny Particles Experiencing Binary Chemical Reaction and Thermal Radiation: Comparative Analysis" Mathematical and Computational Applications 27, no. 3: 52. https://doi.org/10.3390/mca27030052

APA Style

Eswaramoorthi, S., Thamaraiselvi, S., & Loganathan, K. (2022). Exploration of Darcy–Forchheimer Flows of Non-Newtonian Casson and Williamson Conveying Tiny Particles Experiencing Binary Chemical Reaction and Thermal Radiation: Comparative Analysis. Mathematical and Computational Applications, 27(3), 52. https://doi.org/10.3390/mca27030052

Article Metrics

Back to TopTop