删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

FEM for Blood-Based SWCNTs Flow Through a Circular Cylinder in a Porous Medium with Electromagnetic

本站小编 Free考研考试/2022-01-02

Mohamed R. Eid,1,2,?, A. F. Al-Hossainy3,4, M. Sh. Zoromba5,6 Department of Mathematics, Faculty of Science, Northern Border University, Arar, 1321, Saudi Arabia
Department of Mathematics, Faculty of Science, New Valley University, Al-Kharga, Al-Wadi Al-Gadid, 72511, Egypt
Department of Chemistry, Faculty of Science, New Valley University, Al-Kharga, Al-Wadi Al-Gadid, 72511, Egypt
Department of Chemistry, Faculty of Science, Northern Border University, Arar, 1321, Saudi Arabia
Department of Chemical and Materials Engineering, King Abdulaziz University, Rabigh 21911, Saudi Arabia
Department of Chemistry, Faculty of Science, Port Said University, 42521 Port-Said, Egypt

Corresponding authors: ?E-mail:m_r_eid@yahoo.com

Received:2019-05-22Online:2019-12-1
Fund supported:*Supported by Deanship of Scientific Research, Northern Border University under Gran.No.SCI-2018-3-9-F-7614.KSA


Abstract
This work aims to study magnetohydrodynamic flow through a circular cylinder in a horizontal position of SWCNTs in blood as a base fluid in the existence of non-linear thermal radiation and heat source/sink. Three kinds of nanoparticles shapes are considered. The study is employed the finite element technique to explore and enhance the influences of essential parameters on temperature profiles and is debated the heat transport within blood injects with SWCNTs and exposes to electromagnetic radiation. The treatment with thermal analysis and heat transfer rate being a better substitute more than surgery and chemotherapy for cancer therapy. Utilizing of nanoparticles thermal features is a mounting area of nanomedicine field because of the probable for purposeful demolition of cancer cells. This remedy is relied on many parameters, including nanofluid thermal conductivity, nanoparticles volume fraction, thermal radiation and power and heat source. The numerical solutions for flow and heat transfer features are assessed for diverse governing parameters values. The obtained results are substantiated against the relevant numerical results in the published researches. Results show that both flow velocity and temperature increase for larger values of thermal radiation, heat source and SWCNTs volume fraction with lamina and cylinder shapes. Also, spherical shape of SWCNTs occurs high disturbances in velocity and temperature distribution in the case of cooled cylinder.
Keywords: finite element method;heat transfer;blood;SWCNTs;electromagnetic radiation


PDF (5156KB)MetadataMetricsRelated articlesExportEndNote|Ris|BibtexFavorite
Cite this article
Mohamed R. Eid, A. F. Al-Hossainy, M. Sh. Zoromba. FEM for Blood-Based SWCNTs Flow Through a Circular Cylinder in a Porous Medium with Electromagnetic Radiation*. [J], 2019, 71(12): 1425-1434 doi:10.1088/0253-6102/71/12/1425

Nomenclature



Greek symbols



Subscripts



Abbreviations



1 Introduction

The model of nano-liquid was a new category of liquids with different heat transference features that can be architected by suspended mineral nanoparticles in classic heat transference fluids. The successor fluid acknowledged as nanofluid is potential to demonstrate superscript thermal conductivities versus to those of existing utilized heat transference fluids. The nanoparticles employed in nanofluids are manufactured of metals, like alumina, copper, metal oxides or non-metals (CNTs, graphite) and commonly ethylene glycol, water and blood are employed as basis fluid. CNTs have convened a fundamental role in nanostructures because of their singular mechanical structural and electronical features. CNTs have exceptional conductivity that helps them to constitute a grid of conductive tubes. CNTs have been depreciated as thermal boundary substances for thermal maintain. Highly conductivity CNTs have attracted considerable interest from scientist. CNTs utilized in nanofluids are generally of two kinds, exclusively single-walls CNTs (SWCNTs) and multiple-walls CNTs (MWCNTs). The CNTs width varies from 1 nm to 100 nm and CNTs length is in micrometer.

The research topic of nano-liquids has receipted considerable interest universality. Originally high thermal conductivity is a basic constraint in enhancement energy-competence heat transfer liquids that are desired for ultrahigh-performance heating. A very tiny magnitude of guest nanoparticles, when suspended stably in host liquids, can offer enhancements in the thermal characterizations of the host liquids. Buongiorno[1] suggested a form for nanofluids convective transport integrating the impacts of Brownian motion and thermophoresis diffusion. This form was implemented to the problem of Rayleigh-Benard by Refs. [2--3] and they demonstrated that both Brownian motion and thermophoresis influences perform to transit-diffusion conditions that are in certain methods corresponding to the commonplace Soret and Dufour diffusion that established with a dual liquid. Nanofluids in one and two-phase models in diverse situations with various parameter impacts have been examined by Refs. [4--16].

Heat transference performing of CNT nanofluids flow across a horizontal pipe was achieved by Ref. [17]. They existed significant enhancement of heat transference by convection and remarked that the betterment based on Reynolds number and CNTS volume fraction. Meyer $et al.$[18] provided experimentally the heat transfer improvements of MWCNTs flux via a straight horizontal pipe. The thermal conductivity promotion of engine oil and ethylene glycol with MWCNTs is indicated by Ref. [19]. They existed that ethylene glycol-CNTs pendants have the largest thermal conductivities when it compares with ethylene glycol deprived of CNTs. Wang $et al.$[20] discussed the flow and heat transfer of nanofluids with pressure drop in a horizontal circular tubing including CNTs experimentally. They revealed that nanofluids at lowland concentration reinforce the heat transfer with bit additional sanction in pumping power. The experimental examination of convective heat transfer properties of byproduct refrigerants-depend on CNTs nanofluid in a ducted heat exchanger considered by Ref. [21].

