Heat Transfer Analysis of a Magneto-Bio-Fluid Transport with Variable Thermal Viscosity Through a Vertical Ciliated Channel

: We communicate the responses of various physiological ﬂuids containing hemoglobin and other ionic constituents when they propagate in the presence of an electromagnetic body force ﬁeld with the mechanisms of heat generation and conduction. A fully developed mixed convective ﬂow of a Newtonian ﬂuid takes place through a 2D vertical channel in the presence of an external magnetic ﬁeld acting in the direction normal to the ﬂow. The inner surface of the channel is carpeted with a thick mat of cilia, which propagates a sinusoidal metachronal wave travelling in the direction of ﬂow. Coupled, nonlinear governing Naiver-Stokes and temperature equations are simpliﬁed by utilizing the creeping ﬂow and long wavelength approximations. This enables us to formulate the exact analytical solution of the temperature distribution; whereas, the velocity distribution is evaluated from the momentum equations by using the Adomian decomposition method. In order to determine the pumping characteristics, the formulae of volume ﬂow rate and the pressure rise are also obtained. Trapping due to the ciliary system is highlighted by graphing the stream function. The ﬁndings of the present model have signiﬁcant outputs, which can be applicable in the physiological transport of human semen through the male reproduction system. that the present results are reasonably utilized in investigating the physiological ﬂows inside the human vessels and as well in contriving the magnetically actuated artiﬁcial cilia in lab-on-chip devices.


Introduction
Hair-like microscopic structures, which are recognized as cilia have featured roles in the transport processes of various physiological systems in the human bodies. The mechanical activity of cilia provides means for the locomotion of sperms and many other unicellular organisms and for the propulsion of materials and objects across the ciliated surface in multicellular animals. The role of cilia as a sensory organelle is recently admitted in addition to their capabilities in transportation, alimentation, and reproductions [1][2][3][4][5]. In the human body, cilia are found, for example, in lung trachea and bronchial tubes, where they remove mucus from the respiratory tract and the lungs; in fallopian tubes, where they support the movements of cervical mucus and the gamete cells in the oviduct; in the male reproductive system, where the cilia of ductile afferents propel sperm towards the epididymis; and in the brain as well, where the ependymal cilia facilitate the transport of cerebrospinal fluid at

Statement of the Problem
In order to describe the physical model mathematically, we designate a vertical channel of uniform width and coordinate system X, Y , in which the X-axis forms the central axis of the channel and the Y-axis is perpendicular to it, as shown Figure 1. A fully developed mixed convective flow of an electrically conducting Newtonian fluid is occurred in this channel in the presence of an electromagnetic body force. The fluid viscosity is thermally variable and obeys the Reynolds' law of viscosity. A ciliated epithelium covers the inner surface of the channel whose collective beatings structure a travelling wave known as metachronal wave over the ciliated surface. We delete the effects of the induced magnetic field due to a negligibly small magnetic Reynolds number. The effects of electric field are also ignored in this study, i.e., E=0. Utilizing these constraints, the body force (Lorentz force) due to the magnetic field, J × B, will take the form σ(V × B × B) = −σB 2 0 U, 0, 0 . The current density vector, J, is defined as J = σ(V × B) and the magnetic flux density vector B is given as B = (0, B 0 , 0), where V = U, V, 0 is the velocity vector and σ is the electrical conductivity of the fluid. The channel boundary is kept at uniform temperature T 1 , while at the centerline we use the well-known symmetry condition for temperature.

