Heat and mass transfer characteristics in flow of bi-viscosity fluid through a curved channel with contracting and expanding walls: A finite difference approach

Abstract

This article investigates the heat and mass transfer in flow of bi-viscosity fluid through a porous saturated curved channel with sinusoidally deformed walls. The magnetic field and Joule heating effects are also taken into account. The equations describing the flow and heat/mass transfer are developed using curvilinear coordinates. A reduction of these equations is made based on lubrication approximation. The reduced linear ordinary differential equations are integrated numerically using an implicit finite difference scheme. It is observed that, the bi-viscosity fluid parameter, permeability parameter, and Hartmann number have analogous effects on the longitudinal velocity. Moreover, temperature of the fluid, heat coefficient, and mass concentration increase by increasing bi-viscosity fluid parameter, Brinkmann number, and Hartmann number. Further, mass concentration increases by increasing the rate of chemical reaction and bi-viscosity parameter. The size of circulating roll in lower half of the channel boosts up with larger variation of bi-viscosity parameter and permeability parameter. The flow patterns in the channel illustrating the effects of bi-viscosity parameter, permeability parameter, and Hartmann number are also displayed.

Keywords

Bi-viscosity fluid curved channel concentration heat and mass transfer finite difference technique

Introduction

Flow due to sinusoidally deformed vessel walls is commonly known as peristaltic flow. Peristaltic flows are caused by rhythmic contractions and relaxations of smooth muscles of the vessel. Peristalsis is a major mechanism found in the human body for physiological fluid transport. Examples of physiological fluids transported by peristaltic activity are urine, blood, chyme, spermatic fluid, etc. Modern machines which operate on the principle of peristalsis are heart-lung machine, diabetic pumps, roller and finger pumps, pharmacological delivery systems, etc.

The fluid dynamics of peristalsis are extensively investigated in the literature in past few decades. The first theoretical and experimental investigation on peristaltic flows was presented by Shapiro et al.¹ based on the absence of curvature and inertial features. This prestigious investigation has been further extended for nonlinear fluid models,^2–9 curvature consequences,^10–14 magnetohydrodynamic prospective,^15,16 electro-osmotic effects,¹⁷ and slip effects.^18–20 The important phenomena of heat and mass transfer have also been investigated by researchers for peristaltic transport as it involves motivated significance in area of laser therapy, hemodialysis, chemotherapy, and cryosurgery.^21–27 More recently peristaltic flow of nanofluids is an active area of research because of novel applications in micropumps, drug delivery systems, and pharmacological engineering.²⁸

Due to multidisciplinary applications of non-Newtonian materials, researchers have concentrated on explanation of complex physical properties in recent days. It is well justified fact that all physical properties of non-Newtonian liquids cannot be completely predicted by intruding a single mathematical relation. Following such distinguish rheological features, bi-viscosity fluid characterizes to non-Newtonian fluids because of complex rheological properties. The relation for shear stress and deformation rate is extensively different from traditional viscous fluid model. The unique feature associated with this fluid model behaves like solid when yield stress is larger as compared to shear stress and deforms continuously when amount of stress becomes dominant with respect to yield stress. The motivated application encountered by bi-viscosity liquid reports physical applications in food processing, metallurgy and drilling processes. Biviscosity fluid materials include priceless applications in many biological fluids, blood, synovial liquids, plasma, lubricants, many pharmaceutical products, proteins, sewage sludge, china clay, and many more. Number of researchers investigated rheological impact of this model in various flow scenarios. For instance Dash et al.²⁹ described the flow characteristics of bi-viscosity fluid through a homogeneous porous medium bounded by a circular tube. The mixed convection and variable thermal features for dissipative flow of bi-viscosity nanofluid have been intended by Hayat et al.³⁰ The effects of chemical reaction on flow of biviscosity fluid are thoroughly studied by Shehzad et al.³¹ and Reddy et al.³² Further, studies regarding flows of bi-viscosity fluids can be found in references.^33–35

The porous medium is a stuff that encompasses spaces between solid areas through which fluid can pass. In recent years, considerable interest has been noticed to examine the flow through saturated porous space due to its fundamental applications in petroleum engineering, soil sickness, geophysical processes, and medical applications. Some fundamental applications for saturated porous media flow include filtration and extraction of oil from wells, water draining applications for irrigation, oil flow in porous rocks, energy extraction associated with geo-thermal systems, the fluid flow in ion-exchange beds, and oil-spills cleaning. The foremost contribution on flows through porous media was made by Brinkman,³⁶ who calculated the viscous force exerted by a ﬂowing ﬂuid through a porous media. Peristaltic flows through porous medium are also an active area of research. Several authors contributed to this area. Mention may be made to the work of El Shehawey,³⁷ Mekheimer,³⁸ Kothandapani and Srinivas,³⁹ and Hayat et al.⁴⁰ in various scenarios. Javed et al.⁴¹ studied the effect of mixed convection through micropolar nanofluid under the impact of magnetohydrodynamics. In a separate study Javed et al.⁴² examined the effect of heat transfer in water–Cu nanofluid saturated porous medium through trapezoidal cavities under the influence of magnetic field. Javed et al.⁴³ investigated effect of natural convection in a ferro-fluid filled square cavity in the presence of uniform magnetic field. Mehmood et al.⁴⁴ investigated the combined effect of MHD-Mixed convection flow in a lid driven trapezoidal cavity in presence of heated bottom wall. In a separate study Mehmood⁴⁵ simulated numerically the effect of mixed thermomagnetic convection through two lid-driven entrapped trapezoidal cavities enclosing ferrofluid saturated porous medium. Mehmood and Javed⁴⁶ considered the impact of MHD on energy transmission through ferrofluid saturating porous medium contained in a lid-driven trapezoidal container with triangular notched heater configuration.

The phenomenon of heat transfer in flowing materials involves a variety of applications in many chemical and mechanical industries, thermal engineering, electronic devices, nuclear thermal-hydraulic processes, etc. Owing to such interesting applications, many authors analyzed the heat transfer characteristics in flow of different fluids. For instance, Belhocine⁴⁷ examined the heat transfer analysis in flow developed flow configured by a circular tube. The work of Belhocine and Omar⁴⁸ deals with study of convective heat transfer due to circular tube where analytical solution were computed for formulated flow problem. In another useful contribution, Belhocine and Abdullah⁴⁹ deals the flow and heat transfer of viscous fluid by using constant wall temperature assumptions, usually called the Graetz problem. The thermal boundary layer analysis for modified Levèque problem has been performed by Belhocine and Abdullah.⁵⁰ Joule heating appears due to the resistance offered to electric current passing through some conductive material. There are number of systems in which Joule heating effect has vital role such as dielectrophoretic trapping, electric fuses, PCR reactors, hot plate, microvalves for fluid control, electric heaters and stoves, thermistors and soldering irons, etc. The effect of Joule heating in peristalsis has also been investigated. For detail see references.^51–55

