Articles

Dufour and Soret effects on pulsatile hydromagnetic flow of Casson fluid in a vertical non-Darcian porous space

Suripeddi Srinivasa
VIT-AP University, India
bChalla Kalyan Kumar
Vellore Institute of Technology, India
cAnala Subramanyam Reddyc
Vellore Institute of Technology, India

Dufour and Soret effects on pulsatile hydromagnetic flow of Casson fluid in a vertical non-Darcian porous space

Nonlinear Analysis: Modelling and Control, vol. 27, núm. 4, pp. 669-683, 2022

Vilniaus Universitetas

Recepción: 26 Julio 2021

Revisado: 31 Enero 2022

DOI: https://doi.org/0000-0003-4667-8338

Abstract: This article aims to inspect the pulsating hydromagnetic slip flow of Casson fluid in a vertical porous channel with heat and mass transfer. The fluid is injected into the channel from the left wall and removed at the opposite wall with the same velocity. The impact of non-Darcy, Soret, and Dufour effects are taken under consideration. The governing partial differential equations (PDEs) are converted to ordinary differential equations (ODEs) using perturbation method and solved by utilizing 4th-order Runge–Kutta (R–K) technique together with shooting method. The impact of dissimilar parameters on flow, heat and mass transfer characteristics are displayed and discussed.

Keywords: Casson fluid, slip parameter, pulsatile flow, convective boundary, Dufour and Soret effects.

1 Introduction

Studies pertaining to the MHD flows of non-Newtonian fluids in a porous medium are important because of its applications in irrigation problems, process of petroleum, heat- storage beds, paper, textile, and polymer composite industries. The most famous among these fluids is Casson fluid. Casson fluid model was introduced by Casson [4] for predic- tion of the flow behaviour of pigment-oil suspensions in lithographic varnishes. We can characterize Casson fluid as a shear thinning liquid which is considered to have an infinite viscosity at zero rate of shear, a yield stress below which no flow occurs, and a zero viscos- ity at an infinite rate of shear. The examples of Casson fluid are concentrated fruit juices, honey, tomato sauce, jelly, and human blood [8, 10, 14, 17]. Chamkha [7] investigated the hydromagnetic fully developed laminar mixed convection flow in a vertical channel with symmetric and asymmetric wall heating conditions in the presence or absence of heat generation or absorption effects. Loganathan and Deepa [21] researched the EMHD flow of the Casson fluid on a permeable Riga-plate. Pantokratoras [25] investigated the natural convection of non-Newtonian power-law fluids over a vertical plate. Walawender et al. [33] employed Casson model for portraying blood flow curves.

Pulsatile flow in a porous pipe or a channel is an important study due to its application in both engineering systems (microelectromechanical systems or MEMS, pulse combus- tors, IC engines, and reciprocating pumps) and natural systems (vascular diseases, respi- ratory systems, circulatory systems). Especially, the pulsatile flow in a porous channel is significant in the dialysis of blood in artificial kidneys [1, 3, 16, 19, 26]. Chamkha [6] ex- amined the problem of flow and heat transfer of two electrically conducting and heat gen- erating or absorbing immiscible fluids in a vertical infinitely long channel in the presence or absence of a porous medium and applied magnetic field. Kumar et al. [20] analyzed the problem of fully-developed convective flow of a micropolar and viscous fluid between parallel plate vertical channels with asymmetric wall temperature distribution. Celli and Kuznetsov [5] investigated the pulsatile viscous flow inside a horizontal infinitely wide channel. Haddad et al. [12] examined the pulsating laminar and incompressible fully developed pipe flows. Srinivas et al. [30] studied the pulsatile flow of hydromagnetic Casson fluid in a porous channel. Wang [34] illustrated the pulsating flow in a porous channel.

The fluids displaying boundary slip are significant in technological applications like inertial cavities and polishing of artificial heart valves [22,27,35]. The study of convective boundary condition is of extraordinary significance due to its application in engineering and industrial processes like material drying and transpiration cooling process [9, 31]. Malathy et al. [23] have elucidated the influences of slip and thermal radiation on MHD pulsatile flow of an Oldroyd-B fluid in a porous space with convective boundary condition. Sayed et al. [28] studied the peristaltic transport of nanofluid in an inclined asymmetric channel in the presence of slip and convective boundary conditions.

