# Radiation-Conduction Interaction of Steady Streamwise Surface Temperature Variations on Vertical Free Convection

Authors: Md. Kutub Uddin; Rabindra Nath Mondal
Journal: ISRN Applied Mathematics
 (2011)

---

## Abstract

The combined effects of the steady free convective boundary-layer flow induced by a vertical heated surface in the presence of sinusoidal surface temperature variations about a constant mean value with the effect of radiation are examined. The problem is studied using fully numerical techniques. The surface rate of heat transfer eventually alternates in sign with distance from the leading edge, but no separation occurs unless the amplitude of the thermal modulation is sufficiently high. Numerical results are obtained for different values of the physical parameters, the radiation parameterRd, the Prandtl number Pr, and the surface temperature wave amplitude a. It is found that both the local shear stress and the rate of heat transfer decrease when values of Rd increase.


---

## Body

## 1. Introduction

It is well known that power-law surface temperature distributions (and also power-law surface heat fluxes) give rise to self-similar boundary layer flows [1, 2].  Various papers have been published which deal with the effects of surface variations; for example, Yao [3] and Moulic and Yao [4, 5] have sought to investigate the effects of streamwise surface undulations of free and mixed convection from vertical surfaces held at uniform temperatures. Chiu and Chou [6], Hossain et al. [7], and Kim [8] extended these analyses to micropolar fluids, magnetohydrodynamic convection, and non-Newtonian convection, respectively. In a series of papers, Rees and Pop [9–12] and Rees [13] have also considered a large variety of analogous flows in porous media. Of these, only Rees [13] has been concerned with the effect of sinusoidal surface temperature variations, although in that case the surface variations were spanwise, thereby giving rise to a three-dimensional flow field. In the present paper, we are interested in another form of surface variation, namely, sinusoidal variations about a mean temperature which is held above the ambient temperature of the fluid. As by Rees [14], this type of surface distribution may be taken as a simplified model of the effects of a periodical array of heaters behind or within the heated surface, an accurate analysis of such a configuration requires a detailed investigation of the effects of solid conduction within the heated surface, but the aim of the present work is to simplify the problem by imposing a surface temperature distribution. In this way, we can determine a large amount of information about the resulting flow using the numerical techniques.Radiative convective flows are encountered in many industrial and environmental processes, for example, heating and cooling chambers, fossil fuel combustion energy processes, evaporation from large open water reservoirs, astrophysical flows, solar power technology, and space vehicle re-entry. Mathematically, the equations for radiative heat transfer with absorption, scattering, and emission can be generated by one of two approaches, namely, the continuum model or the spectral radiative treatment of a single particle. Details of the derivation of the general equation of radiative heat transfer are provided in the classic monograph by Chandrasekhar [15].Little is currently known about the boundary layer flows of radiating fluids. The inclusion of conduction-radiation effects in the energy equation, however, leads to a more highly nonlinear partial differential equation. The majority of studies concerned with the interaction of thermal radiation and natural convection were made by Sparrow and Cess [2], Cess [16], Arpaci [17], Cheng and Ozisik [18], Hasegawa et al. [19, 20], and Bankston et al. [21] for the case of a vertical semi-infinite plate. Recently, the problem of natural convection-radiation interaction on boundary layer flows with the Rosseland diffusion approximation has been studied by Hossain and Alim [22] and Hossain et al. [7]. Hossain and Rees [23] investigated the effect of radiation-conduction interaction in the mixed convective flow along a slender impermeable vertical cylinder. The inclusion of radiation terms is complicated and the resulting equations are very difficult to solve.  Grief et al.  [24] showed that in the optically thin limit, the physical situation can be simplified and they derived an exact solution of the problem of fully developed radiating laminar convection flow in an infinite vertical heated channel.  Their analysis was based on the work by Cogley et al. [25]. In the optically thin limit, the fluid does not absorb its own emitted radiation but the fluid absorbs the radiation emitted by the boundaries. It was shown by Cogley et al. [25] that in optically thin limit for a gray-gas near equilibrium, the following relation holds:
(1.1)∂qr∂y=4(T-Tw)∫0∞κλw(∂ebλ∂T)wdλ=4(T-Tw)I,whereI=∫0∞κλw(∂ebλ∂T)wdλ.
Here qr, is the radiative flux, κλ is the absorption coefficient, ebλ is the Planck function, and the subscript w represents the value of a quantity at the wall. Further simplification may be made concerning the spectral properties of radiating gases, but this is not essential for the present analysis. It should be mentioned that Soundalgekar [26]   considered the radiative-free convective flow of an optically thin grey-gas past a semi-infinite vertical plate. However, the Rosseland model is valid for isotropic local intensity, and high optical density of the medium and the radiative heat flux is given by [27](1.2)qr=-16σ∇T33(aR+σs)∇T,
where T denotes the temperature, aR is the Rosseland absorption coefficient, σs is the scattering coefficient, and σ is the Stefan-Boltzmann constant. The thermal boundary-layer equation can be written as(1.3)ρCp(u∂T∂x+v∂T∂y)=∂∂y[(16σT33aR+κ)∂T∂y].
The present paper considers, in detail, how the combined effects of surface radiation and sinusoidal surface temperature profiles in the streamwise direction modify the otherwise self-similar boundary-layer flow. Solutions are presented in terms of the surface rate of heat transfer and shear stress and detailed isotherms are also given. An important feature of the flow is that a near-wall layer develops at large distances downstream of the leading edge. The numerical evidence suggests that this inner layer decreases in thickness with distance downstream. A finite difference method is employed in obtaining the numerical solutions. In this paper, we investigate the combined effects of surface temperature variations and radiation on the steady boundary-layer flow of a Newtonian fluid from a heated vertical surface. The effects of varying different physical parameters on the local skin-friction and local rate of heat transfer are also presented.