Khan $et al.$[22] mentioned heat transfer features of CNTs over a flattened plate with boundary conditions of uniformly heat flux and Navier slip numerically, they utilized single and multi-walls CNTs and MWCNTs with three various norm fluids. Al-Haq $et al.$[23] examined water-based flow with SWCNTs and MWCNTs, they studied a superior heat transfer rate for single and multiple walls CNTs. Kandasamy $et al.$[24] evidenced numerically the chemical react influence on copper, Al$_{2}$O$_{3}$, and SWNTs nanofluid flowing with the conditions of slip by utilizing RKF technique based on shooting approach and the chemical reaction effect on SWCNTs, Al$_{2}$O$_{3}$ and Cu-nanoparticles on convective mass transfer in having of water as a base fluid by utilizing numerical regime also deliberated by Ref. [25]. Lately, analytical solution of heat transfer improvement in Maxwell nanofluids flow by natural convection with the two types of CNTs past a vertical plate utilizing Laplace transform method provided by Ref. [26]. Casson model[27] is determined as the majority favored rheological form for characterizing blood flow of human body.\wen{ 28-29} Eid $et al.$[11] elaborated thermal radiation impact on Sisko-model flow of gold nanoparticles in human blood as base fluid.

So far there is no research in the literature on the effect of electromagnetic radiation on SWCNTs-blood nanofluid mixed convection embedded in a porous material with heat source/sink. The presence of electromagnetic radiation on nanofluid usually has the effect of inhibiting the development of destabilization. This paper designed to study the blood flow involving SWCNTs is inspected in the existence of magnetohydrodynamic, thermal radiation, heat source/sink and porosity parameter. The NPs shapes are taken in account in this study. Human blood is adopted as base fluid with SWCNTs insider.

2 Flow Analysis

Consider the issue of steady boundary-layer flow of mixed convection over a cooled/heated circular cylinder placed horizontally of radius $r$ in a porous material saturated with a nanofluid in the existence of applied magnetic field. It is supposed that the external flow velocity is ${\bar{u}}_e(\bar{x})$, where $\bar{x}$ is the coordinate scaled together the cylinder surface starting from the bottom stagnation-point. Also assumed that the constant temperature of the surface of the cylinder is $T_w$, whereas the value of constant temperature ambient is $T_{\infty{}}$, where $T_w>T_{\infty{}}$ for the assisting flow (a heated cylinder) and $T_w<T_{\infty{}}$ for the opposing flow (a cooled cylinder), respectively. In conserving with the Darcy law, and taking the nanofluid model, the Boussinesq approximation is utilized as homogeneous and thermally equilibrium in the porous material is supposed. The main equations can be scripted in Cartesian axis\textit{x} and \textit{y} as (see Refs. [30--31]):

$$ \frac{\partial{}\bar{u}}{\partial{}\bar{x}}+\frac{\partial{}\bar{v}}{\partial{}\bar{y}}=0, $$
$$ \frac{{\mu{}}_{nf}}{K}\bar{u}\!=\!-\!\frac{\partial{}\bar{p}}{\partial{}\bar{x}}\!+\!{\left(\rho{}\beta{}\right)}_{nf}g\left(T\!-\!T_{\infty{}}\right)\sin{\left(\frac{\bar{x}}{r}\right)\!-\!{\sigma{}}_{nf}B_0^2}\bar{u}, $$
$$ \frac{{\mu{}}_{nf}}{K}\bar{v}=-\frac{\partial{}\bar{p}}{\partial{}\bar{y}}-{\left(\rho{}\beta{}\right)}_{nf}g\left(T-T_{\infty{}}\right)\cos{\left(\frac{\bar{x}}{r}\right)}, $$
$$ {\left(\rho{}c_p\right)}_{nf}\left(\bar{u}\frac{\partial{}T}{\partial{}\bar{x}}+\bar{v}\frac{\partial{}T}{\partial{}\bar{y}}\right)=k_{nf}\left(\frac{{\partial{}}^2T}{\partial{}{\bar{x}}^2}+\frac{{\partial{}}^2T}{\partial{}{\bar{y}}^2}\right)\nonumber\\ \quad +Q_0(T-T_{\infty{}})-\frac{\partial{}q_r}{\partial{}\bar{y}}, $$

New window|Download| PPT slide
Fig. 1(Color online) Flow configuration model.



with the associated boundary conditions:

$$ \bar{v}=0,\quad T=T_w\ ,\quad {\rm at}\ \bar{y}=0,\quad 0\leq{}\bar{x}\leq{}\pi{},\ \\ \bar{u}={\bar{u}}_e,\quad T=T_{\infty{}}\ ,\quad {\rm as}\ \bar{y}\rightarrow{}\infty{},\quad 0\leq{}\bar{x}\leq{}\pi{}, $$
where $\bar{u}$ and $\bar{v}$ are the components of velocity in the $\bar{x}$ and $\bar{y}$-axis, respectively, where $\bar{y}$ is the axis scaled in the orientation normal to the cylinder surface, $T$ is the temperature of nanofluid,$p$ is pressure, $g$ is the gravity acceleration volume, ${\sigma{}}_{nf}$ is the electric conductivity, and $B_0$ is the uniform magnetic field strength. The associated boundary conditions (5) characterize flow velocity, temperature, and NPs concentration values of at the boundary of the inner and the outer circular cylinder boundary. Properties ${\mu{}}_{nf}$ and ${\rho{}}_{nf}$ are the nanofluid viscosity and density, and $K$ is the modified porous medium permeability. Also, ${\beta{}}_{nf}$ and ${\alpha{}}_{nf}$ are the thermal expansion and thermal diffusivity coefficients of the saturation porous medium, respectively, which are identified as:[7]

$$ {\rho{}}_{nf}=\left(1-\varphi{}\right){\rho{}}_b+\varphi{}{\rho{}}_{\rm CNT},\nonumber\\ {\mu{}}_{nf}=\frac{{\mu{}}_b}{{\left(1-\varphi{}\right)}^{2.5}},\quad {\alpha{}}_{nf}=\frac{k_{nf}}{{\left(\rho{}c_p\right)}_{nf}},\quad\nonumber\\ {\left(\rho{}c_p\right)}_{nf}=\left(1-\varphi{}\right){\left(\rho{}c_p\right)}_b+{\varphi{}\left(\rho{}c_p\right)}_{\rm CNT},\quad\nonumber\\ {\sigma{}}_{nf}=\left(1-\varphi{}\right){\sigma{}}_b+{\varphi{}\sigma{}}_{\rm CNT},\quad\nonumber\\ {\left(\rho\beta{}\right)}_{nf}=\left(1-\varphi{}\right){\left(\rho{}\beta{}\right)}_b+{\varphi{}\left(\rho{}\beta{}\right)}_{\rm CNT},\quad\nonumber\\ \frac{k_{nf}}{k_b}=\frac{k_{\rm CNT}+\left(m-1\right)k_b-\left(m-1\right)\varphi{}(k_b-k_{\rm CNT})}{k_{\rm CNT}+\left(m-1\right)k_b+\varphi{}(k_b-k_{\rm CNT})}, $$
where $\varphi{}$ is nanoparticles volume fraction, $m$ is nanoparticles shape parameter, ${\mu{}}_b$ is the blood dynamic viscosity, ${\rho{}}_b$ and ${\rho{}}_{\rm CNT}$ are the densities of the blood and SWCNTs nanoparticle, ${\sigma{}}_b$ and ${\sigma{}}_{\rm CNT}$ are electric conductivities of blood and SWCNTs, $k_b$ and $k_{\rm CNT}$ are thermal conductivities of the blood and SWCNTs, ${\beta{}}_b$ and ${\beta{}}_{\rm CNT}$ arethermal expansion coefficients of the blood and nanoparticle, ${\left(\rho{}c_p\right)}_b$ and ${\left(\rho{}c_p\right)}_{\rm CNT}$ are the specific heat factors of the blood and SWNTs, respectively.