From literature cited above, some investigations regarding the peristaltic flow of bi-viscosity model encountered by a porous saturated straight vessel have been performed but no continuation has been utilized to explore the peristaltic phenomenon in a curved porous-saturated channel. Such motivated applications insist us to report the heat and mass transportation in magnetized bi-viscosity fluid configured by a curved channel in presence of porous medium. The objective here is to report the consequences of formulated dimensionless parameters namely Brinkmann number, bi-viscosity fluid constant, radius of curvature, Hartmann number, and porous medium permeability on combined heat/mass transportation. It is remarked that although variety of studies are available in the literature regarding the peristaltic flow of non-Newtonian fluids in curved channel but analysis regarding peristaltic transport of magnetized bi-viscosity fluid in presence of Joule heating, porous medium and heat and mass transfer phenomenon has not be investigated yet. The aim of current contribution is to fulfill this research gap. The governing equations for the concerned flow model are constituted under the assumption of dominant viscous effects as compared to the inertial features. The solution for formulated flow problem has been suggested by employing implicit finite difference scheme with desirable accuracy. Later on, a comprehensive graphical analysis for each flow parameters is reported with relevant physical consequences.

Modeling of flow problem

We assume a curved porous channel having width $2 w$ coiled a circle with center $O$ along with radius $R_{0} .$ The curved channel is filled with incompressible bi-viscosity liquid. The flow has been induced sinusoidal type channel walls deformation where impact of magnetic force is employed radially. Let $T_{0}$ and $C_{0}$ are being upper wall temperature and concentration respectively. The upper wall maintained constant temperature and concentration represented by $T_{1}$ and $C_{1},$ respectively. The flow pattern has been followed by using curvilinear coordinate system $(R, χ, Z)$ in which $R$ is specified in radial direction, $χ$ is taken the direction fluid flow while $Z$ is taken perpendicular plane associated with $R$ and $χ$ The physical problem is described by sketching Figure 1. Equations representing both walls shape are illustrated as follows^12,22,23,25:

H_{1} (χ, t) = w + a \sin ((\frac{2 π}{λ^{*}}) (χ - ct)), Upper wall,

(1)

H_{2} (χ, t) = - w - b \sin ((\frac{2 π}{λ^{*}}) (χ - ct)), Lower wall,

(2)

where $t, c, λ^{*}$ and $(a, b)$ are being time, wave speed, wavelength and amplitudes, respectively. For bi-viscosity fluid, the governing flow equations regarding heat/mass transfer phenomenon are given by^12,23,25:

\nabla \cdot U = 0, (C o n t i n u i t y E q u a t i o n)

(3)

ρ \frac{d U}{d t} = \nabla \cdot τ - \frac{μ}{k^{*}} U + J \times B, (M o m e n t u m E q u a t i o n)

(4)

ρ c_{p} \frac{d T}{d t} = k \nabla^{2} T + μ Φ + \frac{σ B^{* 2} {\tilde{R}}^{2}}{{(R + \tilde{R})}^{2}} U^{2}, (E n e r g y E q u a t i o n)

(5)

\begin{matrix} \frac{d C}{dt} = D \nabla^{2} C + \frac{D K_{T}}{T_{m}} \nabla^{2} T - k_{1} C, \\ (Mass Concentration Equation) \end{matrix}

(6)

where $τ$ notify the Cauchy stress tensor, $μ$ determine the viscosity, $β$ symbolized the bi-viscosity fluid parameter, T is for temperature, $k$ is the thermal conductivity, $c_{p}$ noted specific heat at constant pressure, $U$ is the velocity, $k^{*}$ is the permeability of porous medium, $C$ is the mass concentration, $ρ$ expresses the fluid density, $D$ mass diffusivity coefficient, $K_{T}$ is thermal diffusivity, $T_{m}$ depicts the mean temperature, $Φ$ represents the relation for dissipation function, $σ$ is the Stefan-Boltzman constant while $k_{1}$ is the rate of chemical reaction. The viscosity of a fluid is a measure of its resistance to deformation at a given rate while porosity of the medium is a measure of the void spaces in a material.

Figure 1.

Geometry of the problem.

For bi-viscosity fluid model²⁹