The energy flux caused by a concentration gradient was found in 1873 by Dufour and was correspondingly named the Dufour effect. It is additionally known as the diffusion- thermo effect. Then again, mass flux is able to be made by a temperature gradient as was recognized by Soret. This is the thermal-diffusion effect. Dufour and Soret effects assume a significant role when large density contrasts exist in the flow regime. Radiation, Dufour, and Soret effects on MHD flows emerge in applied physics and numerous areas of engineering like catalytic reactors, heat insulation, geothermal systems, MHD generators, and drying technology [13, 24, 29]. Dzulkifli et al. [11] numerically discussed the Dufour and Soret parameters on the boundary layer flow in nanofluid through shrinking/stretching sheet. Khan et al. [15] discussed the cross diffusion effects on Carreau–Yasuda fluid flow over a porous stretchable surface. Umavathi and Chamkha [32] examined the stability analysis of cross diffusion when a nanofluid saturated with porous space was filled in a horizontal channel. Recently, Kumar and Srinivas [18] performed the simulation for pulsation flow of hydromagnetic Casson fluid in a vertical channel with Dufour and Soret effects.

The literature reveals that no study related to hydromagnetic pulsatile Casson fluid in a vertical permeable channel has been explored so far. Inspired by the past investigations [1, 23, 26, 34] and keeping in the perspective on wide applications, we made an endeavour to portray the cross diffusion impacts on pulsating hydromagnetic slip flow of Casson fluid in a vertical porous channel with convective boundary. The coupled PDEs are converted to ODEs using perturbation method and solved by utilizing 4th-order R–K technique together with shooting method.

2 Formulation of the model

Consider the laminar and incompressible pulsating flow of Casson fluid between two vertical parallel walls at a distance The strength .. of a uniform magnetic field is applied opposite to the flow direction. We assume that the plates are very wide and very long, so that the flow is essentially axial. So that only component ofof of the velocity does not vanish. The condition of fully developed flow implies that . Since the velocity is solenoidal, we obtain . As a consequence, the velocity component is constant in any channel section and is equal to zero at the channel walls, so must be vanishing at any position. The -momentum balance equation can be expressed as

(see [2, 12, 19]). The slip parameter, Joule heating, convective boundary, Dufour and Soret effects are considered. A Cartesian coordinate system is taken so that the .axis is taken along the flow direction (vertical), and -axis is orthogonal to the walls (see Fig. 1). The channel walls possess the characteristics of convective-type boundary condition. The left and right walls maintain temperatures are and con- centrations are , respectively. The convective boundary conditions at the left wall and the right wall are respectively. The fluid is injected into the channel from the left wall with a velocity and removed at the opposite wall with the same velocity. The stress and strain relationship is designed as (Kumar et al. [19], Loganathan and Deepa [21]):

Figure 1. Schematic diagram of vertical porous channel.

[Figure 1. Schematic diagram of vertical porous channel.]

The component of the deformation rate and shear stress tensor are respectively. is yield stress, is the product of shear rate components,. is the critical value of is the plastic dynamic viscisity. The governing equations are (Adesanya et al. [1], Radhakrishnamacharya and Maiti [26])

(1)

(2)

(3)

where is the thermal conductivity, . is the dynamic viscosity,. is the electrical conductivity,. is the frequency,. is the kinematic viscosity,. is the density of the fluid, are the temperature and concentration of the fluid are the porosity and permeability of porous medium, is the coefficient of concentration expansion, is the coefficient of mass diffusivity, is the 1st-order chemical reaction rate, is time is the acceleration due to gravity ,

is the Casson fluid parameter, Tm is the mean temperature of the fluid, kTis the thermaldiffusion ratio, Cs is the concentration susceptibility, qr is the radiative heat flux, Q0 is thecoefficient of heat source/sink, Cbis the form of drag coefficient, Cpis the specific heat atconstant pressure. Thermal radiation is simulated utilizing the Rosseland approximation(Makinde et al. [22]), and as per this, qris specified by