Statement of the Problem
In order to describe the physical model mathematically, we designate a vertical channel of uniform width and coordinate system ( ) , X Y   , in which the X  -axis forms the central axis of the channel and the Y  -axis is perpendicular to it, as shown Figure 1. A fully developed mixed convective flow of an electrically conducting Newtonian fluid is occurred in this channel in the presence of an electromagnetic body force. The fluid viscosity is thermally variable and obeys the Reynolds' law of viscosity. A ciliated epithelium covers the inner surface of the channel whose collective beatings structure a travelling wave known as metachronal wave over the ciliated surface. We delete the effects of the induced magnetic field due to a negligibly small magnetic Reynolds number. The effects of electric field are also ignored in this study, i.e., E = 0 . Utilizing these constraints, the body force (Lorentz force) due to the magnetic field, × J B , will take the form ( . The current density vector, J, is defined as σ and the magnetic flux density vector B is given as is the velocity vector and σ is the electrical conductivity of the fluid. The channel boundary is kept at uniform temperature 1 T , while at the centerline we use the well-known symmetry condition for temperature. The mathematical equations for the present flow are defined as follows: The mathematical equations for the present flow are defined as follows: In the above mathematical model, we have employed the following two assumptions: (i) in the momentum equations, we have used the well-known Boussinesq approximation so that the density variation is detained only in the buoyancy term; and (ii) in the energy equation, we have ignored the effects of viscous dissipation. Here, Q 0 represents the source term, µ(T) is the variable viscosity of the fluid and the remaining terms have their well-known meanings. It is pointed out in many studies [3,9] that the beating patterns of cilia are elliptical in nature whose parametric equations are expressed mathematically as In the above equations X and Y, respectively, describe the positions of the ciliary paths in the horizontal and vertical directions, d is the mean distance from the centerline of the channel, ε is the dimensionless measurement of cilia lengths, δ is the eccentricity measurement of the elliptical paths, χ 0 is the reference position. Furthermore, λ is the wavelength and c is the wave velocity.
With the help of Equations (5) and (6), we can determine the velocities of cilia tips along the horizontal and vertical directions as Using Equations (5) and (6) in Equations (7) and (8) and then solving for U and V, we arrive at The flow behaves as a steady state when we observe it in the wave frame (x, y) attached to the metachronal wave of cilia, whereas it is unsteady in the fixed frame ( X, Y). These two frames are For 2D steady flow, the velocity, pressure, and the temperature fields have the following patterns: Also, we redefine the boundary conditions as ∂u ∂y = 0, v = 0 and ∂T ∂y = 0 at y = 0, In order to normalize the governing equations, we allocate the following variables: Here, Re (modified Reynolds number), β (wave number), and ε (dimensionless measure with respect to the cilia length). We will now employ the above transformations along with Equations (12) and (13) into Equations (1)-(6) and then exploit the creeping flow (i.e., Re → 0) and the long wavelength (i.e., β 1) assumptions. Finally, the simplified system of equations we need to solve is written-after dropping the asterisks-as ∂u ∂x ∂p ∂y = 0 (18) The boundary conditions in dimensionless variables are manipulated as ∂u ∂y = 0, v = 0 and ∂θ ∂y = 0 at y = 0 In the above equations, g t (Grashof number), M (Hartmann number), and S t (constant heat source/sink parameter). In the present investigation, we define µ(θ) = e 1−αθ , µ(θ) = 1 − αθ for α 1 (22) which reduces to the constant viscosity model when α = 0. Here, α represents the viscosity variation parameter. Note that this choice of µ(θ) has been justified physiologically by several authors in their recent articles [25][26][27].

