MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi Publishing Corporation 946495 10.1155/2014/946495 946495 Research Article Oscillatory Free Convection about a Horizontal Circular Cylinder in the Presence of Heat Generation http://orcid.org/0000-0002-3377-0488 Zainuddin Najahulfazliah 1 Ismoen Muhaimin 1 Roslan Rozaini 1 Hashim Ishak 2,3 Abd El Aziz Mohamed 1 Centre for Research in Computational Mathematics Faculty of Science, Technology and Human Development Universiti Tun Hussein Onn 86400 Batu Pahat, Johor Malaysia uthm.edu.my 2 Centre for Modelling & Data Analysis School of Mathematical Sciences Universiti Kebangsaan Malaysia (UKM) 43600 Bangi Selangor Malaysia ukm.my 3 Department of Mathematics Faculty of Science King Abdulaziz University P.O. Box 80257 Jeddah 21589 Saudi Arabia kau.edu.sa 2014 2242014 2014 04 02 2014 31 03 2014 22 4 2014 2014 Copyright © 2014 Najahulfazliah Zainuddin et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

The oscillatory free convection about a horizontal circular cylinder in a Newtonian fluid in the presence of heat generation is investigated numerically by using the finite difference method. The surface temperature of the cylinder oscillates harmonically about the temperature of the surrounding medium. Heat is generated internally within the Newtonian fluid at a rate proportional to a power of the temperature difference. It is found that the presence of heat generation significantly increases the temperature and velocity distribution. The effects of the heat generation parameter and the Prandtl number on the surface rate of heat transfer, in terms of the local Nusselt number, and the surface shear stress, in terms of the local skin friction, are shown graphically from the stagnation point of the circular cylinder.

1. Introduction