## 2. Governing Equations and Boundary-Layer Analysis

A steady two-dimensional laminar-free convection flow of the boundary layer induced by a heated semi-infinite surface immersed in an incompressible Newtonian fluid is considered. In particular, the heated surface is maintained at the steady temperature, fluid having temperature,T, from a vertical cone. The physical coordinates (x,y) are chosen such that x is measured from the leading edge in the stream-wise direction and y is measured normal to the surface of the cone. The coordinate system, velocity direction, and the gravity orientation are shown in Figure 1.Figure 1
Physical model and the coordinate system.The boundary layer form of the equations for flow is(2.1)T=T∞+(Tw-T∞)(1-asin⁡(πx̂d)),
where T∞ is the ambient fluid temperature, Tw is the mean-surface temperature which is such that Tw>T∞, a is the relative amplitude of the surface temperature variations, and 2d is the wavelength of the variations. After a suitable nondimensionalization, the steady two-dimensional equations of motion are given by(2.2)ux+vy=0,(2.3)uux+vuy=-px+Gr-1/2(uxx+uyy)+θ,(2.4)uvx+vvy=-py+Gr-1/2(vxx+vyy)+θ,(2.5)uθx+vθy=αGr-1/2[θyy+16σ3κ(a+σs){θ3θy}y],θ=T-T∞Tw-T∞,
where Gr is the Grashof number and σ is the Prandtl number. In the derivation of (2.2), the Boussinesq approximation has been assumed. We note that the Grashof number has been based on d, half the dimensional wavelength of the thermal waves. In the equations, u and v are, respectively, the velocity components in the x and y directions, T is the fluid temperature, νis the kinematic viscosity, β is the thermal expansion coefficient, α is the thermal diffusivity, κ is the thermal conductivity, a is the Rosseland mean absorption coefficient, σ is the Stephan-Boltzman constant, and σs is the scattering coefficient.When the surface temperature is uniform and the Grashof number is very large, the resulting boundary-layer flow is self-similar. But the presence of sinusoidal surface temperature distributions, such as that given by (2.1), renders the boundary-layer flow nonsimilar. The boundary-layer equations are obtained by introducing the scaling:
(2.6)u=u*,v=Gr-1/4v*,x=x*,y=Gr-1/4y*,p=Gr-1/2p*,θ=θ*,
into (2.2), formally letting Gr become asymptotically large and retaining only the leading order terms. Thus, we obtain (2.7)ux+vy=0,(2.8)uux+vuy=uyy+θ,(2.9)uvx+vvy=-py+vyy,(2.10)uθx+vθy=αGr-1/2[θyy+16σ3κ(a+σs){θ3θy}y],
where the asterisk superscripts have been omitted for clarity of presentation. Equation (2.9) serves to define the pressure field in terms of the two velocity components and is decoupled from the other three equations. Therefore, we should not consider it further. As the equations are two-dimensional, we define a stream function ψ in the usual way: (2.11)u=ψy,v=-ψx,
and, therefore, (2.7) is satisfied automatically. Guided by the familiar self-similar form corresponding to a uniform surface temperature, we use the substitution
(2.12)ψ=x3/4f(η,x),θ=g(η,x),
where
(2.13)η=yx1/4
is the pseudosimilarity variable. Equations (2.8) and (2.9) reduce to (2.14)f′′′+g+34ff′′-12f′f′+x(fxf′′-fxf′)=0,(2.15)1Pr[{1+43Rd(1+(θw-1)g)3}g′]′+34fg′+x(fxg′-f′gx)=0,
and the boundary conditions are
(2.16)f=0,f′=0,g=1+asin⁡πx,atη=0,f′g⟶0asη⟶∞.
In (2.14)–(2.16), primes denote derivatives with respect to η.


