Abstract

In this study, the competency of numerical computational framework based on Lobatto IIIA technique is utilized for dynamical analysis of the Darcy-Forchheimer flow of Sisko nanomaterial with nonlinear thermal radiation. The resultant PDEs of the Sisko fluid model expressions are transformed into system of nonlinear ODEs by exploiting the similarity variables. Graphical representations and numerical illustrations are used to envisage the characteristics of various physical parameters of interest on velocity profile, nanoparticles concentration, and temperature distribution of Sisko fluidic system. In addition, skin friction and Nusselt number are numerically examined with observation that the material parameter of Sisko fluid increases the velocity profile as well as Nusselt number while decreasing temperature, concentration profiles, and skin friction coefficient.

1. Introduction

Material of the type having pores is termed as porous medium, in fluid dynamics; usually it is considered that the porous medium is filled with nanofluid. The study of nanofluid flow through porous medium has wide range of applications due to the fact that it has become a field of high interest for researchers, mathematicians, engineers, geologists, and architectures. Some of the considerable applications are flow of liquid in reservoirs, catalytic reactors, production of oil, heat exchanger, etc. Transfer of heat in porous medium has extensively been used in the heat pipe technology, nonwoven ingredients, production of papers, electronic based technology, storage of energy, etc. [14]. Henry Darc, in 1856, a French engineer, pioneered the flow of homogeneous liquids through porous medium during his valuable work on liquid flow through the sand beds. The law of Darcy is not effective for the situation when boundary and inertial effects appears at higher rate of fluid flow. After that the work was extended by Dutch scientist named Philippes Forchheimer in 1901. He extended the velocity equation by adding the square of velocity term in the equation of momentum of the problem to predict the behavior of boundary layer flow and inertia [5]. The term Forchheimer was later introduced by Muskat [6]. Pal and Mondal [7] constructed a fluid model with Darcy-Forchheimer effects over stretching sheet and he examined that with the decrease of the nanoparticle concentration profile, the value of electric field parameter increased. Ganesh et al. [8] studied the Darcy-Forchheimer MHD nanoliquid.

Thermal radiation, either linear or nonlinear, has a very significant role in different high temperature processes in the field of fluid dynamics. It has diverse uses in several types of power plants and industries. For final possessions of desired qualities, the control of heat transfer is very significant and necessary in industries. Modern system linked to astrophysical flow, reactors cooling, electrical power generators, vehicles in the space, and plasmas are working by useful applications of radiation either linear or nonlinear. Latif et al. [9] studied the nonuniform viscosity fluid model with nonlinear radiation. Unsteady flow in existence of source/sink and radiation effects is modeled by Hayat et al. [10]. Pal [11] studied thermal radiation in a channel the peristaltic transport by using convective conditions. Pandey and kumar [12] studied the thermal radiation as well as mass flux in porous medium on natural convection. The effect of thermal radiation natural convection has been investigated by Huang [13]. Kothandapani and Prakash [14] worked out upon the influence of radiation over the peristalsis of magneto hydrodynamics in a channel. Some renowned researchers also worked on thermal radiation [1523]. In these studies, the radiation with nonlinearity was introduced by the researchers with in fluid flow models. In these investigations extra variable contains difference between temperature of wall and ambient temperature. Due to high temperature variation, the transfer of radiation amongst the two objects is fortified. More details regarding nonlinear radiation are given [2433]. Some recent work related to Non-Newtonian fluids, radiation, porous medium, and nanofluid is provided here [3437]. These are the motivating factors for the authors to investigate in the model and analysis of nonlinear radiated nanofluid flow of Sisko fluid. The noticeable features of the study are listed as follows.(i)An innovative design of Sisko nanofluid with nonlinear radiation and Darcy-Forchheimer flow involving porous medium has been modeled and PDEs of the model are transformed into an equivalent system of ODEs using suitable similarity transformations.(ii)The concept of boundary layer is incorporated to develop the fluidic problem. The energy expression of the problem is developed by adding nonlinear radiation term.(iii)Numerical analysis of the Sisko fluid model involving nanoparticles is presented by invoking the Lobatto IIIA method in order to investigate the effects of involved physical variables on velocity, concentration, and temperature by variation of physical quantities.(iv)The application of Lobatto IIIA method to this model is an innovative work with these effects. Lobatto IIIA is a version of bvp4c method, which is based on finite difference method. Other methods can also solve this problem, i.e., NDSolve, RK-4, HAM, etc. this method has a strength to solve highly nonlinear systems of ODEs.(v)Results presented in numerical and graphical representations endorsed the worth of the methodology as viable and effective computing solver for visualizing the dynamics of Sisko nanofluidic model.