where χ and σˆ are the Rosseland mean absorption coefficient and Stefan–Boltzmann constant. Assuming an adequately small temperature difference in the flow and expanding by Taylor’s series about T0, we get (higher-order terms are neglected). In view of qr and T, Eq. (2) becomes

(4)

The boundary conditions are (Malathy et al. [23], Xinhui et al. [35])

where are slip coefficient at the surface of the porous walls and heat coefficient, respectively.

The following dimensionless quantities are invoked:

(5)

Using Eq. (5), Eqs. (1), (4), and (3) transformed to

(6)

(7)

(8)

where U is the characteristic velocity, are dimensionless temperature and concentration, is solutal Grashof number, is the radiation parameter, is the Hartmann number, is the Forchheimernumber, is the Eckert number,is cross flow Reynolds number,is the Soretnumber , is the chemical reaction parameter, is the Schmidtnumber,is the Dufour number, and

The corresponding boundary conditions are

where is the heat transfer Biot number, and is the slip parameter.

3 Solution of the problem

To acquire the solution of Eqs. (6)–(8), a perturbative solution has been assumed in thefollowing form:

and neglecting higher orders. Here ϵ is the suitably chosen positive quantity, φ1(y) is unsteady concentration profile, θ1(y) is unsteady temperature profile, u1(y) is unsteady velocity profile, φ0(y) is steady concentration profile, θ0(y) is steady temperature profile, u0(y) is steady velocity profile, u is nondimensional velocity, and A0, A1 are positive constants.

Substituting Eq. (9) into Eqs. (6)–(8) and comparing the coefficients of various powers of , we obtain

(10)

(11)

(12)

(13)

(14)

(15)

The corresponding boundary conditions are

(16)

Further, the dimensionless Nusselt and Sherwood numbers at the walls are given by

It is noted that the system of ODEs (10)–(15) along with associated boundary conditions(16) is nonlinear and coupled. We have employed the 4th-order Runge–Kutta techniquetogether with shooting method for finding the numerical solution. Throughout thecalculations, the employed parametric values unless otherwise stated.

4 Results and discussion

The influences of various physical parameters on velocity, temperature, concentration profiles are elucidated graphically in Figs. 2–5. Figures 2(a)–2(f) describe the impact of the slip parameter (L), Darcy number of the porous media (Da), Casson fluid parameter (β), Hartmann number (M ), Forchheimer number (Fs), and Grashof number (Gr ) on the velocity profile. Figure 2(a) reveals that an increase in L results in rise of velocity profile. Figure 2(b) elucidates that the velocity increases with an enhancing Darcy number. Figures 2(c) and 2(f) depict that there is an enhancement in velocity with an enhancing Casson fluid parameter and Grashof number. Figure 2(d) depicts that for a rise in M , there is a decrease in velocity. This can be because of Lorentz forces created by the applied magnetic field act as resistive drag forces opposite to the flow direction. Hence there is a decrease in velocity. Figure 2(e) delineates that a rise in Forchheimer number creates a resistance in fluid flow which results an abatement in velocity.

The variation of temperature distribution for various values of Bi, Rd , Du, and Sr are shown in Figs. 3(a)–3(d). Figure 3(a) illustrates the influence of B. on .. It is observed that for a given rise Biot number, there is a decrease in the temperature. A superior value

Figure 2. Velocity distributions.
Figure 2. Velocity distributions.

observe

of B. includes a higher degree of convective cooling at the channel walls, subsequently inflicting lower temperature at the channel walls and additionally within the bulk fluid. It is predictable that as. , the convective boundary conditions will turn into the prescribed wall temperatures. Figure 3(b) depicts the effect of thermal radiation parameter on .. It is observed that. increases with an enhancement of Rd . This phenomenon can be ascribed to the physical fact that the thermal boundary layer thickness rises with an enhancing Rd . The variation of temperature with respect to the Dufour number is shown in Fig. 3(c). It is seen that there is an enhancement in temperature with an enhancing of Du. Actually, Du is associated in energy flux caused by a concentration gradient.

Figure 3. Temperature distributions.
Figure 3. Temperature distributions.