\bar{τ} = {\begin{matrix} (μ_{b} + \frac{P_{y}}{\sqrt{2 Π}}) {\bar{A}}_{1}, when Π > Π_{C} \\ (μ_{b} + \frac{P_{y}}{\sqrt{2 Π}}) {\bar{A}}_{1}, when Π < Π_{C} \end{matrix}

(7)

where $P_{y}$ is the yield stress of the fluid given by

P_{y} = \frac{μ_{b} \sqrt{2 Π}}{β},

(8)

in which $μ_{b}$ is the plastic dynamic viscosity, $Π$ is second invariant of ${\bar{A}}_{1}$ , $Π_{C}$ is the critical value based on non-Newtonian model. For bi-viscosity fluid flow, where $Π > Π_{C}$ , it is possible to write

μ = μ_{b} + \frac{P_{y}}{\sqrt{2 Π}} .

(9)

Using equations (8) and (9), we get stress tensor for bi-viscosity fluid flow

\bar{τ} = μ_{b} (1 + \frac{1}{β}) {\bar{A}}_{1}, when Π > Π_{C} .

(10)

In view of equation (10), equation (4) becomes²⁹

ρ \frac{d U}{dt} = μ_{b} (1 + \frac{1}{β}) \nabla \cdot {\bar{A}}_{1} - \frac{μ}{k^{*}} U + J \times B .

(11)

The form of the radial magnetic field B is given by²²

B = (\frac{B^{*} \tilde{R}}{R + \tilde{R}}) e_{R},

(12)

where $B^{*}$ is the characteristic magnetic induction in the limit $\tilde{R} \to \infty$ and $e_{R}$ is the unit vector along radial direction. It is remarked that expression for magnetic force defined in equation (12) is solenoidal and satisfies Maxwell’s equations.

Using equation (12), the term $J \times B$ in equation (4) is given by¹³

J \times B = - (\frac{σ {B^{*}}^{2} U_{2} {\tilde{R}}^{2}}{{(R + \tilde{R})}^{2}}) e_{X},

(13)

where $e_{X}$ azimuthal direction unit vector.

Assuming $U = [U_{1} (χ, R, t), U_{2} (χ, R, t), 0], T = T (χ, R, t), C = C (χ, R, t),$ equations (3) to (6) yield^12,22,23,25

\frac{\partial}{\partial R} {(R + \tilde{R}) U_{1}} + \tilde{R} \frac{\partial U_{2}}{\partial χ} = 0,

(14)

\begin{matrix} \frac{\partial U_{1}}{\partial t} + (U_{1} . \nabla) U_{1} - \frac{{U_{2}}^{2}}{R + \tilde{R}} = - \frac{1}{ρ} \frac{\partial P}{\partial R} + \frac{μ_{b}}{ρ} (1 + \frac{1}{β}) \\ [\nabla^{2} U_{1} - \frac{U_{1}}{{(R + \tilde{R})}^{2}} - \frac{2 \tilde{R}}{{(R + \tilde{R})}^{2}} \frac{\partial U_{2}}{\partial χ}] - \frac{ν}{k^{*}} U_{1}, \end{matrix}

(15)

\begin{matrix} \frac{\partial U_{2}}{\partial t} + (U_{1} . \nabla) U_{2} + \frac{U_{1} U_{2}}{R + \tilde{R}} = - \frac{\tilde{R}}{ρ (R + \tilde{R})} \frac{\partial P}{\partial χ} + \frac{μ_{b}}{ρ} (1 + \frac{1}{β}) \\ [\nabla^{2} U_{2} - \frac{U_{2}}{{(R + \tilde{R})}^{2}} + \frac{2 \tilde{R}}{{(R + \tilde{R})}^{2}} \frac{\partial U_{1}}{\partial χ}] - \frac{ν}{k^{*}} U_{2} \\ - \frac{σ {B^{*}}^{2} {\tilde{R}}^{2}}{ρ {(R + \tilde{R})}^{2}} U_{2}, \end{matrix}

(16)

ρ c_{p} [\frac{\partial T}{\partial t} + (U_{1} . \nabla) T] = k \nabla^{2} T + μ Φ + \frac{σ {B^{*}}^{2} {\tilde{R}}^{2} U_{2}^{2}}{{(R + \tilde{R})}^{2}},

(17)

\begin{matrix} \frac{\partial C}{\partial t} + (U_{1} . \nabla) C = D \nabla^{2} C + \frac{D K_{T}}{T_{m}} \\ (\frac{\partial^{2} T}{\partial R^{2}} + \frac{\tilde{R}}{R + \tilde{R}} \frac{\partial T}{\partial R} + {(\frac{\tilde{R}}{R + \tilde{R}})}^{2} \frac{\partial^{2} T}{\partial χ^{2}}) - k_{1} C, \end{matrix}

(18)

where

U_{1} . \nabla = U_{1} \frac{\partial}{\partial R} + \frac{\tilde{R} U_{2}}{R + \tilde{R}} \frac{\partial}{\partial χ},

(19)

\nabla^{2} = \frac{1}{R + \tilde{R}} \frac{\partial}{\partial R} {(R + \tilde{R}) \frac{\partial}{\partial R}} + {(\frac{\tilde{R}}{R + \tilde{R}})}^{2} \frac{\partial^{2}}{\partial χ^{2}},

(20)

\begin{array}{l} Φ = [2 {(\frac{\partial U_{1}}{\partial R})}^{2} + (\frac{\tilde{R}}{R + \tilde{R}} \frac{\partial U_{1}}{\partial χ} - \frac{U_{2}}{R + \tilde{R}}) \\ (\frac{\partial U_{2}}{\partial R} + \frac{\tilde{R}}{R + \tilde{R}} \frac{\partial U_{1}}{\partial χ} - \frac{U_{2}}{R + \tilde{R}}) + \frac{\partial U_{2}}{\partial R} \\ [\frac{\partial U_{2}}{\partial R} + \frac{\tilde{R}}{R + \tilde{R}} \frac{\partial U_{1}}{\partial χ} - \frac{U_{2}}{R + \tilde{R}}] + 2 (\frac{\tilde{R}}{R + \tilde{R}} \frac{\partial U_{2}}{\partial χ} + \frac{U_{1}}{R + \tilde{R}})] . \end{array}

(21)

The boundary conditions associated with equations (14) to (18) are²³

U_{2} = 0, U_{1} = \frac{\partial H_{1}}{\partial t}, T = T_{o}, C = C_{o} at R = H_{1} (χ, t),

(22)

U_{2} = 0, U_{1} = \frac{\partial H_{2}}{\partial t}, T = T_{1}, C = C_{1} at R = H_{2} (χ, t) .

(23)

In order to shift from fixed frame $(R, χ)$ to wave frame $(r, x)$ , we employ the transformations^22,23

\begin{matrix} x = χ - c t, r = R, p = P, u_{1} = U_{1}, u_{2} \\ = U_{2} - c, T = T . \end{matrix}

(24)

The governing equations are transformed to wave frame by using above transformations. These equations after defining the dimensionless variables^22,23

$\bar{x} = \frac{2 π}{λ^{*}} x, η = \frac{r}{a_{1}},$ ${\bar{u}}_{1} = \frac{u_{1}}{c_{1}},$ ${\bar{u}}_{2} = \frac{u_{2}}{c_{1}}, \bar{p} = \frac{2 π {a_{1}}^{2}}{λ^{*} μ_{b} c} p,$ $θ = \frac{T - T_{1}}{T_{0} - T_{1}},$ $ϕ = \frac{C - C_{1}}{C_{0} - C_{1}},$ $δ = \frac{2 π a_{1}}{λ^{*}},$ $Re = \frac{ρ c a_{1}}{μ_{b}},$ and invoking lubrication approximations ( $δ \approx 0, Re \approx 0$ ) reduce to:

\frac{\partial p}{\partial η} = 0,

(25)

\begin{matrix} - \frac{\partial p}{\partial x} + \frac{1}{γ} {- \frac{\partial}{\partial η} ((η + γ) \frac{\partial^{2} ψ}{\partial η^{2}}) - \frac{1}{η + γ} (1 - \frac{\partial ψ}{\partial η})} \\ - \frac{η + γ}{K^{*} γ} (1 - \frac{\partial ψ}{\partial η}) \\ - \frac{γ H a^{2}}{η + γ} (1 - \frac{\partial ψ}{\partial η}) = 0, \end{matrix}

(26)

\begin{matrix} \frac{\partial^{2} θ}{\partial η^{2}} + \frac{1}{(η + γ)} \frac{\partial θ}{\partial η} \\ + Br {(- \frac{1}{η + γ} (1 - \frac{\partial ψ}{\partial η}) - \frac{\partial^{2} ψ}{\partial η^{2}})}^{2} \\ + \frac{BrH a^{2} γ^{2}}{{(η + γ)}^{2}} {(1 - \frac{\partial ψ}{\partial η})}^{2} = 0, \end{matrix}

(27)

\begin{matrix} \frac{\partial^{2} ϕ}{\partial η^{2}} + \frac{1}{(γ + η)} \frac{\partial ϕ}{\partial η} - R_{c} ϕ \\ = - SrSc (\frac{\partial^{2} θ}{\partial η^{2}} + \frac{1}{(γ + η)} \frac{\partial θ}{\partial η}) + R_{c} . \end{matrix}

(28)

In above equations Re, δ, $γ$ , $K^{*}$ , and $R_{c}$ represent the Reynolds number, the wave number, the dimensionless radius of curvature, dimensionless permeability parameter and dimensionless rate of chemical reaction respectively. The stream function $ψ$ and velocity components $u_{1}$ and $u_{2}$ are related through the expressions²²

u_{1} = δ \frac{γ}{η + γ} \frac{\partial ψ}{\partial x}, u_{2} = - \frac{\partial ψ}{\partial η} .

(29)

After eliminating pressure between equations (25) and (26), we get

\begin{matrix} \frac{\partial}{\partial η} [\frac{1}{γ} (\frac{\partial}{\partial η} {(η + γ) \frac{\partial^{2} ψ}{\partial η^{2}}} + \frac{1}{η + γ} (1 - \frac{\partial ψ}{\partial η})) \\ + \frac{η + γ}{K^{*}} γ (1 - \frac{\partial ψ}{\partial η}) + \frac{γ {Ha}^{2}}{η + γ} (1 - \frac{\partial ψ}{\partial η})] = 0 . \end{matrix}

(30)

The transmuted boundary conditions (22) and (23) are

\begin{matrix} ψ = - \frac{q}{2}, \frac{\partial ψ}{\partial η} = 1, θ = 0, ϕ = 0, at η = h_{1} \\ = 1 + λ \sin x, \end{matrix}

(31)

\begin{matrix} ψ = \frac{q}{2}, \frac{\partial ψ}{\partial η} = 1, θ = 1, ϕ = 1, at η \\ = h_{2} = - 1 - λ' \sin x, \end{matrix}

(32)

where $λ = a / w$ and $λ' = b / w$ stands for amplitude ratios. We need to compute the solution of equations (27), (28), and (30) by using boundary conditions in equations (31) and (32). Following expressions are initiated for pressure per wavelength $Δ p,$ local Nusselt number at both the wall $z_{i} (i = 1, 2)$ and Sherwood number $S h_{i} (i = 1, 2)$ ^15,22,25

Δ p = \int_{0}^{2 π} \frac{dp}{dx} dx,

(33)

z_{i} = \frac{\partial h_{i}}{\partial x} {\frac{\partial θ}{\partial η} |}_{η = h_{i}}, i = 1, 2 .

(34)

Sh = \frac{\partial h_{i}}{\partial x} {\frac{\partial ϕ}{\partial η} |}_{η = h_{i}}, i = 1, 2 .

(35)

Finite difference scheme

The set of equations (27), (28), and (30) subject to boundary conditions given in equations (31) and (32) are quite complex and highly nonlinear for which analytical solution is not possible. On this end, we follow finite difference approach to simulate the numerical solution of such highly complex equations. This method successfully tacked nonlinear differential equations with excellent accuracy. This method is also sufficient to compute the solution of problem associated with curved boundaries and complex boundary conditions. The numerical computations are performed by using FORTRAN software. This method is based on following four steps:

(i) First step is to construct an iterative technique to allow the new nonlinear boundary value problem (BVP) to be transformed into a linear one at the (m + 1)th iterative stage. For this specific problem, the subsequent iterative technique is suggested:

\begin{matrix} \frac{\partial^{4} ψ^{(m + 1)}}{\partial η^{4}} + \frac{2}{(η + γ)} \frac{\partial^{3} ψ^{(m + 1)}}{\partial η^{3}} \\ - {\frac{1}{{(η + γ)}^{2}} + \frac{γ^{2} H a^{2}}{{(η + γ)}^{2}} + \frac{1}{K^{*}}} \frac{\partial^{2} ψ^{(m + 1)}}{\partial η^{2}} \end{matrix}

\begin{matrix} + {\frac{1}{{(η + γ)}^{3}} - \frac{γ^{2} H a^{2}}{{(η + γ)}^{3}} - \frac{1}{K^{*} (η + γ)}} \\ \frac{\partial ψ^{(m + 1)}}{\partial η} - \frac{1}{{(η + γ)}^{3}} + \frac{γ^{2} H a^{2}}{{(η + γ)}^{3}} + \frac{1}{K^{*} (η + γ)} = 0 . \end{matrix}

(36)

\begin{matrix} \frac{\partial^{2} θ^{(m + 1)}}{\partial η^{2}} + \frac{1}{γ + η} \frac{\partial θ^{(m + 1)}}{\partial η} \\ = - Br {(- \frac{1}{γ + η} (1 - \frac{\partial ψ^{(m)}}{\partial η}) - \frac{\partial^{2} ψ^{(m)}}{\partial η^{2}})}^{2} \\ - \frac{BrH a^{2} γ^{2}}{{(γ + η)}^{2}} {(1 - \frac{\partial ψ^{(m)}}{\partial η})}^{2} \end{matrix}

(37)

\begin{matrix} \frac{\partial^{2} ϕ^{(m + 1)}}{\partial η^{2}} + \frac{1}{η + γ} \frac{\partial ϕ^{(m + 1)}}{\partial η} - R_{c} ϕ^{(m + 1)} \\ = - SrSc (\frac{\partial^{2} θ^{(m)}}{\partial η^{2}} + \frac{1}{(γ + η)} \frac{\partial θ^{(m)}}{\partial η}) + R_{c}, \end{matrix}

(38)

\begin{matrix} ψ^{m + 1} = - \frac{q}{2}, \frac{\partial ψ^{m + 1}}{\partial η} = 1, θ^{(m + 1)} = 0, ϕ^{(m + 1)} = 0, \\ at η = h_{1}, \end{matrix}

(39)

\begin{matrix} ψ^{(m + 1)} = \frac{q}{2}, \frac{\partial ψ^{(m + 1)}}{\partial η} = 1, θ^{(m + 1)} = 1, ϕ^{(m + 1)} = 1, \\ at η = h_{2} . \end{matrix}

(40)

Here, the index (m) shows the iterative step. It is now clear that above BVP is linear in $ψ^{(m + 1)}$ .

(ii) In this step, we insert finite difference approximations of $ψ^{(m + 1)}$ , $θ^{(m + 1)}$ , $ϕ^{(m + 1)}$ and their derivatives into equations (36) to (38). In this way, we get a system of linear algebraic equations at each iterative step.

(iii) In third step, the system of algebraic equations obtained in previous step are solved at each cross-section to get numerical results of $ψ^{(m + 1)}$ , $θ^{(m + 1)}$ , and $ϕ^{(m + 1)}$ . Evidently, suitable initial guesses are required for $ψ^{(m)}$ , $θ^{(m)}$ , and $ϕ^{(m)}$ at each cross-section to start the iterative procedure. The iterative procedure at cross-section is carried out until a convergent solution is reached.

(iv) In this step, we employed method of successive under-relaxation to get rapid convergent solution. In this method the values of ${\tilde{ψ}}^{(m + 1)},$ ${\tilde{θ}}^{(m + 1)}$ , and ${\tilde{ϕ}}^{(m + 1)},$ at (m+1)th iterative step are refined as

ψ^{(m + 1)} = ψ^{(m)} + τ ({\tilde{ψ}}^{(m + 1)} - ψ^{(m)}),

(41)

θ^{(m + 1)} = θ^{(m)} + τ ({\tilde{θ}}^{(m + 1)} - θ^{(m)}),

(42)

ϕ^{(m + 1)} = ϕ^{(m)} + τ ({\tilde{ϕ}}^{(m + 1)} - ϕ^{(m)}),

(43)

where $τ$ is under relaxation parameter.

Usually, $τ$ is chosen small for rapid convergence. In present computation, the iterative procedure is terminated after achieving the values of $ψ, θ$ , and $ϕ$ convergent to $10^{- 8} .$ Some authors have already been used for instance, Wang et al.,¹⁵ Hayat et al.,¹⁹ Ali et al.,²³ and Abbas et al.²⁷ The complete iterative process in represented in Figure 2.

Figure 2.

The flow chart of the numerical method.

Verification of solution

The numerical solution computed in previous section is compared for velocity profile with work of Ali et al.²² as a limiting case $(K^{*} \to \infty, Ha = 0, γ = 2)$ . On this end, Figure 3 is prepared for which superimposed line show the solution reported by Ali et al.²² while solid line represents the current numerical solution. An excellent accuracy between solutions is noted.

Figure 3.

Comparison of solution for velocity: superimposed line denotes the work of Ali et al.¹¹ while solid line represents numerical solution.

Results and discussion

In this section, it is intended to analyze the important features of peristaltic motion such as velocity, pressure rise, temperature distribution, mass concentration and trapping for various values of the curvature parameter $(γ)$ , bi-viscosity parameter $(β),$ Brinkman number $(Br)$ , rate of chemical action ( $R_{c}$ ), permeability parameter $(K^{*})$ , and Hartmann number $(H a) .$ The heat transfer coefficient at upper wall is also discussed. To this end, Figures 4 –28 are plotted.

Figure 4.

$u_{2} (η)$ for $β$ with $γ$ = 2.5, $K^{*}$ = 0.2, Ha = 0.5, $λ = 0.4,$ and $Θ$ = 1.5.

Figure 5.

$u_{2} (η)$ for $K^{*}$ with $γ$ = 2.5, $β$ = 0.2, Ha = 0.2, $λ = 0.4,$ and $Θ$ = 1.5.

Figure 6.

$u_{2} (η)$ for $Ha$ with $γ$ = 2.5, $K^{*}$ = 0.2, $β$ = 0.2, $λ = 0.4,$ and $Θ$ = 1.5.

Figure 7.

$u_{2} (η)$ for $γ$ with $β$ = 0.2, $K^{*}$ = 0.2, Ha = 0.5, $λ = 0.4,$ and $Θ$ = 1.5.

Figure 8.

$Δ p$ for $β$ with $γ$ = 2.5, $K^{*}$ = 0.2, Ha =0.5, $λ = 0.4,$ and $Θ$ = 1.5.

Figure 9.

$Δ p$ for $K^{*}$ with $γ$ = 2.5, Ha = 0.5, $β$ = 0.2, $λ = 0.4,$ and $Θ$ = 1.5.

Figure 10.

$Δ p$ for $Ha$ with $γ$ = 2.5, $K^{*}$ = 0.2, $β$ = 0.2, $λ = 0.4,$ and $Θ$ = 1.5.

Figure 11.

$Δ p$ for $γ$ with $K^{*}$ = 0.2, $β$ = 0.2, Ha = 0.5, $λ = 0.4,$ and $Θ$ = 1.5.

Figure 12.

$θ (η)$ for $β$ with Br = 0.2, $γ$ = 2, $K^{*}$ = 1, Ha = 0.5, $λ = 0.4,$ and $Θ$ = 1.5.

Figure 13.

$θ (η)$ for Br with $β$ =0.2, $γ$ = 2, $K^{*}$ = 1.5, Ha = 0.5, $λ = 0.4,$ and $Θ$ = 1.5.

Figure 14.

$θ (η)$ for $K^{*}$ with Br = 0.2, $γ$ = 2, Ha = 0.5, $β$ = 0.2, $λ = 0.4,$ and $Θ$ = 1.5.

Figure 15.

$θ (η)$ for Ha with Br = 0.2, $γ$ = 2, $K^{*}$ = 1.5, $β$ = 0.2, $λ = 0.4,$ and $Θ$ = 1.5.

Figure 16.

z at upper wall for $β$ with Br = 0.2, $γ$ = 2, $K^{*}$ = 1, Ha = 0.5, $λ = 0.4,$ and $Θ$ = 1.5.

Figure 17.

z at upper wall for Br with $β$ = 0.2, $γ$ = 2, $K^{*}$ = 1.5, Ha = 0.5, $λ = 0.4,$ and $Θ$ = 1.5.

Figure 18.

z at upper wall $K^{*}$ with Br = 0.2, $γ$ = 2, Ha = 0.5, $β$ = 0.2, $λ = 0.4,$ and $Θ$ = 1.5.

Figure 19.

z at upper wall for $Ha$ with Br = 0.2, $γ$ = 2, $K^{*}$ = 1.5, $β$ = 0.2, $λ = 0.4,$ and $Θ$ = 1.5.

Figure 20.

ϕ $(η)$ for $β$ with Br = 2, $γ$ = 2, $R_{c} = 0.2,$ Sr = 1.5, Sc = 1.2, $λ = 0.4,$ and $Θ$ = 1.5.

Figure 21.

ϕ $(η)$ for Br with $β$ = 0.4, $γ$ = 2, $R_{c} = 0.2,$ Sr = 1.5, Sc = 1.2, $λ = 0.4,$ and $Θ$ = 1.5.

Figure 22.

ϕ $(η)$ for $R_{c}$ with $β$ = 0.4, Br = 2, $γ$ = 2, Sr = 1.5, Sc = 1.2, $λ = 0.4,$ and $Θ$ = 1.5.

Figure 23.

ϕ $(η)$ for Ha with $β$ = 0.4, Br = 2, $γ$ = 2, $R_{c}$ = 0.2, Sr = 1.5, Sc = 1.2, $λ = 0.4,$ and $Θ$ = 1.5.

Figure 24.

ϕ $(η)$ for Sc with $β$ = 0.4, Br = 2, $γ$ = 2, $R_{c} = 0.2,$ Sr = 1.5, $λ = 0.4,$ and $Θ$ = 1.5.

Figure 25.

ϕ $(η)$ for Sr with $β$ = 0.4, Br = 2, $γ$ = 2, $R_{c} = 0.2,$ Sc = 1.2, $λ = 0.4,$ and $Θ$ = 1.5.

Figure 26.

Streamlines in wave frame for (a) $β = 0.1,$ (b) $β = 0.5,$ (c) $β = 1$ and (d) $β = 1.5$ . The other parameters chosen are $γ$ = 2, $K^{*}$ = 1, Ha = 0.5, $λ = 0.4,$ and $Θ$ = 1.5.

Figure 27.

Streamlines in wave frame for (a) $K^{*} = 0.1,$ (b) $K^{*} = 0.15,$ (c) $K^{*} = 0.2,$ and (d) $K^{*} = 0.3 .$ The other parameters chosen are $γ$ = 2, Ha = 0.5, $β = 0.5,$ $λ = 0.4,$ and $Θ$ = 1.5.

Figure 28.

Streamlines in wave frame for (a) $Ha = 1,$ (b) $Ha = 5,$ (c) $Ha = 10,$ and (d) $Ha = 15 .$ The other parameters chosen are $γ$ = 2, $K^{*}$ = 0.2, $β = 0.5,$ $λ = 0.4,$ and $θ$ = 1.5.

The axial velocity distribution for some specific values of bi-viscosity parameter $(β),$ permeability parameter $(K^{*})$ , Hartmann number $(H a)$ , and curvature parameter $(γ)$ is shown in Figures 4 –7, respectively. Figure 4 depicts that velocity profile increases by increasing bi-viscosity parameter. The truncation in velocity component $u_{2}$ for permeability parameter $(K^{*})$ is verified by illustrating Figure 5. The weaker porous medium consequences are deliberated due to lower $(K^{*})$ values in contrast to higher values which denote the case when resistance corresponds to porous matrix is more dominant. Moreover, the amplitude of velocity due to variation of $(K^{*})$ slow down the velocity amplitude the velocity is shift from lower channel wall to central line $(η = 0) .$ Figure 6 exposes the consequences of applied magnetic field on the flow velocity. Here it is quite obvious that flow velocity exhibits boundary layer character for larger values of Hartmann number. In fact for large values of Hartmann number, the disturbance in flow velocity is confined in thin layer near both upper and lower walls. The fluid outside the boundary layer moves with a velocity which varies linearly with radial coordinate $η$ . In contrast, the fluid outside the boundary layers moves with constant velocity in a straight channel. The consequences of radius of curvature $(γ)$ on velocity have been predicted in Figure 7. The velocity distribution is observed to be asymmetric about $η = 0$ associated with lower values of $γ .$ It is further noticed that for larger values of $γ,$ velocity recovers its symmetry about $η = 0 .$ This is expected because for large values of $γ$ channel becomes straight. The graphical visualization of pressure rise per wavelength against diverse flow parameters has been reported in Figures 8 –11. Figure 8 shows the effects of bi-viscosity parameter $β$ on $Δ p$ . It is observed that $Δ p$ is decreasing by increasing $β$ in the peristaltic pumping region. The profile of pressure rise per wavelength for different values of $K^{*}$ (permeability parameter) is shown in Figure 9. A declining trend in $Δ p$ due to enhanced porous medium permeability has been seen in the region of peristaltic pumping $(Θ > 0, Δ p > 0) .$ Thus, the normal fluid flow inside the channel impedes due to porous medium and peristalsis play a valuable role against the pressure rise to maintain the same flux as in the case of clear medium inside the channel. This fact reduced the pumping efficiency subsequently. Figures 10 and 11 are plotted to show the variation of $Δ p$ against dimensionless mean ﬂow rate $Θ$ for various values of Ha and $γ$ , respectively. The maximum variation in Ha and $γ$ , an improved pressure rise per wavelength has been resulted in the pumping region. However, for case of free pumping ( $Δ p$ = 0) and co-pumping case $(Θ > 0, Δ p < 0),$ this situation is quite different. A lower trend in $Δ p$ is observed due to Ha and $γ$ . Figure 11 revealed that change in $Δ p$ is maximum in case of straight channel as compared to curved configuration up to a certain values of $Θ$ . However, opposite observations are referred above this critical numerical value. Similar reverse styles encountered in region of co-pumping and free pumping. The alteration in radial temperature distribution inside the channel against flow parameters like bi-viscosity parameter $(β),$ permeability parameter $(K^{*})$ Hartmann number (Ha) and Brinkmann number $(Br)$ are executed in Figures 12 –15. Figures 12, 13, and 15, respectively, report that $θ (η)$ increases by increasing $β, Br$ and Ha inside the channel. In contrast, Figure 14 shows a decrease in temperature with increasing permeability parameter. It is due to the fact that heat transportation appeared from the walls to the fluid is reduced for clear medium and as a result, a decrease in temperature is noted inside the channel. Figures 16 –19 indicate change in heat transfer coefficient $z$ at the upper channel wall for $K^{*}, β, Br,$ and Ha. The distribution of $z$ oscillates periodically which is associated with the oscillatory nature of wall surfaces. The results reveal that a relatively improved oscillatory profile is referred to the maximum values of $β, Br$ and Ha. However, a declining periodically oscillating phenomenon is evaluated with enlarge values of $K^{*} .$

The mass concentration variation inside the channel for several values of bi-viscosity parameter, Brinkmann number, rate of chemical reaction, Hartmann number, Schmidt number and Soret number are shown in Figures 20 –25. It is observed that mass concentration inside the channel increases with increasing each of these parameters. The streamlines of flow inside the channel for different values of bi-viscosity parameter $(β),$ permeability parameter $(K^{*})$ , Hartmann number (Ha) and curvature parameter $(γ)$ are shown in Figures 26 –29. The aim of these curves is to highlight the trapping phenomenon where fluid volume known as bolus is trapped within closed streamlines. Figure 26 shows two nearly symmetric circulating rolls in lower and upper halves of the channel. However, lower roll increases in size while the upper one shrinks with increasing $β .$ Figure 27 illustrates the effects of $K^{*}$ on trapping phenomenon. For smaller values of $K^{*}$ i.e. for strong permeability effects the circulating roll in lower half contains two small eddies. The size of lower half is greater than the upper one. With increasing $K^{*},$ the two eddies inside the lower roll merge into a single cell. Moreover, the size of upper roll reduces with increasing $K^{*}$ . The formation of two circulating eddies in the lower roll for small values of $K^{*}$ is never reported in the earlier available literature on the peristaltic flows through porous-saturated curved channel. Figure 28 demonstrates the effects of applied magnetic field on trapping phenomenon. It is observed that the role of applied magnetic field is to increase the size of circulating roll. The streamlines highlighting the effects of dimensionless radius of curvature on trapping has been conveyed through Figure 29. The graphical analysis revealed that both circulating roll become in size with increasing $γ .$

Figure 29.

Streamlines in wave frame for (a) $γ = 2,$ (b) $γ = 3.5$ and (c) $γ \to \infty .$ The other parameters chosen are Ha = 0.5, $K^{*}$ = 0.2, $β = 0.5,$ $λ = 0.4,$ and $θ$ = 1.5.

Concluding remarks

The peristaltic flow of bi-viscosity fluid in a curved porous channel with extended and contracting walls has been studied under the assumptions of long wavelength. The channel with sinusoidally deformed walls is assumed to be porous. Rheological non-Newtonian material characteristics are evaluated by using bi-viscosity fluid model. The heat and mass transfer features are examined in presence of Joule heating effects. The numerical simulations based on finite difference scheme are performed with desirable accuracy. The novel characteristics of pressure, velocity, temperature, concentration and streamlines are graphically underlined. This analysis predicts following interesting observations:

The flow velocity is observed up to maximum level with variation of $β .$

Pressure rise per wavelength decreases with increasing $β$ and $K^{*} .$

A declining temperature variation has been observed due to impact of $K^{*} .$ However, it increases with increasing $β .$

The mass concentration increases with increasing $R_{c}$ and $β .$

The size of circulating roll in lower half of the channel boosts up with larger variation of $β$ and $K^{*} .$

The symmetry in velocity and streamlines pattern is observed when $γ \to \infty$ .

The results for Newtonian fluid are achieved for larger values of $β .$

Supplemental Material

Appendix – Supplemental material for Heat and mass transfer characteristics in flow of bi-viscosity fluid through a curved channel with contracting and expanding walls: A finite difference approach

Supplemental material, Appendix for Heat and mass transfer characteristics in flow of bi-viscosity fluid through a curved channel with contracting and expanding walls: A finite difference approach by Raheel Ahmed, Nasir Ali, Sami Ullah Khan, Ali Chamkha and Iskander Tlili in Advances in Mechanical Engineering

Footnotes

Appendix

Handling Editor: James Baldwin

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) received no financial support for the research, authorship, and/or publication of this article.