Solution Methodology
Equation (18) points out that p = p(x), so we can replace ∂p ∂x with dp dx in Equation (17). The temperature field of the fluid is obtained by solving the energy Equation (19), subjected to the relevant boundary conditions on temperature as Upon making use of Equations (22) and (23), the momentum Equation (17) will take the form Obviously, Equation (24) does not admit an exact closed form solution. In order to solve it, we use the well-known Adomian decomposition method. Accordingly, we rewrite it in the operator form in the following way [27][28][29][30]: where L yy = ∂ 2 ∂y 2 . Since L yy is a second order differential operator, L −1 yy will be a twofold integral operator defined by Applying the inverse operator L −1 yy , Equation (25) can be written as In the above equation, the functions C(x) and D(x) can be formulated from the boundary conditions (20) and (21). According to the well-known Adomian decomposition method, we can write Now, we substitute Equation (28) into Equation (27) and then decompose u as and By utilizing the recursive relation given in Equation (30), we will compute the components u 1 , u 2 . Consequently, the three terms ADM solution for the velocity field is found to be u = a 11 y 6 + a 12 + a 13 dp dx y 4 + a 14 + a 15 dp dx We define the volume flow rate q in the wave frame as [9,12], In the fixed frame, we construct which can be used to define Q (time-mean volume flow rate) in the following way [9,12], Upon making use of Equation (31) in Equation (32) and then solving the resultant integral, we obtain the expression for the pressure gradient as Integrating Equation (35) over the interval [0, 1], we get which yields the expression for the pressure rise. Upon solving it for Q, we obtain where In order to visualize the stream line patterns of ciliary induced flow, we compute ψ (the stream function) by introducing into Equations (16) and (31). Here, we have used the symbolic software Maple to create the graphical representations of the streamlines, which would be helpful in analyzing the trapping phenomena in the present model.

Results and Discussions
This section provides physical interpretations of the graphical results computed for the velocity (u), the pressure gradient (dp/dx), pressure difference across one wavelength (∆p), volume flow rate (Q) and the 2D stream function (ψ) under changing values of the emerging parameters. These parameters are the Hartmann number (M), which represents the effect of magnetic field, the viscosity variation parameter due to temperature (α), which represents the fluid viscosity, the Grashof number (g t ), which represents the buoyancy effects, and the heat source (sink) parameter (S t ), which represents the heat transfer.

Velocity Distribution
The graphs in Figure 2a represent that the fluid velocity is diminished in the central region (i.e., in the interval [−0.4, 0.4]) and is amplified near the channel walls. Physically, this phenomenon is associated with the magnetic body force i.e., the Lorentz force, which provides a resistance to the fluid motion [12,24]. Figure 2b shows that the fluid velocity increases with an increase in α in the center of the channel, whereas it decreases slightly near the boundaries. Figure 2c,d explain that the behavior of the axial velocity is not uniform in whole region under the impacts of g t and S t . When we increase the values of the Grashof number (g t is taken as negative, zero and positive) and the heat source sink parameter (S t is taken as negative, zero and positive), the velocity rises in the center of the channel, whereas it decreases near the walls. The negative values in the axial velocity are noted for the negative values of the Grashof number and the heat sink parameter, which indicates the existence of flow reversal or backward flow in the human vessels.
representations of the streamlines, which would be helpful in analyzing the trapping phenomena in the present model.

Results and Discussions
This section provides physical interpretations of the graphical results computed for the velocity ( u), the pressure gradient ( d p d x ), pressure difference across one wavelength ( p Δ ), volume flow rate ( Q ) and the 2D stream function (ψ ) under changing values of the emerging parameters. These parameters are the Hartmann number ( M ), which represents the effect of magnetic field, the viscosity variation parameter due to temperature ( α ), which represents the fluid viscosity, the Grashof number ( t g ), which represents the buoyancy effects, and the heat source (sink) parameter ( t S ), which represents the heat transfer.

Velocity Distribution
The graphs in Figure 2a  ) and is amplified near the channel walls. Physically, this phenomenon is associated with the magnetic body force i.e., the Lorentz force, which provides a resistance to the fluid motion [12,24]. Figure 2b shows that the fluid velocity increases with an increase in α in the center of the channel, whereas it decreases slightly near the boundaries.   Figure 3b describes that the pressure will show non-uniform behavior over the entire region when the viscosity of the fluid depends upon its temperature, however, it may be not very significant. The graphs in Figure 3c,d indicate that there is a significant rise in when we increase the values of t g and t S . Thus, the heat transfer effect supports the development of the larger pressure gradient in the channel. These graphs clearly indicate a sinusoidal response of the pressure gradient, which endorses the correctness of our analysis.