By introducing the following non-dimensional variables:

$$ x=\frac{\bar{x}}{r}\ ,\quad\ y=\frac{\bar{y}}{r}\ Pe^{{1}/{2}},\quad u=\frac{\bar{u}}{U_{\infty{}}}\ ,\\ v=\frac{\bar{v}}{U_{\infty{}}}\ Pe^{{1}/{2}},\quad\ \theta{}=\frac{T-T_{\infty{}}}{T_w-T_{\infty{}}}\ ,\quad u_e=\frac{{\bar{u}}_e}{U_{\infty{}}}\ ,\ \nonumber\\ p=\frac{\bar{p}-p_{\infty{}}}{\left({\rho{}}_{nf}U_{\infty{}}^2\right)}\ ,\quad Pe=\frac{U_{\infty{}}r}{{\alpha{}}_f}. $$
Excluding $\bar{p}$ from Eqs. (2)--(3) by cross derivation, replacing Eqs. (6) and (7) into the result equation and deeming the boundary-layer calculation $(Pe\rightarrow{}\infty{})$, where $Pe$ is the Péclet number. To get the subsequent boundary-layer equations for the current problem:

$$ \frac{\partial{}u}{\partial{}x}+\frac{\partial{}v}{\partial{}y}=0,$$
$$ \Big[\frac{1}{{\left(1-\varphi{}\right)}^{2.5}}-A_2M\Big]\frac{\partial{}u}{\partial{}y}=A_1\lambda{}\frac{\partial{}\theta{}}{\partial{}y}\ {\rm sin} x,$$
or this equation can be written in the form

$ \Big[\frac{1}{{\left(1-\varphi{}\right)}^{2.5}}-A_2M\Big]\left(u-u_e\right)=A_1\lambda{}\theta{}\ {\rm sin} x,$
where

$ A_1=\left(1-\varphi{}\right)+\varphi{}\frac{{\left(\rho{}\beta{}\right)}_{\rm CNT}}{{\left(\rho{}\beta{}\right)}_b},\\ A_2=\left(1-\varphi{}\right)+\varphi{}\frac{{\sigma{}}_{\rm CNT}}{{\sigma{}}_b},$

$M={{\sigma{}}_b{KB}_0^2}/{{\mu{}}_b}$ is the magnetic field parameter.

$\lambda{}={Ra_r}/{Pe}$ is the parameter of mixed convection.

$Ra_r={Kgr{\beta{}}_b\left(T_w-T_{\infty{}}\right)}/{{\alpha{}}_b{\upsilon{}}_b}$ is the Rayleigh number for the porous medium.

$ \because{}\;\;{\left(\rho{}c_p\right)}_{nf}\left(\bar{u}\frac{\partial{}T}{\partial{}\bar{x}}+\bar{v}\frac{\partial{}T}{\partial{}\bar{y}}\right)=k_{nf}\left(\frac{{\partial{}}^2T}{\partial{}{\bar{x}}^2}+\frac{{\partial{}}^2T}{\partial{}{\bar{y}}^2}\right)\\ \quad +\,Q_0(T-T_{\infty{}})-\frac{\partial{}q_r}{\partial{}\bar{y}},$

where $q_r=-({4{\sigma{}}^*}/{3k^*})({\partial{}T^4}/{\partial{}y})$ is the radiative heat flux.

${\sigma{}}^*$ and $k^*$ are the Stefan-Boltzmann factor and the average absorption coefficient, respectively.

$$ T^4\cong{}4T_{\infty{}}^3T-3T_{\infty{}}^4,\\ u\frac{\partial{}\theta{}}{\partial{}x}+v\frac{\partial{}\theta{}}{\partial{}y}=\frac{{\alpha{}}_{nf}}{{\alpha{}}_b}\left(1+\frac{4}{3}R_d\right)\frac{{\partial{}}^2\theta{}}{\partial{}y^2}+s\theta{}. $$
From the aforementioned equation it is viewed that the radiation effect is to reinforce the thermal diffusivity.

where $R_d={4{\sigma{}}^*T_{\infty{}}^3}/{k^*k_{nf}}$ is the radiation parameter. It is worth mentioning here that the normal solve for heat equation,Eq. (10) without non-linear radiation impact can be gained from the foregoing equation which diminishes to $u({{\partial{}\theta{}}/{\partial{}x}})+v({{\partial{}\theta{}}/{\partial{}y}})=({{{\alpha{}}_{nf}}/{{\alpha{}}_b}})({{{\partial{}}^2\theta{}}/{\partial{}y^2}})$ and eliminates heat source.

The subject boundary conditions (5) become:

$$ v=0; \quad \theta{}=1\ {\rm at}\ y=0; \quad 0\leq{}x\leq{}\pi{},\ \ \ \ \ \ \\ u\rightarrow{}u_e; \quad\ \theta{}\rightarrow{}0\ {\rm as}\ y\rightarrow{}\infty{}; \quad\ 0\leq{}x\leq{}\pi{}. $$
Here we assume that $u_e=\sin x$ (see Refs. [30--31]). Incorporating Eq. (9a) and utilizing boundary-conditions (11), Eq. (9b) can be written as:

$$ \Big[\frac{1}{{\left(1-\varphi{}\right)}^{2.5}}-A_2M\Big]\left(u-\sin x\right)=A_1\lambda{}\theta{}\sin x. $$
To solve Eqs. (8)--(10) and (12) with the boundary-conditions (11) introduce the following variables:

$$ \psi{}=xf\left(x,\eta{}\right),\quad \ \theta{}=\theta{}\left(x,\eta{}\right),\quad\ \eta{}=y. $$
The equation of continuity is automatically achieved by identifying a stream function $\psi{}(x,y)$, defines as: $u={\partial{}\psi{}}/{\partial{}y}$ and $v=-({{\partial{}\psi{}}/{\partial{}x}})$. Replacing the variables (13) into Eqs. (10) and (12), we obtain:

$ \frac{k_{nf}/k_b}{A_3}{\left(1+\frac{4}{3}R_d\right)\theta{}}{'{}'}+f{\theta{}}{'}+s\theta{}=x\left(f'\frac{\partial{}\theta{}}{\partial{}x}-{\theta{}}{'}\frac{\partial{}f}{\partial{}x}\right),$
$$ \Big[\frac{1}{{\left(1-\varphi{}\right)}^{2.5}}-A_2M\Big]f'=\Big[\frac{1}{{\left(1-\varphi{}\right)}^{2.5}}-A_2M{+A}_1\lambda{}\theta{}\Big]\frac{\sin x}{x}, $$
where

$ A_3=\left(1-\varphi{}\right)+\varphi{}\Big[{{\left(\rho{}c_p\right)}_{\rm CNT}}/{{\left(\rho{}c_p\right)}_b}\Big]\,,$

with the boundary conditions:

$$ f\left(0\right)=0, \quad\ \theta{}\left(0\right)=1,\\ f^{\prime}\left(\infty{}\right)\rightarrow{}\frac{u_e}{x}=\frac{\sin x}{x},\quad\ \theta{}(\infty{})\rightarrow{}0. $$
At $x\approx{}0$, i.e.\ nearby the bottom stagnation-point of circular cylinder surface, Eqs. (14) and (15) diminish to the subsequent ordinary DEs:

$$ \left[\frac{1}{{\left(1-\varphi{}\right)}^{2.5}}\!-\!A_2M\right]f'\!=\!\frac{1}{{\left(1-\varphi{}\right)}^{2.5}}\!-\!A_2M{\!+\!A}_1\lambda{}\theta{}, $$
$ \frac{k_{nf}/k_b}{A_3}\left(1+\frac{4}{3}R_d\right){\theta{}}^{'{}'}+f{\theta{}}^{'}+s\theta{}=0,$ (18)}

with the boundary conditions:

$$ f\left(0\right)=0,\quad\ \theta{}\left(0\right)=1,\ \ \ \ \\ f^{'}\left(\infty{}\right)\rightarrow{}1,\quad\\theta{}(\infty{})\rightarrow{}0, $$
where primes indicate to $\eta{}$-differentiation. The practical interest quantities, in the current investigation are the coefficient of skin friction $C_{fx}$ and Nusselt number $Nu_x$ which are defined as:

$$ C_{fx}=\frac{{\tau{}}_w}{{\rho{}}_bU_{\infty{}}^2}=\frac{{\mu{}}_{nf}{\left(\frac{\partial{}\bar{u}}{\partial{}\bar{y}}\right)}_{\bar{y}=0}}{{\rho{}}_bU_{\infty{}}^2}, $$
$$ Nu_x=\frac{q_w\bar{x}}{k_b\left(T_w-T_{\infty{}}\right)}=\frac{{\left.-k_{nf}\frac{\partial{}T}{\partial{}\bar{y}}\right\vert{}}_{y=0}\bar{x}}{k_b\left(T_w-T_{\infty{}}\right)}, $$
where $\tau_w$ and $q_w$ denote the surface shear stress and the local heat transfer coefficient. Substituting from Eqs. (7) and (13) into Eqs. (20) and (21), to gain the skin friction and reduced Nusselt number as:

$$ PrPe^{1/2}C_{fx}=\frac{x}{{\left(1-\varphi{}\right)}^{2.5}}f^{'{}'}(x,0), $$
$$ Pe^{1/2}Nu_x=-\frac{k_{nf}}{k_b}{\theta{}}^{'}(x,0). $$

3 Solution Methodology

3.1 Finite Element Technique

To resolve the problem formulated above by the DEs (17)--(18), with the boundary condition (19), we suppose that:

$ f^{\prime}=h.$

The system of Eqs. (17) and (18) reduce to

$$ \frac{k_{nf}/k_b}{A_3}\left(1+\frac{4}{3}R_d\right){\theta{}}^{'{}'}+f{\theta{}}^{'}+s\theta{}=0,$$ $$
$\left[\frac{1}{{\left(1-\varphi{}\right)}^{2.5}}-A_2M\right]h'=A_1\lambda{}\theta{}',$ $$
subject to

$$ f\left(0\right)=0,\quad\ \theta{}\left(0\right)=1,\ \ \ \ \\ h\left(\infty{}\right)\rightarrow{}1,\quad\\theta{}(\infty{})\rightarrow{}0. $$

3.2 Formulation of Variation

The variation form coupled with Eqs. (24)--(25) over a standard linear two-nodded $({\eta{}}_e,{\eta{}}_{e+1})$ is specified as:

$$ \int_{{\eta{}}_e}^{{\eta{}}_{e+1}}w_1\left(f^{'}-h\right){\rm d}\eta{}=0, $$
$\int_{{\eta{}}_e}^{{\eta{}}_{e+1}}w_2\left\{\left[\frac{1}{{\left(1-\varphi{}\right)}^{2.5}}\!-\!A_2M\right]h'\!-\!A_1\lambda{}\theta{}'\right\}{\rm d}\eta{}\!=\!0,$
$$ \int_{{\eta{}}_e}^{{\eta{}}_{e+1}}w_3\left\{\frac{k_{nf}/k_b}{A_3}\left(1\!+\!\frac{4}{3}R_d\right){\theta{}}^{'{}'}\!+\!f{\theta{}}^{'}\!+\!s\theta{}\right\}{\rm d}\eta{}\!=\!0, $$
where $w_1,\ w_2$, and $w_3$ are arbitrary functions in $f$, $h$, and $\theta$.

3.3 Formulas of Finite Element

The model of finite element may be gained from Eqs. (27)--(29) by replacement approximations of finite element in the form:

$ f=\sum_{j=1}^2f_j{\psi{}}_j,\quad\ h=\sum_{j=1}^2h_j{\psi{}}_j,\quad\ \theta{}=\sum_{j=1}^2{\theta{}}_j{\psi{}}_j,$

with $w_1=w_2=w_3={\psi{}}_i\ (i=1,2)$, where ${\psi{}}_i$'s are the functions of shape for a squaredelement $({\eta{}}_e,{\eta{}}_{e+1})$ and are defined as:

$ {\psi{}}_1^e=\frac{{\eta{}}_{e+1}-\eta{}}{{\eta{}}_{e+1}-{\eta{}}_e},\quad\ {\psi{}}_2^e=\frac{\eta{}-{\eta{}}_e}{{\eta{}}_{e+1}-{\eta{}}_e},\quad\ {\eta{}}_e\leq{}\eta{}\leq{}{\eta{}}_{e+1}.$

The finite element system of equations is:

$ \left[\begin{array}{ccc} \left[K^{11}\right] \ \ \left[K^{12}\right] \ \ \left[K^{13}\right] \\ \left[K^{21}\right] \ \ \left[K^{22}\right] & \left[K^{23}\right] \\ \left[K^{31}\right] \ \ \left[K^{32}\right] \ \ \left[K^{33}\right]\end{array}\right]\left[\begin{array}{ccc}\left\{f\right\} \\ \left\{h\right\} \\ \left\{\theta{}\right\}\end{array}\right]=\left[\begin{array}{ccc}\left\{b^1\right\} \\ \left\{b^2\right\} \\\left\{b^3\right\} \end{array}\right],$

where $\left[K^{mn}\right]$ and $\left[b^m\right] (m,n=1,2,3)$ are defined as:

$ K_{ij}^{11}=\int_{{\eta{}}_e}^{{\eta{}}_{e+1}}{\psi{}}_i\frac{{\rm d}{\psi{}}_j}{{\rm d}\eta{}}{\rm d}\eta{},\quad\ {\ K}_{ij}^{12}=-\int_{{\eta{}}_e}^{{\eta{}}_{e+1}}{\psi{}}_i{\psi{}}_j{\rm d}\eta{},\quad\{\ K}_{ij}^{13}=0,\quad\ {\ K}_{ij}^{21}=0,\\ K_{ij}^{22}=\int_{{\eta{}}_e}^{{\eta{}}_{e+1}}\left\{\left[\frac{1}{{\left(1-\varphi{}\right)}^{2.5}}-A_2M\right]{\psi{}}_i\frac{{\rm d}{\psi{}}_j}{{\rm d}\eta{}}-A_1\lambda{}{\psi{}}_i\frac{{\rm d}{\psi{}}_j}{{\rm d}\eta{}}\right\}{\rm d}\eta{},\quad\ {\ K}_{ij}^{23}={\ K}_{ij}^{31}={\ K}_{ij}^{32}=0,\\ K_{ij}^{33}=\int_{{\eta{}}_e}^{{\eta{}}_{e+1}}\left\{\frac{1}{A_3}\left(\frac{k_{nf}}{k_b}\right)\frac{{\rm d}{\psi{}}_i}{{\rm d}\eta{}}\frac{{\rm d}{\psi{}}_j}{{\rm d}\eta{}}\left(1+\frac{4}{3}R_d\right)+\bar{f}{\psi{}}_i\frac{{\rm d}{\psi{}}_j}{{\rm d}\eta{}}+s{\psi{}}_i{\psi{}}_j\right\}{\rm d}\eta{},\\ b_i^1=0,\quad\ b_i^2=-{\left({\psi{}}_i\frac{{\rm d}h}{{\rm d}\eta{}}\right)}_{{\eta{}}_e}^{{\eta{}}_{e+1}},\quad\b_i^3=-{\left({\psi{}}_i\frac{{\rm d}\theta{}}{{\rm d}\eta{}}\right)}_{{\eta{}}_e}^{{\eta{}}_{e+1}},$

where,

$ \bar{f}=\sum_{j=1}^2\bar{f_i}{\psi{}}_i,\quad\ \bar{h}=\sum_{j=1}^2\bar{h_i}{\psi{}}_i.$

The same system after grouping of all the elemental equations is non-linear therefore an iterative technique is used in the solution, which is resolved by utilizing the Gaussian elimination technique with preserving an accuracy of ${10}^{-5}$.

4 Discussion

CNTs are cylindrical projectiles manufactured, in notion, by rolling graphene plates into a smooth cylinder. CNTs found in two types SWNTs or MWNTs. The SWNT composed of one graphene plate, whichis a flat benzene molecules matrix,including only hexagonal loops with single and double C-C bonding. The axis rolling selection relative to the graphene plate hexagonal grid and the closing cylinder radius provides for diverse kinds of SWNTs, which vary from insulating to conducting. Figure 2 depicts three various types of SWNTs:[32]

New window|Download| PPT slide
Fig. 2(Color online) Schematic SWNTs structures.



The impacts of NPs shapes in the existence of heat source/sink and radiation on the problem of MHD blood based on SWCNTs flow through a circular cylinder in horizontal position are analyzed and examined. The SWNTs shapes in SWNTs-blood play an important influence on temperature profiles. These three types of NPs shapes are illustrated in Fig. 3. The numerical calculations of the thermo-physical characteristics of SWCNTs and blood are utilized from Table 1. In order to admit that our method is acceptable, it is agreed from Table 2 that the numerical result of the heat transfer rate $-{\theta{}}^{\prime}(0)$ for diverse parameter of the $\lambda{}$ are associated with the results have been obtained by of Nazar $et al.$[30-31] and Cheng.[33] An ideal accord can be remarked among them. Table 3 introduces the calculated values of heat transfer rates $-\theta^{\prime}(0)$ for diverse $\varphi$, $\lambda$, $s$, $M$ and $R_d$ values with various NPs shapes. It is illustrated that the spherical NPs have a considerable effect in heat transfer rate generally, whereas the effects of governing parameters increase or decrease according to each parameter on $-\theta^{\prime}(0)$. This will be discussed in the subsequent section.

New window|Download| PPT slide
Fig. 3(Color online) NPs diverse shapes.[34]




Table 1
Table 1The fluid and NPs thermo-physical properties.[35]

New window|CSV


Table 2
Table 2Comparison values of $-\theta{'}(0)$ when $\varphi{}=0 \ $(pure blood) and $M=R_d=s=0$ with $\lambda$.

New window|CSV

Figure 4 illustrated the NPs volume fraction effects of $\varphi$ ($0.05\leq{}\varphi{}\leq{}0.2$) on the velocity $f'(\eta{})$ and temperature $\theta{}(\eta{})$ distributions for nanofluid (blood based on SWNTs) when $M\!=\!0.1,\ \lambda{}\!=\!R_d=1.0$, and $\ s\!=\!0.1$. It is remarkedthat for both cylindrical and lamina NPs shapes, as the $\varphi{}$ raises, the flow velocity and the nanofluid temperature diminish. The opposite trend takes place in the temperature profile with $\varphi{}$ in the case of spherical NPs. It is also observed that variations in $\varphi{}$ points to adaptations in nanofluid flow velocity and temperature which plays an important role in the cooling and heating procedures like, therapy of hematology diseases related to high body temperature. The effect of magnetic field parameter $M$ on velocity and temperature distributions is illustrated in Fig. 5 when $\varphi{}\!=\!0.05,\ \lambda{}\!=\!R_d\!=\!1.0$, and $\ s\!=\!0.5$. It is noticed that for growing the $M$ value, the velocity exceeds, on contrary the temperature. The magnetic field $M$ has powerful impacts on the flow velocity $f'$ than on the temperature $\theta{}$. This result to a physical internal structure characterization of the nanofluid in motion with NPs shape, which is at the own time conformity with the Lorentz force. Whereas $M$ depends on porous parameter $K$, the porous medium resistance is scaled down which hikes the momentum enhancement of the nanofluid flow regime definitively develops the velocity field for both SWNTs shapes. Figure 6 detects the influence of the heat source parameter $s$ on flow velocity and profile of temperature when $\varphi{}=0.05,\ \lambda{}=R_d=1.0,\ M=0.3$. The nanofluid velocity and temperature profiles exceed (for all types of NPs shapes) with mounting values $s$ resulting to an increment in the boundary layer thickness of the velocity. Physically, the increasing temperature results to heat source, which is boosted NPs movement and therefore enhances the velocity of the nanofluid. In Fig. 7, the distributions of velocity and temperature for diverse amounts of thermal radiation parameter $R_d$ are drawn when $\varphi{}=0.1,\ \lambda{}=1.0,\ s=0.5$ and $M=0.3$. It is noted that for considerable values of $R_d$, flow velocity and temperature increase which interpret the physical status that as $R_d$ raises. The velocity and boundary-layer thickness of temperature tend to become more stable with increasing values of $R_d$.


Table 3
Table 3Values of $-\theta{}{'}(0)$ for diverse$\varphi{}$, $\lambda{}$, $s$, $M$ and $R_d$ with different nanoparticle shapes

New window|CSV

New window|Download| PPT slide
Fig. 4(Color online) $f{'}\left(\eta{}\right)$ and$\ \theta{}\left(\eta{}\right)\ $via$\ \varphi{}$\ for diverse NPs shape.



New window|Download| PPT slide
Fig. 5(Color online) $f{'}\left(\eta{}\right)$ and $\ \theta{}\left(\eta{}\right)$ via $\ M$\ for diverse NPs shape.



New window|Download| PPT slide
Fig. 6(Color online) $f{'}\left(\eta{}\right)$ \ and\ $\theta{}\left(\eta{}\right)$\ via\ $s$\ for diverse NPs shape.



New window|Download| PPT slide
Fig. 7(Color online) $f{'}\left(\eta{}\right)$ and\ $ \theta{}\left(\eta{}\right)$\ via\ $R_d$\ fordiverse NPs shape.



Figure 8 exposes the effects of mixed convection parameter $\lambda{}$ on the fields of velocity and temperature at $\varphi{}=0.1,\ R_d=3.0,\ s=1.0$, and $M=0.3$. The nanofluid velocity raises steadily with rising values of $\lambda{}$ for diverse SWNTs shapes. The opposite direction happens in temperature profile with increasing $\lambda{}$. Fundamentally, this is the result to the causing that a rise in $\lambda{}$ indicators to an increment in the porosity of the medium $K$ and other the fluid viscosity ${\upsilon{}}_b$ diminishes. In cooled cylinder ($\lambda{}<0$), SWNTs spherical shape causes a disturbance in the flow velocity and temperature. This means that the flow is in the opposite trend or there is no flux when $\lambda{}<0$ and SWNTs with spherical shape.

New window|Download| PPT slide
Fig. 8(Color online) $f{'}\left(\eta{}\right)$ and\ $\theta{}\left(\eta{}\right)$\ via\ $\lambda{}$\ for diverse NPs shape.



Figure 9 demonstrates the heat transfer rate difference for SWNTs shapes for miscellaneous NPs volume fraction $\varphi{}$ with $\lambda{}$ with $R_d=2.0,\ s=0.5$, and $M=0.1$. It is monitored that with the rising of NPs promote heat transfer rate, except in the case of lamina shape. Depending on $k_{nf}$, the blood has declined heat transfer rate as compared with the blood based SWNTs.Lamina shape has a great impact when applied case of negative flow (i.e. opposing flow $\lambda{}<0$). Other, the growth values of $\lambda{}$increase the rate of heat transfer. Likewise, Fig. 10 describes the heat transfer differ with magnetic field $M$ for diverse shapes of NPs at $\varphi{}=0.05,\ R_d=1.0$, and $s=0.5$. It is shown that heat transfer rate hikes result to enforcement of applied transverse magnetic field. This because magnetic field yields electric current in nanofluid which yields heat and hence the magnetic field enhances phenomena. Figure 11 the radiation effect $R_d$ on heat transfer rate $-\theta{}'(0)$ when $\varphi{}=0.1,\ M=0.3\ $ and $\ s=0.5$. It is found displays that the rising $R_d$ raises transfer rate in opposing flow ($\lambda{}<0$) and reduces transfer rate in adding flow ($\lambda{}>0$) case. The impact of heat source or sink $s$ on heat transfer rate is drawn in Fig. 12. $-{\theta{}}^{'}\left(0\right)$ diminishes with the rising values of $s\ $ at $\varphi{}=0.1,\ R_d=1.0$,\ and $\ M=0.5$. The lamina shape is the significant affect in the case of opposing flow, but in positive flow case ($\lambda{}>0$) the spherical shape is more efficient than others.

New window|Download| PPT slide
Fig. 9(Color online) $-\theta{}'\left(0\right)$ via\ $ \lambda{}$\ with $\varphi{}$\ for diverse NPs shape.



New window|Download| PPT slide
Fig. 10(Color online) $-\theta{}'\left(0\right)$ via$\ \lambda{} $\ with $M$\ for diverse NPs shape.



New window|Download| PPT slide
Fig. 11(Color online) $-\theta{}'\left(0\right)$ via$\ \lambda{}$\ with $R_d$ for diverse NPs shape.



New window|Download| PPT slide
Fig. 12(Color online) $-\theta{}'\left(0\right)$ via$\ \lambda{}$\ with $s$\ for diverse NPs shape.



5 Closing Remarks

The impacts of NPs shapes in the existence of radiation and heat source/sink on the problem of MHD. The SWNTs shapes in SWNTs -- blood play an important influence on temperature profiles. The problem of MHD mixed convective nanofluid in the case blood based on SWCNTs flow through a circular cylinder are investigated numerically. Spherical, cylindrical and laminar NPs shapes in the presence of radiation and heat source or sink influences are analyzed. The flow velocity, temperature and rate of heat transfer inside the nanofluid with diverse NPs shapes are pointed in figures and tables terms with various parameters impacts and the computational outcomes are incorporated as follows:

$\bullet\ \ $The lamina shape ($m=16.1576$) of NPs in SWNTs -- blood plays an overriding role on flow velocity and temperature profiles with diverse governing parameters in the case of adding flow, but the spherical shape is dominant in opposing flow.

$\bullet\ \ $The flow velocity increases with the increase of NPs volume fraction, magnetic field, heat source/sink and radiation.

$\bullet\ \ $The spherical shape of SWNTs gives rise to opposite effects on temperature profile with rise of NPs volume fraction, which boosts the temperature for all the other shapes.

$\bullet\ \ $The increasing of mixed parameter and magnetic field lead to an increasing in temperature, but the contrary happens with growing values of radiation and heat source.

$\bullet\ \ $The NPs spherical shape has dynamic accomplishments in adding flow (heated cylinder) enhancement over a heat transfer rate.

$\bullet\ \ $The existence of magnetic field and NPs volume fraction has more advisable enhancement on heat transfer as compared with the heat source and radiation for diverse NPs shapes.

$\bullet\ \ $In opposing flow (cooled cylinder), heat transfer declines significantly,especially with the spherical NPs shape and diminishes increasingly with heat source increment, which makes it an important application in the cooling processes.

Conflicts of interest

The authors announce that there is no conflict of interests about the publication of this article.

Reference By original order
By published year
By cited within times
By Impact factor

J. Buongiorno , J. Heat Transf .128(2006) 240.
DOI:10.1115/1.2150834URL [Cited within: 1]

D. Nield and A. Kuznetsov , Int. J. Heat Mass Transf. 52(2009) 5792.
DOI:10.1016/j.ijheatmasstransfer.2009.07.024URL [Cited within: 1]

D. Nieldand A. Kuznetsov , Int. J. Heat Mass Transf. 52(2009) 5796.
DOI:10.1016/j.ijheatmasstransfer.2009.07.023URL [Cited within: 1]

F. Hady, F. Ibrahim, S. Abdel-Gaied, M. Eid , Int. Commun. Heat Mass Transf. 38(2011) 1414.
DOI:10.1016/j.icheatmasstransfer.2011.07.008URL [Cited within: 1]

F. Hady, F. Ibrahim, S. Abdel-Gaied, M. R. Eid, J. Mech. Sci . Technol. 25(2011) 2043.


F. Hady, F. Ibrahim, S. Abdel-Gaied, M. R. Eid , Appl. Math. Mech-Engl. 32(2011) 1577.
DOI:10.1007/s10483-011-1524-7URL
The free convective heat transfer to the power-law non-Newtonian flow from a vertical plate in a porous medium saturated with nanofluid under laminar conditions is investigated. It is considered that the non-Newtonian nanofluid obeys the mathematical model of power-law. The model used for the nanofluid incorporates the effects of Brownian motion and thermophoresis. The partial differential system governing the problem is transformed into an ordinary system via a usual similarity transformation. The numerical solutions of the resulting ordinary system are obtained. These solutions depend on the power-law index n, Lewis number Le, buoyancy-ratio number Nr, Brownian motion number Nb, and thermophoresis number Nt. For various values of n and Le, the effects of the influence parameters on the fluid behavior as well as the reduced Nusselt number are presented and discussed.

F. M. Hady, F. S. Ibrahim, S.M Abdel-Gaied ., and M.R. Eid, Nanoscale Res. Lett. 7(2012) 229.
DOI:10.1186/1556-276X-7-229URLPMID:22520273 [Cited within: 1]
In this work, we study the flow and heat transfer characteristics of a viscous nanofluid over a nonlinearly stretching sheet in the presence of thermal radiation, included in the energy equation, and variable wall temperature. A similarity transformation was used to transform the governing partial differential equations to a system of nonlinear ordinary differential equations. An efficient numerical shooting technique with a fourth-order Runge-Kutta scheme was used to obtain the solution of the boundary value problem. The variations of dimensionless surface temperature, as well as flow and heat-transfer characteristics with the governing dimensionless parameters of the problem, which include the nanoparticle volume fraction ?, the nonlinearly stretching sheet parameter n, the thermal radiation parameter NR, and the viscous dissipation parameter Ec, were graphed and tabulated. Excellent validation of the present numerical results has been achieved with the earlier nonlinearly stretching sheet problem of Cortell for local Nusselt number without taking the effect of nanoparticles.

M. R. Eid, J. Mol. Liq . 220(2016) 718.


M. R. Eid , J. Nanofluids 6 ( 2017) 550.


M.R Eid .and K. L. Mahny, Adv. Powder Technol. 28(2017) 3063.
DOI:10.1016/j.apt.2017.09.021URL

M. R. Eid, A. Alsaedi, T. Muhammad, T. Hayat , Results Phys. 7(2017) 4388.
DOI:10.1166/jnn.2009.m65URLPMID:19916462 [Cited within: 1]
Recent experiments [F. E. Pinkerton, M. S. Meyer, G. P. Meisner, M. P. Balogh, and J. J. Vajo, J. Phys. Chem. C 111, 12881 (2007) and J. J. Vajo and G. L. Olson, Scripta Mater. 56, 829 (2007)] demonstrated that the recycling of hydrogen in the coupled LiBH4/MgH2 system is fully reversible. The rehydrogenation of MgB2 is an important step toward the reversibility. By using ab initio density functional theory calculations, we found that the activation barrier for the dissociation of H2 are 0.49 and 0.58 eV for the B and Mg-terminated MgB2(0001) surface, respectively. This implies that the dissociation kinetics of H2 on a MgB2(0001) surface should be greatly improved compared to that in pure Mg materials. Additionally, the diffusion of dissociated H atom on the Mg-terminated MgB2(0001) surface is almost barrier-less. Our results shed light on the experimentally-observed reversibility and improved kinetics for the coupled LiBH4/MgH2 system.

M.R Eid K. . L.Mahny, Heat Transf. — Asian Res. 47(2018) 54.


M. R. Eid, K. L. Mahny, T. Muhammad, M. Sheikholeslami , Results Phys. 8(2018) 1185.
DOI:10.1016/j.rinp.2018.01.070URL

A.Wakif, Z. Boulahia, F. Ali , et al., Appl. Comput.Math. 4(2018) 81.


M. R . Eid and O. Makinde, Int.
[J]. Chem. React. Eng. 16( 2018), DOI: https//doi.org/10.1515/ijcre-2017-0212.



T. Muhammad, D. C. Lu, B. Mahanthesh , et al., Commun. Theor.Phys. 70(2018) 361.
[Cited within: 1]

Y. Ding, H. Alias, D. Wen, R. A. Williams , Int. J. Heat Mass Transf. 49(2006) 240.
DOI:10.1016/j.ijheatmasstransfer.2005.07.009URL [Cited within: 1]

J. P. Meyer, T. McKrell, K. Grote , Int. J. Heat Mass Transf. 58(2013) 597.
DOI:10.1016/j.ijheatmasstransfer.2012.11.074URL [Cited within: 1]

M. S. Liu, M. C. C. Lin, H. I. Te, C. C. Wang , Int. Commun. Heat Mass Transf. 32( 9) ( 2005) 1202-1210.
[Cited within: 1]

J. Wang, J. Zhu, X. Zhang, Y. Chen , Exp. Therm. Fluid Sci. 44(2013) 716.
DOI:10.1016/j.expthermflusci.2012.09.013URL [Cited within: 1]
The heat transfer and pressure drop of nanofluids containing carbon nanotubes in a horizontal circular tube are experimentally investigated. The friction factor of the dilute nanofluids shows a good agreement with the prediction from the Hagen-Poiseuille flow theory. A considerable enhancement in the average convective heat transfer is also observed compared with the distilled water. For the nanofluids with volumetric concentration of 0.05% and 0.24%, the heat transfer enhancement are 70% and 190% at Reynolds number of about 120 respectively, while the enhancement of thermal conductivity is less than 10%, therefore, the large heat transfer increase cannot be solely attributed to the enhanced thermal conductivity. By measuring the pump power supply and the thermal conductance of the test tube, our results suggest that the nanofluids at low concentration enhance the heat transfer with little extra penalty in pump power, thus have great potential for applications in the heat transfer systems. (C) 2012 Elsevier Inc.

V. Kumaresan, R. Velraj, S. K. Das , Int. J. Refrig. 35(2012) 2287.
DOI:10.1016/j.ijrefrig.2012.08.009URL [Cited within: 1]
This paper reports an experimental study on the convective heat transfer characteristics of a secondary refrigerant based CNT nanofluids in a tubular heat exchanger. The nanofluid is prepared by dispersing the multiwalled carbon nanotubes (MWCNT) in the water ethylene glycol mixture (70:30 by volume). The results show that the convective heat transfer coefficient is enhanced to a maximum of similar to 160% for the nanofluid containing 0.45 vol. % MWCNT, which could not be attributed solely by enhanced thermal conductivity of the nanofluids. Particle rearrangement, the very high aspect ratio and delaying the boundary layer development due the random movement of the carbon nanotubes are the possible mechanisms. Further, there is an appreciable decrease in Reynolds number for a given velocity for all the nanofluids. The enhancement in the friction factor is negligible at higher velocity and higher temperature for the nanofluids with 0.15 vol. % MWCNT. (C) 2012 Elsevier Ltd and IIR.

W. A. Khan, Z. H. Khan, M. Rahi , App. Nanosci. 4(2014) 633.
[Cited within: 1]

R. U. Haq, S. Nadeem, Z. H. Khan, N. F. M. Noor , Physica B 457 ( 2015) 40.
[Cited within: 1]

R. Kandasamy, R. Mohamad, M. Ismoen , Eng. Sci. Technol. 19(2016) 700.
[Cited within: 1]

R. Kandasamy, I. Muhaimin, R. Mohammad , Alex. Eng. 55(2016) 275.
[Cited within: 1]

S. Aman, I. Khan, Z. Ismail , et al., Sci. Rep.-UK 7 ( 2017) 2445.
[Cited within: 1]

R.C Chaudhary .and A. Jain, Rom.J. Phys. 52(2007) 505.
[Cited within: 1]

S. Mukhopadhyay, K. Bhattacharyya, T. Hayat , Chin. Phys. B 22 ( 2013) 114701.


S. Mukhopadhyay , Chin. Phys. B 22 ( 2013) 074701.


R. Nazar, N. Amin, I. Pop , Int. J. Heat Mass Transf. 46(2003) 3167.
DOI:10.1016/S0017-9310(03)00122-4URL [Cited within: 3]

R. Nazar, L. Tham, I. Pop, D. B. Ingham , Transport Porous Med. 86(2011) 517.
DOI:10.1007/s11242-010-9637-1URL [Cited within: 3]

M. Zhang and J. Li , Mater. Today 12 ( 2009) 12.
[Cited within: 1]

P. Cheng , Int. J. Heat Mass Transf. 20(1977) 807.
DOI:10.1016/0017-9310(77)90110-7URL [Cited within: 1]

R. Kandasamy, R. Mohammad, N. A. B. M. Zailani, N. F. B . Jaafar, J. Mol. Liq. 233(2017) 156.
[Cited within: 1]

A. Khalid, I. Khan, A. Khan , et al., Case Studies Therm. Eng. 12(2018) 374.
DOI:10.1016/j.csite.2018.04.004URL [Cited within: 1]

相关话题/BloodBased SWCNTs Through