## 3. Numerical Solutions

The parabolic system of (2.14)–(2.16) is nonsimilar and its numerical solution must be obtained using a marching method. The results presented here were obtained using the Keller-box method. After reducing (2.14)-(2.15) to first-order form in η, the subsequent second-order accurate discretisation based halfway between the grid points in both the η- and x-directions yields a set of nonlinear difference equations which are solved using a multidimensional Newton-Raphson iteration scheme. The results presented in Figures 2–5 are based on uniform grids in both coordinate directions. There were 201 grid points lying between η=0 and η=20 and 401 between x=0 and x=20. We restrict the presentation of our results to the two values of the Prandtl number, Pr = 0.7 (air) and Pr = 7.0 (water).(a) Skin friction and (b) rate of heat transfer againstξ for a = 0.0, 0.2, 0.4, 0.6, 0.8, 1.0 at Pr = 0.7.
(a)(b)(a) Skin friction and (b) rate of heat transfer againstξ for a = 0.0, 0.2, 0.4, 0.6, 0.8, 1.0 for Pr = 7.0.
(a)(b)(a) Skin friction and (b) rate of heat transfer againstξfor Rd=0.0, 1.0, 5.0, 10 for Pr = 0.7.
(a)(b)(a) Skin friction and (b) rate of heat transfer againstξ for Rd=0.0, 1.0, 5.0, 10 for Pr = 7.0.
(a)(b)Figure2(a) shows the evolution of f′′(η=0) with x, a scaled surface shear stress, for various values of the temperature wave amplitude, a, and the constant radiation parameter Rd  for Pr = 0.7. The corresponding rates of surface heat transfer are shown in Figure 2(b). Some aspects of the overall behavior of these curves may be explained by observing that the boundary layer is thinner when the surface temperature is relatively high and thicker when it is low. This arises because relatively high surface temperatures induce relatively large upward fluid velocities with the consequent increase in the rate of entertainment into the boundary layer. This causes, in turn, a thinning of the boundary layer. Thus, we should expect high shear stresses and rates of heat transfer at, or perhaps just beyond, where the surface temperature attains its maximum values. There is an obvious qualitative difference between the curves shown in Figure 2(a) and those in Figure 2(b). It is seen that as x increases, the amplitude of oscillation of the shear stress curves decays slowly, whereas the amplitude of heat transfer curves increases with x. Indeed, the curves in Figure 2(b) suggest that, whatever the value of a is, there will always be a value of x beyond which some part of the rate of the heat transfer curve between successive surface temperature maxima will be positive. This somewhat unusual phenomenon for boundary layer flows may be explained by noting that when relatively hot fluid encounters a relatively cold part of the heated surface, the overall heat transfer will occur from the fluid into the surface rather than the other way around.Figure3(a) shows the surface shear stress for various values of the temperature wave amplitude a and constant radiation parameter Rd for Pr = 7.0. The corresponding rates of surface heat transfer are shown in Figure 3(b). The same result, as shown above, is also happening in this case. Besides this, shear stress and rates of surface heat transfer are decreasing at Pr = 7.0 than that at Pr = 0.7. Figure 4(a) presents surface shear stress for various values of radiation parameter Rd and constant temperature wave amplitude a for Pr = 0.7. The corresponding rates of surface heat transfer are shown in Figure 4(b). This figure describes that when the radiation parameter is increasing both the shear stress and rates of surface heat transfer are decreasing. It might be mentioned here that whenever the value of radiation parameter Rd=0.0,  that is, when there is no effect of radiation, then these results have an excellent agreement with the results of Rees [14]. Figure 5(a) shows surface shear stress for various values of radiation parameter Rd and constant temperature wave amplitude a for Pr = 7.0. The corresponding rates of surface heat transfer are shown in Figure 5(b). Here, we observe that both the shear stress and rates of surface heat transfer are decreasing significantly for Pr = 7.0 than that for Pr = 0.7.In Figure6, we show the isotherms for Pr = 0.7 and a=0.2 for Rd=0.0,5.0, and 10.0. Here, we see that the boundary layer maintains its overall thickness in terms of ηwhen x is large, although variations in thickness are clearly presented when x is small. The thickness of the region, in which strong surface, induced temperature variations are present, reduces slowly in size as x increases. It is also clear that isolines are dense in the areas of high values (e.g., high temperature) and sparse in the areas with low intensity.Isotherms for Pr = 0.7 anda = 0.2. (a) Rd=0.0, (b) Rd=10.0.
(a)(b)In Table1, we show shear stress and surface rate of heat transfer for different values of Pr when Rd=0.0 and 5.0. As seen in Table 1, both the shear stress and the rate of surface heat transfer are decreasing due to increasing the value of Pr as well as radiation parameter Rd. Moreover, when Rd=0.0, that is, when there is no radiation, these results have an excellent agreement with the result found by Rees [14].Table 1
Shear stress and surface rate of heat transfer for different values of Pr whenRd=0.0 and 5.0.
Prf”g’Rd=0.0Rd=5.0Rd = 0.0Rd=5.00.011.394760.12113−0.05752−4.887600.021.355740.11860−0.07897−5.132630.051.284380.11860−0.12000−5.135330.101.215040.11861−0.16274−5.136260.201.132870.11862−0.21772−5.136350.501.008550.11863−0.31195−5.136040.700.960120.11863−0.35320−5.135531.000.908190.11864−0.40103−5.134962.000.807890.11865−0.50661−5.134383.000.750770.11865−0.57667−5.133864.000.711290.11866−0.63042−5.133405.000.681360.11866−0.67458−5.133016.000.657400.11866−0.71234−5.132697.000.637510.11866−0.74551−5.132458.000.620560.11867−0.77519−5.1322610.00.592840.11867−0.82684−5.1321420.00.511900.11867−1.00519−5.1321540.00.439210.11866−1.21461−5.1324770.00.386650.11865−1.41044−5.13332100.00.355970.11862−1.54955−5.13484


## 4. Conclusions

In this paper, the combined effects of radiaton-conduction interaction with steady streamwise surface temperature variations on vertical free convection have been investigated numerically by using a finite difference method. The effects of variations in the surface temperature wave amplitude and the Prandtl number on the shear stress and rate of surface heat transfer have been given.We calculated the surface shear stress (skin friction) and rate of heat transfer for various values of the temperature wave amplitude at a constant radiation parameter for the Prandtl numbers Pr = 0.7 and  Pr = 7.0. It is seen that asx increases, the amplitude of oscillation of the shear stress curves decays slowly, whereas the amplitude of heat transfer curves increases with x. We then calculated the surface shear stress (skin friction) and the rate of heat transfer for various values of radiation parameter at a constant temperature wave amplitude for Pr = 0.7 and Pr = 0.7 and it is found that when the radiation parameter is increasing both the shear stress and rates of surface heat transfer are decreasing, and at Rd=0.0, that is, when there is no effect of radiation, the results have an excellent agreement with the results of Rees [14]. It is also found that both the shear stress and rates of surface heat transfer are decreasing significantly for Pr = 0.7 than that at Pr = 0.7.



---

*Source: 927416-2011-12-14.xml*