Entropy generation in MHD Casson fluid flow with variable heat conductance and thermal conductivity over non-linear bi-directional stretching surface.

This consideration highlights the belongings of momentum, entropy generation, species and thermal dissemination on boundary layer flow (BLF) of Casson liquid over a linearly elongating surface considering radiation and Joule heating effects significant. Transportation of thermal and species are offered by using the temperature-dependent models of thermal conductivity and mass diffusion coefficient. Arising problem appear in the form of nonlinear partial differential equations (NPDEs) against the conservation laws of mass, momentum, thermal and species transportation. Appropriate renovation transfigures the demonstrated problem into ordinary differential equations. Numerical solutions of renovated boundary layer ordinary differential equations (ODEs) are attained by a proficient and reliable technique namely optimal homotopy analysis method (OHAM). A graphical and tabular interpretation is given for convergence of analytic solutions through error table and flow behavior of convoluted physical parameters on calculated solutions are presented and explicated in this examination. Reliability and effectiveness of the anticipated algorithm is established by comparing the results of present contemplation as a limiting case of available work, and it is found to be in excellent settlement. Decline in fluid velocity and enhancement in thermal and species transportation is recorded against the fluctuating values of Hartman number. Also reverse comportment of Prandtl number and radiation parameter is portrayed. Moreover, it is conveyed that supplementing values of the magnetic parameter condenses the fluid velocity and upsurges the thermal and concentration distributions. Negative impact of elevating Joule heating phenomenon is noted on the molecular stability of the system via Brinkman number [Formula: see text] Furthermore, the system's stability at a molecular level is controlled by diminishing values of radiation [Formula: see text] temperature difference [Formula: see text] concentration difference [Formula: see text] diffusion parameters [Formula: see text] and Brinkman number [Formula: see text].

Introduction