Supplemental material

Supplemental material for this article is available online.

References

Shapiro

Jaffrin

Weinberg

SL.

Peristaltic pumping with long wavelength at low Reynolds number. J Fluid Mech 1969; 37: 799–825.

Raju

Devanathan

Peristaltic motion of non-Newtonian fluid. Acta 1972; 11: 170–178.

Böhme

Friedrich

Peristaltic flow of viscoelastic liquids. J Fluid Mech 1983; 128: 109–122.

Srivastava

VP.

Interaction of peristaltic flow with pulsatile flow in a circular cylindrical tube. J Biomech 1985; 18: 247–253.

Siddiqui

Provost

Schwarz

WH.

Peristaltic pumping of second order fluid in a planar channel. Rheol Acta 1991; 30: 249–262.

Hayat

Wang

Siddiqui

, et al. Peristaltic transport of third order fluid in a cylindrical tube. Math Models Methods Appl Sci 2002; 12: 169–176.

Hayat

Wang

Siddiqui

, et al. Peristaltic motion of Johnson-Segalman fluid in a planar channel. Math Problems Eng 2003; 1: 1–23.

Haroun

MH.

Nonlinear peristaltic flow of a fourth grade fluid in an inclined asymmetric channel. Comput Mater Sci 2007; 39: 324–333.