Figure 4 illustrates the influences of Sc, Sr , and on concentration distribution. Figure 4(a) exhibits that the concentration profile decreases for given rise in Schmidt number. This is observed due to a rise in Sc that it turn makes the concentration boundary layer thinner than momentum boundary layer. Figure 4(b) depicts the influence of Soret number on. One can infer that φ is enhanced with increasing Sr . Figure 4(c) depicts effect on concentration. It is seen that the concentration falls with a rise in the destructive chemical reaction. The contrary pattern can be noticed for the case of generative chemical reaction. The effect of Du and Sr on the temperature and concentration distributions is shown in Figs. 5(a)–5(b). Figure 5(a) elucidates that a rise in Du with a decrease in Sr rises the thermal boundary layer growth. It is noticed from Fig. 5(b) that a decrease in Sr with a rise in Du has the tendency to decrease the concentration distribution. Physically, Soret effect reports that mass flux is made once a system is underneath a temperature gradient. Further, a rise in Du improves the convention velocity over combined influences of thermal and solutal buoyancy forces that leads to enhance the heat transfer but fall the mass transfer of the fluid.

The present numerical values corresponding to the Nusselt number compared with the previously published numerical results of Kumar et al. [19] are shown in Table 1. This comparison shows that the present results in limiting case are in good agrement with the

Figure 4. Concentration distributions.
Figure 4. Concentration distributions.

Figure 5. Effect of Du and Sr on the temperature and concentration distributions.
Figure 5. Effect of Du and Sr on the temperature and concentration distributions.

published results. Table 2 compares the findings obtained by an analytical method (double perturbation) with the results obtained by a numerical method (4th-order R–K method together with shooting technique) to verify the validity of the current model. Table 3 shows the variations in Nu and Sh for different values of. The Nusselt and Sherwood numbers at the left and right walls are denoted by, and, respectively. It is noticed that Nu falls with a rise in Forchheimer number at the left wall, while it rises at the right wall. The Nusselt number rises at both walls by

Table 1. Comparison of results of Nusselt number for several values of Casson fluid parameter in the absence of Brownian motion parameter, thermophoresis parameter, Lewis number, Dufour number, cross flow Reynolds number, Forchheimer number, slip parameter, and the heat transfer Biot number.
Nu0Nu1Nu0Nu1
2 36.00734 5.37199−5.80169 −5.139256.00734 5.37199−5.80169 −5.13925

Table 2. Comparison of the analytical results with the numerical results by taking R =Fs = Du = Bi = Ec = 0.
u(y)θ(y)φ(y)u(y)θ(y)φ(y)
10.10.20960.06480.11120.20950.06480.1112
20.30.40450.20190.33170.40450.20190.3317
30.50.53660.36240.54560.53670.36240.5456
40.70.55620.56470.74650.55610.56470.7465
50.90.39010.83240.92430.39010.83240.9243

Table 3. Effect of Fs, β, Bi, Du, Sr , Sc on Nusselt and Sherwood number distributions when A0 = 1, A1 = 1, γ = 1, Gr = 7, Rd = 2, Q = 1, Da = 0.5, L = 0.06, R = 1, t = π/4, H = 2, ε = 0.01, Gc = 7, Re = 3, Pr = 21, Ec = 0.5, K1 = 0.001, M = 2.
Parameter Values Nu0Nu1Sh0Sh1Fsβ 0 2 4 6 1 2 3 4 0.75329 0.69289 0.64031 0.59371 0.65735 0.73733 0.77105 0.79022 −0.83727 −0.82247 −0.81183 −0.80410 −0.75266 −0.83309 −0.88217 −0.91441 0.43837 0.41899 0.40492 0.39459 0.33749 0.43291 0.49110 0.52927 3.42651 3.30946 3.21333 3.13252 3.08491 3.39487 3.56435 3.67176 Bi 1 2 3 4 −1.29919 −0.89213 −0.79586 −0.75594 −1.24512 −1.37839 −1.63896 −1.87780 0.87477 0.78766 0.73488 0.69815 1.78748 2.36379 2.75861 3.07305 Du 0.2 0.3 0.4 0.6 2.01968 1.13861 −0.47402 −1.94072 0.35629 0.00978 −0.45283 −1.16345 −1.05216 −0.61865 0.77237 1.98978 1.74097 0.63832 0.74398 Sr 0 1 2 3 −1.41501 −1.36021 −1.29919 −1.23228 −1.28365 −1.26557 −1.24512 −1.22238 0.63868 0.77031 0.87477 0.94845 1.56538 1.66325 1.78748 1.94132 Sc 0.22 −1.36315 −1.26158 0.91487 1.22006