Stretched flows have much applications in different industrial phenomenon. Time dependent boundary layer flow streaming past a shrinking surface with influence of transverse magnetic field was premeditated by Merkin and Kumaran[1]. They showed that the fluid seemed to be effected by the range of magnetic parameter and dimensionless time parameter, such that for magnetic parameter less than 1, steady state was achieved but for times very large but for magnetic parameter taken as 1, boundary layer flow was achieved for all values of time and the wideness of boundary layer was observed to increase with time. Shehzad et al.[2] premeditated the contribution of temperature dependent thermal conductivity on Oldroyd-B model past over a stretched surface. Thermal distribution with heat generation on viscoelastic model produced due to the stretching of the sheet was reported by Khan et al.[3]. They have offered the demeanor of several embryonic parameters on velocity field and temperature profile through graphs and tables. They presented the comparative study too in order to discourse the authenticity of the archived solutions. Moreover, they recorded that positive mounting values of heat generation factor serves to enhance the temperature distribution. Khan et al.[4] dissected the heat and mass transference in a viscous fluid past over a nonlinear lengthening sheet. For comparative study, they offered the numerical and analytical solutions both in their contribution. They found excellent settlement in both the computed solutions. Natural convection flow of copper–water based nanofluid through oddly shaped geometry was studied by Parvin and Chamkha[5] using penalty finite element method with Galerkin weighted residual technique. The study highlighted the importance of Rayleigh number as with larger value of Rayleigh number the contribution of heat transfers in entropy generation rose while the contribution of viscous dissipation dropped. Hence, with the right range of Rayleigh number, it was believed that stability could be maximized for this particular physical system. Time dependent MHD Casson fluid free flow over an oscillating vertical plate covered with porous media was analyzed by Khalid et al.[6] using Laplace transform. The results shined light on reduced velocity and risen skin friction on the surface with elevating magnetic parameter. Rauf et al.[7] investigated boundary layer flow (BLF) of Casson nanofluid (CN) on a sheet elongating in two directions exposed to transverse magnetic field and thermal radiation with mixed convection conditions using Range-Kutta-Fehlberg (RKF45) technique. It was noticed that the system improved thermally with higher radiation. Zaib et al.[8] concentrated on Casson fluid with effects of viscous dissipation flowing on an exponentially shrinking sheet. The mathematical system was solved numerically using shooting method and dual solutions were obtained for velocity and temperature. The analysis showed a decline in system’s temperature with rising Prandtl number whereas system heated up for enhanced effects of viscous dissipation. Raju et al.[9] compared Casson fluid with Newtonian fluid over an exponentially broadening surface under the effects of thermal radiation, viscous dissipation and magnetism. Using MATLAB bvp4c package, the results showed an improvement in the system thermally for higher viscous dissipation effects and that heat transfer rates were better for Casson fluid compared to Newtonian fluid. An almost similar study done by Soluchana et al[10] in which three dimensional Casson nanofluid is compared to Newtonian fluid using RK shooting technique. The study incorporated MHD effects and the fluid is passed over a stretching plane. Both mass and heat transfer rates boosted for larger stretching ratio parameter and here again Casson fluid had better heat, and mass, transfer rates compared to Newtonian fluid. Kumaran and Sandeep[11] worked on the comparison of MHD Casson and Williamson fluid streaming on top of an upper paraboloid of revolution taking into account the thermophoresis and Brownian motion impacts. Using the RK method with shooting technique, Casson fluid turned out to be better than Williamson fluid in terms of heat and mass transfer. MHD Casson fluid streaming past a wedge with influence of binary chemical reactions and activation energy was premeditated by Zaib et al.[12]. With the use of modified Arrhenius function to represent activation energy the model was developed and solved using keller box method. The results highlighted decaying temperature and concentration boundary layers and enhanced fluid stream velocity for elevating Casson fluid parameter. Irfan et al.[13] analyzed the impact of stretching rate ratio parameter on the three dimensional forced convection Carreau nanofluid flow over a surface stretching in two directions keeping variable thermal conduction and heat generation/absorption as significant effects. After solving the equations using MATLAB bvp4c package, it was seen that elevated stretching rate ratio caused a drop in skin friction coefficient in both and directions for shear thickening as well as thinning cases of the Carreau fluid. MHD mixed convection Casson fluid flow with effects of double stratification and heat cohort/immersion was considered by Rehman et al.[14]. The study proceeded with the use of RK method with shooting technique and the results showed a deescalating velocity in regard with rising Casson fluid parameter. Merkin et al.[15] worked on stagnation point flow over an exponentially elongating/dwindling cylinder using shooting method. It was seen that the solution came out to be unique for the stretching case, whereas dual solutions appeared for the shrinking case further developing the study based on the critical values defined through patterns observed. Stagnation point nanofluid flow streaming on top of a stretching surface was investigated by Jalilpour et al.[16] using the RK-4 numerical technique considering the impact of thermal radiation. Higher thermophoresis brought about a negative change in the Nusselt number while a rise in the mass transfer rate, while radiation boosted the Sherwood number. Sohail et al.[17] studied the second grade fluid on a stretching surface using HAM. Considerations of generalized Fourier’s and Fick’s laws were taken to incorporate the Cattaneo-Christov heat and mass flux in the physical system. It was seen that compared to non-Fick’s and non-Fourier’s laws, temperature and concentrations was higher for classical Fick’s and Fourier’s laws. Zaib et al.[18] emphasized on stagnation point flow of Casson nanofluid on top of a plate immersed in Darcy-Brinkman porous medium. Non-linear radiation and activation energy was incorporated using modified Arrhenius function and impact of binary chemical reaction was also considered. With the use of shooting method nanoparticles concentration became dense for rising influence of activation energy while opposite happened for larger reaction rate. Abrar et al.[19] studied the entropy generation in a water based titanium dioxide nanofluid that was transported through cilia. Viscous dissipation, radiation and MHD effects were considered significant. Similarly, Rashidi et al.[20] considered MHD third grade fluid streaming past a sheet stretching in a linear manner and worked on the system using optimal homotopy analysis method (OHAM). Both the studies showed magnetic field and total entropy generation were in direct relation. Raju and Sandeep[21] worked on Casson nanofluid considered to be flowing on a rotating cone present in a rotating frame using the Range-Kutta (RK) shooting technique. Their results highlighted escalating drag force on the surface and heat transfer rates for higher volume fraction of nanoparticles present in the fluid. Non-Newtonian fluids have been a strong focus of researchers due to their vast properties compared to simple Newtonian fluids. These fluids have a non-linear relation of shear stress to shear strain. Examples of non-Newtonian fluids are numerous, some including: ketchup, honey, blood, molten polymers, clays, paints, etc. Some materialistic characteristics are reported in[22]. Casson fluid is a non-Newtonian fluid that best describes fluids like blood, honey, jelly, etc., where fluids behave like elastic solids and have molecular chains connecting the particles within. This fluid model is better at describing the rheological properties of such fluids that cannot be described using viscoelastic fluid models. Casson fluid is a shear thinning fluid liquid with three assumptions: (i) at zero shear rate the viscosity is infinite, (ii) no flow exists below the yield stress and (iii) at infinite shear the viscosity is neglible as explicated by Reddy et al.[23]. Unsteady Casson fluid streaming on top of a enlarging sheet was investigated by Khan et al.[24] using HAM. They noticed that the thickness of the boundary layer for unsteady case was very less compared to the steady case. Sheikholeslami et al.[25] worked on an unsteady MHD radiative nanofluid flowing over an oscillating vertical plate with effects of heat generation and absorption taken into account. Flow stream velocity was seen to decline with heightening magnetic encouragement. Zaib et al.[26] emphasized on the stability of stagnation point flow of Williamson nanofluid flowing on top of a moving plate through entropy generation. Impacts of activation energy and binary chemical reactions were taken into consideration, where activation energy was incorporated using modified Arrhenius function. The mathematical system was analyzed numerically using shooting method and the results portrayed heightened concentration boundary layer for large activation energy whereas the opposite was noted for increased reaction rates. Moreover, entropy generation for Williamson fluid was greater in comparison to Newtonian fluid. Usman et al.[27] premeditated the impression of non-linear thermal radiation and thermal conductivity depending on time on the rotating flow of hybrid copper–aluminum oxide nanofluid with water as base fluid. The fluid flowed on three dimensional stretching surface. Using least square method (LSM) to solve the mathematical system it was found that the hybrid nanofluid had lower thermal advantages compared to individual copper and aluminum oxide nanofluids. Steady axisymmetric rotational stagnation point flow was analyzed by Lok et al.[28] using the bvp4c MATLAB package over a rotating permeable shrinking/stretching disk. The study concluded the existence of three solutions depending on the stretching rate, the rotation rate and the fluid withdrawal/injection. Reddy et al.[29] compared the Casson and Maxwell MHD fluids over a stretched surface using the numerical scheme RK Fehlberg. The result highlighted higher mass transfer rates but reduced heat transfer rates for Maxwell fluid in comparison to Casson fluid. Kumar et al.[30] considered MHD Casson and Maxwell fluid with effects of chemical reaction and heat generation/absorption. With the use of RK method along with shooting technique comparison between both the fluids was studied. The results highlighted reduced heat and mass transfer rates for heightened magnetic influence for both Casson as well as Maxwell fluid. Mehryan et al.[31] considered using periodic magnetic field on the natural convection flow of -water based nanofluid within a square enclosed system to study the entropy generation. With the help of Galerkin finite element method, for values of Hartmann number less than 10, the presence of nanoparticles seemed to boost entropy generation, while for values greater than 10, nanoparticles’ influence positively impacted the stability of the molecules within the system. Hamid et al.[32] studied the characteristics of variable thermal conductivity on mixed convective radiative flow of rotating nanofluid past over an elongating sheet which is stretched horizontally. They presented the finite element analysis of the modelled physical system. They have mentioned that brick shape nanoparticles can work as a coolant for the thermal system. Legendre wavelets approach was presented by Soomro et al.[33] to investigate the influence of thermal and momentum slips on mixed convective Williamson model. They concluded that augmenting values of Prandtl number increases the heat transfer rate and mass transport rate is assisted due to increment in Lewis number. Shah et al.[34] emphasized on the Casson micropolar ferrofluid flowing on a linearly stretching sheet taking transverse magnetohydrodynamic (MHD) and thermal radiation effects into contemplation. Homotopy analysis method (HAM) was used to solve the mathematical system. The results concluded a significant rise in the surface drag force with large magnetic influence and higher heat transfer rates with upsurge in Prandtl values. Hamid et al.[35] considered a stretching sheet covered with MHD Casson fluid influenced by linearly acting thermal radiation. Dual solutions were calculated and was seen that both upper and lower profiles of the temperature were an increasing function of radiation. MHD nanofluid flowing through an enclosed area filled with minor pores was investigated by Sheikholeslami[36]. A new numerical technique called the control volume finite element process was used to solve the mathematical system and the results shows a downfall in the heat transfer rates with rising magnetic influence. Ullah et al.[37] studied the MHD squeezing viscous fluid flowing between two parallel plates flooded with minor pores. The focus was on the comparison of numerous schemes like optimal homotopy analysis method (OHAM), homotopy perturbation (HP), differential transform (DT), Daftardar Jafari (DJ) and Adomian decomposition method (ADM). In view of the results, they recorded that OHAM and HPM were more accurate. Heat transfer rates were studied for natural convection flow of Casson fluid flowing through a partially heated trapezoidal cavity using Galerkin finite element method (GFEM) by Hamid et al.[38]. Through results it was seen by elevating Casson fluid parameter enhanced heat transfer rates in the middle of the cavity. Sheikholeslami et al.[39] analyzed the Nusselt number for alumina nanofluid streaming between two parallel plates, top plate being permeable, using neural networks attaining the numerical data by RK method. Nusselt number was positively influenced by the nanofluid concentration. Hosseini et al.[40] premeditated the belongings of using aluminum oxide water based nanofluid paired with symmetric heating and MHD on the entropy generation in the porous horizontal channel. Irreversibility of fluid heat transfer under the influence of MHD were less than when MHD impact was considered negligible. Hamid et al.[41] used the finite difference approach to captures the features of slip flow of Prandtl model with mass and heat transportation under thermos-diffusion effects. They commented that unsteadiness parameter upsurges the velocity, mass and temperature distribution profiles. MHD flow in a cavity with radiative heat transportation and mass transfer was explored by Usman et al.[42]. They presented that radiation parameter serves as a coolant for the thermal system which is much needed because higher warming effects can cause a failure for the thermal system. Numerically the flow of Carreau fluid with heat and mass transportation induced due to the stretching of the sheet was surveyed by Ali et al.[43]. They have shown that Lorentz force cause to boost the temperature field and lessen the velocity profile. Ali et al.[44] discussed the features of chemically reactive species with thermal distribution in a Jeffery liquid past over a stretching cylinder. Generalized model for heat flux is presented to model the energy equation. Contribution of heat and mass in a yield exhibiting fluid model in a cylinder was explored by Ali et al.[45]. Revised definition for heat flux model is used with heat generation. They have shown that curvature parameter and Lorentz force shows opposite impact of flow velocity. Mixed convective flow of Newtonian liquid in a stretched cylinder with heat and mass transfer was studied by Ali et al.[46]. Significant support of joule heating and heat cohort are encompassed in their study. Ali et al.[47] discussed the flow past over a rotating disk with heat transportation. They presented the involvement of homogeneous-heterogeneous reaction in their worthy investigation.

