Sage Journals: Discover world-class research

Abstract

This paper employs to discuss the cilia induced flow of a compressible Jeffrey fluid with variable characteristic of thermal conductivity. Thermal conductivity is thermally variable. Effects of thermal radiation and MHD are also explored. The velocity slip at the wall is also considered. Flow equations have been modeled mathematically and tackled analytically using regular perturbation method. A proper attention is given to plot and discuss the behavior of quantities of interest against the relevant parameters. We explore growth of temperature with Mach number but this growth is strongly visible in supersonic regime as compared to the subsonic regime. The discoveries of the current model have noteworthy yields, which can be pertinent in various areas of aircraft industry, geophysics, biological transport processes, and in the operation of other mechanical devices.

Keywords

Cilia compressible fluid peristalsis variable thermal conductivity radiation

Introduction

Cilia are small, slim, and tiny hair-like organelles exists at the outline of certain cells. The significant capacities performed by cilia include movement, alimentation, respiration, and tactile capacities. They adopt a significant part in cell cycle and furthermore in the progress of humans and animals. Two kinds of cilia comprises motile and non-motile cilia. Motile cilia in human respiratory parcel clears the organic fluid and bring out residue from lungs. The latter one is essential cilia and act as sensory organelles. They participate in signaling between cells. For instance, they help in legitimate urine stream by flagging the kidney cells. Additionally they empower the exchange of significant particles from one side of the light-sensitive cells to another in the retina. Ciliary-prompted peristalsis shows up altogether in different organic transport cycles, for example; in biomedicine, physiology and atomic reactors. Theoretical study on ciliated structures can be found in the literature.^1,2 Theory of propulsion of micro-organism by cilia is directed by Blake.³ A suitable modeling approach for multi-cilia configurations is analyzed by Gueron and Liron.⁴ They presented a general method for the investigation of metachronism; a process in which swing of cilia tips form a travelling waves. The respiratory tract and the esophagus are genuine examples of peristaltic stream in a ciliated tube. Radhakarishnamacharya and Sharma.⁵ deliberated the peristaltic transport of a self-assertive micro-organism settle in a channel containing a viscous fluid. Maqbool et al.⁶ explored MHD Jeffrey fluid flow persuaded by cilia in an inclined tube. They adopted lubrication approach and observed an increment in the amplitude of pressure gradient with cilia length. MHD effect on metachronal beating of cilia for Casson fluid model is examined by Akbar and Khan.⁷ A direct relation of velocity profile with Hartman number is perceived near the channel walls while reverse relation is perceived at the center of channel. Krishna Kumari et al.⁸ carried out combined impact of peristalsis and cilia on motion of MHD Micropolar fluid. Very recently Hasen and Abdulhadi⁹ studied MHD effect on peristaltic transport for Rabinowitsch fluid in cilia channel with porous medium. An oscillatory nature is investigated for pressure gradient as well as heat transfer. While velocity profile is parabolic in nature with greatest incentive at the center of channel.

Compressible fluids have caught the eyes of researchers as a result of their common applications in the development of planes, airplane channels, vehicles at high speed, and supersonic air streams. Existing literature shows that only few attempts explore the impact of MHD on peristaltic flows of compressible fluids. Basic studies on the solution of equations of compressible viscous fluids are addressed in Maccormack¹⁰ and Guerra and Gustafsson.¹¹ In both studies solutions are obtained numerically. Mekhmier et al.¹² analyzed the impression of MHD on compressible Maxwell fluid prompted by a progressive wave in a porous microchannel and examined a lower stream rate for slip fluid comparative with the non-slip fluid. Abdelsalam and Vafai¹³ explored the consequence of magnetic field on compressible Jeffrey fluid flowing in a microfluidic channel under peristalsis. Some more attempts regarding compressible fluid flow have been provided in the studies.^14–16 Physiological progression of biomedical compressible fluids in a symmetric ciliated channel is examined by Saleem et al.¹⁷ They explored a declination in Jeffrey fluid fluctuation as it switches from hydrodynamic to hydromagnetic fluid therefore retardation time diminishes.

Thermal radiation is one of the imperative variables controlling the heat transfer. Heat transfer alongside radiation has gotten a subject of various exploration because of its noteworthy contributions in atomic plants, aircrafts, gas turbines and so forth. The thermal radiation ends up being more noteworthy whenever high temperature zone is taken into account. Very recently Khan and Rafaqat¹⁸ explored the peristaltic stream of compressible Jeffrey fluid with the impact of MHD and radiation. They observed that Reynolds number and radiation number are responsible for increment of temperature of the fluid while Prandtl number indicates a decrement in temperature. We can find applications of this problem in fluid transport mechanism specifically in biological systems. Since nature of biological fluids is electrically conducting and cilia are attached on the lining of lungs which protects the lungs from harmful substances (dust particles, bacteria etc.) by a highly viscous fluid containing various proteins showing the relevance of cilia prompted flow of MHD compressible Jeffrey fluid in the respiratory tract of humans.

All these examinations are done assessing consistent physical properties of the fluid but variable attributes of physical properties commanded for practical conditions. The study of heat transfer with variable thermal conductivity is significant in dealing with physical problems. Hussain et al.¹⁹ deliberated the heat transfer of MHD Jeffrey fluid with peristalsis and variable thermal conductivity. Hayat et al.²⁰ considered 3D stretched stream of Jeffrey fluid with radiation and varying thermal conductivity. They investigated reduction in temperature with increasing Prandtl number. Some ongoing researches concerning the effect of heat transfer along with the variable thermal conductivity are quoted in Refs.^21–24 In all the above mentioned attempts thermal conductivity is presumed to be linearly varying function of temperature. Khan et al.²⁵ inspects heat and mass transfer on a third-grade MHD fluid with chemical reaction and variable reactive index. Analysis of heat transfer of squeezing unsteady nanofluid with the impact of inclined magnetic field and varying thermal conductivity is explored by Lahmar et al.²⁶ They realized a decrement in heat transfer due to the existence of varying thermal conductivity. Further exploration about heat transfer with varying thermal conductivity could be understood through.^27–29

As far as author could possibly know thermal effects on compressible fluid flow with variable fluid properties have not been done in the existing literature. In this paper transport characteristics of 2D flow of MHD compressible Jeffrey fluid inside a ciliated channel is deliberated. The main novelty of this work is to intend the combined influence of heat transfer and thermal radiation on cilia driven transport when thermal conductivity is thermally variable. This study is sorted out as follows: the physical detailing of this problem is presented in section 2. The perturbation technique used to comprehend our system of equations is given in section 3. The numerical results are interpreted via graphs in section 4. In the last section, all results are summarized.

Problem formulation

In order to designate the physical model of present study mathematically, we consider a symmetric channel of uniform width 2d, with $(x, y)$ as coordinate axis, where x-axis shapes the central axis of the channel and y-axis transverse to it as appeared in Figure 1. A 2D flow of an electrically conducting compressible Jeffrey fluid is occurred in the ciliated channel. Cilia are available in the physical model as we considered a ciliated channel and are hinged at the inner walls of channel and their swing caused the fluid motion. Channel walls are ciliated with metachronal wave patterns; produced as the result of mutual pulsing of cilia with constant speed $c$ along the channel walls. The schematic illustration of ciliated channel and metachronal wave pattern has been given in Figure 1(a) and (b). The consequences of induced magnetic field are ignored while external magnetic field of strength $B_{0}$ is applied in $y$ - direction. The bottom and top walls of the channel are kept at temperature $T_{0}$ and $T_{1}$ , respectively. Thermal radiation and viscous dissipation in heat transfer are considered. Also, thermal conductivity is thermally variable. Geometry of flow domain and its physical parameters has been shown in Figure 1(c).