Pressure Distribution
Figure 3a-d depict the behaviors of dp/dx (the axial pressure gradient) versus x under the effects of parameters M, α, g t , and S t . Figure 3a indicates that the pressure gradient is reduced when we increase the values of M, i.e., the amount of the pressure gradient produced for MHD Newtonian fluid is smaller as compared to the non MHD Newtonian fluid in the channel, where they are moving. Figure 3b describes that the pressure will show non-uniform behavior over the entire region when the viscosity of the fluid depends upon its temperature, however, it may be not very significant. The graphs in Figure 3c,d indicate that there is a significant rise in when we increase the values of g t and S t . Thus, the heat transfer effect supports the development of the larger pressure gradient in the channel. These graphs clearly indicate a sinusoidal response of the pressure gradient, which endorses the correctness of our analysis. Figure 4a-d point out the variations in ∆p (the pressure difference across one wavelength) versus Q (the volume flow rate) for various values of M, α, g t , and S t . In Figure 4a, we see that the pressure difference declines significantly for higher values of the Hartmann number. This behavior highlights the pumping characteristics for MHD fluid, i.e., the pumping rate for MHD Newtonian fluid is smaller as compared to the non MHD Newtonian fluid. Figure 4b describes the effects of α on ∆p versus Q. We note that the impacts of α on ∆p substantially depend upon the amount of the flow rate. When Q < 1, the effect of α on the pressure difference is not much indicative. However, when Q ≥ 1, α exerts potent effects on the pumping rate. Figure 4c,d explain that there is substantial growth in ∆p when we increase the values of g t and S t . Thus, the heat transfer effects promote the pumping rate.

Volume Flow Rate
The graphs in Figure 5a depict the impacts of the Hartmann number M on Q versus ε , the dimensionless cilia height estimation. Volume flow rate decreases substantially due to increasing values of the Hartmann number. From these observations, we can conclude that the presence of a magnetic field creates MHD drag force that will cause the velocity field across the channel length to decrease, and hence reduce the flow rate. This characteristic of the magnetic field has many applications in physiological flows; one of these applications is the use of magnetic field in controlling the excessive bleeding during critical surgeries. From Figure 5b, we notice that the volume flow rate for the high viscosity (temperature dependent) fluids is larger as compared to small viscosity fluids. Moreover, these graphs indicate that the cilia heights hold down the flow rate in the channel flow, i.e., they act as the controlling parameter for the fluid motion. Figure 5c,d depict that the increasing

Volume Flow Rate
The graphs in Figure 5a depict the impacts of the Hartmann number M on Q versus ε, the dimensionless cilia height estimation. Volume flow rate decreases substantially due to increasing values of the Hartmann number. From these observations, we can conclude that the presence of a magnetic field creates MHD drag force that will cause the velocity field across the channel length to decrease, and hence reduce the flow rate. This characteristic of the magnetic field has many applications in physiological flows; one of these applications is the use of magnetic field in controlling the excessive bleeding during critical surgeries. From Figure 5b, we notice that the volume flow rate for the high viscosity (temperature dependent) fluids is larger as compared to small viscosity fluids. Moreover, these graphs indicate that the cilia heights hold down the flow rate in the channel flow, i.e., they act as the controlling parameter for the fluid motion. Figure 5c,d depict that the increasing values of the thermal Grashof number and the heat source (sink) parameters produce accelerating effects on the volume flow rate. These observations endorse the actual physical situations and are according to the results of previous studies mentioned in the literature [24][25][26][27][28][29][30]. Furthermore, it can be seen that the variation in the flow rate is not same for negative values of g t and S t as compared to their positive values, i.e., the flow rate increases with the cilia heights when we take negative values of these parameters. Thus, the heat transfer effects support the volumetric flux in the physiological conduits. Figure 6a,b describe the impacts of M (Hartmann number) and ε(the cilia length parameter) on Q (the volume flow rate) versus δ (the eccentricity parameter). We note that the volume flow rate decreases with the Hartmann number and the cilia heights, whereas it increases with the eccentricity parameter. values of the thermal Grashof number and the heat source (sink) parameters produce accelerating effects on the volume flow rate. These observations endorse the actual physical situations and are according to the results of previous studies mentioned in the literature [24][25][26][27][28][29][30]. Furthermore, it can be seen that the variation in the flow rate is not same for negative values of t g and t S as compared to their positive values, i.e., the flow rate increases with the cilia heights when we take negative values of these parameters. Thus, the heat transfer effects support the volumetric flux in the physiological conduits. Figure 6a,b describe the impacts of M (Hartmann number) and ε (the cilia length parameter) on Q (the volume flow rate) versus δ (the eccentricity parameter). We note that the volume flow rate decreases with the Hartmann number and the cilia heights, whereas it increases with the eccentricity parameter.