The study of oscillatory free convective flow is very important in the engineering field of fluid mechanics. This is because such flows often occur in practice, while the heat generation is important for physical problems. Much information has been obtained about free convection over a circular cylinder, but less has been found about when this takes place in the presence of heat generation. Besides that, there are many related problems with the presence of heat generation but where the surface temperature is uniformly constant or steady. Vajravelu and Hadjinicolaou  presented heat transfer characteristics in a laminar boundary layer flow of viscous fluid over a linearly stretching continuous surface with viscous dissipation or frictional heating and internal heat generation. In this study, they considered the volumetric rate of heat generation, q [ W / m 3 ] , to be (1) q ′′′ = { Q 0 ( T - T ) , for T T , 0 , for T < T , where Q 0 is the heat generation constant. The above relation is valid for the state of some exothermic processes having T as the onset temperature. This study was extended by Vajravelu  to natural convection on a heated semi-infinite vertical plate with internal heat generation. The effects of heat generation over a plate, a cylinder, and a sphere have been investigated by researchers. The work conducted by Vajravelu and Hadjinicolaou  has been continued by other researchers. Chamkha and Issa  studied the effects of heat generation or absorption on a hydromagnetic flow over a flat surface, while Molla et al.  studied the natural convection flow along a vertical wavy surface with a uniform surface temperature in the presence of heat generation or absorption. Meanwhile, Ferdousi and Alim  considered the effect of heat generation on natural convection flow from a porous vertical plate.

Some works have also appeared on natural convection in the presence of heat generation over a horizontal circular cylinder. Molla et al.  studied an isothermal horizontal circular cylinder, and this work was extended in  to the existence of uniform heat flux. Molla et al.’s work [6, 7] was reviewed by Cheng  in relation to the study of a natural convection boundary layer on a horizontal elliptical cylinder with constant heat flux and internal heat generation, while the effect of radiation on a natural convection flow over a sphere in the presence of heat generation was investigated by Miraj et al. .

Unsteady free convection about a circular cylinder where the temperature of the surface either oscillates or is constant was studied in application technology. The oscillatory circular cylinder problem can be solved either analytically or numerically. Merkin  was the first who studied this related problem in oscillatory free convection for an infinite horizontal cylinder. He solved it analytically by using a matching technique between the inner and outer boundary layers. Continuing the work of Merkin , Chatterjee and Debnath  solved the problem by using the technique of separating the flow field into a steady and an unsteady part from the beginning in the governing equations. Ganesan and Rani  presented the problem of transient natural convection along a vertical cylinder with heat and mass transfer and solved it by using the implicit finite difference Crank-Nicolson scheme that had been used by Rani and Devaraj  in their study of the unsteady flow past a vertical cylinder with temperature oscillations. Jaman and Hossain  investigated the effect of fluctuating surface temperature on natural convection flow over cylinders of elliptic cross section and solved the problem by the most efficient and accurate implicit finite difference method together with the Keller-box elimination technique (also known as the Keller-box method).

Some work on free convection over other types of geometric shapes was carried out by Verma and Singh , who studied the effect of plate oscillations on horizontal free convection flow. Meanwhile, Hossain et al.  presented the heat transfer response of MHD free convection flow from a vertical heated plate to an oscillating surface heat flux. The solutions for general frequencies were obtained using the Keller-box method. The partial differential equations were reduced to a pair of linear ordinary differential equations and the solutions were obtained using a straightforward shooting method. A study of unsteady free convection heat and mass transfer in Walters-B viscoelastic flow past a semi-infinite vertical plate has been carried out by Prasad et al.  using an efficient, accurate, and implicit finite difference scheme of the Crank-Nicolson type. This method was originally developed by Prasad et al.  in their study of radiation and mass transfer effects on two-dimensional flow past an impulsively started infinite vertical plate. Udin and Kumar  studied unsteady free convection in a fluid past an inclined plate immersed in a porous medium numerically by using explicit finite difference, as demonstrated in Carnahan et al. . Recently, Kasim et al.  considered the unsteady free convection flow over a three-dimensional stagnation point with internal heat generation or absorption; the governing equations were solved by an implicit finite difference method.

Oscillatory free convection about a circular cylinder in the presence of heat generation is new and has not yet been studied. The governing equations and boundary equations are first transformed into a dimensionless form. The resulting system of equations is solved by a finite difference Crank-Nicolson scheme.

2. Formulation of Problem

Two-dimensional oscillatory free convection on a circular cylinder of radius a in presence of heat generation γ which is aligned in an ambient fluid of temperature T -    is considered. The physical configuration considered is as shown in Figure 1. The surface temperature of the cylinder is assumed to be T - s    where T - s > T - .. Extending the work of Molla et al.  and Merkin , who considered the effect of heat generation on free convection, and oscillatory free convection on a cylinder, respectively, the equations for mass continuity, momentum, and energy are governed by (2) u x + v y = 0 , u t - + u u x + v u y = v ( 2 u y 2 ) + g β ( T - - T - ) sin ( x a ) , T - t - + u T - x + v T - y = k ρ C ρ ( 2 T - y 2 ) + Q 0 ρ C ρ ( T - - T - ) , where x is measured along the surface of the circular cylinder and y is measured normal to the surface while t - is the time and u and v denote the velocity components in the x and y directions, respectively. The gravitational acceleration g acts downward and β is the coefficient of thermal expansion. ( x / a ) is the angle of the y -axis. T - is the temperature of fluid and T - is the temperature that is far away from the cylinder. k is the thermal conductivity, ρ is the fluid density, C ρ is the specific heat at constant pressure, and Q 0 is a heat generation constant. The boundary conditions are defined as follows: (3) u = v = 0 , T - = ( T - s - T - ) cos ω t - + T - s at    y = 0 , u 0 , T - T - at    y . Here T - s is the temperature of the cylinder. ω is the frequency parameter. The above equations become dimensionless by substituting the following nondimensional variables: (4) X = x a , Y = y a Gr 1 / 4 , t - = v a 2 Gr 1 / 2 t , T = ( T - - T - ) ( T - s - T - ) , U = a υ Gr 1 / 2 u , V = a υ Gr - 1 / 4 v , Gr = g β ( T - s - T - ) υ 2 a 3 , υ = μ ρ , Pr = μ C ρ k , γ = Q 0 a 2 μ C ρ Gr 1 / 2 , where υ is the reference kinematic viscosity, Gr is the Grashof number, Pr is the Prandtl number, and γ is the heat generation parameter. By substituting the variables (4) into (2), the equations are reduced to the following form: (5) U X + V Y = 0 , (6) U t + U U X + V U Y = 2 U Y 2 + T sin X , (7) T t + U T X + V T Y = 1 Pr 2 T Y 2 + γ T . By assuming τ = ω t - , the corresponding initial and boundary conditions in dimensionless form are as follows: (8) U = V = 0 , T = 1 + cos τ at    Y = 0 , U 0 ,    T T 0 as    Y .    By using the nondimensional quantities (4), the surface rate of heat transfer in terms of the local Nusselt number,and the surface shear stress in terms of the local skin friction are as follows: (9) C f Gr 1 / 4 = X ( U Y ) Y = 0 , Nu Gr - 1 / 4 = - ( T Y ) Y = 0 .

3. Solution Procedure

The nonlinear partial differential equations (5) and (6) under the initial and boundary conditions (8) are solved numerically by an implicit finite difference Crank-Nicolson scheme. The region of integration is considered where an axis along the surface X max = π and an axis normal to the surface Y max = 30.0 corresponding to    Y . There are M and N grid spacing in the X and Y    directions, respectively. The mesh points in the ( X , Y )    plane by X i and Y j , where i = 0,1 , 2,3 , , M and j = 0,1 , 2,3 , , N .    The mesh sizes have been fixed as Δ X = π / 150 , Δ Y = 0.12 with time step Δ τ = 0.0 1 . The iterative procedure is stopped when the time steps reach τ max = 3 0 . The numerical values of the variables such as U , V and the T temperature are obtained at specific points. The finite difference equations corresponding to (5), (6) and (7) are as follows: (10) 1 2 ( V i , j k + 1 - V i , j - 1 k + 1 Δ Y + V i , j k - V i , j - 1 k Δ Y ) = 1 2 ( U i , j k + 1 - U i - 1 , j k + 1 Δ X + U i , j k - U i - 1 , j k Δ X ) , (11) U i , j k + 1 - U i , j k Δ t + 1 2 [ U i , j k ( U i , j k + 1 - U i - 1 , j k + 1 + U i , j k - U i - 1 , j k Δ X ) v v v v v v v v v v v v v v v v + V i , j k ( U i , j k + 1 - U i , j - 1 k + 1 + U i , j k - U i , j - 1 k Δ Y ) ] = 1 2 [ ( U i , j + 1 k + 1 - 2 U i , j k + 1 + U i , j - 1 k + 1 + U i , j + 1 k - 2 U T i , j k + U i , j - 1 k Δ Y 2 ) v v v v v v v v v v v - γ ( T i , j k + 1 + T i , j k ) sin X i ( U i , j + 1 k + 1 - 2 U i , j k + 1 + U i , j - 1 k + 1 + U i , j + 1 k - 2 U T i , j k + U i , j - 1 k Y 2 ) ] , (12) T i , j k + 1 - T i , j k Δ t + 1 2 [ U i , j k ( T i , j k + 1 - T i - 1 , j k + 1 + T i , j k - T i - 1 , j k Δ X ) v v v v v v v v v v v v v v v v v v + V i , j k ( T i , j k + 1 - T i , j - 1 k + 1 + T i , j k - T i , j - 1 k Δ Y ) ] = 1 2 [ 1 Pr ( T i , j + 1 k + 1 - 2 T i , j k + 1 + T i , j - 1 k + 1 + T i , j + 1 k - 2 T i , j k + T i , j - 1 k Δ Y 2 ) v v v v v v v v v v + γ ( T i , j k + 1 + T i , j k ) ( T i , j + 1 k + 1 - 2 T i , j k + 1 + T i , j - 1 k + 1 + T i , j + 1 k - 2 T i , j k + T i , j - 1 k Y 2 ) ] . The coefficients U i , j k and V i , j k in (11) and (12) are treated as constant at any one time step. At time τ = 0 , the values of T i , j k , U i , j k , and V i , j k are known while the values of T i , j k + 1 , U i , j k + 1 , and V i , j k + 1 for the next time step are unknown and calculated by using Crank-Nicolson method.

4. Results and Discussion

The effect of heat generation on free convection on a circular cylinder with oscillating surface temperature was investigated. Solutions were obtained for several Prandtl numbers, Pr = 0.7 , 2.0 , 3.0 , 7.0 and different values of the heat generation parameter γ = 0.0 , 0.2 , 0.4 , 0.6 , 0.8 .

A grid independence test was carried out when τ = 6.00 with γ = 0.2 and Pr = 0.7 by using sets of grids in the range of 50 × 150 , 100 × 200 , and 150 × 250 mesh points. Table 1 shows that there are insignificant differences between sets of grids with 150 × 250 mesh points and above. Therefore, for all computations of similar problems in this section a 150 × 250 grid will be used.

Temperature values when τ = 6.00 at different M × N with γ = 0.2 and Pr = 0.7 .

 Y = 0.12 X 50 × 150 100 × 200 150 × 250 0 0.3713 0.3985 0.4159 π/6 0.5041 0.5429 0.5652 π/3 0.5432 0.5764 0.5992 π/2 0.6011 0.6361 0.6571 2π/3 0.6863 0.7266 0.7447 5π/6 0.8229 0.8468 0.8688 π 0.9910 1.0224 1.0397

For the steady state, the comparison of the numerical results of the present study with those of Molla et al.  is discussed in Table 2 without the effects of heat generation ( γ = 0.0 ) and with Prandtl number Pr = 1.0 . The results obtained start from the lower stagnation point of the circular cylinder X = 0 and go up to the upper stagnation point X i = π on the surface of the circular cylinder. The results of the present study and those of Molla et al.  are found to be in excellent agreement.

Comparison of the numerical results in the present study for the rate of heat transfer Nu Gr - 1 / 4 and the local skin friction coefficients C f Gr 1 / 4 with those of Molla et al.  for γ = 0.0 and Pr = 1.0 .

X Nu Gr - 1 / 4 C f Gr 1 / 4
Present study Molla et al.  Present study Molla et al. 
0 0.4221 0.4241 0.0000 0.0000
π/6 0.4193 0.4161 0.4064 0.4145
π/3 0.4050 0.4005 0.7392 0.7539
π/2 0.3810 0.3740 0.9562 0.9541
2π/3 0.3306 0.3355 0.9668 0.9696
5π/6 0.2481 0.2812 0.7523 0.7739
π 0.1975 0.1917 0.3245 0.3264

The temperature and velocity profiles are shown against the axial coordinate X at different values of the radial coordinate Y . Here, the numerical results are those obtained at time τ = 6.00    with γ = 0.2 and Pr = 0.7 . In Figure 2, the temperature profile increases with X . It is proved that the increasing values of X from the lower stagnation point of the circular cylinder X = 0 to the upper stagnation point X i = π lead to an increase in the temperature profile. In addition, the temperature profiles near the surface are higher, and they decrease as the distance from the surface increases. The velocity profiles shown in Figure 3 increase from the surface of the cylinder and decrease for some value of Y .

Temperature profile against X with γ = 0.2 , Pr = 0.7 , and τ = 6.00 .

Velocity profile against X with γ = 0.2 , Pr = 0.7 , and τ = 6.00 .

The temperature profiles against the time τ at various positions of the radial coordinate Y have been plotted in Figure 4 with γ = 0.2 and at X = π / 6 . The Prandtl number Pr = 0.7 is chosen. From Figure 4, it can be observed that the temperature profile oscillates harmonically near the surface of the cylinder and that there are no more oscillations as we move away from the surface. This is because the oscillating surface temperature makes the temperature of the fluid near the surface oscillate too. At a distance far away from the surface, the fluid temperature does not oscillate frequently because of the influence of the ambient temperature.

Temperature profile against time, τ , with γ = 0.2 and Pr = 0.7 .

The effect of heat generation γ = 0.0 , 0.2 , 0.4 , 0.6 , 0.8 on the temperature and velocity profiles with Pr = 0 . 7 at X = π / 6    and τ = 6.00 is illustrated in Figures 5 and 6. Both Figures 5 and 6 show that the increasing values of γ = 0.0 , 0.2 , 0.4 , 0.6 , 0.8 lead to increases in the temperature and the velocity profiles. The increase in the heat generation parameter causes the thermal boundary layer near the surface, causing the fluid temperature near the surface to increase. As γ increases, the increase in the temperature near the surface leads to an increase in the velocity gradient and enhances the fluid velocity. For γ = 0.0 , 0.2 , 0.4 , 0.6 , 0.8 , the maximum values of the velocity are shown to be 0.2720, 0.3125, 0.3532, 0.4060, and 0.4714, and these occur at Y = 0 . 84 with the corresponding temperatures of 0.8652, 0.9723, 1.1188, 1.3277, and 1.6241. Thus, it can be observed that the temperature at Y = 0.84 increases by 75.89% and the velocity increases by 19.94% as γ increases from 0.0 to 0.8. At Y = 1.80 , at which the velocity distribution tends to decrease are recorded to be 0.1548, 0.1867, 0.2197, 0.2566, and 0.2910 corresponding to temperatures of 0.2451, 0.3125, 0.4136, 0.5586, and 0.7443, so that as γ increases from 0.0 to 0.8, the velocity increases by 13.62% and the temperature by 49.92%.

Temperature profile against Y with Pr = 0.7 and τ = 6.00 at X = π / 6 .

Velocity profile against Y with Pr = 0.7 and τ = 6.00 at X = π / 6 .

The local values of the skin friction and the Nusselt number are obtained for different values of the heat generation parameter γ on the surface of the circular cylinder against the axial coordinate value X from ( X = 0 ) to ( X = π ) at time τ = 6.28 . The numerical results have been illustrated in Figures 7 and 8, respectively. It can be seen that an increase in the heat generation parameter γ leads to a decrease in the local Nusselt number NuGr - 1 / 4 which becomes negative for large values of γ . This is to be expected, since the heat generation mechanism creates a layer of hot fluid near the surface and hence the rate of heat transfer from the surface decreases. This phenomenon was discussed earlier in the section on Figures 5 and 6 as γ increases the fluid temperature near the surface increases. Since the fluid temperature near the surface increases, the viscosity of the fluid increases, leading the local skin friction    C f Gr 1 / 4 to increase.

Nusselt number for different values of γ with Pr = 0.7 at τ = 6.28 .

Skin friction for different values of γ with Pr = 0.7 at τ = 6.28 .

Figures 9 and 10 present the effect of different values of the Prandtl number Pr = 0.7 , 2.0 , 3.0 , 7.0 on the temperature and velocity profiles at X = π / 6 and τ = 6.00 with heat generation γ = 0.2 . These figures show that the fluid temperature and velocity distribution tend to be lower for higher Prandtl numbers. This is because, when the value of Pr is higher, the thermal boundary layer becomes thinner, leading to a decrease in the fluid temperature and in the velocity distribution. The reason for this is that Pr controls the relative thickness of the momentum and thermal boundary layers. As Pr decreases, the heat diffuses more quickly because the thickness of the thermal boundary layer is greater. This means that for gases, for which the Pr range is between 0.7 and 1.0, the thickness of the thermal boundary layer is much greater compared to water, for which the range of Pr is around 7.0 to 8.0. Furthermore, the temperature and velocity profiles decrease and approach zero as the coordinate of Y increases or the distance from the cylinder becomes greater.

Temperature profile against Y with γ = 0.2 and τ = 6.00 at X = π / 6 .

Velocity profile against Y with γ = 0.2 and τ = 6.00 at X = π / 6 .

Figure 11 shows the effects of the Prandtl number on the local skin friction C f Gr 1 / 4 and the heat transfer rate NuGr - 1 / 4 when γ = 0.2 at τ = 6.28 . The local skin friction decreases as Pr increases from 0.7 to 7.0 while the local heat transfer rate increases. Increasing values of Pr will decrease the shear stress because Pr gives thicker velocity profiles; this can be seen in Figure 8. This is the reason why the skin friction decreases. On the other hand, an increased value of Pr will increase the local heat transfer rate because when Pr becomes larger, the thermal boundary layer becomes thinner and the fluid temperature decreases.

Local skin friction, Nusselt number of Pr with γ = 0.2 at τ = 6.28 .

5. Conclusion

Oscillatory free convection on a circular cylinder in presence of heat generation has been investigated in this paper with the different values of Prandtl number Pr = 0.7,2.0,3.0,7.0    and  different heat generation parameters γ = 0.2 , 0.4 , 0.6 , 0.8. The governing boundary layer equations were transformed into a nondimensional form, and the resulting nonlinear systems of partial differential equations were solved numerically by an implicit finite difference scheme of Crank-Nicolson type. The solutions are obtained for the entire region of the circular cylinder, from the lower stagnation point ( X = 0 ) to upper stagnation point ( X = π ) . The effect of heat generation with different values of Pr on free convection about a circular cylinder with oscillating surface temperature has been analyzed. From the investigation, the interesting conclusion is summarized as follows.

The effects of heat generation parameter γ on velocity and temperature profiles as well as on the skin friction and the heat transfer rate have been investigated. An increase in the values of heat generation parameter γ leads to increase both velocity and temperature profiles as well as the skin frictions coefficient but the local rate of heat transfer decreases for Pr = 0.7 .

The increase in Prandtl number, Pr , with a fixed heat generation parameter leads to a decrease in all the velocity profile, the temperature profile, and the local skin friction coefficient, C f , but increases the local rate of heat transfer, Nu .

The steady flow existing around the circular cylinder has been found. The temperature profile oscillates harmonically near the surface of the cylinder and there are no more oscillations as we move away from the surface.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors would like to acknowledge the financial support received from the University Tun Hussein Onn Malaysia (UTHM) (Fundamental Research Grant Scheme FRGS/2/2010/SG/UTHM/02/4/0761) and the Universiti Kebangsaan Malaysia Grant no. DIP-2012-12.

Vajravelu K. Hadjinicolaou A. Heat transfer in a viscous fluid over a stretching sheet with viscous dissipation and internal heat generation International Communications in Heat and Mass Transfer 1993 20 3 417 430 2-s2.0-0027595994 10.1016/0735-1933(93)90026-R Vajravelu K. Natural convection at a heated semiinfinite vertical plate with internal heat generation Acta Mechanica 1979 34 1-2 153 159 2-s2.0-33745160513 10.1007/BF01176264 ZBL0437.76076 Chamkha A. J. Issa C. Effects of heat generation/absorption and thermophoresis on hydromagnetic flow with heat and mass transfer over a flat surface International Journal of Numerical Methods for Heat and Fluid Flow 2000 10 4 432 448 2-s2.0-0039768471 10.1108/09615530010327404 ZBL0948.76636 Molla M. M. Hossain M. A. Yao L. S. Natural convection flow along a vertical wavy surface with uniform surface temperature in presence of heat generation/absorption International Journal of Thermal Sciences 2004 43 2 157 163 2-s2.0-0742306425 10.1016/j.ijthermalsci.2003.04.001 Ferdousi A. Alim M. A. Natural convection flow from a porous vertical plate in presence of heat generation International University Journal Science Technology 2010 5 1 73 80 Molla M. M. Hossain M. A. Paul M. C. Natural convection flow from an isothermal horizontal circular cylinder in presence of heat generation International Journal of Engineering Science 2006 44 13-14 949 958 2-s2.0-33747787310 10.1016/j.ijengsci.2006.05.002 ZBL1213.76176 Mamun Molla M. Paul S. C. Anwar Hossain M. Natural convection flow from a horizontal circular cylinder with uniform heat flux in presence of heat generation Applied Mathematical Modelling 2009 33 7 3226 3236 2-s2.0-61449098515 10.1016/j.apm.2008.10.039 ZBL1205.76246 Cheng C.-Y. Natural convection heat transfer from a horizontal isothermal elliptical cylinder with internal heat generation International Communications in Heat and Mass Transfer 2009 36 4 346 350 2-s2.0-62549145676 10.1016/j.icheatmasstransfer.2009.01.004 Miraj M. Alim M. A. Mamun M. A. H. Effect of radiation on natural convection flow on a sphere in presence of heat generation International Communications in Heat and Mass Transfer 2010 37 6 660 665 2-s2.0-77953536804 10.1016/j.icheatmasstransfer.2010.01.013 Merkin J. H. Oscillatory free convection from an infinite horizontal cylinder Journal of Fluid Mechanics 1967 30 561 575 10.1017/S0022112067001612 ZBL0158.44702 Chatterjee A. K. Debnath L. Double boundary layers in oscillatory convective flow Il Nuovo Cimento B Series 11 1979 52 1 29 44 2-s2.0-51649169477 10.1007/BF02743567 Ganesan P. Rani H. P. Transient natural convection along vertical cylinder with Heat and Mass transfer Heat and Mass Transfer/Waerme- und Stoffuebertragung 1998 33 5-6 449 455 2-s2.0-0032048878 10.1007/s002310050214 Rani H. P. Devaraj R. Numerical solution of unsteady flow past a vertical cylinder with temperature oscillations Forschung im Ingenieurwesen/Engineering Research 2003 68 2 75 78 2-s2.0-0346342553 10.1007/s10010-003-0108-5 Jaman M. K. Hossain M. A. Effect of fluctuating surface temperature on natural convection flow over cylinders of elliptic cross section Open Transport Phenomena Journal 2010 2 1 35 47 2-s2.0-84859491341 10.2174/1877729501002010035 Verma R. L. Singh P. The effect of plate oscillations on horizontal free convection flow Australian Journal of Physics 1977 30 3 335 345 2-s2.0-0017501821 10.1071/PH770335 Hossain M. A. Das S. K. Pop I. Heat transfer response of MHD free convection flow along a vertical plate to surface temperature oscillations International Journal of Non-Linear Mechanics 1998 33 3 541 553 2-s2.0-0006589112 10.1016/S0020-7462(96)00151-5 ZBL0910.76080 Prasad V. R. Vasu B. Anwar O. Parsha R. Unsteady free convection heat and mass transfer in a Walters-B viscoelastic flow past a semi-infinite vertical plate: A Numerical Study International Journal of Thermal Sciences 2011 15 291 305 Prasad V. R. Reddy N. B. Muthucumaraswamy R. Radiation and mass transfer effects on two-dimensional flow past an impulsively started infinite vertical plate International Journal of Thermal Sciences 2007 46 12 1251 1258 2-s2.0-35148838172 10.1016/j.ijthermalsci.2007.01.004 Udin Z. Kumar M. Unsteady free convection in a fluid past an inclined plate immersed in a porous medium Computer Modelling and New Technologies 2010 14 41 47 Carnahan B. Luther H. A. Wilkes J. O. Applied Numerical Methods 1969 New York, NY, USA John Wiley & Sons Kasim A. R. M. Admon M. A. Shafie S. Free convection boundary layer flow of a viscoelastic fluid in the presence of heat generation World Academy of Science, Engineering and Technology 2011 75 492 499 2-s2.0-79953651964