Figure 1

(a) Ciliated channel, (b) metachronal wave pattern, and (c) illustrative showing flow domain and its physical parameters.

The density in terms of pressure is given as,

\frac{1}{ρ} \frac{\partial ρ}{\partial p} = k_{c},

(1)

with the solution,

ρ = ρ_{0} \exp (k_{c} (p - p_{c})),

(2)

where the reference pressure and reference density are $p_{c}$ and $ρ_{0}$ , respectively. $k_{c}$ is the fluid compressibility.

The governing equations of the flow are¹⁸

ρ (\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y}) + \frac{\partial ρ}{\partial t} + u \frac{\partial ρ}{\partial x} + v \frac{\partial ρ}{\partial y} = 0,

(3)

(1 + λ_{1} \frac{\partial}{\partial t}) [ρ (\frac{\partial u}{\partial t} + u \frac{\partial u}{\partial x} + v \frac{\partial u}{\partial y})] = - (1 + λ_{1} \frac{\partial}{\partial t}) \frac{\partial p}{\partial x} + μ (1 + λ_{2} \frac{\partial}{\partial t}) [\nabla^{2} u + \frac{1}{3} \frac{\partial}{\partial x} (\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y})] - σ {B_{0}}^{2} (1 + λ_{1} \frac{\partial}{\partial t}) u,

(4)

(1 + λ_{1} \frac{\partial}{\partial t}) [ρ (\frac{\partial v}{\partial t} + u \frac{\partial v}{\partial x} + v \frac{\partial v}{\partial y})] = - (1 + λ_{1} \frac{\partial}{\partial t}) \frac{\partial p}{\partial y} + μ (1 + λ_{2} \frac{\partial}{\partial t}) [\nabla^{2} v + \frac{1}{3} \frac{\partial}{\partial y} (\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y})],

(5)

ρ C_{p} (\frac{\partial T}{\partial t} + u \frac{\partial T}{\partial x} + v \frac{\partial T}{\partial y}) = \frac{\partial p}{\partial t} + u \frac{\partial p}{\partial x} + v \frac{\partial p}{\partial y} + \frac{\partial}{\partial x} (k (T) \frac{\partial T}{\partial x}) + \frac{\partial}{\partial y} (k (T) \frac{\partial T}{\partial y}) + \frac{16 δ^{*} {T_{\infty}}^{3}}{3 k^{*}} (\frac{\partial^{2} T}{\partial x \partial y} + \frac{\partial^{2} T}{\partial y^{2}}) + φ,

(6)

where viscous dissipation term is given by

(1 + λ_{1} \frac{\partial}{\partial t}) φ = μ (1 + λ_{2} \frac{\partial}{\partial t}) [2 {(\frac{\partial u}{\partial x})}^{2} + 2 {(\frac{\partial v}{\partial y})}^{2} + {(\frac{\partial v}{\partial x} + \frac{\partial u}{\partial y})}^{2} - \frac{2}{3} {(\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y})}^{2}] .

(7)

Thermal conductivity as a linear function of temperature is given as²⁷

k (T) = k_{0} [1 - β' (T - T_{0})] .

(8)

Many studies explored the nature of different pattern of cilia tips. Here, we assume the elliptic nature of cilia tips. The parametric equations for the beating pattern of cilia as proposed by Sleigh¹ and Lardner and Shack² can be expressed mathematically as

y = h (x, t) = \pm b ϕ \cos [\frac{2 π}{λ} (x - ct)] \pm d,

(9)

where $x$ and y describes the paths of cilia tips in horizontal and vertical position, respectively. $ϵ$ and $α$ signify measurement of cilia tips and elliptical paths, respectively. $b$ and $ϕ$ together makes the amplitude of metachronal wave. Wavelength and speed of metachronal wave are $λ$ and $c$ , respectively. Also, $x_{0}$ is the reference position of the particle. Inspection of various pattern of cilia in Sleigh¹ and Lardner and Shack² motivate us to consider the cilia to move in the elliptical paths so that implicit form of horizontal orientation of cilia tips can be written as

x = f (x, x_{0}, t) = b α ϕ \sin [\frac{2 π}{λ} (x - ct)] + x_{0} .

(10)

The velocity of cilia tips are

u = \frac{\partial x}{\partial t} = (\frac{\partial}{\partial t} + u \frac{\partial}{\partial x}) f,

(11)

v = \frac{\partial y}{\partial t} = (\frac{\partial}{\partial t} + u \frac{\partial}{\partial x}) h .

(12)

The suitable boundary conditions are

u (x, \pm h, t) = \mp A \frac{\partial u}{\partial y} + (\frac{- \frac{2 π}{λ} b ϕ α c \cos [\frac{2 π}{λ} (x - ct)]}{1 - \frac{2 π}{λ} b ϕ α \cos [\frac{2 π}{λ} (x - ct)]}),

(13)

v (x, \pm h, t) = \pm \frac{\partial η (x, t)}{\partial t} \pm (\frac{\frac{2 π}{λ} b ϕ c \sin [\frac{2 π}{λ} (x - ct)]}{1 - \frac{2 π}{λ} b ϕ α \cos [\frac{2 π}{λ} (x - ct)]}),

(14)

T (x, h, t) = T_{1}, and T (x, - h, t) = T_{0} .

(15)

We allocate the following non- dimensional variables as

\begin{matrix} x' = \frac{x}{d}, h' = \frac{h}{d}, u^{'} = \frac{u}{c}, ρ' = \frac{ρ}{ρ_{0}}, t' = \frac{ct}{d}, p' = \frac{p}{ρ_{0} c^{2}}, ρ_{c}^{'} = \frac{p_{c}}{ρ_{0} c^{2}}, λ_{1}^{'} = \frac{c λ_{1}}{d}, λ_{2}^{'} = \\ \frac{c λ_{2}}{d}, θ = \frac{T - T_{0}}{T_{1} - T_{0}}, L = - α ϕ, α = \frac{2 π d}{λ}, β = β' (T_{1} - T_{0}), R = \frac{ρ_{0} cd}{μ}, M = \frac{σ d}{ρ_{0} c} {B_{0}}^{2}, χ = \\ k_{c} ρ_{0} c^{2}, K_{n} = \frac{A}{d}, Rd = \frac{16 δ^{*}}{3 k^{*} k_{0}} {T_{\infty}}^{3}, Ma = \frac{c}{a^{*}}, v' = \frac{v}{c}, y' = \frac{y}{d}, ϵ = \frac{b ϕ}{d}, \Pr = \frac{μ C_{p}}{k_{0}}, \end{matrix}

(16)

where $ϵ, α, β, R, M, \Pr, χ, K_{n}, Rd$ , and $Ma$ are the amplitude ratio, the wave number, the variable thermal conductivity parameter, Reynolds number, the magnetic number, Prandtl number, the compressibility parameter, slip parameter, Radiation number, and Mach number, respectively.

In term of equation (16), equations (2–10) and equations (13–15) reduce to

ρ = \exp (χ (p - p_{c}),

(17)

ρ (\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y}) + \frac{\partial ρ}{\partial t} + u \frac{\partial ρ}{\partial x} + v \frac{\partial ρ}{\partial y} = 0,

(18)

(1 + λ_{1} \frac{\partial}{\partial t}) [ρ (\frac{\partial u}{\partial t} + u \frac{\partial u}{\partial x} + v \frac{\partial u}{\partial y})] = \frac{1}{R} (1 + λ_{2} \frac{\partial}{\partial t}) [\nabla^{2} u + \frac{1}{3} \frac{\partial}{\partial x} (\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y})] - (1 + λ_{1} \frac{\partial}{\partial t}) (\frac{\partial p}{\partial x} + Mu),

(19)

(1 + λ_{1} \frac{\partial}{\partial t}) [ρ (\frac{\partial v}{\partial t} + u \frac{\partial v}{\partial x} + v \frac{\partial v}{\partial y})] = \frac{1}{R} (1 + λ_{2} \frac{\partial}{\partial t}) [\nabla^{2} v + \frac{1}{3} \frac{\partial}{\partial y} (\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y})] - (1 + λ_{1} \frac{\partial}{\partial t}) \frac{\partial p}{\partial y},

(20)

\begin{matrix} (1 + λ_{1} \frac{\partial}{\partial t}) [ρ (u \frac{\partial θ}{\partial x} + v \frac{\partial θ}{\partial y} + \frac{\partial θ}{\partial t})] = (1 + λ_{1} \frac{\partial}{\partial t}) \\ [(γ^{*} - 1) M a^{2} (\frac{\partial p}{\partial t} + u \frac{\partial p}{\partial x} + v \frac{\partial p}{\partial y}) + \frac{β}{\Pr \cdot R} [{(\frac{\partial θ}{\partial x})}^{2} + {(\frac{\partial θ}{\partial y})}^{2}] + \frac{(1 + β θ)}{\Pr . R} \nabla^{2} θ + \frac{Rd}{\Pr \cdot R} \frac{\partial}{\partial y} (\frac{\partial θ}{\partial x} + \frac{\partial θ}{\partial y})] \\ + (1 + λ_{2} \frac{\partial}{\partial t}) (\frac{1}{R} (γ^{*} - 1) {Ma}^{2}) [2 {(\frac{\partial u}{\partial x})}^{2} + 2 {(\frac{\partial v}{\partial y})}^{2} + {(\frac{\partial v}{\partial x} + \frac{\partial u}{\partial y})}^{2} - \frac{2}{3} {(\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y})}^{2}], \end{matrix}