rising Sr , whereas it falls at both walls by rising Du. One can observe that at the left wall, Nu enhances with a rise in the, while it falls at the right wall. One can infer that at both walls, Sh rises by rising . and Sr , while it falls at both walls by increasing Fs. It can be noticed that Sh decreases with an increase in B. and Sc at the left wall, whereas it rises at the right wall. The reverse trend can be seen for the case of Dufour number.

5 Conclusions

An analysis is made on pulsating MHD slip flow of Casson fluid in a vertical non- Darcian porous space with convective boundary condition. The considered investigation is significant as flow of Casson fluids (drilling muds, greases, clay coating, certain oils, numerous impulsions, and blood) in porous channel are utilized in modelling biological and science research. The analytical and numerical solutions are constructed for flow variables. The salient points of this investigation are as follows:




The velocity increases with an increasing slip and Casson fluid parameter, while itfalls with an enhancing Forchheimer number.
It is noticed that the temperature rises with an enhancing Dufour number and radiation parameter, while it falls with an enhancing Soret number.
Temperature enhances for a given rise in Dufour number with a decrease in Soretnumber.
Concentration falls by enhancing Schmidt number and rises by rising the Soret number.
Nusselt number escalates with the increase in Casson fluid parameter and heattransfer Biot number at the left wall, while it is a decreasing function at right wall.
Sherwood number rises at the both walls for a given increase in Casson fluid parameterand Soret number.

Fuente: The velocity increases with an increasing slip and Casson fluid parameter, while itfalls with an enhancing Forchheimer number. It is noticed that the temperature rises with an enhancing Dufour number and radiationparameter, while it falls with an enhancing Soret number. Temperature enhances for a given rise in Dufour number with a decrease in Soretnumber. Concentration falls by enhancing Schmidt number and rises by rising the Soretnumber. Nusselt number escalates with the increase in Casson fluid parameter and heattransfer Biot number at the left wall, while it is a decreasing function at right wall. Sherwood number rises at the both walls for a given increase in Casson fluid parameterand Soret number.

References

1 S.O. Adesanya, J.A. Falade, O.D. Makinde, Pulsating flow through vertical porous channel with viscous dissipation effect, U.P.B. Sci. Bull., Ser. D, 77(1):25–36, 2015, https://www. scientificbulletin.upb.ro/rev_docs_arhiva/fulldb6_966044.

2 S.O. Ajadi, A note on the unsteady flow of dusty viscous fluid between two parallel plates, J. Appl. Math. Comput.,18:393–403, 2005, https://doi.org/10.1007/ BF02936582.

3 R. Bhargava, H.S. Takhar, S. Rawat, T.A. Beg, O.A. Beg, Finite element solutions for non-Newtonian pulsatile flow in a non-darcian porous medium conduit, Nonlinear Anal. Model. Control, 12(3):317–327, 2007, https://doi.org/10.15388/NA.2007.12. 3.14690 .

4 N. Casson, A flow equation for the pigment oil suspensions of the printing ink type, in C.C. Mill (Ed.), Rheology of Disperse Systems, Pergamon, New York, 1959, pp. 84–102.

5 M. Celli, A.V. Kuznetsov, A new mechanism for buoyancy driven convection in pulsating viscous flows: A theoretical study, Int. J. Heat Mass Transfer, 118:340–348, 2018, https://doi.org/10.1016/j.ijheatmasstransfer.2017.10.112.

6 A.J. Chamkha, Flow of two-immiscible fluids in porous and nonporous channels, J. Fluids Eng., 122(1):117–124, 2000, https://doi.org/10.1115/1.483233.

7 A.J. Chamkha, On laminar hydromagnetic mixed convection flow in a vertical channel with symmetric and asymmetric wall heating conditions, Int. J. Heat Mass Transfer, 45:2509–2525, 2002, https://doi.org/10.1016/S0017-9310(01)00342-8.

