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 Br.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left( {Br} \right).$$\end{document} Furthermore, the system’s stability at a molecular level is controlled by diminishing values of radiation R,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left( R \right),$$\end{document} temperature difference ∈1,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left( { \in_{1} } \right),$$\end{document} concentration difference ∈2,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left( { \in_{2} } \right),$$\end{document} diffusion parameters ∈3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left( { \in_{3} } \right)$$\end{document} and Brinkman number Br.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left( {Br} \right).$$\end{document}

www.nature.com/scientificreports/ 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 x and y 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  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  www.nature.com/scientificreports/ 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 Fe 3 O 4 -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  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][49][50][51][52][53][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".
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 B = [0, 0, B N ] 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 www.nature.com/scientificreports/

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  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.

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 K T * 1 = 1 + ε a θ the temperature dependent thermal conductivity, temperature dependent diffusion coefficient D T * 1 = 1 + ε b θ, the magnetic parameter (Ha), the Prandtl number (Pr), the Schmidt number (Sc), the radiation parameter (R), the Brinkman number (Br), temperature difference parameter (∈ 1 ), concentration  Velocity components x 1 −,x 2 − and x 3 − Space coordinates Fluid temperature and concentration c p Specific heat (K m ) Mean absorption coefficient Brinkman number (∈ 1 ) Temperature difference parameter Radiation parameter (∈ 2 ) and (∈ 3 ) Concentration difference and diffusion parameter www.nature.com/scientificreports/ difference parameter (∈ 2 ) and the diffusion parameter (∈ 3 ) are noted on momentum, thermal, mass boundary layers and entropy generation through graphical representations. Figures 2 and 3 highlight slowed stream flow in x 1 and x 2 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 (Ha) on the velocity, thermal and mass behaviors. (Ha) 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 (R) on the temperature of the system. This happens because higher (R) 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 (ε a ). This is because rising (ε a ) 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 (Pr) 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 (Sc) to the mass diffusivity, the concentration boundary layer becomes thinner for incrementing values of (Sc) 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 (ε b ) on mass boundary layer as with increasing (ε b ) the mass diffusivity of the fluid improves and hence mass boundary layer thickens. Dimensionless     . Since (Br) and (R) are contributions of Joule heating and thermal radiation, respectively, hence their rise causes more instability in the system. Rise in temperature difference parameter (∈ 1 ) , concentration difference parameter (∈ 2 ) and the diffusion parameter (∈ 3 ) heighten the phenomenon of concentration irreversibility, which too is a direct contributing    Table 3 is constructed to validate the applied scheme in comparison with already published work for − 1 + 1 β f ′′ (0) and − 1 + 1 β g ′′ (0) taking β = ∞, Ha = 0. 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

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 x 1 and x 2 directions reduces when values of Casson fluid (β) and magnetic (Ha) parameters are enhanced. Figure 19. Influence of (∈ 3 ) on Eg(ξ ). Table 3. Analysis of − 1 + 1 β f ′′ 0) and − 1 + 1 β g ′′ (0) for β = ∞, Ha = 0.  Table 4. Comparative studies for − 1 + 1 β f ′′ 0) and − 1 + 1 β g ′′ (0) when β = ∞, Pr = 0.7, Sc = 0.5, Ha = 0.  www.nature.com/scientificreports/ • The system heats up with elevating Casson fluid parameter (β), magnetic parameter (Ha), small scalar parameter (ε a ) and radiation (R) parameters while opposite impact is noticed for the mounting values of Prandtl number (Pr). • Thicker boundary layer is recorded against (ε a ) and (R).
• 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 (Ha) and small scalar parameter (ε b ) parameters while thinned out for Schmidt number (Sc). • The system's stability at a molecular level is controlled by diminishing values of radiation (R), temperature difference (∈ 1 ), concentration difference (∈ 2 ), diffusion parameters (∈ 3 ) and Brinkman number (Br). • Thermal transmission is the growing function of (β), magnetic parameter (Ha), small scalar parameter (ε a ) and radiation (R).