(21)

h (x, t) = \pm 1 \pm η (x, t) = \pm 1 \pm ϵ \cos (α (x - t)),

(22)

u (x, y, t) = \mp K_{n} \frac{\partial u}{\partial y} + (\frac{- α^{2} ϵ \cos (α (x - t))}{1 - α^{2} ϵ \cos (α (x - t))}), at y = \pm 1 \pm η

(23)

v (x, y, t) = \pm \frac{\partial η (x, t)}{\partial t} \pm (\frac{α^{2} ϵ \sin (α (x - t))}{1 - α^{2} ϵ \cos (α (x - t))}), at y = \pm 1 \pm η

(24)

θ (x, y, t) = 1, at y = 1 + η

(25)

θ (x, y, t) = 0 . at y = - 1 - η

(26)

in which primes have been dropped.

Solution of the problem

For steady parallel flow case in which $u = u_{0} (y)$ , $v = 0$ and $\frac{dp}{dx} = \frac{d p_{0}}{dx} =$ constant, the exact solution is given as follows

u_{0} (y) = \frac{R}{ζ^{2}} \frac{d p_{0}}{dx} \frac{(\cosh ζ y - \cosh ζ - K_{n} ζ \sinh ζ)}{(\cosh ζ + K_{n} ζ \sinh ζ)} + (\frac{- α^{2} ϵ \cos (α (x - t))}{1 - α^{2} ϵ \cos (α (x - t))}) (\frac{\cosh ζ y}{\cosh ζ + K_{n} ζ \sinh ζ y}),

(27)

where $ζ = \sqrt{MR}$ .

When $K_{n} \to 0$ that is, in case of non-slip fluid equation (27) becomes

u_{0} (y) = \frac{R}{ζ^{2}} \frac{d p_{0}}{dx} \frac{(\cosh ζ y - \cosh ζ)}{\cosh ζ} + (\frac{- α^{2} ϵ \cos (α (x - t))}{1 - α^{2} ϵ \cos (α (x - t))}) (\frac{\cosh ζ y}{\cosh ζ}) .

(28)

Pure peristalsis¹⁸ is considered here in which $\frac{\partial p_{0}}{\partial x} = 0$ .

Now consider $θ = θ_{0} (y), u = u_{0} (y)$ , $v = 0$ and $\frac{dp}{dx} = \frac{d p_{0}}{dx} = constant$ . The solution of the equation (21) along with equations (25) and (26) will be

θ_{0} (y) = - \frac{(1 + Rd)}{β} + \sqrt{{(\frac{1 + Rd}{β})}^{2} + \frac{1}{2} (1 + 2 \frac{(1 + Rd)}{β}) (1 + y)} .

(29)

We use regular perturbation method to solve the above non-linear equations. The perturbed solution assumptions of the velocity components, pressure, temperature, the viscous dissipation and density are as follows:

\begin{matrix} u = ϵ u_{1} (x, y, t) + ϵ^{2} u_{2} (x, y, t) + . . ., \\ v = ϵ v_{1} (x, y, t) + ϵ^{2} v_{2} (x, y, t) + . . ., \\ p = ϵ p_{1} (x, y, t) + ϵ^{2} p_{2} (x, y, t) + . . ., \\ θ = θ_{0} (y) + ϵ θ_{1} (x, y, t) + ϵ^{2} θ_{2} (x, y, t) + . . ., \\ φ = ϵ φ_{1} (x, y, t) + ϵ^{2} φ_{2} (x, y, t) + \dots, \\ ρ = 1 + ϵ ρ_{1} (x, y, t) + ϵ^{2} ρ_{2} (x, y, t) + . \dots \end{matrix}

(30)

The following system was obtained after incorporating expansion (30) into (17–26) and equating $O (ϵ)$ and $O (ϵ^{2})$ ,we get

for $O (ϵ)$ ;

ρ_{1} = χ p_{1}, \frac{\partial ρ_{1}}{\partial t} + \frac{\partial v_{1}}{\partial y} + \frac{\partial u_{1}}{\partial x} = 0,

(1 + λ_{1} \frac{\partial}{\partial t}) \frac{\partial u_{1}}{\partial t} = - (1 + λ_{1} \frac{\partial}{\partial t}) \frac{\partial p_{1}}{\partial x} - M (1 + λ_{1} \frac{\partial}{\partial t}) u_{1} + \frac{1}{R} (1 + λ_{2} \frac{\partial}{\partial t}) [\nabla^{2} u_{1} + \frac{1}{3} \frac{\partial}{\partial x} (\frac{\partial u_{1}}{\partial x} + \frac{\partial v_{1}}{\partial y})],

(1 + λ_{1} \frac{\partial}{\partial t}) \frac{\partial v_{1}}{\partial t} = - (1 + λ_{1} \frac{\partial}{\partial t}) \frac{\partial p_{1}}{\partial y} + \frac{1}{R} (1 + λ_{2} \frac{\partial}{\partial t}) [\frac{1}{3} \frac{\partial}{\partial y} (\frac{\partial u_{1}}{\partial x} + \frac{\partial v_{1}}{\partial y}) + \nabla^{2} v_{1}],

\frac{\partial θ_{1}}{\partial t} + v_{1} \frac{\partial θ_{0}}{\partial y} = (γ^{*} - 1) M a^{2} \frac{\partial p_{1}}{\partial t} + \frac{1}{\Pr \cdot R} [(1 + β θ_{0}) \nabla^{2} θ_{1} + β (2 \frac{\partial θ_{0}}{\partial y} \frac{\partial θ_{1}}{\partial y} + θ_{1} \frac{\partial^{2} θ_{0}}{\partial y^{2}}) + Rd \frac{\partial}{\partial y} (\frac{\partial θ_{1}}{\partial x} + \frac{\partial θ_{1}}{\partial y})],