Böhme

Müller

Analysis of non-Newtonian effects in Peristaltic pumping. J Non-Newtonian Fluid Mech 2013; 201: 107–119.

10.

Sato

Kawai

Fujita

, et al. Two dimensional peristaltic flow in curved channels. Jpn Soc Mech Eng 2000; 66: 679–685.

11.

Ali

Sajid

Hayat

Long wavelength ﬂow analysis in a curved channel. Z Naturforsch 2010; 65a: 191–196.

12.

Ali

Javid

Sajid

, et al. Numerical simulation of peristaltic flow of a biorheological fluid with shear-dependent viscosity in a curved channel. Comput Methods Biomech Biomed Eng 2016; 19: 614–627.

13.

Kalantari

Sadeghy

Sadeqi

Peristaltic flow of non-Newtonian fluids through curved channels: a numerical study. Annu Trans Nordic Rheol Soc 2013; 21: 163–170.

14.

Ramanamurthy

Prasad

Narla

VK.

Unsteady peristaltic transport in curved channels. Phys Fluids 2013; 25: 091903–091920.

15.

Wang

Hayat

Ali

, et al. Magnetohydrodynamic peristaltic motion of a Sisko fluid in a symmetric or asymmetric channel. Physica A 2008; 387: 347–362.

16.