Comprehensive cited literature survey fails to report the irreversibility analysis of Casson fluid obeying temperature dependent conductance and thermal conductivity. Over the past years, MHD Casson fluid under the inspiration of thermal radiation, temperature dependent conductivity and temperature dependent mass diffusivity over a linearly elongating surface has not been considered. Moreover, investigation of entropy generation for such scenario is new as well which plays a vital role in numerous industrial disciplines. This study aims to focus on the considered phenomenon with the help of optimal homotopy analysis method (OHAM)[48-54] and obtain graphical representations for the velocity, temperature, concentration changes and molecular stability of the system. Comparative study is also presented in order to verify the effectiveness of the anticipated computational tool. This exploration is organized as follows: literature survey is presented in “Introduction”, physical happening is listed in “Mathematical drafting (MD) for the physical happening via boundary layer theory (BLT)”, irreversibility phenomenon is listed in “Entropy generation computation (EGC)”, used procedure is mentioned in “Solution via OHAM”, physical interpretation is organized in “Results and discussion” and the last section is “Conclusion”.

Mathematical drafting (MD) for the physical happening via boundary layer theory (BLT)

Magnetohydrodynamic boundary layer flow (MBLF) of yield stress demonstrating Casson liquid[6–12, 14, 18, 21, 23, 24, 30, 32, 35, 38, 45, 51] over a linear stretchable surface with the transportation of heat and mass is conveyed in current work. Schematic diagram for the deliberated physical happening is shown through Fig. 1.