8 A.J. Chamkha, T. Grosan, I. Pop, Fully developed free convection of a micropolar fluid in a vertical channel, Int. Commun. Heat Mass Transfer, 29(8):1119–1127, 2002, https://doi.org/10.1016/S0735-1933(02)00440-2.

9 T. Chinyoka, O.D. Makinde, Analysis of transient generalized couette flow of a reactive variable viscosity third-grade liquid with asymmetric convective cooling, Math. Comput. Modelling, 54(1–2):160–174, 2011, https://doi.org/10.1016/j.mcm.2011.01. 047.

10 R.K. Dash, K.N. Mehta, G. Jayaraman, Casson fluid flow in a pipe filled with a homogeneous porous medium, Int. J. Eng. Sci., 34(10):1145–1156, 1996, https://doi.org/10. 1016/0020-7225(96)00012-2.

11 N.F. Dzulkifli, N. Bachok, I. Pop, N.A. Yacob, N.Md. Arifin, H. Rosali, Soret and Dufour effects on unsteady boundary layer flow and heat transfer of nanofluid over a stretching/shrinking sheet: A stability analysis, J. Chem. Eng. Process Technol., .(3):1–9, 2017, https://doi.org/10.4172/2157-7048.1000336.

12 K. Haddad, O. Ertunc, M. Mishra, A. Delgado, Pulsating laminar fully developed channel and pipe flows, Phys. Rev. E, 81(1):016303, 2010, https://doi.org/10.1103/ PhysRevE.81.016303.

13 A.S. Idowu, B.O. Falodun, Soret–Dufour effects on MHD heat and mass transfer of Walter’s-B viscoelastic fluid over a semi-infinite vertical plate: spectral relaxation analysis, J. Taibah Univ. Sci., 13(1):49–62, 2018, https://doi.org/10.1080/16583655.2018.1523527.

14 A. Khalid, I. Khan, A. Khan, S. Shafie, Unsteady MHD free convection flow of Casson fluid past over an oscillating vertical plate embedded in a porous medium, Eng. Sci. Technol. Int. J., 18:309–317, 2015, https://doi.org/10.1016/j.jestch.2014.12.006.

15 M.I. Khan, T. Hayat, S. Afzal, M.I. Khan, A. Alsaedi, Theoretical and numerical investigation of Carreau-Yasuda fluid flow subject to Soret and Dufour effects, Comput. Methods Programs Biomed., 186:105145, 2020, https://doi.org/10.1016/j.cmpb.2019.105145.

16 C.K. Kumar, S. Srinivas, Simultaneous effects of thermal radiation and chemical reaction on hydromagnetic pulsatile flow of a Casson fluid in a porous space, Eng. Trans., 65(3):461–481, 2017, https://et.ippt.gov.pl/index.php/et/article/view/400.

17 C.K. Kumar, S. Srinivas, Influence of joule heating and thermal radiation on unsteady hydromagnetic flow of chemically reacting Casson fluid over an inclined porous stretching sheet, Spec. Top. Rev. Porous Media, 10(4):385–400, 2019, https://doi.org/10. 1615/SpecialTopicsRevPorousMedia.2019026908.

18 C.K. Kumar, S. Srinivas, Pulsating hydromagnetic flow of Casson fluid in a vertical channel filled with non-Darcian porous medium, Heat Transfer, 50(6):5225–5239, 2021, https://doi.org/10.1002/htj.22121.

19 C.K. Kumar, S. Srinivas, A.S. Reddy, MHD pulsating flow of Casson nanofluid in a vertical porous space with thermal radiation and joule heating, J. Mech., 36(4):535–549, 2020, https://doi.org/10.1017/jmech.2020.5.

20 J.P. Kumar, J.C. Umavathi, A.J. Chamkha, I. Pop, Fully-developed free-convective flow of micropolar and viscous fluids in a vertical channel, Appl. Math. Modelling, 34(5):1175–1186, 2010, https://doi.org/10.1016/j.apm.2009.08.007.

21 P. Loganathan, K. Deepa, Computational exploration of Casson fluid flow over a Riga-plate with variable chemical reaction and linear stratification, Nonlinear Anal. Model. Control,,25(3):443–460, 2020, https://doi.org/10.15388/namc.2020.25.16659.

