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 kalyankumar.challa@gmail.com

Vellore Institute of Technology, India

cAnala Subramanyam Reddyc

Vellore Institute of Technology, India

Esta obra está bajo una Licencia Creative Commons Atribución-NoComercial-SinDerivar 4.0 Internacional.

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.

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, k_Tis the thermaldiffusion ratio, C_s is the concentration susceptibility, q_r is the radiative heat flux, Q₀ is thecoefficient of heat source/sink, C_bis the form of drag coefficient, C_pis the specific heat atconstant pressure. Thermal radiation is simulated utilizing the Rosseland approximation(Makinde et al. [22]), and as per this, q_ris 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 T₀, we get (higher-order terms are neglected). In view of q_r 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.

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.

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.

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

Table 2. Comparison of the analytical results with the numerical results by taking R =Fs = Du = B_i = Ec = 0.

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.

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.

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.