Tripathi

Bég

OA.

A study of unsteady physiological magneto-fluid flow and heat transfer through a finite length channel by peristaltic pumping. Proc Inst Mech Eng Part H J Eng Med 2012; 226: 631–644.

17.

Bandopadhyay

Tripathi

Chakraborty

Electroosmosis-modulated peristaltic transport in microfluidic channel. Phys Fluids 2016; 28: 052002.

18.

Ali

Wang

Hayat

, et al. Slip effects on the peristaltic flow of a third grade fluid in a circular cylindrical tube. J Appl Mech 2009; 76: 1–10.

19.

Hayat

Hina

Ali

Simultaneous effects of slip and heat transfer on the peristaltic flow. Commun Nonlinear Sci Numer Simulat 2010; 15: 1526–1537.

20.

Hayat

Hina

The influence of wall properties on the MHD peristaltic flow of a Maxwell fluid with heat and mass transfer. Nonlinear Anal Real World Appl 2010; 11: 3155–3169.

21.

Mekheimer

Abd Elmaboud

The influence of heat transfer and magnetic field on peristaltic transport of a Newtonian fluid in a vertical annulus: application of endoscope. Phys Lett A 2008; 372: 1657–1665.

22.

Ali

Sajid

Abbas

, et al. Non-Newtonian fluid-flow induced by peristaltic waves in a curved channel. Eur J Mech B Fluids 2010; 29: 384–397.