22 O.D. Makinde, M.G. Reddy, K.V. Reddy, Effects of thermal radiation on MHD peristaltic motion of Walters-B fluid with heat source and slip conditions, J. Appl. Fluid Mech., 10(4):1105–1112, 2017, https://doi.org/10.18869/ACADPUB.JAFM.73.241. 27082.

23 T. Malathy, S. Srinivas, A.S. Reddy, Chemical reaction and thermal radiation effects on MHD pulsatile flow of an Oldroyd-B fluid in a porous medium with slip and convective boundary conditions, J. Porous Media, 20(4):287–301, 2017, https://doi.org/10. 1615/JPorMedia.v20.i4.10.

24 D. Pal, G. Mandal, K. Vajravelu, Soret and Dufour effects on MHD convective-radiative heat and mass transfer of nanofluids over a vertical non-linear stretching/shrinking sheet, Appl. Math. Comput.,287-288:184–200, 2016, https://doi.org/10.1016/j.amc.2016. 04.037.

25 A. Pantokratoras, Natural convection over a vertical isothermal plate in a non-Newtonian power-law fluid: New results, Adv. Mech. Eng., .(5):1–12, 2016, https://doi.org/ 10.1177/1687814016644111.

26 G. Radhakrishnamacharya, M.K. Maiti, Heat transfer to pulsatile flow in a porous channel, Int. J. Heat Mass Transfer, 20(2):171–173, 1977, https://doi.org/10.1016/0017- 9310(77)90009-6.

27 M. Sajid, M. Awais, S. Nadeem, T. Hayat, The influence of slip condition on thin film flow of a fourth grade fluid by the homotopy analysis method, Comput. Math. Appl., 56:2019–2026, 2008, https://doi.org/10.1016/j.camwa.2008.04.022.

28 H.M. Sayed, E.H. Aly, K. Vajravelu, Influence of slip and convective boundary conditions on peristaltic transport of non-Newtonian nanofluids in an inclined asymmetric channel, Alexandria Eng. J., 55:2209–2220, 2016, https://doi.org/10.1016/j.aej.2016. 04.041.

29 S. Shateyi, S.S. Motsa, P. Sibanda, The effects of thermal radiation, hall currents, Soret, and Dufour on MHD flow by mixed convection over a vertical surface in porous media, Math. Probl. Eng., 2010:Article ID 627475, 2010, https://doi.org/10.1155/2010/ 627475.

30 S. Srinivas, C.K. Kumar, A.S. Reddy, Pulsating flow of Casson fluid in a porous channel with thermal radiation, chemical reaction and applied magnetic field, Nonlinear Anal. Model. Control, 23(2):213–233, 2018, https://doi.org/10.15388/NA.2018.2.5.

31 S. Srinivas, T. Malathy, A.S. Reddy, A note on thermal-diffusion and chemical reaction effects on MHD pulsating flow in a porous channel with slip and convective boundary conditions, J. King Saud Univ., Eng. Sci., 28(2):213–221, 2016, https://doi.org/10.1016/j. jksues.2014.03.011.

32 J.C. Umavathi, A.J. Chamkha, Stability analysis of cross diffusion for the walters B fluid model saturated with permeable nanofluid, J. Therm. Sci. Eng. Appl., 11(4):041014, 2019, https://doi.org/10.1115/1.4044078.

33 W.P. Walawender, T.Y. Chen, D.F. Cala, An approximate Casson fluid model for tube flow of blood, Biorheology, 12(2):111–119, 1975, https://doi.org/10.3233/bir-1975- 2202.

34 C.Y. Wang, Pulsatile flow in a porous channel, J. Appl. Mech., 38(2):553–555, 1971, https://doi.org/10.1115/1.3408822.

35 S. Xinhui, Z. Liancun, C. Xuehui, Z. Xinxin, C. Limer, L. Min, The effects of slip velocity on a micropolar fluid through a porous channel with expanding or contracting walls, Comput. Methods Biomech. Biomed. Eng., 17(4):423–432, 2014, https://doi.org/10.1080/10255842.2012.688108.

HTML generado a partir de XML-JATS4R por