Outline of the rest of the paper is given as follows: mathematical model of the problem is represented in Section 2, in Section 3, results and discussion of numerical experiments are provided, while the conclusion and future research work are presented in the last section.

2. Mathematical Development

Consider two-dimensional Sisko nanofluid flow through a porous medium and nonlinear stretching sheet as shown in Figure 1. The coordinate y is taken as normal to the sheet while x-axis is along the sheet. A hot fluid is taken for heating or cooling the surface of sheet connectively. It is considered that uniformly concentrated fluid with nanoparticles is taken at the surface of the stretching sheet. The two dimensional governing equations with boundary layer are given below [38, 39]:with the boundary conditionswhere the inertia coefficient F is defined byRosseland radiation approximation is given asTerm in the above equation is Stefan-Boltzmann parameter while is mean absorption parameter. After simplifying (7), we haveIn order to obtain ODEs from the above PDEs, we use the similarity transformationsContinuity equation is identically satisfied, while (2)–(4) yieldwith boundary conditionswhere parameters appearing in (9)-(16) are defined aswhere A is material parameter of Sisko fluid. The local Nusselt number is and is the heat flux, given by

3. Results and Discussion

The transformed fluid flow problem defined in (10)–(13) has been solved numerically by exploiting the strength of Labotto IIIA method and invoking MATLAB. routine “bvp4c” for solution of boundary value problem of differential equation. We refer the interested reader to the article [40] for the necessary additional details of the Labotto III solver in [40] and citations therein. The solutions of (10)–(13) are also calculated with explicit Runge-Kutta (RK) method by using the Mathematica routine “NDSolve” for numerical solution of differential equations. The results of Lobatto III A are compared with explicit RK method and error analysis shows the matching of the order 10−08 to 10−09 exist between both methodologies which validate the correctness of the proposed numerical procedure. Therefore, in the rest of the study numerical procedure of Labotto IIIA is used for numerical solution of the problem.