23.

Ali

Sajid

Javed

, et al. Heat transfer analysis of peristaltic flow in a curved channel. Int J Heat Mass Transf 2010; 53: 3319–3325.

24.

Srinivas

Gayathri

Kothandapani

Mixed convective heat and mass transfer in an asymmetric channel with peristalsis. Commun Nonlinear Sci Numer Simul 2011; 16: 1845–1862.

25.

Hayat

Hina

Hendi

, et al. Effect of wall properties on the peristaltic ﬂow of a third grade ﬂuid in a curved channel with heat and mass transfer. International Journal of Heat and Mass Transf 2011; 54: 5126–5136.

26.

Tripathi

Bég

OA.

Mathematical modelling of heat transfer effects on swallowing dynamics of viscoelastic flood bolus through the human esophagus. Int J Thermal Sci 2013; 70: 41–53.

27.

Abbas

Wang

Hayat

Slip effects and heat transfer analysis in a viscous fluid over an oscillatory stretching surface. Int J Numer Methods Fluids 2009; 59: 443–458.

28.

Tripathi

Bég

OA.

A study on peristaltic flow of nanofluids: application in drug delivery systems. Int J Heat Mass Transf 2014; 70: 61–70.

29.

Dash

Mehta

Jayaraman

Effect of yield stress on the flow of a Casson fluid in a homogeneous porous medium bounded by a circular tube. Appl Sci Res 1996; 57: 133–149.