u_{1} (x, \pm 1, t) = \mp K_{n} u_{1 y} (x, \pm 1, t) - \frac{α^{2}}{2} (e^{i α (x - t)} - e^{- i α (x - t)}),

v_{1} (x, \pm 1, t) = \mp \frac{i α}{2} (e^{i α (x - t)} - e^{- i α (x - t)}) \pm \frac{α}{2} (e^{i α (x - t)} - e^{- i α (x - t)}),

θ_{1} (x, \pm 1, t) = \mp \frac{1}{2} (e^{i α (x - t)} + e^{- i α (x - t)}) θ_{0 y} (\pm 1) .

(31)

For $O (ϵ^{2});$

ρ_{2} = χ p_{2} + \frac{1}{2} χ^{2} {p_{1}}^{2},

\frac{\partial ρ_{2}}{\partial t} + u_{1} \frac{\partial ρ_{1}}{\partial x} + v_{1} \frac{\partial ρ_{1}}{\partial y} + (\frac{\partial u_{2}}{\partial x} + \frac{\partial v_{2}}{\partial y}) + ρ_{1} (\frac{\partial u_{1}}{\partial x} + \frac{\partial v_{1}}{\partial y}) = 0,

\begin{matrix} (1 + λ_{1} \frac{\partial}{\partial t}) [\frac{\partial u_{2}}{\partial t} + u_{1} \frac{\partial u_{1}}{\partial x} + v_{1} \frac{\partial v_{1}}{\partial y} + ρ_{1} \frac{\partial u_{1}}{\partial t}] = - (1 + λ_{1} \frac{\partial}{\partial t}) \frac{\partial p_{2}}{\partial x} + \frac{1}{R} (1 + λ_{2} \frac{\partial}{\partial t}) \\ [\nabla^{2} u_{2} + \frac{1}{3} \frac{\partial}{\partial x} (\frac{\partial u_{2}}{\partial x} + \frac{\partial v_{2}}{\partial y})] - M (1 + λ_{1} \frac{\partial}{\partial t}) u_{2}, \end{matrix}

(1 + λ_{1} \frac{\partial}{\partial t}) [\frac{\partial v_{2}}{\partial t} + u_{1} \frac{\partial v_{1}}{\partial x} + v_{1} \frac{\partial v_{1}}{\partial y} + ρ_{1} \frac{\partial v_{1}}{\partial t}] = - (1 + λ_{1} \frac{\partial}{\partial t}) \frac{\partial p_{2}}{\partial y} + \frac{1}{R} (1 + λ_{2} \frac{\partial}{\partial t}) [\nabla^{2} v_{2} + \frac{1}{3} \frac{\partial}{\partial y} (\frac{\partial u_{2}}{\partial x} + \frac{\partial v_{2}}{\partial y})],

\begin{matrix} (1 + λ_{1} \frac{\partial}{\partial t}) [\frac{\partial θ_{2}}{\partial t} + u_{1} \frac{\partial θ_{1}}{\partial x} + v_{1} \frac{\partial θ_{1}}{\partial y} + ρ_{1} \frac{\partial θ_{1}}{\partial t} + ρ_{1} v_{1} \frac{\partial θ_{0}}{\partial y} + v_{2} \frac{\partial θ_{0}}{\partial y}] = (1 + λ_{1} \frac{\partial}{\partial t}) \\ [\begin{matrix} (γ^{*} - 1) M a^{2} (\frac{\partial p_{2}}{\partial t} + u_{1} \frac{\partial p_{1}}{\partial x} + v_{1} \frac{\partial p_{1}}{\partial y}) \\ + \frac{1}{\Pr \cdot R} [(1 + β θ_{0}) \nabla^{2} θ_{2} + β (θ_{2} \frac{\partial^{2} θ_{0}}{\partial y^{2}} + θ_{1} (\nabla^{2} θ_{1}) + {(\frac{\partial θ_{1}}{\partial x})}^{2} + {(\frac{\partial θ_{1}}{\partial y})}^{2} + 2 \frac{\partial θ_{0}}{\partial y} \frac{\partial θ_{2}}{\partial y})] \\ + \frac{Rd}{\Pr \cdot R} \frac{\partial}{\partial y} (\frac{\partial θ_{2}}{\partial x} + \frac{\partial θ_{2}}{\partial y}) \end{matrix}] \\ + (1 + λ_{2} \frac{\partial}{\partial t}) (\frac{1}{R} (γ^{*} - 1) M a^{2}) [2 {(\frac{\partial u_{1}}{\partial x})}^{2} + 2 {(\frac{\partial v_{1}}{\partial y})}^{2} + {(\frac{\partial v_{1}}{\partial x} + \frac{\partial u_{1}}{\partial y})}^{2} - \frac{2}{3} {(\frac{\partial u_{1}}{\partial x} + \frac{\partial v_{1}}{\partial y})}^{2}], \end{matrix}

\begin{matrix} u_{2} (x, \pm 1, t) \pm \frac{1}{2} (e^{i α (x - t)} + e^{- i α (x - t)}) u_{1 y} (x, \pm 1, t) = \mp K_{n} \\ [u_{2 y} (x, \pm 1, t) \pm \frac{1}{2} (e^{i α (x - t)} + e^{- i α (x - t)}) u_{1 yy} (x, \pm 1, t)] - \frac{α^{4}}{4} {(e^{i α (x - t)} + e^{- i α (x - t)})}^{2}, \end{matrix}

v_{2} (x, \pm 1, t) \pm \frac{1}{2} (e^{i α (x - t)} + e^{- i α (x - t)}) v_{1 y} (x, \pm 1, t) = \pm \frac{α^{3}}{4} (e^{2 i α (x - t)} - e^{- 2 i α (x - t)}),

θ_{2} (x, \pm 1, t) \pm \frac{1}{2} (e^{i α (x - t)} + e^{- i α (x - t)}) θ_{1 y} (x, \pm 1, t) \pm \frac{1}{2} {(\frac{e^{i α (x - t)} + e^{- i α (x - t)}}{2})}^{2} θ_{0 yy} (\pm 1) = 0 .

(32)

Here we obtained the boundary condition by following the same strategy as proposed by Mekheimer et al.¹² and Khan and Rafaqat¹⁸

We adopt the solution in the form

\begin{matrix} u_{1} (x, y, t) = U_{1} (y) e^{i α (x - t)} + \bar{U_{1}} (y) e^{- i α (x - t)}, \\ v_{1} (x, y, t) = V_{1} (y) e^{i α (x - t)} + \bar{V_{1}} (y) e^{- i α (x - t)}, \\ p_{1} (x, y, t) = P_{1} (y) e^{i α (x - t)} + \bar{P_{1}} (y) e^{- i α (x - t)}, \\ ρ_{1} (x, y, t) = χ P_{1} (y) e^{i α (x - t)} + χ \bar{P_{1}} (y) e^{- i α (x - t)}, \\ θ_{1} (x, y, t) = θ_{11} (y) e^{i α (x - t)} + \bar{θ_{11}} (y) e^{- i α (x - t)} . \end{matrix}

(33)

and

\begin{matrix} u_{2} (x, y, t) = U_{20} (y) + U_{2} (y) e^{2 i α (x - t)} + \bar{U_{2}} (y) e^{- 2 i α (x - t)}, \\ v_{2} (x, y, t) = V_{20} (y) + V_{2} (y) e^{2 i α (x - t)} + \bar{V_{2}} (y) e^{- 2 i α (x - t)}, \\ p_{2} (x, y, t) = P_{20} (y) + P_{2} (y) e^{2 i α (x - t)} + \bar{P_{2}} (y) e^{- 2 i α (x - t)}, \\ ρ_{2} (x, y, t) = D_{20} (y) + D_{2} (y) e^{2 i α (x - t)} + \bar{D_{2}} (y) e^{- 2 i α (x - t)}, \\ θ_{2} (x, y, t) = θ_{20} (y) + θ_{22} (y) e^{2 i α (x - t)} + \bar{θ_{22}} (y) e^{- 2 i α (x - t)} . \end{matrix}