Trapping Phenomenon
In ciliary-induced flows, trapping will exist due to metachronal waves of cilia and is recognized as a natural mechanism of physiological flows in the human body. It arises due to the circulation of streamlines at larger amount of volume flow rate. The streamlines are split to trap a bolus of the flowing material that will move as a whole with the wave speed along its direction. Figures 7-10 explain the effects of key parameters M , α , t g , and t S on the 2D stream function ψ . These graphs will explain the trapping behavior and circulation of the flow characteristics. The impacts of the magnetic field on the streamlines are displayed in Figure 7, which shows that the circulation as well as the size of trapped boluses gradually decrease due to increasing values of M (the Hartmann number). From Figure 8, one can observe that the rise in the viscosity variation parameter markedly boosts the bolus growth and circulations in the streamlines. Figure 9 explains that the trapped boluses increase in size as well as in numbers with the increasing magnitude of t g , which indicates that the buoyancy force parameter promotes the formation of boluses in the cilia-driven flows. The streamlines in Figure 10 explain that the size as well as the numbers of the trapped boluses increase with the increasing magnitude of t S (heat source parameter). This confirms the previous observations that the occurrence of trapping will be enhanced when we increase the effects of heat transfer [24][25][26][27][28][29].

Trapping Phenomenon
In ciliary-induced flows, trapping will exist due to metachronal waves of cilia and is recognized as a natural mechanism of physiological flows in the human body. It arises due to the circulation of streamlines at larger amount of volume flow rate. The streamlines are split to trap a bolus of the flowing material that will move as a whole with the wave speed along its direction. Figures 7-10 explain the effects of key parameters M, α, g t , and S t on the 2D stream function ψ. These graphs will explain the trapping behavior and circulation of the flow characteristics. The impacts of the magnetic field on the streamlines are displayed in Figure 7, which shows that the circulation as well as the size of trapped boluses gradually decrease due to increasing values of M (the Hartmann number). From Figure 8, one can observe that the rise in the viscosity variation parameter markedly boosts the bolus growth and circulations in the streamlines. Figure 9 explains that the trapped boluses increase in size as well as in numbers with the increasing magnitude of g t , which indicates that the buoyancy force parameter promotes the formation of boluses in the cilia-driven flows. The streamlines in Figure 10 explain that the size as well as the numbers of the trapped boluses increase with the increasing magnitude of S t (heat source parameter). This confirms the previous observations that the occurrence of trapping will be enhanced when we increase the effects of heat transfer [24][25][26][27][28][29]. (a) (b) Figure 9. Streamlines for (a) g t = 0, (b) g t = 3, (c) g t = 5. The other parameters are fixed as δ = 0.2,

Application: Physiological Transport of Semen
The present study has a close relevance to the physiological flow of human semen through the efferent ducts of the male reproductive tract. It is now widely accepted that ciliary movements generating the metachronal waves have a significant role in the propulsion of semen through the