Characteristics of some influential parameters like heat Prandtl number , the Froude number Fr, Brownian motion Nb parameter, thermophoresis number Nt, Lewis number Le, fitted rate constant n, the material parameter of Sisko fluid A, the radiation parameter Rd, the porosity permeability parameter , the heat generation/absorption variable , the Biot number , and the temperature ratio parameter on the velocity , temperature , and concentration are shown. These numerical outcomes have been illustrated in Figures 210 and Tables 17. Numerical computation has been performed by choosing the default estimations of involved variables in the problem as , , , , , , , , and . Figures 2(a), 2(b), and 2(c) show the influence of Fr on , , and . For increasing values of Fr, and decrease while increases. Figure 2(d) discloses the effect of n on . By increasing n, shows decreasing behavior. Figures 3(a), 3(b), and 3(c) depict the decreasing behavior of , , and , respectively, for increasing values of n. In Figure 3(d) increasing quality is illustrated for for higher values of Prandtl number Pr. The concentration profile shows decreasing characteristic for increasing Prandtl number Pr, Lewis number Le, and Brownian number Nb, which are shown in Figures 4(a), 4(b), and 4(d) respectively. Increase in Lewis Number Le has decreasing effect because the larger Lewis number makes diffusivity of molecules smaller and due to this reason, the concentration field is decreased. The temperature profile is enhanced in Figure 4(c) for higher Nb. In Figures 5(a) and 5(b), it is clearly seen that increasing Nt causes increasing temperature and concentration profiles. Enhancement in temperature profile due to increase in thermophoresis number is because of nanoparticles. The nanoparticles which are closer to hot fluid shift towards cold fluid due to increase in thermophoresis parameter, which makes the temperature boundary layer thick. It is obvious in Figures 5(c) and 5(d) that increase in A causes increase in and , respectively. From Figures 6(a) and 6(b), it is clear that temperature and concentration , respectively, both are decreasing for higher values of A. Also, Figures 6(c) and 6(d) both represent the decreasing behavior of and , respectively, for enhanced values of parameter s. Enhancement is seen in and for increase in s, which is shown in Figures 7(a) and 7(b), respectively. An increase in results in decrease in and , which is visualized in Figures 7(c) and 7(d), respectively, and increasing effect of is shown for and for increasing values of in Figures 8(a) and 8(b), respectively. It is found in Figures 8(c) and 8(d) that and are enhanced, respectively, for higher estimations of . Figures 9(a), 9(b), 9(c), and 9(d) visualize the effect of on , , , and , respectively. Higher estimations of cause decrease in and , while increase in and dual behavior in concentration profile . In the domain from 0 to 2, decreases while in the domain from 2 to 7 it shows opposite behavior. Figure 10(a) shows that increase in Rd causes decrease in and dual behavior of is seen for higher estimations of Rd in Figure 10(b). In the domain from 0 to 2, increases while in the domain from 2 to 7, decreases. Figure 10(c) discloses the increasing behavior of for higher values of the temperature ratio parameter on temperature distribution . This is because greater ratio in the temperature increases the temperature boundary layer. Figure 10(d) implies the dual characteristic of for higher values of . In the domain from 0 to 2, diminishes while in domain from 2 to 7, enhances. Table 1 shows the variation of values of involved parameters in different scenarios. Scenario 1 shows the variation of Fr; scenario 2, variation of n; scenario 3, variation of Pr; scenario 4, variation of Le; scenario 5, variation of Nb; scenario 6, variation of Nt; scenario 7, variation of A; Scenario 8, variation of s; scenario 9, variation of ; scenario 10, variation of ; scenario 11, variation of ; scenario 12, variation of ; and scenario 13, variation of . Table 2 shows skin friction coefficients for various scenarios. Table 3 shows numerical values for surface heat transfer rate in different scenarios.

Error and complexity analysis of Lobatto IIIA solver is presented in Tables 47. Table 4 explains the number of meshes for each case of different scenarios. Table 5 provides maximum error for each case of different scenarios, which show that reasonable accuracy attained by the scheme. Table 6 depicts the number of ODEs evaluated for each case of various scenarios. Table 7 explains the number of boundary conditions (BCs) evaluated during each case for different scenarios. The presented complexity indices can be used for comparison of the given scheme with future similar studies for the Sisko nanofluidic model.

4. Conclusions

Characteristics of radiated Sisko nanofluid are studied. Main features of the present analysis are listed below:(i)The parameter A increases and while decreasing and .(ii)Velocity profile decreases while concentration profile increases with enhancement in s.(iii)An increase in shows decrease in and while showing increase in and .(iv)Enhancement is observed in and for higher .(v)A decrease in and while increasing behavior is observed for higher values of .(vi)Skin friction enhances for higher values of , , , and while diminishing for higher values of and .(vii)Higher values of , , , , , , , , and decay the heat transfer rate while reverse behavior is observed for n, A, , and .

Nomenclature

:Porosity parameter
:Thermal conductivity
F:Inertia coefficient
:Stefan-Boltzmann parameter
:Mean absorption parameter
:Material parameter
:Heat flux
Pr:Prandtl number
:Radiation parameter
Fr:Froude number
Le:Lewis number
Nt:Thermophoresis number
Nb:Brownian motion parameter
:Velocity profile parameter
:Temperature profile parameter
:Concentration profile parameter
:Permeability parameter
:Heat generation/absorption variable
s:Stretching rate of surface
n:Fitted rate constant
:Temperature ratio parameter
:Biot number
ODE:Ordinary differential equation
PDE:Partial differential equation
MHD:MagnetoHydroDynamics.

Data Availability

All the data used in the research is available within the article.

Ethical Approval

All the authors of the manuscript declared that there is no research involving human participants and/or animal.

All the authors of the manuscript declared that there is no material that required informed consent.

Disclosure

Authors are responsible for any fee if any associated with this manuscript.

Conflicts of Interest

All the authors of the manuscript declared that there are no potential conflicts of interest.