Figure 1

Geometry of the modeled physical problem with Coordinate system.

Flow is manufactured due to the enlarging of the bi-directional strained surface[4, 49–51] under constant pressure. In current assessment, Cartesian coordinates system is used to model the conservation laws. Magnetic field is tarnished normal to stretchable surface. Attractive Reynolds number is anticipated to be less, so that the predominant magnetic field is deserted associated to pragmatic magnetic field. Heat and mass repositioning contrivances are recognized by exhausting the temperature dependent thermal conductivity[27, 32] and diffusion coefficient. Moreover, impact of radiation[32, 43] and Joule heating[46] is taken in this exploration. With such implication predominant boundary layer equations that entitles the flow state are enunciated as

Following boundary conditions appears against the modeled situation recorded in Eqs. (2–5) has been developed via no slip perception

Entropy generation computation (EGC)

Entropy phenomenon[5, 12, 19, 20, 26, 31, 40, 49, 52, 53] for Casson fluid past over a bidirectional strained surface with heat and mass conveyance having variable thermal conductivity and temperature dependent mass diffusion coefficient is communicated as

Equation (7) comprehends concentration, Joule heating and thermal irreversibility and its dimensionless relation is articulated as

Essential transmutations are

Exploiting above declared renovations Eqs. (2–3), (10, 11) and (14) takes the form