(34)

where overbar denotes the complex conjugate.

Incorporating equations (33) and (34) into equations (31) and (32), we get these sets of following first-order and second-order equations, respectively:

V_{1}' + i α U_{1} = i α χ P_{1},

U_{1}'' - (\frac{γ_{1} R}{γ_{2}} (M - i α) + \frac{4 α^{2}}{3}) U_{1} - i α \frac{γ_{1}}{γ_{2}} {RP}_{1} + \frac{i α}{3} V_{1}' = 0,

V_{1}'' - \frac{3}{4} (α^{2} - i α \frac{γ_{1} R}{γ_{2}}) V_{1} - \frac{3}{4} \frac{γ_{1} R}{γ_{2}} P_{1}' + \frac{i α}{4} U_{1}' = 0,

θ_{11}'' + β_{2} θ_{11}' + β_{3} θ_{11} - \frac{\Pr . R}{(1 + Rd + β N_{0})} [N_{1} V_{1} + i α (γ^{*} - 1) {Ma}^{2} P_{1}] = 0,

U_{1} (\pm 1) = \mp K_{n} U_{1}' (\pm 1) - \frac{α^{2}}{2}, V_{1} (\pm 1) = \mp \frac{i α}{2} \pm \frac{α}{2},

θ_{11} (1) = - \frac{1}{8} \frac{(1 + 2 \frac{(1 + Rd)}{β})}{(1 + \frac{(1 + Rd)}{β})}, θ_{11} (- 1) = \frac{1}{4} (1 + \frac{β}{2 (1 + Rd)}) .

(35)

D_{20} = χ P_{20} + χ^{2} \bar{P_{1}} P_{1},

- P_{20}' + \frac{4}{3 R} V_{20}'' = F,

U_{20}'' - MR U_{20} = RL,

V_{20}' = - χ E,

