CC BY
83
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ФИЗИЧЕСКИХ ПРОЦЕССОВ
DOI: 10.5862/JPM.253.3 UDC 536.24
N.V. Ganesh1, B. Ganga2, A.K.A. Hakeem3, S. Saranya3, R. Kalaivanan3
1 Karpagam Academy of Higher Education, Coimbatore, Tamil Nadu, India
2 Providence College for Women, Coimbatore, Tamil Nadu, India
3 Sri Ramakrishna Mission Vidyalaya College of Arts and Science,
Coimbatore, Tamil Nadu, India
HYDROMAGNETIC AXISYMMETRIC SLIP FLOW ALONG A VERTICAL STRETCHING CYLINDER wITH A CONVECTIVE BOuNDARY CONDITION
Steady axisymmetric laminar boundary-layer slip flow of a viscous incompressible fluid and heat transfer towards a vertical stretching cylinder in the presence of a uniform magnetic field is investigated. It is assumed that the left surface of the cylinder is heated by a hot convective flow. Using a similarity transformation, the governing system of partial differential equations is first transformed into a system of coupled nonlinear ordinary differential equations. The resulting intricate nonlinear boundary value problem is solved numerically by the fourth-order Runge — Kutta method with the shooting iteration technique. The analytical solutions are presented for a special case. The effects of various physical parameters on the velocity and temperature profiles are discussed through graphs. The values of the skin friction coefficient and the Nusselt number are tabulated and examined. It is found that the thermal boundary layer thickness increases with an increase in the velocity slip, the magnetic field, the surface convection parameter and the curvature parameter and decreases with the Prandtl number.
AXISYMMETRIC FLOW, MHD, SLIP FLOW, STRETCHING CYLINDER, CONVECTIVE BOUNDARY.
1. Introduction
Flow and heat transfer due to a stretching cylinder or a flat plate in a quiescent or moving fluid is important in a number of industrial manufacturing processes that include manufacturing both metal and polymer sheets. Flow over a cylinder is considered to be two-dimensional if the body's radius is large and comparable to the boundary-layer thickness. On the other hand, for a thin or a slender cylinder, the radius of the cylinder can be of the same order as the boundary-layer thickness. Therefore, the flow may be considered axisymmetric instead of two-dimensional along a vertical or a horizontal cylinder [1 — 6].
Magnetohydrodynamic (MHD) flow and heat transfer for fluid have broad applications in various engineering problems such as MHD power generators, petroleum industries, plasma studies, geothermal energy extractions, boundary layer control in the field of aerodynamics, and many others [7—27]. MHD flow along a stretching cylinder has been considered by the authors of Refs. [28-31].
The nonadherence of the fluid to a solid boundary, known as velocity slip, is a phenomenon that has been observed under certain circumstances. When fluid is encountered in microelectromechanical systems, the no-slip condition at the solid-fluid
interface is abandoned in favor of a slip flow model which represents the nonequilibrium region near the interface more accurately. In all of the above-mentioned investigations, the no-slip condition at the boundary was assumed. Even in literature, the slip flow over a flat plate has not been studied sufficiently. Zheng, et al. [32] investigated the flow and the radiation heat transfer of nanofluid over a stretching sheet with the velocity slip and the temperature jump in a porous medium. Mukhopadhyay [33] analyzed the slip effects on MHD flow along a horizontal stretching cylinder. Very recently, Abbas, et al. [34] studied the slip effects on the flow over an unsteady stretching/shrinking cylinder in the presence of suction.
Heat transfer characteristics of the stretching sheet problem have been restricted to two boundary conditions of either prescribed temperatures or heat flux on the wall in the published papers. Most recently, the heat transfer problems for a boundary layer flow concerning a convective boundary condition was investigated by Aziz for Blasius flow [35]. Following Aziz, many researchers investigated boundary layer flows with the convective boundary condition [36—40].
All the above-mentioned studies on the flow along a stretching cylinder have been carried out with a prescribed surface temperature (PST) and prescribed surface heat flux (PHF) boundary conditions only. This motivates us to do the present work. We investigated the slip effect on MHD flow along a vertical stretching cylinder with a convective boundary condition. Using a similarity transformation, the governing system of partial differential equations was first transformed into coupled nonlinear ordinary differential equations and then solved through the numerical shooting method.
2. Equations of motion
Let us consider the steady axisymmetric slip flow of incompressible fluid along a vertical stretching cylinder in the presence of a uniform magnetic field. The x-axis is measured along the tube axis and the r-axis is measured in the radial direction. It is assumed that a uniform magnetic field of intensity B0 acts in the radial direction. The magnetic Reynolds number
Fig. 1. The coordinate system and the physical model
is assumed to be small so that the induced magnetic field is negligible in comparison with the applied magnetic field. It is also assumed that the left side of the cylinder is heated by convection from the hot fluid at temperature Tf which provides a heat transfer coefficient h, and T is the ambient fluid temperature (see
Fig. 1).
The continuity, the momentum and the energy equations governing such type of flow [331 are as follows:
dx
du du u — + v — dx dr
d(ru) d(rv) _ n ■ + u,
dr
vd_ r dr
du dr
J J
+
+ gB(T - r )
dT dT
u — + v-
dx dr
a d ( ( dT^
r dr
r
\ \
dr
(1)
(2)
(3)
where u and v are the velocity components, respectively in x and r directions; v = ^ / p is the kinematic viscosity (^ is the dynamic viscosity of the fluid); p is the fluid density; g is the acceleration due to gravity; B is the
volumetric coefficient of thermal expansion; o is the electric conductivity of the medium, B0 is the uniform magnetic field, a is the thermal diffusivity of the fluid; T is the fluid temperature.
2.1. Boundary conditions
The appropriate boundary conditions for the problem are given by the following formulae:
at r = R;
u = U(x) + A1v —, v = 0, dr
- k dT = hf ( Tf - T ) dr
u ^ 0, T ^ T» as r ^ ».
(4)
(5)
Here U(x) = U0x/L is the stretching velocity (U0 is the reference velocity, L is the characteristic length), A1 is the velocity slip, k is the thermal conductivity, Tf is the hot fluid temperature, hf is the convective heat transfer coefficient, R is the cylinder radius.
2.2. Method of solution
The continuity equation is automatically satisfied by the introduction of a stream function y as
1 -1 u =--—, v =--—.
r dr r dx
Introducing the similarity variable as
r - R ( U
0.5
П = ' ,
2R V vx
y = (U vx)05 RF (n),
0(n) = , AT = Tf - T»
AT
(6)
and upon substitution of (6) in Eqs. (2) — (5) the governing equations and the boundary conditions are reduced to
(1 + 2Cn)F'" + 2CF" + FF" -- F'2 - M2F' +10 = 0,
(7)
(1 + 2Cn)0" + 2C0 + Pr(F0' - F'0) = 0, (8)
F = 0, F' = 1 + AF", 0' = -Y(l - 0) at n = o,
F' ^ 0,
0 as n ^ »,
(9) (10)
where the prime denotes differentiation with respect to n; A is the slip parameter, M2 is the magnetic parameter,
A = A
U0 v
M 2 = gB2 L
L ) ' pUo '
C is the curvature parameter, 1 is the mixed convection parameter,
g3A70L2 .
C =
/ \0.5
^ vL Л
V U0R у
; i =
u2
Pr is the Prandtl number, y is the surface convection parameter,
Pr =
mC
k
Y =
hf ( Lv^
k
V U0 у
We assume that AT = AT0 x. Other important characteristics of the present investigation are the skin friction coefficient C and the Nusselt number N which
f u
are defined as follows:
'duN
2m
Cf =
dr
r =R
x
N.. = —
PU2 dT dr
= 2F"(0)(Rex)
0.5.
_ - = -(Rex)05 0'(O).
w «
Here Rex is the local Reynolds number, tj Ux
Re x =— .
V
2.4. Analytical solutions for a special case
It is significant th=t tf C =0 (in. R^ <»), the problem under consideration (with t = M = A = 0 and y ^ <» ) is reduced to the boundary layer flow along a stretching at the plate. The closed form of the solution of Eq. (7), with the boundary conditions given in Eqs. (9) and (10), in the absence of 1 and C is obtained as
f ( n ) =
en Л
where
ß = -(3 A)-1 -
1 - e
Aß2 +ß
21/3ß1
(11)
3 A(ß2 + V4ßT+ßf )1/3
(р2 W4 в3 + Р2 )1/3
3 А21/3
with
p1 = -(1 + 3 M2 A2), p2 = A2(27 + 18M2) - 2.
The analytical solution of Eq. (8) with the boundary conditions in Eqs. (9) and (10) is obtained by the power series method in terms of Kummer's hypergeometric function [41] as
0(n) = c1e-a°enM\ a0 -1, a0 +1,
- Pr X
-Pn
P2
•(12)
The heat transfer rate at the surface is given
by
0'(O) = c1
-a0 pM a0 - 1, a0 + 1,
- Pr X
P2
a0 - 1V Pr X
V a0 + 1
a
M
2 + ~ Pr X
a0 , 2 + a0 , q2
(13)
where
c1 = -Y
Pr X
Pr X
x M \ a0 - 1, a0 + 1,
- Pr X
a0 - 1 a0 + 1
X M \ a0,2 + a0,-
Pr X
YM \ a0 - 1, a0 + 1,
- Pr X
a0 = Pr2X, X = (ЛР + 1)-1.
2.5. Numerical methodforsolution
The nonlinear coupled differential Eqs. (7) and (8) along with the boundary conditions (9) and (10) form a two-point boundary value problem. This problem is solved using the shooting technique together with the fourth-order Runge-Kutta integration scheme by converting it into an initial value problem. In terms of this method we have to choose a suitable finite value of for instance n„.
We set the following first-order system:
У1 = У y' = У
Уз' = (1 + Сn)-1(+ M2У2 -
- 2Суз - y Уз - Ay,),
У 4 = У5>
Уз' = -(1 + Сn)-1[2Су5 + Pr(У1 y5 - y2У4)] with the following boundary conditions: y^O) = 0, У2(0) = 1 + АузЗ),
У4(0) = 1 + ' У5(0)
(14)
(15)
To solve Eq. (14) with conditions (15) as an initial value problem, the values fory3(0), i.e., F"(0), and y5(0), i.e. 0'(0), are required, but such values are not given. In order to obtain the solution, we chose the initial guess values for F"(0) and 0'(0) and applied the fourth-order Runge-Kutta integration scheme. Then we compared the calculated values of F'(n) and 0(n) at n„ with the given boundary conditions and adjusted the values F"(0) and 0'(0) using the shooting technique to achieve the better approximation for the solution. The process was repeated until we obtained the correct results up to the desired accuracy of 10-9 level, which fulfilled the convergence criterion.
3. Results and discussion
For the verification of the present numerical results, we have compared F "(0) values with those obtained by Mukhopadhyay [11] in the absence of a slip and mixed convection parameters which are presented in Table 1. The comparison results are found to be good.
In order to analyze the obtained results, the numerical computations were carried out for various values of the curvature parameter С, the mixed convection parameter 1, the magnetic parameter M, the Prandtl number Pr and the surface convection parameter y. The numerical values are plotted in Figs. 2—6 for illustration of the results. The values of the skin friction coefficient and the Nusselt number have been calculated and are presented in Table 2.
The effects of С and А on the dimensionless horizontal velocity profile are presented
+
Tab le 1
Comparison between the results obtained for the local skin friction coefficient F ''(0)
Magnetic parameter Mukhopadhyay [11] (without slip parameters) Present results (without the slip and 1)
Without curvature
0.5 -1.11803400 -1.11803399
1.0 -1.41421350 -1.41421356
1.5 -1.80277564 -1.80277564
Fig. 2. Effects of slip and curvature parameters on the velocity (a) and the temperature (b) profiles
with M = 1, l = 1, Pr = 8, y = 0.6;
a — solid lines correspond to C variation with A = 0.2 (variation I), dashed ones to A variation with C = 0.1 (variation II); b — solid lines correspond to variation II, dashed ones to variation I
а)
Пч)
b)
0(л)
0.5 1.0 1.5 2.0 2.5 3.0 3.5
Fig. 3. Effects of magnetic and mixed convection parameters on the velocity (a) and the temperature (b)
profiles with C = 0.1, A = 0.2, Pr = 8, y = 0.6. Solid lines correspond to l variation with M = 1, dashed ones do to M variation with X = 1
in Fig. 2, a. It is clear that the curvature parameter decreases the velocity profile near the wall and increases it far away from the wall. When a slip occurs, the flow velocity near the stretching wall is no longer equal to the stretching velocity of the wall. With an increase in A such
a slip velocity increases and consequently the fluid velocity decreases because under the slip condition, the pulling of the stretching wall can only be partly transmitted to the fluid. Thus, the increasing values of the slip parameter decelerate the dimensionless velocity profile.
е(л>
0.5 1.0 1.5 2.0 Л
Fig. 4. Effects of the Prandtl number and the surface convection parameter on the temperature profile with C = 0.1, A = 0.2, M = l = 1;
Y variation with Pr = 8 (solid lines) and Pr variation with y = 0.6 (dashed ones) are presented
F"(0) Or
Fig. 5. Effects of curvature, slip, mixed convection and magnetic parameters on F "(0). C variation with M = 1, l = 1, Pr = 8, A = 0.2 and y = 0.6 (solid lines); X variation with M = 1, C = 0.1, Pr = 8, A = 0.2 and y = 0.6 (dashed ones); A variation with M = 1, C = 0.1, Pr = 8, l = 1, and y = 0.6
(dotted lines with stars) are presented
Fig. 3, a depicts the effects of M and X on the horizontal velocity profile. It can be seen that increasing values of the magnetic parameter reduce the horizontal velocity profile. The figure clearly shows that the transverse magnetic field
opposes the transport phenomena. This is due to the fact that M variation leads to the variation of the Lorentz force induced by the magnetic field, and this force produces more resistance to transport phenomena. With an increase in 1,
the horizontal velocity is found to increase the buoyancy-induced flow (1 > 0) but decreases the buoyancy-opposed one (1 < 0). For 1 > 0, there is a favorable pressure gradient caused by the buoyancy forces, which results in the flow being accelerated.
The effects of C and A on the dimensionless temperature profile are shown in fig. 2, b. It can be seen that increasing the values of the
curvature parameter enhances the thermal boundary-layer thickness. The temperature profile increases with an increase in the slip parameter. This is because an increase in the slip parameter leads to the thickening of the thermal boundary layer.
fig. 3, b demonstrates the effects of 1 and M on the temperature profile. with increasing 1, the temperature is found to decrease in the
e'(o)
0.51 0.50 0.49 0.48 0.47 0.46 0.45
-e'(0)r
.....*..
...................
2.8 M
0.50
0.45
0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 M
fig. 6. Effects of curvature, slip, mixed convection, surface convection and magnetic parameters on 0'(0). Here we present: C variation with M = 1, l = 1, Pr = 8, A = 0.2 and y = 0.6 (solid lines); l variation with M = 1, C = 0.1, Pr = 8, A = 0.2 and y = 0.6 (dashed ones); A variation with M = 1, C = 0.1, Pr = 8, l = 1, and y = 0.6 (dotted lines with stars); y variation with M = 1, C = 0.1, Pr = 8, l = 1, A = 0.2 (dotted ones with circles)
buoyancy-induced flow and to increase with 1 in the buoyancy-opposed one. The thermal boundary-layer thickness decreases with an increase in the values of the mixed convection parameter. The temperature profile enhances with an increase in the magnetic parameter. Thus, the presence of the magnetic field leads to increasing the thermal boundary-layer thickness.
The effects of Pr and y on the dimensionless temperature profile are elucidated in Fig. 4. It is clear that the increasing values of Pr reduce the temperature profile. The fluids with lower Prandtl numbers have higher thermal conductivities (and thicker thermal boundary-layer structures) so that heat can diffuse from the sheet faster than for fluids with higher Pr numbers (that have thinner boundary-layers). Hence the Prandtl number can be used to increase the rate of cooling in conductive fluids. It is also seen that the temperature profile
enhances with the surface convection parameter. The y parameter is directly proportional to the heat transfer coefficient associated with the hot fluid h. The thermal resistance on the hot fluid side is inversely proportional to h. Thus, as y increases, the hot fluid side convection resistance decreases and consequently, the temperature profile increases.
Physically, the negative sign of F"(0) implies that the surface exerts a dragging force on the fluid and the positive sign implies the opposite effect. This is consistent with the present case, as a stretching cylinder which induces the flow is considered here. Fig. 5 shows the variation of the skin friction coefficient. It is clear that this coefficient F "(0) grows with increasing 1 and A, and decreases with increasing C and M. The variation of the Nusselt number is shown in Fig. 6. It is observed that the Nusselt number grows in value with increasing values of C, 1 and y, and decreases with M and A.
Table 2
Various values of governing parameters and /""(0), 8'(0) values for them
Governing parameter Value
Governing parameter —F" (0) -0'(0)
C 0.1 0.2 0.3 1.06754708 1.09310236 1.11775841 0.18479261 0.18493968 0.18508489
M 1 2 3 1.06754708 1.52007910 1.91790594 0.18479261 0.18220507 0.17847387
l 1 2 3 1.06754708 1.05348361 1.03963697 0.18479261 0.18484647 0.18489862
A 0.2 0.3 0.4 1.06754708 0.94852567 0.85471337 0.18479261 0.18407353 0.18341212
Pr 6 7 8 1.06754708 1.06950873 1.07099257 0.18479261 0.18595638 0.18689031
Y 0.1 0.2 0.3 1.07438124 1.06754708 1.06125922 0.09604074 0.18479261 0.26708712
Note: While studying the effect of individual parameters the following values were assumed: C = 0.1, M = 1, l = 1, A = 0.2, Pr = 6, y = 0.2.
4. Conclusion
The steady axisymmetric laminar boundary-layer slip flow of viscous incompressible fluid and heat transfer towards a vertical stretching cylinder in the presence of the uniform magnetic field has been investigated. Using a similarity transformation, the governing system of partial differential equations was first transformed into coupled nonlinear ordinary differential equations and then was solved with the help of the numerical shooting method. The parameters involved in this study significantly affect the flow and heat transfer. The following conclusions can be drawn as a result of the computations:
(i) The thickness of the momentum boundary layer decreases in the presence of the velocity slip and the magnetic field and increases with the mixed convection and curvature parameters;
(ii) The thermal boundary-layer thickness increases with an increase in the velocity slip,
the magnetic field, the surface convection parameter and the curvature parameter, and decreases with the Prandtl number;
(iii) The skin friction coefficient enhances with the increasing values of mixed convection slip and surface convection parameters and decreases with the curvature parameter, the magnetic parameter and the Prandtl number. The Nusselt number increases with the curvature parameter, the mixed convection parameter, the Prandtl number and the surface convection parameter, and decreases with magnetic and slip parameters.
we expect that the physics of flow over a stretching cylinder can be utilized as a basis for many engineering and scientific applications especially in petroleum engineering with the help of the present model. The results pertaining to the present study may be useful for investigating different models. The findings of the present problem are also of great interest in different areas of science and technology, where the surface layers are being stretched.
REFERENCES
[1] M. Kumari, G. Nath, Mixed convection boundary layer flow over a thin vertical cylinder with localized injection/suction and cooling/heating, Int. J. Heat Mass Transfer. 47(5) (2004) 969-976.
[2] E.M.E. Elbarbary, N.S. Elgazery, Flow and heat transfer of a micropolar fluid in an axi symmetric stagnation flow on a cylinder with variable properties and suction, Acta Mech. 176(3-4) (2005) 213-229.
[3] P. Datta, D.A. Kumar, S. Roy, N.C. Mahanti, Effect of non-uniform slot injection (suction) on a forced flow over a slender cylinder, Int. J. Heat Mass Transfer. 49 (13-14) (2006) 2366-2371.
[4] A. Ishak, R. Nazar, I. Pop, Uniform suction /blowing effect on flow and heat transfer due to a stretching cylinder, Appl. Math. Model. 32(10) (2008) 2059-2066.
[5] A. Ishak, R. Nazar, Laminar boundary layer flow along a stretching cylinder, Eur. J. Sci. Res. 36(1) (2009) 22-29.
[6] S. Mukhopadhyay, Mixed convection boundary layer flow along a stretching cylinder in porous medium, Journal of Petroleum Science and Engineering. 96-97 (2012) 73-78.
[7] M. Turkyilmazoglu, Analytic heat and mass transfer of the mixed hydrodynamic/thermal slip MHD viscous flow over a stretching sheet, Int. J. Mech. Sci. 53(10) (2011) 886-896.
[8] L. Zheng, J. Niu, X. Zhang, Y. Gao, MHD
flow and heat transfer over a porous shrinking surface with velocity slip and temperature jump, Math. Comput. Model. 56 (5- 6) (2012) 133-144.
[9] M. Turkyilmazoglu, MHD fluid flow and heat transfer due to a stretching rotating disk, Int. J. Therm. Sci. 51 (2012) 195-201.
[10] M. Turkyilmazoglu, Multiple analytic solutions of heat and mass transfer of magneto-hydrodynamic slip flow for two types of viscoelastic fluids over a stretching surface, J. Heat Transfer. 134(7) (2012) doi:10.1115/1.4006165.
[11] M.S. Kandelousi, Effect of spatially variable magnetic field on ferrofluid flow and heat transfer considering constant heat flux boundary condition, Eur. Phys. J. Plus. 129 (11) (2014) 248.
[12] M. Sheikholeslami, S. Abelman, D.D. Ganji, Numerical simulation of MHD nanofluid flow and heat transfer considering viscous dissipation, Int. J. Heat Mass Transfer. 79 (2014) 212-222.
[13] M. Turkyilmazoglu, An analytical treatment for the exact solutions of MHD flow and heat over two-three dimensional deforming bodies, Int. J. Heat Mass Transfer. 90 (2015) 781-789.
[14] M. Sheikholeslami, M.M. Rashidi, D.D. Ganji, Numerical investigation of magnetic nanofluid forced convective heat transfer in existence of variable magnetic field using two phase model , J. Mol. Liq. 212 (2015) 117-126.
[15] M. Sheikholeslami, R. Ellahi, Three-
dimensional mesoscopic simulation of magnetic field effect on natural convection of nanofluid, Int. J. Heat Mass Trans. 89 (2015) 799-808.
[16] M. Sheikholeslami, M.M. Rashidi, Effect of space dependent magnetic field on free convection of Fe304-water nanofluid, J. Taiwan Inst. Chem. Eng. 56 (2015) 6-15.
[17] M. Sheikholeslami, S. Abelman, Two-phase simulation of nanofluid low and heat transfer in an annulus in the presence of an axial magnetic field, IEEE Trans. Nanotechnol. 14 (3) (2015) 561-569.
[18] M. Sheikholeslami, M.M. Rashidi, D.D. Ganji, Effect of non-uniform magnetic field on forced convection heat transfer of Fe304-water nanofluid, Comput. Methods Appl. Mech. Eng. 294 (2015) 299-312.
[19] M. Sheikholeslami, D.D. Ganji, M.Y. Javed, R. Ellahi, Effect of thermal radiation on magnetohydrodynamics nanofluid flow and heat transfer by means of two phase model, J. Magn. Magn. Mater. 374 (2015) 36-43.
[20] C. Zhang, L. Zheng, X. Zhang, G. Chen, MHD flow and radiation heat transfer of nanofluids in porous media with variable surface heat flux and chemical reaction, Appl. Math. Model. 39 (1) (2015) 165-181.
[21] Y. Lin, L. Zheng, X. Zhang, L. Ma, G. Chen, MHD pseudo-plastic nanofluid unsteady flow and heat transfer in a finite thin film over stretching surface with internal heat generation, Int. J. Heat Mass Transfer. 84 (2015) 903-911.
[22] M. Sheikholeslami, K. Vajravelu, M.M. Rashidi, Forced convection heat transfer in a semi annulus under the influence of a variable magnetic field, Int. J. Heat Mass Transfer. 92 (2016) 339-348.
[23] M. Sheikholeslami, M.M. Rashidi, T. Hayat, D.D. Ganji, Free convection of magnetic nanofluid considering MFD viscosity effect, J. Mol. Liq. 218 (2016) 393-399.
[24] M. Sheikholeslami, T. Hayat, A. Alsaedi, MHD free convection of Al203-water nanofluid considering thermal radiation: A numerical study, Int. J. Heat Mass Transfer. 96 (2016) 513-524.
[25] M. Sheikholeslami, H.R. Ashorynejad, P. Rana, Lattice Boltzmann simulation of nanofluid heat transfer enhancement and entropy generation, J. Mol. Liq. 214 (2016) 86-95.
[26] M. Sheikholeslami, S. Soleimani, D.D. Ganji, Effect of electric field on hydrothermal behavior of nanofluid in a complex geometry, J. Mol. Liq. 213 (2016) 153-161.
[27] M. Sheikholeslami, M.M. Rashidi, Non-uniform magnetic field effect on nanofluid hydrothermal treatment considering Brownian motion and thermophoresis effects, J. Braz. Soc. Mech. Sci. Eng. 38 (4) (2016) 1171-1184.
[28] I. Anuar, N. Roslinda, I. Pop, Magneto-hydrodynamic (MHD) flow and heat transfer due to a stretching cylinder, Energy Conversion and Management. 49 (11) (2008) 3265-3269.
[29] K. Vajravelu, K.V. Prasad, S.R. Santhi, Axisymmetric magneto-hydrodynamic (MHD) flow and heat transfer at a non-isothermal stretching cylinder, Applied Mathematics and Computation. 219(8) (2012) 3993-4005.
[30] A.S. Butt, A. Ali, Entropy analysis of magnetohydrodynamic flow and heat transfer due to a stretching cylinder, Journal of the Taiwan Institute of Chemical Engineers. 45 (3) (2014) 780-786.
[31] T. Hayat, A. Shafiq, A. Alsaedi, MHD axisymmetric flow of third grade fluid by a stretching cylinder, Alexandria Engineering Journal. 54(2) (2015) 205-212.
[32] L. Zheng, C. Zhang, X. Zhang, J. Zhang, Flow and radiation heat transfer of a Nanofluid over a stretching sheet with velocity slip and temperature jump in porous medium, J. Franklin Inst. 350 (5) (2013) 990-1007.
[33] S. Mukhopadhyay, MHD boundary layer slip flow along a stretching cylinder, Ain Shams Eng. J. 4 (2) (2013) 317-324.
[34] Z. Abbas, S. Rasool, M.M. Rashidi, Heat transfer analysis due to an unsteady stretching/ shrinking cylinder with partial slip condition and suction, Ain Shams Eng. J. 6(3) (2015) 939-945.
[35] A. Aziz, A similarity solution for laminar thermal boundary layer over a flat plate with a convective surface boundary condition, Commun. Nonlinear Sci Numer. Simul. 14(4) (2009) 1064-1068.
[36] O.D. Makinde, Similarity solution of hydromagnetic heat and mass transfer over a vertical plate with a convective surface boundary condition, Int. J. Phys. Sci. 5 (6) (2010) 700-710.
[37] A. Ishak, Nor A. Yacob, N. Bachok, Radiation effects on the thermal boundary layer flow over a moving plate with convective boundary condition, Meccanica. 46 (4) (2011) 795-801.
[38] T. Hayat, Z. Iqbal, M. Qasim, S. Obaidat, Steady flow of an Eyring Powell fluid over a moving surface with convective boundary conditions, Int. J. Heat Mass Transfer. 55 (7-8) (2012) 1817-1822.
[39] N.S. Akbar, S. Nadeem, R. Ul Haq, Z.H. Khan, Radiation effects on MHD stagnation point flow of nano fluid towards a stretching surface with convective boundary condition, Chinese Journal of Aeronautics. 26(6) (2013) 1389-1397.
[40] S. Nadeem, R. Mehmood, N.S. Akbar, Optimized analytical solution for oblique flow of a Casson nano-fluid with convective boundary conditions, Int. J. Thermal Sciences. 78 (2014) 90-100.
[41] S.P. Anjali Devi, B. Ganga, Effects of surface embedded in a porous medium, Nonlinear viscous and Joules dissipation on MHD flow, Analysis: Modelling and Control. 14 (3) (2009) heat and mass transfer past a stretching porous 14303-14314.
THE AuTHORS
GANESH N. Vishnu
Karpagam Academy of Higher Education
Pollachi Main Road, Eachanari, Coimbatore, Tamil Nadu 641021, India GANGA B.
Providence College for Women
Spring Field, Coonoor, Bandishola, Tamil Nadu 643104, India gangabhose@gmail.com
HAKEEM Abdul A.K.
Sri Ramakrishna Mission Vidyalaya College of Arts and Science
Periyanaickenpalayam, Sri Ramakrishna Vidyalaya Post, Coimbatore, Tamil Nadu 641020, India abdulhakeem6@gmail.com
SARANYA Sarani
Sri Ramakrishna Mission Vidyalaya College of Arts and Science
Periyanaickenpalayam, Sri Ramakrishna Vidyalaya Post, Coimbatore, Tamil Nadu 641020, India
KALAIVANAN Raja
Sri Ramakrishna Mission Vidyalaya College of Arts and Science
Periyanaickenpalayam, Sri Ramakrishna Vidyalaya Post, Coimbatore, Tamil Nadu 641020, India
Ганеш Н.В., Ганга Б., Хаким А.К.А., Саранья С., Калайванан Р. ГИДРОМАГНИТНЫЙ ОСЕСИММЕТРИЧНЫЙ ОБТЕКАЮЩИЙ ПОТОК ВДОЛЬ ВЕРТИКАЛЬНОГО РАСТЯНУТОГО ЦИЛИНДРА С КОНВЕКТИВНЫМ ГРАНИЧНЫМ УСЛОВИЕМ.
Исследован стационарный осесимметричный ламинарный обтекающий поток вязкой несжимаемой жидкости в пограничном слое, а также его теплопередача по направлению к вертикальному растянутому цилиндру в однородном магнитном поле. Сделано предположение, что левая поверхность цилиндра нагревается горячим конвективным потоком. В качестве первого шага мы преобразовали определяющую систему дифференциальных уравнений с частными производными в систему парных нелинейных обыкновенных дифференциальных уравнений, используя преобразование подобия. Полученная сложная нелинейная краевая задача была решена численно методом Рунге-Кутты 4-го порядка с помощью итерационной процедуры пристрелки. Для специального случая нами представлены аналитические решения. Влияние различных физических параметров на профили скорости и температуры обсуждается с привлечением графического материала. Значения коэффициента поверхностного трения и числа Нуссельта сведены в таблицы и проанализированы. Установлено, что толщина термического граничного слоя увеличивается с ростом скоростного скольжения, величины магнитного поля, параметров поверхностной конвекции и кривизны и уменьшается с ростом числа Прандтля.
ОСЕСИММЕТРИЧНЫй ПОТОК, МАГНИТОГИДРОДИНАМИКА, ОБТЕКАЮЩИЙ ПОТОК, РАСТЯНУТЫЙ ЦИЛИНДР, КОНВЕКТИВНАЯ ГРАНИЦА.
СПИСОК ЛИТЕРАТУРЫ
[1] Kumari M., Nath G. Mixed convection boundary layer flow over a thin vertical cylinder with localized injection/suction and cooling/heating // Int. J. Heat Mass Transfer. 2004. Vol. 47. No. 5.
Pp. 969-976.
[2] Elbarbary E.M.E., Elgazery N.S. Flow and heat transfer of a micropolar fluid in an axi symmetric stagnation flow on a cylinder with variable properties
and suction // Acta Mech. 2005. Vol. 176. No. 3-4. Pp. 213-229.
[3] Datta P., Kumar D.A., Roy S., Mahanti N.C. Effect of non-uniform slot injection (suction) on a forced flow over a slender cylinder // Int. J. Heat Mass Transfer. 2006. Vol. 49. No. 13-14. Pp. 2366-2371.
[4] Ishak A., Nazar R., Pop I. Uniform suction /blowing effect on flow and heat transfer due to a stretching cylinder // Appl. Math. Model. 2008. Vol. 32. No. 10. Pp. 2059-2066.
[5] Ishak A., Nazar R. Laminar boundary layer flow along a stretching cylinder // Eur. J. Sci. Res. 2009. Vol. 36. No. 1. Pp. 22-29.
[6] Mukhopadhyay S. Mixed convection boundary layer flow along a stretching cylinder in porous medium // Journal of Petroleum Science and Engineering. 2012. Vol. 96-97 (October). Pp. 73-78.
[7] Turkyilmazoglu M. Analytic heat and mass transfer of the mixed hydrodynamic/thermal slip MHD viscous flow over a stretching sheet // Int. J. Mech. Sci. 2011. Vol. 53. No. 10. Pp. 886-896.
[8] Zheng L., Niu J., Zhang X., Gao Y. MHD flow and heat transfer over a porous shrinking surface with velocity slip and temperature jump // Math. Comput. Model. 2012. Vol. 56. No. 5-6. Pp. 133-144.
[9] Turkyilmazoglu M. MHD fluid flow and heat transfer due to a stretching rotating disk // Int. J. Therm. Sci. 2012. Vol. 51 (January). Pp. 195-201.
[10] Turkyilmazoglu M. Multiple analytic solutions of heat and mass transfer of magnetohydrody-namic slip flow for two types of viscoelastic fluids over a stretching surface // J. Heat Transfer. 2012. Vol. 134. No. 7. doi:10.1115/1.4006165.
[11] Kandelousi M.S. Effect of spatially variable magnetic field on ferrofluid flow and heat transfer considering constant heat flux boundary condition // Eur. Phys. J. Plus. 2014. Vol. 129. No. 11. P. 248.
[12] Sheikholeslami M., Abelman S., Ganji D.D.
Numerical simulation of MHD nanofluid flow and heat transfer considering viscous dissipation // Int. J. Heat Mass Transfer. 2014. Vol. 79 (December). Pp. 212-222.
[13] Turkyilmazoglu M. An analytical treatment for the exact solutions of MHD flow and heat over two-three dimensional deforming bodies // Int. J. Heat Mass Transfer. 2015. Vol. 90 (November). Pp. 781-789.
[14] Sheikholeslami M., Rashidi M.M., Ganji D.D. Numerical investigation of magnetic nanofluid forced convective heat transfer in existence of variable magnetic field using two phase model // J. Mol. Liq. 2015. Vol. 212 (December). Pp. 117-126.
[15] Sheikholeslami M., Ellahi R. Three di-
mensional mesoscopic simulation of magnetic field effect on natural convection of nanofluid // Int. J. Heat Mass Transfer. 2015. Vol. 89 (October). Pp. 799-808.
[16] Sheikholeslami M., Rashidi M.M. Effect of space dependent magnetic field on free convection of Fe3O4-water nanofluid // J. Taiwan Inst. Chem. Eng. 2015. Vol. 56 (November). Pp. 6-15.
[17] Sheikholeslami M., Abelman S. Two-phase simulation of nanofluid low and heat transfer in an annulus in the presence of an axial magnetic field // IEEE Trans. Nanotechnol. 2015. Vol. 14. No. 3. Pp. 561-569.
[18] Sheikholeslami M., Rashidi M.M., Ganji D.D. Effect of non-uniform magnetic field on forced convection heat transfer of Fe3O4-water nanofluid // Comput. Methods Appl. Mech. Eng. 2015. Vol. 294 (September). Pp. 299-312.
[19] Sheikholeslami M., Ganji D.D., Javed M.Y., Ellahi R. Effect of thermal radiation on magnetohydrodynamics nanofluid flow and heat transfer by means of two phase model // J. Magn. Magn. Mater. 2015. Vol. 374 (15 January). Pp. 36-43.
[20] Zhang C., Zheng L., Zhang X., Chen G.
MHD flow and radiation heat transfer of nanofluids in porous media with variable surface heat flux and chemical reaction // Appl. Math. Model. 2015. Vol. 39. No. 1. Pp. 165-181.
[21] Lin Y., Zheng L., Zhang X., Ma L., Chen G. MHD pseudo-plastic nanofluid unsteady flow and heat transfer in a finite thin film over stretching surface with internal heat generation // Int. J. Heat Mass Transfer. 2015. Vol. 84 (May). Pp. 903-911.
[22] Sheikholeslami M., Vajravelu K., Rashidi M.M. Forced convection heat transfer in a semi annulus under the influence of a variable magnetic field // Int. J. Heat Mass Transfer. 2016. Vol. 92 (January). Pp. 339-348.
[23] Sheikholeslami M., Rashidi M.M., Hayat T., Ganji D.D. Free convection of magnetic nanofluid considering MFD viscosity effect // J. Mol. Liq. 2016. Vol. 218 (June). Pp. 393-399.
[24] Sheikholeslami M., Hayat T., Alsaedi A. MHD free convection of Al2o3-water nanofluid considering thermal radiation: A numerical study // Int. J. Heat Mass Transfer. 2016. Vol. 96 (May). Pp. 513-524.
[25] Sheikholeslami M., Ashorynejad H.R., Rana P. Lattice Boltzmann simulation of nanofluid heat transfer enhancement and entropy generation // J. Mol. Liq. 2016. Vol. 214 (February). Pp. 86-95.
[26] Sheikholeslami M., Soleimani S., Ganji D.D. Effect of electric field on hydrothermal behavior of nanofluid in a complex geometry // J. Mol. Liq. 2016. Vol. 213 (January). Pp. 153-161.
[27] Sheikholeslami M., Rashidi M.M. Nonuniform magnetic field effect on nanofluid hydrothermal treatment considering Brownian motion and thermophoresis effects // J. Braz. Soc. Mech. Sci. Eng. 2016. Vol. 38. No. 4. Pp. 1171-1184.
[28] Anuar I., Roslinda N., Pop I. Magneto-hydrodynamic (MHD) flow and heat transfer due to a stretching cylinder // Energy Conversion and Management. 2008. Vol. 49. No. 11. Pp. 3265-3269.
[29] Vajravelu K., Prasad K.V., Santhi S.R. Axisymmetric magneto-hydrodynamic (MHD) flow and heat transfer at a non-isothermal stretching cylinder // Applied Mathematics and computation. 2012. Vol. 219. No. 8. Pp. 3993-4005.
[30] Butt A.S., Ali A. Entropy analysis of magnetohydrodynamic flow and heat transfer due to a stretching cylinder // Journal of the Taiwan Institute of Chemical Engineers. 2014. Vol. 45. No. 3. Pp. 780-786.
[31] Hayat T., Shafiq A., Alsaedi A. MHD axisymmetric flow of third grade fluid by a stretching cylinder // Alexandria Engineering Journal. 2015. Vol. 54. No. 2. Pp. 205-212.
[32] Zheng L., Zhang C., Zhang X., Zhang J. Flow and radiation heat transfer of a nanofluid over a stretching sheet with velocity slip and temperature jump in porous medium // J. Franklin Inst. 2013. Vol. 350. No. 5. Pp. 990-1007.
[33] Mukhopadhyay S. MHD boundary layer slip flow along a stretching cylinder // Ain Shams Eng. J. 2013. Vol. 4. No. 2. Pp. 317-324.
[34] Abbas Z., Rasool S., Rashidi M.M. Heat transfer analysis due to an unsteady stretching/ shrinking cylinder with partial slip condition and suction // Ain Shams Eng. J. 2015. Vol. 6. No. 3.
Pp. 939-945.
[35] Aziz A. A similarity solution for laminar thermal boundary layer over a flat plate with a convective surface boundary condition // Commun. Nonlinear Sci. Numer. Simul. 2009. Vol. 14. No. 4. Pp. 1064-1068.
[36] Makinde O.D. Similarity solution of hydromagnetic heat and mass transfer over a vertical plate with a convective surface boundary condition // Int. J. Phys. Sci. 2010. Vol. 5. No. 6. Pp. 700-710.
[37] Ishak A., Nor Yacob A., Bachok N. Radiation effects on the thermal boundary layer flow over a moving plate with convective boundary condition // Meccanica. 2011. Vol. 46. No. 4. Pp. 795-801.
[38] Hayat T., Iqbal Z., Qasim M., Obaidat S. Steady flow of an Eyring Powell fluid over a moving surface with convective boundary conditions // Int. J. Heat Mass Transfer. 2012. Vol. 55. No. 7-8. Pp. 1817-1822.
[39] Akbar N.S., Nadeem S., Ul Haq R., Khan Z.H. Radiation effects on MHD stagnation point flow of nanofluid towards a stretching surface with convective boundary condition // Chinese Journal of Aeronautics. 2013. Vol. 26. No. 6. Pp. 1389-1397.
[40] Nadeem S., Rashid Mehmood, Akbar N.S. Optimized analytical solution for oblique flow of a Casson-nano fluid with convective boundary conditions // Int. J. Thermal Sciences. 2014. Vol. 78. Pp. 90-100.
[41] Anjali Devi S.P., Ganga B. Effects of viscous and Joules dissipation on MHD flow, heat and mass transfer past a stretching porous surface embedded in a porous medium // Nonlinear Analysis: Modelling and Control. 2009. Vol. 14. No. 3. Pp. 14303-14314.
СВЕДЕНИЯ ОБ АВТОРАХ
ГАНЕш Н. вишну — сотрудник Академии высшего образования г. Карпагама, штат Тамилнад, Индия.
Pollachi Main Road, Eachanari, Coimbatore, Tamil Nadu 641021, India
Ганга Б. — Ph.D., доцент кафедры математики Провиденс-колледжа для женщин, деревня Бан-дишола, штат Тамилнад, Индия.
Spring Field, Coonoor, Bandishola, Tamil Nadu 643104, India
gangabhose@gmail.com
ХАким Абдул к.А. — Ph.D., доцент кафедры математики Колледжа искусства и науки Миссии Шри Рамакришны, штат Тамилнад, Индия.
Periyanaickenpalayam, Sri Ramakrishna Vidyalaya Post, Coimbatore, Tamil Nadu 641020, India
abdulhakeem6@gmail.com
САРАНьЯ Сарани — сотрудник Колледжа искусства и науки Миссии Шри Рамакришны, г. Коим-батур, штат Тамилнад, Индия.
Periyanaickenpalayam, Sri Ramakrishna Vidyalaya Post, Coimbatore, Tamil Nadu 641020, India
кАЛАИвАНАН Раджа M.Sc., сотрудник Колледжа искусства и науки Миссии Шри Рамакришны, г. Коимбатур, штат Тамилнад, Индия.
Periyanaickenpalayam, Sri Ramakrishna Vidyalaya Post, Coimbatore, Tamil Nadu 641020, India
© Санкт-Петербургский политехнический университет Петра Великого, 2016