Renovated boundary conditions are

All dimensionalized variables and parameters are shown in Table 1.

Table 1

Nomenclature.

u1∗,u2∗ and u3∗	Velocity components	x_1-,x_2- and x_3-	Space coordinates
β=μ_δ2π_δp_δ	FLUID parameter	ν_N	Kinematic viscosity
σ_N	ELECTRICAL conductivity	B_N	Magnetic field strength
ρ_N	Fluid density	T1∗ and C1∗	Fluid temperature and concentration
c_p	Specific heat	K_m	Mean absorption coefficient
KT1∗=1+ε_aθ	Temperature dependent thermal conductivity	DT1∗=1+ε_bθ	Temperature dependent diffusion coefficient
Ha^2=σ_NB_N^2ρ_Na1∗	Magnetic parameter	Pr=ν_Nδ^∗	Prandtl number
Sc=D_βδ^∗	Schmidt number	ξ=a1∗ν_Nx_3	Dimensionless independent variable
α=b1∗a1∗	Ratio parameter	Ec=Uw2μ_FT_w-T_∞	Eckert number
Br	Brinkman number	∈_1	Temperature difference parameter
R=43kK_mσ_NT∞3	Radiation parameter	∈_2 and ∈_3	Concentration difference and diffusion parameter

Solution via OHAM

Due to a couple nonlinear system for controlling differential equations system it is recommended to optimize homotopy scheme for calculating the solutions. The proposed scheme is parameter free or large and does not have to be set aside. This method has no stability issues as seen in numerical schemes. It can be applied to infinite domains and can address linear, non-linear, homogeneous and inhuman problems in a similar way. This scheme requires selecting a linear operator and initial guess. An initial guess is chosen in a way that satisfies the given boundary conditions. Correct selection of initial guesses guarantees convergence of the homotopic procedure.

Initial guesses with operators are listed below these operators imitate following characteristic where are indiscriminate constants.

Consequences and argument

This subclass contains graphical results for the convergence of homotopic solutions, dimensionless velocity, temperature and concentration. For this exam, we have set some muddled parameters and change one parameter and then figured out its bearing.

Convergence analysis

This section contains tabular results in the form of error analysis for the convergence of expected homotopic solutions for dimensionless velocities, temperature and concentration solutions. For this scrutiny we have set some thin parameters and change one parameter and then pressed its behavior. The fixed values are and Convergence of desired solution is depicted in Table 2 which confirms that by increasing order of approximations, error diminishes which collaterals the convergence of suggested scheme. This table presents the direct relation between the errors and approximation order. Higher order of approximation reduces the errors in obtained solutions.

Table 2

Convergence examination through error reduction via OHAM.