{θ_{20}}^{″} + \frac{2 N_{1} β}{Rd} θ_{20}^{'} + \frac{N_{2} β}{Rd} θ_{20} = \Pr . R . D^{*},

U_{20} (\pm 1) \pm \frac{1}{2} (\bar{U_{1}'} (\pm 1) + U_{1}' (\pm 1)) = \pm \frac{1}{2} (\bar{U_{1} ″} (\pm 1) + U_{1} ″ (\pm 1)) \mp K_{n} ({U_{2}}_{0}' (\pm 1) - \frac{α^{4}}{2},

V_{20} (\pm 1) \pm \frac{1}{2} (\bar{V_{1}'} (\pm 1) + V_{1}' (\pm 1)) = 0,

θ_{20} (\pm 1) \pm \frac{1}{2} (\bar{θ_{11}'} (\pm 1) + θ_{11}' (\pm 1))) \pm \frac{1}{4} ({θ_{0}}^{″} (\pm 1)) = 0 .

(36)

where $N_{0} = θ_{0} (y), N_{1} = \frac{d θ_{0}}{dy}, N_{2} = \frac{d^{2} θ_{0}}{d y^{2}} .$

To accomplish the main equations we follow the same approach as given by Mekheimer et al.,¹² and which was trailed by Khan and Rafaqat¹⁸ later. Thus, we neglect the long calculations and locate the first-order solution for the pressure, velocity, and temperature. We stifled here the prolonged expression for second order solution of heat transfer.

V_{1} (y) = C_{1} \sinh a_{1} y + C_{2} \sinh a_{2} y,

(37a)

U_{1} (y) = b_{1} C_{1} \cosh a_{1} y + C_{3} χ + b_{2} {C_{2}}_{2} \cosh a_{2} y,

(37b)

where C_{1} = (\frac{i α}{2}) (\frac{b_{2} g_{2}}{G}) - C_{11}, C_{2} = - \frac{i α}{2} (\frac{1}{\sinh a_{2}}) (1 + \frac{b_{2} g_{2}}{G} \sinh a_{1}) + C_{22},

(37c)

P_{1} (y) = C_{1} (\frac{({a_{1}}^{2} - {β_{1}}^{2})}{γ a_{1}}) \cosh a_{1} y + C_{2} (\frac{({a_{2}}^{2} - {β_{1}}^{2})}{γ a_{2}}) \cosh a_{2} y + C_{22},

(37d)

θ_{11} (y) = C_{4} e^{d_{1} y} + C_{5} e^{d_{2} y} + Q (y) .

(37e)

V_{20} = - χ (P_{1} \bar{V_{1}} + \bar{P_{1}} V_{1}) + D_{1},

U_{20} (y) = D_{2} \cosh ζ y + D_{3} \sinh ζ y + N (y),

P_{20} (y) = D_{4} - \frac{4 χ}{3 R} E (y) - \int_{- 1}^{y} F (r) dr .

The parameters and complex constants are characterized in Appendix.

It may be noted that if we substitute $C_{11} = 0$ and $C_{22} = 0$ in equation (37c) we get same result as reported earlier by Khan and Rafaqat.¹⁸ If we set $α = 0$ in the expression of $D_{2}$ we get solution of $U_{2_{20}} (y)$ same as obtained in Khan and Rafaqat.¹⁸ The velocity perturbation function $G (y)$ and the mean net axial velocity $\bar{u}$ are defined as follow

G (y) = - \frac{200}{α^{2} R^{2}} [N (y) - N (1)],

(38)

\bar{u} = ϵ^{2} U_{20} (y) .

(39)

Results and discussion

This section is dedicated to illustrate the change in perturbation function, mean axial velocity and temperature profile for pertinent parameters.

Figures 2 to 6 depicts the effect of $M, R, λ_{1}, λ_{2}$ and $χ$ on mean velocity perturbation function $G (y)$ . Figures 2 and 3 show the influence of $M$ and $R$ on $G (y)$ . It is investigated that increment in $M$ causes a decrement in $G (y)$ and makes it more flat. Reason behind the flatness of $G (y)$ is the presence of current density in equation of motion in which dominant part is the hindering force with magnitude $σ {B_{0}}^{2} u$ results in decreasing the mean velocity perturbation function. The same results are observed for increasing R. Figures 4 and 5 examine the effect of $G (y)$ on relaxation and retardation time, respectively. Both relaxation and retardation time show decaying impact on $G (y)$ . It is also investigated that $G (y)$ attains its maximum value in the vicinity of center of channel and over some range it shows a constant behavior. Figure 6 addresses the decreasing influence of compressibility parameter on $G (y)$ . In Figures 5 and 6 it is noticeable that large Reynolds number and magnetic field are responsible for the flatness of curve.

Figure 2.

Variation in perturbation function $G (y)$ towards $M$ .

Figure 3.

Variation in perturbation function $G (y)$ towards $R$ .

Figure 4.

Variation in perturbation function $G (y)$ towards $λ_{1}$ .

Figure 5.

Variation in perturbation function $G (y)$ towards $λ_{2}$ .

Figure 6.

Variation in perturbation function $G (y)$ towards $χ$ .

Outcomes of emerging variables involved in present analysis on mean axial velocity can be examined from Figures 7 to 11. Figures 7 and 8 explicate the variation trend of mean axial velocity against $M$ and $R$ . It is witnessed that both $M$ and $R$ depict opposite behavior on axial velocity. Here, $M$ has an inverse impact on axial velocity while R has a direct impact. This is because the resistive nature of magnetic field. Figures 9 and 10 display the effect of $λ_{1}$ and $λ_{2}$ on axial velocity. It is perceived from these figures that enhancement in relaxation and retardation time corresponds to lower velocity. Figure 11 addresses the slowing down of the fluid with compressibility parameter.

Figure 7.

Variation in mean axial velocity towards $M$ .

Figure 8.

Variation in mean axial velocity towards $R$ .

Figure 9.

Variation in mean axial velocity towards $λ_{1}$ .

Figure 10.

Variation in mean axial velocity towards $λ_{1}$ .

Figure 11.

Variation in mean axial velocity towards $χ$ .

The variation trend of temperature against emerging parameters are inferred in Figures 12 to 19. Figure 12 witnessed the growing impact of variable thermal conductivity on temperature. Figures 13 to 15 made to see the effects of $M, R$ , and $Rd$ . An increasing response of temperature towards $R$ is inferred in Figure 14. Since kinetic energy is directly related to fluid temperature and enhancement in $R$ causes an increment in kinetic energy of the fluid particles, which in turn increases fluid temperature. Moreover a similar pattern is noticed in Figures 13 and 15 for rising values of $M$ and $Rd$ . Figure 16 ensures the decay of temperature with increasing $\Pr$ . Since temperature growth or decay related inversely with thermal diffusivity. And higher Prandtl number lessen thermal diffusivity which causes more transfer of heat. Variational trends of temperature for escalating values of Mach number is inferred in Figures 17 to 19. These figures examine the behavior of temperature in subsonic, transonic and supersonic regions, separately. Response of temperature toward Mach number is increasing. Here, strong impact is noticed for supersonic region while a weak impact is noticed for subsonic region. Also, Mach number intensified the heat transfer mechanism in the region $1.5 \leq Ma \leq 5.0$ .

Figure 12.

Variation in temperature towards $β$ .

Figure 13.

Variation in temperature towards $M$ .

Figure 14.

Variation in temperature towards $R$ .

Figure 15.

Variation in temperature towards $Rd$ .

Figure 16.

Variation in temperature towards $\Pr$ .

Figure 17.

Variation in temperature towards $Ma$ (subsonic region).

Figure 18.

Variation in temperature towards $Ma$ (transonic region).

Figure 19.

Variation in temperature towards $Ma$ (supersonic region).

Conclusion

This study presents heat transfer analysis of a compressible Jeffrey fluid with variable thermal conductivity. It also incorporates the impact of MHD and radiation on cilia induced flow. The worth mentioning points of the performed analysis can be concisely written as:

A decaying nature of velocity perturbation function has been reported for all emerging parameters.

The mean axial velocity depicts a direct relation with $R$ while an inverse relation is observed for all other parameters such as $M, λ_{1}, λ_{2}$ , and $χ$ .

Velocity profile diminishes at the center of the channel but increases close to the walls by rising compressibility parameter.

The temperature profile against variable thermal conductivity parameter exhibits a similar behavior to that of mean axial velocity towards $R$ .

It is analyzed that variational trends of temperature with Mach number are direct. Also, these effects are strongly analyzed in supersonic region while weak trends are investigated in subsonic region.

Footnotes

Appendix

The parameters and variables found in section 3 for the first order and second order solution are given as

γ = - (\frac{i α χ}{3}) + (\frac{γ_{1}}{γ_{2}}) R, γ_{1} = (1 - i α λ_{1}), γ_{2} = (1 - i α λ_{2}),

{β_{0}}^{2} = α^{2} - (\frac{γ_{1}}{γ_{2}}) R (i α - M), {β_{1}}^{2} = α^{2} - (\frac{i α γ_{1}}{γ_{2}}) R,

β_{2} = \frac{(i α Rd + 2 β N_{1})}{(1 + Rd + β N_{0})}, β_{3} = \frac{(i α \Pr . R - α^{2} - β (α^{2} N_{0} + N_{2}))}{(1 + Rd + β N_{0})}

a_{1}^{2} = (\frac{({β_{0}}^{2} + v^{2}) + \sqrt{{({β_{0}}^{2} + v^{2})}^{2} - 4 {β_{1}}^{2} {v_{1}}^{2}}}{2}), a_{2}^{2} = (\frac{({β_{0}}^{2} + v^{2}) - \sqrt{{({β_{0}}^{2} + v^{2})}^{2} - 4 {β_{1}}^{2} {v_{1}}^{2}}}{2}),

b_{1} = χ (\frac{({a_{1}}^{2} - {β_{1}}^{2})}{γ a_{1}}) + (\frac{i a_{1}}{α}), b_{2} = χ (\frac{({a_{2}}^{2} - {β_{1}}^{2})}{γ a_{2}}) + (\frac{i a_{2}}{α}),

b_{3} = \frac{i C_{1} ({a_{1}}^{2} Rd + iPr . R . α)}{d_{3}}, b_{4} = \frac{i C_{2} ({a_{2}}^{2} Rd + iPr . R . α)}{d_{4}}, b_{5} = \frac{C_{1} α Rd}{d_{5}}, b_{6} = \frac{C_{2} α Rd}{d_{6}}

d_{1} = (- \frac{i α}{2} - \frac{{(- 1)}^{3 / 4} \sqrt{α (4 P r . R - i α R d)}}{2 \sqrt{R d}}), d_{2} = (- \frac{i α}{2} + \frac{{(- 1)}^{3 / 4} \sqrt{α (4 P r . R - i α R d)}}{2 \sqrt{R d}})

C_{1} = (\frac{i α}{2}) (\frac{b_{2} g_{2}}{G}) - C_{11}, C_{2} = - \frac{i α}{2} (\frac{1}{\sinh a_{2}}) (1 + \frac{b_{2} g_{2}}{G} \sinh a_{1}) + C_{22}, C_{3} = 0 .

C_{11} = \frac{α}{2} \frac{(b_{2} g_{2} + α \sinh a_{2})}{G}, C_{22} = \frac{α}{2} (\frac{1}{\sinh a_{2}}) (1 + \frac{b_{2} g_{2}}{G} \sinh a_{1}) + \frac{α^{2}}{2} \frac{\sinh a_{1}}{G},

\begin{matrix} C_{4} = - \frac{α M a^{2} . Rd Csch (- i α - 2 d_{1}) e^{d_{1}}}{2 γ a_{1} a_{2} d_{5} d_{6}} \\ [\begin{matrix} - a_{2} b_{3} d_{3} d_{6} e^{- 2 d_{2}} ({a_{1}}^{2} - {β_{1}}^{2}) Cosh a_{1} + a_{1} b_{4} d_{4} d_{5} (e^{ι α} + e^{- 2 d_{1}}) ({a_{2}}^{2} - {β_{1}}^{2}) Cosh a_{2} - α Rd \\ [a_{1} a_{2} C_{1} d_{6} e^{- 2 d_{2}} ({a_{1}}^{2} - {β_{1}}^{2}) Sinh a_{1} + a_{1} a_{2} C_{2} d_{5} (e^{ι α} + e^{- 2 d_{1}}) ({a_{2}}^{2} - {β_{1}}^{2}) Sinh a_{2}] \end{matrix}] (γ^{*} - 1) \end{matrix}

\begin{matrix} C_{5} = \frac{α M a^{2} . Rd [Coth (i α + 2 d_{2}) - 1] e^{d_{2}} (e^{- 2 d_{1}} - 1)}{2 γ a_{1} a_{2} d_{5} d_{6}} \\ [\begin{matrix} - a_{2} b_{3} d_{3} d_{6} ({a_{1}}^{2} - {β_{1}}^{2}) Cosh a_{1} + a_{1} b_{4} d_{4} d_{5} ({a_{2}}^{2} - {β_{1}}^{2}) \\ Cosh a_{2} - α Rd [a_{1} a_{2} C_{1} d_{6} ({a_{1}}^{2} - {β_{1}}^{2}) Sinh a_{1} + a_{1} a_{2} C_{2} d_{5} ({a_{2}}^{2} - {β_{1}}^{2}) Sinh a_{2}] \end{matrix}] (γ^{*} - 1) \end{matrix}

G = b_{1} g_{1} \sinh a_{2} - b_{2} g_{2} \sinh a_{1},

where $g_{1} = \cosh a_{1} + K_{n} a_{1} \sinh a_{1}$ , $g_{2} = \cosh a_{2} + K_{n} a_{2} \sinh a_{2}$ .

D_{1} = \frac{α χ}{2} (P_{1} (1) + {\bar{P}}_{1} (1)) + (\frac{i α}{2}) (U_{1} (1) - {\bar{U}}_{1} (1)),

D_{2} = \frac{- 1}{2 (\cosh ζ + K_{n} ζ \sinh ζ)} \times [β_{4} + β_{5} + N (1) + N (- 1) + K_{n} (β_{6} - β_{7} + N^{'} (1) - N' (- 1)) - α^{4}],

D_{3} = \frac{- 1}{2 (\sinh ζ + K_{n} ζ \cosh ζ)} \times [β_{4} - β_{5} + N (1) - N (- 1) + K_{n} (β_{6} + β_{7} + N^{'} (1) + N' (- 1))],

D_{4} = P_{20} (- 1) + \frac{4 χ}{3 R} E (- 1),

where $β_{4} = \frac{1}{2} (U_{1}^{'} (1) + \bar{U_{1}^{'}} (1))$ , $β_{5} = - \frac{1}{2} (U_{1}^{'} (- 1) + \bar{U_{1}^{'}} (- 1))$ ,

β_{6} = \frac{1}{2} (U_{1}^{″} (1) + \bar{U_{1}^{″}} (1)), β_{7} = - \frac{1}{2} (U_{1}^{″} (- 1) + \bar{U_{1}^{″}} (- 1)) .

The functions $F, L, N$ and $Q$ are given as

F = i α χ P_{1} {\bar{V}}_{1} - i α \bar{P_{1}} V_{1} + V_{1}' {\bar{V}}_{1} + V_{1} \bar{V_{1}'} - i α U_{1} {\bar{V}}_{1} + i α \bar{U_{1}} V_{1},

L = \frac{d}{dy} (V_{1} {\bar{U}}_{1} + U_{1} {\bar{V}}_{1}), E = \frac{d}{dy} (P_{1} \bar{V_{1}} + V_{1} {\bar{P}}_{1}),

\begin{matrix} D^{*} = N_{1} V_{20} + χ N_{1} (P_{1} {\bar{V}}_{1} + \bar{P_{1}} V_{1}) + i α (\bar{U_{1}} θ_{11} - U_{1} \bar{θ_{11}}) + (\bar{V_{1}} θ_{11}' + V_{1} \bar{θ_{11}'}) + i α χ (P_{1} \bar{θ_{11}} - \bar{P_{1}} θ_{11}) + \\ (1 - γ^{*}) M a^{2} [i α (U_{1} \bar{P_{1}} - P_{1} \bar{U_{1}}) - (V_{1} \bar{P_{1}'} + P_{1}' \bar{V_{1}}) - P_{20}] \\ - \frac{β}{\Pr . R} (θ_{11} \bar{θ_{11}^{″}} + \bar{θ_{11}} θ_{11}^{″} + 2 θ_{11}' \bar{θ_{11}'}) - α_{1} [(\frac{8}{3}) α^{2} U_{1} \bar{U_{1}} + (\frac{8}{3}) V_{1}' \bar{V_{1}'} + \\ 2 α^{2} V_{1} \bar{V_{1}} + 2 U_{1}' \bar{U_{1}'} + 2 i α V_{1} \bar{U_{1}'} + 2 i α U_{1}' {\bar{V}}_{1} - (\frac{4 i α}{3}) ({\bar{U}}_{1} {V_{1}}^{'} - U_{1} \bar{V_{1}'})], \end{matrix}

where $α_{1} = (γ^{*} - 1) \frac{M a^{2}}{R} (\frac{γ_{2}}{γ_{1}})$ ,

N (y) = R {\begin{matrix} \frac{h_{1}}{{(a_{1} + \bar{a_{1}})}^{2} - ζ^{2}} \cosh (a_{1} + \bar{a_{1}}) y + \frac{h_{2}}{{(a_{1} + \bar{a_{2}})}^{2} - ζ^{2}} \cosh (a_{1} + \bar{a_{2}}) y \\ + \frac{h_{3}}{{(a_{2} + \bar{a_{1}})}^{2} - ζ^{2}} \cosh (a_{2} + \bar{a_{1}}) y + \frac{h_{4}}{{(a_{2} + \bar{a_{2}})}^{2} - ζ^{2}} \cosh (a_{2} + \bar{a_{2}}) y \\ + \frac{h_{11}}{{(a_{1} - \bar{a_{1}})}^{2} - ζ^{2}} \cosh (a_{1} - \bar{a_{1}}) y + \frac{h_{22}}{{(a_{1} - \bar{a_{2}})}^{2} - ζ^{2}} \cosh (a_{1} - \bar{a_{2}}) y \\ + \frac{h_{33}}{{(a_{2} - \bar{a_{1}})}^{2} - ζ^{2}} \cosh (a_{2} - \bar{a_{1}}) y + \frac{h_{44}}{{(a_{2} - \bar{a_{2}})}^{2} - ζ^{2}} \cosh (a_{2} - \bar{a_{2}}) y \end{matrix}},

h_{1} = \frac{1}{2} (\bar{a_{1}} \bar{b_{1}} \bar{C_{1}} C_{1} + a_{1} b_{1} C_{1} \bar{C_{1}} + a_{1} C_{1} \bar{b_{1}} \bar{C_{1}} + \bar{a_{1}} \bar{C_{1}} b_{1} C_{1}),

h_{2} = \frac{1}{2} (\bar{a_{2}} \bar{b_{2}} \bar{C_{2}} C_{1} + a_{1} b_{1} C_{1} \bar{C_{2}} + a_{1} C_{1} \bar{b_{2}} \bar{C_{2}} + \bar{a_{2}} \bar{C_{2}} b_{1} C_{1}),

h_{3} = \frac{1}{2} (\bar{a_{1}} \bar{b_{1}} \bar{C_{1}} C_{2} + a_{2} b_{2} C_{2} \bar{C_{1}} + a_{2} C_{2} \bar{b_{1}} \bar{C_{1}} + \bar{a_{1}} \bar{C_{1}} b_{2} C_{2}),

h_{4} = \frac{1}{2} (\bar{a_{2}} \bar{b_{2}} \bar{C_{2}} C_{2} + a_{2} b_{2} C_{2} \bar{C_{2}} + a_{2} C_{2} \bar{b_{2}} \bar{C_{2}} + \bar{a_{2}} \bar{C_{2}} b_{2} C_{2}),

h_{11} = - \frac{1}{2} (\bar{a_{1}} \bar{b_{1}} \bar{C_{1}} C_{1} + a_{1} b_{1} C_{1} \bar{C_{1}} - a_{1} C_{1} \bar{b_{1}} \bar{C_{1}} - \bar{a_{1}} \bar{C_{1}} b_{1} C_{1}),

h_{22} = - \frac{1}{2} (\bar{a_{2}} \bar{b_{2}} \bar{C_{2}} C_{1} + b_{1} C_{1} \bar{C_{2}} - a_{1} C_{1} \bar{b_{2}} \bar{C_{2}} - \bar{a_{2}} \bar{C_{2}} b_{1} C_{1}),

h_{33} = - \frac{1}{2} (\bar{a_{1}} \bar{b_{1}} \bar{C_{1}} C_{2} + a_{2} b_{2} C_{2} \bar{C_{1}} - a_{2} C_{2} \bar{b_{1}} \bar{C_{1}} - \bar{a_{1}} \bar{C_{1}} b_{2} C_{2}),

h_{44} = - \frac{1}{2} (\bar{a_{2}} \bar{b_{2}} \bar{C_{2}} C_{2} + a_{2} b_{2} C_{2} \bar{C_{2}} - a_{2} C_{2} \bar{b_{2}} \bar{C_{2}} - \bar{a_{2}} \bar{C_{2}} b_{2} C_{2}) .

Q (y) = \frac{(γ^{*} - 1)}{γ} α M a^{2} . Rd [b_{3} ({a_{1}}^{2} - {β_{1}}^{2}) Cosh a_{1} y + b_{4} ({a_{2}}^{2} - {β_{1}}^{2}) Cosh a_{2} y + b_{5} ({a_{1}}^{2} - {β_{1}}^{2}) Sinh a_{1} y + b_{6} ({a_{2}}^{2} - {β_{1}}^{2}) Sinh a_{2} y],

d_{3} = {a_{1}}^{5} R d^{2} - a_{1} P r^{2} R^{2} α^{2} + {a_{1}}^{3} α Rd (2 iPr . R + Rd α),

d_{4} = {a_{2}}^{5} R d^{2} - a_{2} P r^{2} R^{2} α^{2} + {a_{2}}^{3} α Rd (2 iPr . R + Rd α),

d_{5} = {a_{1}}^{4} R d^{2} - P r^{2} R^{2} α^{2} + {a_{1}}^{2} α Rd (2 iPr . R + Rd α),

d_{6} = {a_{2}}^{4} R d^{2} - P r^{2} R^{2} α^{2} + {a_{2}}^{2} α Rd (2 iPr . R + Rd α),

Handling Editor: J Svalastog

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.

ORCID iD

Ambreen Asfar Khan

References

Sleigh

. The biology of Cilia and Flagella: international series of monographs on pure and applied biology: zoology, vol. 12. Elsevier Science, 2016.

Lardner

Shack

. Cilia transport. Bull Math Biophys 1972; 34: 325–335.

Blake

. A model for the micro-structure in ciliated organisms. J Fluid Mech 1972; 55: 1–23.

Gueron

Liron

. Ciliary motion modeling, and dynamic multicilia interactions. Biophys J 1992; 63: 1045–1058.

Radhakrishnamacharya

Sharma

. Motion of a self-propelling micro-organism in a channel under peristalsis: effects of viscosity variation. Nonlinear Anal Model Control 2007; 12: 409–418.

Maqbool

Shaheen

Mann

. Exact solution of cilia induced flow of a Jeffrey fluid in an inclined tube. Springer Plus 2016; 5: 1379–1394.

Akbar

Khan

. Metachronal beating of cilia under the influence of Casson fluid and magnetic field. J Magn Magn Mater 2015; 378: 320–326.

Krishna Kumari

Vernekar

Ravi Kumar

. A study on the effect of peristalsis and cilia of MHD micropolar fluid flow through an inclined porous channel. J Appl Fluid Mech 2018; 11: 1321–1331.

Hasen

Abdulhadi

. MHD effect on peristaltic transport for rabinowitsch fluid through a porous medium in cilia channel. Iraqi J Sci 2020; 61: 1461–1472.

10.

Maccormack

. A numerical method for solving the equations of compressible viscous flow. AIAA J 1982; 20: 1275–1281.

11.

Guerra

Gustafsson

. A numerical method for incompressible and compressible flow problems with smooth solutions. J Comput Phys 1986; 63: 377–397.

12.

Mekheimer

Komy

Abdelsalam

. Simultaneous effects of magnetic field and space porosity on compressible Maxwell fluid transport induced by a surface acoustic wave in a microchannel. Chin Phys B 2013; 22: 124702–124711.

13.

Abdelsalam

Vafai

. Combined effects of magnetic field and rheological properties on the peristaltic flow of a compressible fluid in a microfluidic channel. Eur J Mech B Fluids 2017; 65: 398–411.

14.

Mekheimer

Abdel-Wahab

. Net annulus flow of a compressible viscous liquid with peristalsis. J Aerosp Eng 2012; 25: 660–669.

15.

Mekheimer

Abdel-Wahab

. Effect of wall compliance on compressible fluid transport induced by a surface acoustic wave in a microchannel. Numer Methods Partial Differ Equations 2011; 27: 621–636.

16.

Damianou

Georgiou

Moulitsas

. Combined effects of compressibility and slip in flows of a Herschel–Bulkley fluid. J Non-Newtonian Fluid Mech 2013; 193: 89–102.

17.

Saleem

Qaiser

Nadeem

. Physiological flow of biomedical compressible fluids inside a ciliated symmetric channel. Adv Mech Eng 2020; 12: 1–11.

18.

Khan

Rafaqat

. Effects of radiation and MHD on compressible Jeffrey fluid with peristalsis. J Therm Anal Calorim 2021; 143: 2775–2787.

19.

Hussain

Asghar

Hayat

, et al. Heat transfer analysis in peristaltic flow of MHD Jeffrey fluid with variable thermal conductivity. Appl Math Mech 2015; 36: 499–516.

20.

Hayat

Shehzad

Alsaedi

. Three-dimensional stretched flow of Jeffrey fluid with variable thermal conductivity and thermal radiation. Appl Math Mech 2013; 34: 823–832.

21.

Hayat

Javed

Imtiaz

, et al. Effect of Cattaneo-Christov heat flux on Jeffrey fluid flow with variable thermal conductivity. Results Phys 2018; 8: 341–351.

22.

Hamad

AbdEl-Gaied

Khan

. Thermal jump effects on boundary layer flow of a Jeffrey fluid near the stagnation point on a stretching/shrinking sheet with variable thermal conductivity. J Fluids 2013; 2013: Article ID 749271.

23.

Ahmed

Khan

Ahmad

. MHD swirling flow and heat transfer in Maxwell fluid driven by two coaxially rotating disks with variable thermal conductivity. Chin J Phys 2019; 60: 22–34.

24.

Hamid

Khan

. Unsteady mixed convective flow of Williamson nanofluid with heat transfer in the presence of variable thermal conductivity and magnetic field. J Mol Liq 2018; 260: 436–446.

25.

Khan

Bukhari

Marin

, et al. Effects of chemical reaction on third-grade MHD fluid flow under the influence of heat and mass transfer with variable reactive index. Heat Transfer Res 2019; 50: 1061–1080.

26.

Lahmar

Kezzar

Eid

, et al. Heat transfer of squeezing unsteady nanofluid flow under the effects of an inclined magnetic field and variable thermal conductivity. Physica A 2020; 540: 123138

27.

Bibi

Hayat

Farooq

, et al. Entropy generation analysis in peristaltic motion of Sisko material with variable viscosity and thermal conductivity. J Therm Anal Calorim 2019; 143: 363–375.

28.

Irfan

Khan

. Numerical analysis of unsteady 3D flow of Carreau nanofluid with variable thermal conductivity and heat source/sink. Results Phys 2017; 7: 3315–3324.

29.

Marinca

. Some exact solutions for MHD flow and heat transfer to modified second grade fluid with variable thermal conductivity in the presence of thermal radiation and heat generation/absorption. Comput Math Appl 2018; 76: 1515–1524.

MHD compressible fluid with variable thermal conductivity through a ciliated channel

Abstract

Keywords

Introduction

Problem formulation

Solution of the problem

Results and discussion

Conclusion

Footnotes

Appendix

Declaration of conflicting interests

Funding

ORCID iD

References