30.

Hayat

Khan

Waqas

, et al. Viscous dissipation effect in flow of magneto nanofluid with variable properties. J Mol Liq 2016; 222: 47–54.

31.

Shehzad

Hayat

Qasim

, et al. Effects of mass transfer on MHD flow of Casson fluid with chemical reaction and suction. Braz J Chem Eng 2013; 30: pp. 187–195.

32.

Reddy

Rao

Surender

Soret, Joule heating and Hall effects on free convection in a Casson fluid saturated porous medium in a vertical channel in the presence of viscous dissipation. Proc Eng 2015; 127: 1219–1226.

33.

Makanda

Shaw

Sibanda

Diffusion of chemically reactive species in Casson fluid flow over an unsteady stretching surface in porous medium in the presence of a magnetic field. Math Problems Eng 2015; 2015: 10 pp.

34.

Tripathi

Akbar

Khan

, et al. Peristaltic transport of bi-viscosity fluids through a curved tube: a mathematical model for intestinal flow. Proc Inst Mech Eng Part H J Eng Med 2016; 230: 817–828.

35.

Tamoor

Waqas

Khan

, et al. Magnetohydrodynamic flow of Casson fluid over a stretching cylinder. Results Phys 2017; 7: 498–502.

36.

Brinkman

HC.

A calculation of the viscous force exerted by a ﬂowing ﬂuid on a dense swarm of particles. Appl Sci Res 1949; 1: 27–34.

37.

El Shehawey

. Peristaltic transport through a porous medium. Appl Math Comput 2001. DOI: 10.1016/S0096-3003(01)00236-3

38.

Mekheimer

KS.

Non-linear peristaltic transport through a porous medium in an inclined planar channel. J Porous Media 2003; 6: 189–202.

39.

Kothandapani

Srinivas

On the influence of wall properties in the MHD peristaltic transport with heat transfer and porous medium. Phys Lett A 2008; 372: 4586–4591.

40.

Hayat

Umar Qureshi

Hussain

Effect of heat transfer on peristaltic flow of an electrically conducting fluid in a porous space. Appl Math Model 2009; 33: 1862–1873.

41.

Javed

Mehmood

Siddiqui

AM.

Mixed convection in a triangular cavity permeated with micropolar nanofluid-saturated porous medium under the impact of MHD. Braz Soc Mech Sci Eng 2017; 39: 3897–3909.

42.

Javed

Mehmood

Siddiqui

, et al. Study of heat transfer in water-Cu nanofluid saturated porous medium through two entrapped trapezoidal cavities under the influence of magnetic field. J Mol Liq 2017; 240: 402–411.

43.

Javed

Mehmood

Abbas

Natural convection in square cavity filled with ferrofluid saturated porous medium in the presence of uniform magnetic field. Physica B 2017; 506; 122–132.

44.

Mehmood

Javed

Pop

MHD-mixed convection flow in a lid driven trapezoidal cavity under uniformly/non-uniformly heated bottom wall. Int J Numer Methods Heat Fluid Flow 2017; 27; 1231–1248.

45.

Mehmood

Numerical simulations and linear stability analysis of mixed thermomagnetic convection through two lid-driven entrapped trapezoidal cavities enclosing ferrofluid saturated porous medium. Commun Heat Mass Transfer 2019; 109: 104345.

46.

Mehmood

Javed

Impact of MHD on energy transmission through ferrofluid saturating porous medium contained in a lid-driven trapezoidal container with triangular notched heater configuration. Can J Phys 2019; 97: 1104–1114.

47.

Belhocine

Numerical study of heat transfer in fully developed laminar flow inside a circular tube. Int J Adv Manuf Technol 2016; 85: 2681–2692.

48.

Belhocine

Omar

WZW

. An analytical method for solving exact solutions of the convective heat transfer in fully developed laminar flow through a circular tube. Heat Transfer Asian Res 2017; 46: 1342–1353.

49.

Belhocine

Abdullah

OI.

Numerical simulation of thermally developing turbulent flow through a cylindrical tube. Int J Adv Manuf Technol 2019; 102: 2001–2012.

50.

Belhocine

Abdullah

OI.

Similarity and numerical analysis of the generalized Levèque problem to predict the thermal boundary layer. Int J Interact Design Manuf 2018; 12: 1015–1025.

51.

Hayat

Rafiq

Alsaedi

, et al. Radiative and Joule heating effects on peristaltic transport of dusty fluid in a channel with wall properties. Eur Phys J Plus 2014; 129: 225.

52.

Hayat

Abbasi

Al-Yami

, et al. Slip and Joule heating effects in mixed convection peristaltic transport of nanofluid with Soret and Dufour effects. J Mol Liq 2014; 194: 93–99.

53.

Reddy

KV.

Influence of Joule heating on MHD peristaltic flow of a nanofluid with compliant walls. Process Eng 2015; 127: 1002–1009.

54.

Hayat

Nisar

Ahmad

, et al. Simultaneous effects of slip and wall properties on MHD peristaltic motion of nanofluid with Joule heating. J Magnet Magnet Mater 2015; 395: 48–58.

55.

Abbasi

Hayat

Alsaedi

Effects of inclined magnetic field and Joule heating in mixed convective peristaltic transport of non-Newtonian fluids. Bull Polish Acad Sci Techn Sci 2015; 63: 501–514.

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

0.00 MB

0.07 MB