Application: Physiological Transport of Semen
The present study has a close relevance to the physiological flow of human semen through the efferent ducts of the male reproductive tract. It is now widely accepted that ciliary movements generating the metachronal waves have a significant role in the propulsion of semen through the ductile efferentes. Lardner and Shack [9] conducted experimental studies for various animals and found that the volume flow rate of human semen in the ductile efferent is approximately equal to 0.006 mL/h. In 1972, Lardner and Shack also developed a mathematical model describing the flow of a Newtonian viscous fluid through a 2D symmetric channel whose inner walls are ciliated. They calculated the values of non-dimensional volume flow rate and dimensional volume flow rate of the human semen by considering it to behave as a Newtonian viscous fluid to be 0.022 and 0.00012 mL/h, respectively. While calculating these values, the following data were used: In this study, we aim to investigate the impacts of heat transfer and applied magnetic field on the volume flow rate Q of human semen through a 2D channel. Table 1 describes the numerical values of Q and its dimensional counterpart Q * = (Q × πd 2 c) for various values of the Hartmann number M, the viscosity variation parameter α, the Grashof number (the buoyancy force parameter) g t and the heat source (sink) parameter S t . From Table 1, we note that the flow rate increases with the increasing values of α, g t and S t , whereas it decreases with the increasing values of the Hartmann number. We also observed that the volume flow rate for the conducting fluid is lower in magnitude as compared to the non-conducting fluid and this trend has been noticed in many studies [24][25][26][27][28][29]. Moreover, it is interesting to mention that the results of the present model are similar to those obtained in [9] for M = 0, α = 0, g t = 0 and S t = 0. Table 1. Numerical results of Q and Q * = (Q × 2πd 2 c), when δ = 0.1, ε = 0.1, β = 1.0, and ∆p = 0. These values do not agree with the estimated values of the volume flow rate of human semen (i.e., 0.006 mL/h). This suggests that a considerable improvement is needed in the model of Lardner and Shack [9]. However, when we choose M = 0.1, α = 0.03, g t = 0.9, S t = 0.9, we obtain Q = 1.07876, and Q * = 0.00588. These values are comparatively closer to the experimentally estimated values. We hope that the present results are reasonably utilized in investigating the physiological flows inside the human vessels and as well in contriving the magnetically actuated artificial cilia in lab-on-chip devices.

Conclusions
An analysis of heat transfer effects on the ciliary pumping system in propelling biological fluids like blood, semantic is carried out under the influence of magnetic body force. The heat transfer features are characterized by the adding external heat source (sink) elements. The fluid viscosity depends upon the temperature according to the Reynolds law. In order to provide the generality and correctness of the present model, we would like to say that our model contains the results of Lardner and Shack [9] for M = 0, α = 0, g t = 0, and S t = 0; and the results of Siddiqui et al. [12] for α = 0, g t = 0, and S t = 0. The numerical results of the volume flow rate obtained in this study are in GOOD agreement with the experimental values of the flow rate of the human semen estimated by Lardner and Shack [9]. This study reveals that the applied magnetic field exhibits strong decelerating consequences on the volume flow rate. Thus, the magnetic field can be used to achieve a better control on the flow rate. This phenomenon has useful applications in medical sciences and in the design of artificial cilia for the drug delivery system. Another important observation of this study is that the volume flow rate can be significantly accelerated due to increasing the heat transfer effects (Grashof number and heat source parameter), i.e., heat transfer assists the ciliary induced pumping. We also note that the magnetic field can control the development of trapping. a 21 = −4a 14 a 16 − 7a 12 a 16 a 22 = 5a 14 a 16 − 6a 12 a 23 = 14S t g t a 14 + 28Ma 13 a 24 = 210M 2 a 16 + 128Ma 13 a 24 = 210M 2 a 16 + 128Ma 13