s^∗	E⃛s^∗fx_1	E⃛s^∗g	E⃛s^∗θ	E⃛s^∗φ	CPU time (s)
2	293593×10^-6	2.93593×10^-6	0.000244616	0.0000214515	2.5467
4	1.85543×10^-8	1.85543×10^-8	0.0000187684	1.03026×10^-6	10.0456
6	3.31176 ×10^-10	3.31176 ×10^-10	3.41302 ×10^-6	2.51364 ×10^-8	15.6537
8	8.81889 ×10^-12	8.81889 ×10^-12	5.40688 ×10^-7	1.20006 ×10^-9	28.1432
10	2.86855 ×10^-13	2.86855 ×10^-13	1.14876 ×10^-7	5.96324 ×10^-11	48.9438
12	1.03859 ×10^-14	1.03859 ×10^-14	2.9143 ×10^-8	2.29532 ×10^-12	126.2504
16	1.56777×10^-16	1.56777 ×10^-16	1.72687 ×10^-9	6.03542 ×10^-13	354.5284
20	2.53601×10^-20	2.53601 ×10^-20	1.17993 ×10^-10	2.24072 ×10^-17	528.2963

Results and discussion

This section contains the graphical results for the fluid velocity, temperature and concentration fields against the various influential variables which appears in the transformed system of ordinary differential equations. The impacts of the temperature dependent thermal conductivity, temperature dependent diffusion coefficient the magnetic parameter the Prandtl number the Schmidt number the radiation parameter the Brinkman number temperature difference parameter concentration difference parameter and the diffusion parameter are noted on momentum, thermal, mass boundary layers and entropy generation through graphical representations.

Figures 2 and 3 highlight slowed stream flow in and directions of the stretching surface as Casson fluid parameter is enhanced. This is due to the drop in yield stress at higher values of which causes the fluid to behave more like a Newtonian fluid, hence the velocities are seen to fall. Figures 4 and 5 shine light on the elevated thermal and mass profiles due to high Rise in causes elasticity stress parameter to intensify which in turn effects the temperature and concentration of the system positively. Figures 6, 7, 8 and 9 show the effects of increasing magnetic parameter on the velocity, thermal and mass behaviors. contributing directly to the resistive Lorentz force causes the velocities in both the directions to decay, as seen in Figs. 6 and 7, and due to the fall of fluid speed, the system is noted to heat up and become denser as represented in Figs. 8 and 9, respectively. Decline in velocity happens due to resistance during flow in associated region. Moreover, molecular vibration heightens the thermal and species transmission. Figure 10 focuses on the positive influence of radiation on the temperature of the system. This happens because higher results in the thicker thermal boundary layer which results to enhance temperature. Moreover, involvement of radiation improves conduction phenomenon which upturns the temperature field. Figure 11 shows a rise in temperature profile for incrementing values of This is because rising results in higher thermal conductivity and hence more heat is transferred from plate to liquid. In Fig. 12 the system is seen to cool down for raising values of Prandtl number because of its inverse relation to thermal diffusivity. Against larger values of Prandtl number fluid temperature and associated boundary layer thickness reduces. Inverse bearing of thermal diffusion causes to lessen the thermal distribution. Similarly, due to the inverse relation of Schmidt number to the mass diffusivity, the concentration boundary layer becomes thinner for incrementing values of as noted in Fig. 13. Since Schmidt number has the inverse relation with mass diffusion. Incrimination in Schmidt number is due to lessen in mass diffusion which results to slow down the concentration field. Figure 14 displays the positive effect of on mass boundary layer as with increasing the mass diffusivity of the fluid improves and hence mass boundary layer thickens. Dimensionless entropy generation is studied in Figs. 15, 16, 17, 18 and 19 against various parameters. In Figs. 15 and 16, is noted to rise for increasing radiation parameter and Brinkman number Since and are contributions of Joule heating and thermal radiation, respectively, hence their rise causes more instability in the system. Rise in temperature difference parameter , concentration difference parameter and the diffusion parameter heighten the phenomenon of concentration irreversibility, which too is a direct contributing factor in the process of entropy generation, therefore is noted to grow in Figs. 17, 18 and 19. Table 3 is constructed to validate the applied scheme in comparison with already published work for and taking A admirable arrangement is comprehended for the technique used in this exploration compared with the earlier published studies which use shooting method and MATLAB built in package bvp5c. Moreover, Table 4 is prepared to record the behavior of dimensionless stresses against the mounting values of ratio parameter by setting Results are listed by comparing the obtained values of stresses with three different approaches. The recorded results shown an excellent settlement with those of reported in ref.[2] and ref.[3] in the absence of magnetic effect which retards the flow and control the turbulence. Also these results show an enhancement in dimensionless stresses against the fluctuating values of ratio parameter.

Figure 2

Inspiration of on .

Figure 3

Influence of on .

Figure 4

Comportment of on .

Figure 5

Bearing of on .

Figure 6

Effect of on .

Figure 7

Inspiration of on .

Figure 8

Encouragement of on .

Figure 9

Bearing of on .

Figure 10

Comportment of on .

Figure 11

Involvement of against .

Figure 12

Influence of on .

Figure 13

Inspiration of on .

Figure 14

Influence of on .

Figure 15

Influence of on .

Figure 16

Influence of on .

Figure 17

Influence of on .

Figure 18

Influence of on .

Figure 19

Influence of on .

Table 3

Analysis of and for

n	α	-1+1βf^″0)[4]		-1+1βg^″(0)[4]		Present results using OHAM
		Shooting	Bvp5c	Shooting	Bvp5c	-1+1βf^″0)	-1+1βg^″(0)
1.0	0	1.0	1.0	0.0	0.0	1.0	0.0
-	0.5	1.224745	1.224742	0.612372	0.612371	1.224744	0.612371
-	1.0	1.414214	1.414214	1.414214	1.414214	1.414214	1.414214

Table 4

Comparative studies for and when

α	HPM[2]		Exact[2]		[3]	[3]	Present	Present
	− f^″(0)	g^″(0)	f^″(0)	g^″(0)	f^″(0)	g^″(0)	f^″(0)	g^″(0)
0.0	1.0	0.0	1.0	0.0	1.0	0.0	1.0	0.0
0.1	1.02025	0.06684	1.020259	0.066847	1.02026	0.06685	1.020257	0.066845
0.2	1.03949	0.14873	1.039495	0.148736	1.03949	0.14874	1.039489	0.148738
0.3	1.05795	0.24335	1.057940	0.243359	1.05795	0.24336	1.057938	0.243350
0.4	1.07578	0.34920	1.075788	0.349208	1.07578	0.34921	1.075785	0.349202
0.5	1.09309	0.46520	1.093095	0.465204	1.09309	0.46521	1.093073	0.465201
0.6	1.10994	0.59052	1.109946	0.590528	1.10994	0.59053	1.109942	0.590523
0.7	1.12639	0.72453	1.126397	0.724531	1.12639	0.72453	1.126395	0.724530
0.8	1.14248	0.86668	1.142488	0.866682	1.14249	0.86668	1.142487	0.866680
0.9	1.15825	1.01653	1.158253	1.016538	1.15826	1.016538	1.158250	1.016535
1.0	1.17372	1.17372	1.173720	1.173720	1.173720	1.17372	1.173720	1.173720

Conclusions

Flow of Casson model under boundary layer has been scrutinized over a linearly stretching surface immersed in porous medium. Temperature dependent conductivity and diffusivity has been taken into consideration. Using OHAM analytical solutions for momentum, thermal and mass boundary layers have been obtained and examined against various parameters. Error analysis reveals the effectiveness of utilized algorithm. Proposed procedure has the ability to handle the other nonlinear problems occurring in other discipline of science. Main verdicts of this study are presented as follows:

Fluid velocity in the and directions reduces when values of Casson fluid and magnetic parameters are enhanced.

The system heats up with elevating Casson fluid parameter magnetic parameter small scalar parameter and radiation parameters while opposite impact is noticed for the mounting values of Prandtl number

Thicker boundary layer is recorded against and

Temperature and concentartion fields are enhanced for the mounting values of Hartman number, whereas, decline in the velocity field is observed.

Fluid becomes dense for incrementing Casson fluid magnetic and small scalar parameter parameters while thinned out for Schmidt number

The system’s stability at a molecular level is controlled by diminishing values of radiation temperature difference concentration difference diffusion parameters and Brinkman number

Thermal transmission is the growing function of magnetic parameter small scalar parameter and radiation