Mixed convection cooling of a heated, continuously stretching surface

 An numerical study of the flow and heat transfer characteristics associated with a heated, continuously stretching surface being cooled by a mixed convection flow has been carried out. The relevant heat transfer mechanisms are of interest in a wide variety of practical applications, such as hot rolling, continuous casting, extrusion, and drawing. The surface velocity of the continuously stretching sheet was assumed to vary according to a power-law form, that is, u _w (x)=Cx ^p . Two conditions of surface heating were considered, a variable wall temperature (VWT) in the form T _w (x)−T _∞=Ax ⁿ and a variable surface heat flux (VHF) in the form q _w (x)=Bx ^m . The governing differential equations are transformed by introducing proper nonsimilarity variables and solved numerically using a procedure based on finite difference approximations. Results for the local Nusselt number and the local friction coefficient are obtained for a wide range of governing parameters, such as the surface velocity parameter p, the wall temperature exponent n, the surface heat flux exponent m, the buoyancy force parameters (ξ for the VWT case and χ for the VHF case), and Prandtl number of the fluid. It is found that the local Nusselt number is increased with increasing the velocity exponent parameter p for the VWT case, while the opposite trend is observed for the VHF case. The local friction coefficient is increased for a decelerated stretching surface, while it is decreased for an accelerated stretching surface. Also, appreciable effects of the buoyancy force on the local Nusselt number and the local friction coefficient are observed for both VWT and VHF cases, as expected. Received on 11 January 1999     

Unsteady free convection on a vertical cylinder with variable heat and mass flux

The unsteady natural convection boundary layer flow over a semi-infinite vertical cylinder is considered with combined buoyancy force effects, for the situation in which the surface temperature T ^′ _w(x) and C ^′ _w(x) are subjected to the power-law surface heat and mass flux as K(T ^′/r) = −ax ⁿ and D(C ^′/r) = −bx ^m . The governing equations are solved by an implicit finite difference scheme of Crank-Nicolson method. Numerical results are obtained for different values of Prandtl number, Schmidt number ‘n’ and ‘m’. The velocity, temperature and concentration profiles, local and average skin-friction, Nusselt and Sherwood numbers are shown graphically. The local Nusselt and Sherwood number of the present study are compared with the available result and a good agreement is found to exist. Received on 7 July 1998     

Hydromagnetic flow and heat transfer adjacent to a stretching vertical sheet

An analysis is made for the steady two-dimensional magneto-hydrodynamic flow of an incompressible viscous and electrically conducting fluid over a stretching vertical sheet in its own plane. The stretching velocity, the surface temperature and the transverse magnetic field are assumed to vary in a power-law with the distance from the origin. The transformed boundary layer equations are solved numerically for some values of the involved parameters, namely the magnetic parameter M, the velocity exponent parameter m, the temperature exponent parameter n and the buoyancy parameter λ, while the Prandtl number Pr is fixed, namely Pr = 1, using a finite difference scheme known as the Keller-box method. Similarity solutions are obtained in the presence of the buoyancy force if n = 2m−1. The features of the flow and heat transfer characteristics for different values of the governing parameters are analyzed and discussed. It is found that both the skin friction coefficient and the local Nusselt number decrease as the magnetic parameter M increases for fixed λ and m. For m = 0.2 (i.e. n = −0.6), although the sheet and the fluid are at different temperatures, there is no local heat transfer at the surface of the sheet except at the singular point of the origin (fixed point).     

Laminar Mixed Convection Adjacent to Vertical Continuously Stretching Sheets With Variable Viscosity and Variable Thermal Diffusivity

The paper presents a study of the laminar mixed convection adjacent to vertical continuously stretching sheets, taking into account the effects of variable viscosity and variable thermal diffusivity. The similarity solutions are reported for isothermal sheet moving with a velocity of the form u_w=Bx^0.5 and a continuous linearly stretching sheet with a linear surface temperature distribution. The equations of conservation of mass, momentum and energy, which govern the flow and heat transfer, are solved numerically by using the shooting method. The numerical results obtained for the flow and heat transfer characteristics reveal many interesting behaviors. The numerical results show that, variable viscosity, variable thermal diffusivity, the velocity exponent parameter, the temperature exponent parameter and the buoyancy force parameter have significant influences on the velocity and temperature profiles, shear stress and Nusselt number in two cases air and water.     

Combined heat and mass transfer along a vertical moving cylinder with a free stream

The mixed convection flow over a continuous moving vertical slender cylinder under the combined buoyancy effect of thermal and mass diffusion has been studied. Both uniform wall temperature (concentration) and uniform heat (mass) flux cases are included in the analysis. The problem is formulated in such a manner that when the ratio λ(= u _w/(u _w + u _∞), where u _w and u _∞ are the wall and free stream velocities, is zero, the problem reduces to the flow over a stationary cylinder, and when λ = 1 it reduces to the flow over a moving cylinder in an ambient fluid. The partial differential equations governing the flow have been solved numerically using an implicit finite-difference scheme. We have also obtained the solution using a perturbation technique with Shanks transformation. This transformation has been used to increase the range of the validity of the solution. For some particular cases closed form solutions are obtained. The surface skin friction, heat transfer and mass transfer increase with the buoyancy forces. The buoyancy forces cause considerable overshoot in the velocity profiles. The Prandtl number and the Schmidt number strongly affect the surface heat transfer and the mass transfer, respectively. The surface skin friction decreases as the relative velocity between the surface and free stream decreases. Received on 17 May 1999     

Finite difference analysis of laminar free convection flow past a non isothermal vertical cone

A finite-difference analysis for the transient free convection flow of an incompressible viscous fluid past a vertical cone with variable wall surface temperature T _w^′ (x) = T _∞^′ + a x ⁿ varying as power function of distance from the apex (x = 0) is presented here. The dimensionless governing equations of the flow that are unsteady, coupled and non-linear partial differential equations are solved by an efficient, accurate and unconditionally stable finite difference scheme of Crank-Nicolson type. The velocity and temperature fields have been studied for various parameters such as Prandtl number and n (exponent in power law variation in surface temperature). The local as well as average skin-friction and Nusselt number are also presented and analyzed graphically. The present results are compared with available results in literature and are found to be in good agreement.     

Heat transfer characteristics of boundary-layer flows induced by continuous surfaces stretched with prescribed skin friction

A continuous surface stretched with velocity u _w=u _w (x) and having the temperature distribution T _w=T _w (x) interacts with the viscous fluid in which it is immersed both mechanically and thermally. The thermal interaction is characterized by the surface heat flux q _w=q _w (x) and the mechanical one by the skin friction τ _w=τ _w (x). In the whole previous theoretical research concerned with such processes, either (u _w and T _w) or (u _w and q _w) have been prescribed as known boundary conditions. The goal of the present paper is to initiate the investigation of the boundary layer flows induced by stretching processes for which either (τ _w and T _w ) or (τ _w and q _w) are the prescribed quantities. The case of an isothermal surface stretched with constant skin friction, (τ _w=const., T _w=const. ≠ T _∞) is worked out in detail. The corresponding flow and heat transfer characteristics are compared to those obtained for the (well known) case of a uniformly moving isothermal surface (u _w=const., T _w=const. ≠ T _∞).     

Non-Darcian Forced Convection Flow of Viscous Dissipating Fluid over a Flat Plate Embedded in a Porous Medium

In this study, laminar boundary layer flow over a flat plate embedded in a fluid-saturated porous medium in the presence of viscous dissipation, inertia effect and suction/injection is analyzed using the Keller box finite difference method. The flat plate is assumed to be held at constant temperature. The non-Darcian effects of convection, boundary and inertia are considered. Results for the local heat transfer parameter and the local skin friction parameter as well as the velocity and temperature profiles are presented for various values of the governing parameters. The non-Darcian effects are shown to decrease the velocity and to increase the temperature. It is also shown that the local heat transfer parameter and the local skin friction parameter increase due to suction of fluid while injection reverses this trend. It is disclosed that the effect of the viscous dissipation for negative values of Ec (T _w < T _∞) is to enhance the heat transfer coefficient while the opposite is true for positive values of Ec (T _w > T _∞). The results are compared with those available in the existing literature and an excellent agreement is obtained.     

Laminar mixed convection from a continuously moving vertical surface with suction or injection

The boundary layer flow over a uniformly moving vertical surface with suction or injection is studied when the buoyancy forces assist or oppose the flow. Similarity solutions are obtained for the boundary layer equations subject to power law temperature and velocity boundary conditions. The effect is of various governing parameters, such as Prandtl number Pr, temperature exponent n, injection parameter d, and the mixed convection parameter λ=Gr/Re², which determine the velocity and temperature distributions and the heat transfer coefficient, are studied. The heat transfer coefficient increases as λ assisting the flow for all d at Pr=0.72 however, for n=−1 it decreases sharply with λ. On the other hand, increasing λ has no effect on heat transfer coefficient for Pr=10 at n=0, and 1 for almost all values of d studied. However, for n=−1 it has similar effect as for Pr=0.72. It is also found that Nusselt number increases as n increases for fixed λ and d. Received on 26 March 1997     

Heat transfer in the flow of a viscoelastic fluid over a stretching surface

An analysis is made of heat transfer in the boundary layer of a viscoelastic fluid flowing over a stretching surface. The velocity of the surface varies linearly with the distance x from a fixed point and the surface is held at a uniform temperature T _w higher than the temperature T _∞ of the ambient fluid. An exact analytical solution for the temperature distribution is found by solving the energy equation after taking into account strain energy stored in the fluid (due to its elastic property) and viscous dissipation. It is shown that the temperature profiles are nonsimilar in marked contrast with the case when these profiles are found to be similar in the absence of viscous dissipation and strain energy. It is also found that temperature at a point increases due to the combined influence of these two effects in comparison with its corresponding value in the absence of these two effects. A novel result of this analysis is that for small values of x, heat flows from the surface to the fluid while for moderate and large values of x, heat flows from the fluid to the surface even when T _w >T _∞. Temperature distribution and the surface heat flux are determined for various values of the Prandtl number P, the elastic parameter K ₁ and the viscous dissipation parameter a. Numerical solutions are also obtained through a fourth-order accurate compact finite difference scheme. Received on 14 October 1997     

Mixed convection flow of second grade fluid along a vertical stretching flat surface with variable surface temperature

In the present study we have explored the effects of thermal buoyancy on flow of a viscoelastic second grade fluid past a vertical, continuous stretching sheet of which the velocity and temperature distributions are assumed to vary according to a power-law form. The governing differential equations are transformed into dimensionless form using appropriate transformations and then solved numerically. The methods here employed are (1) the perturbation method together with the Shanks transformation, (2) the local non-similarity method with second level of truncation and (3) the implicit finite difference method for values of ξ ( = Gr _x /Re _x ², defined as local mixed convection parameter) ranging in [0, 10]. The comparison between the solutions obtained by the aforementioned methods found in excellent agreement. Effects of the elasticity parameter λ on the skin-friction and heat transfer coefficients have been shown graphically for the fluids having the values of the Prandtl number equal to 0.72, 7.03 and 15.0. Effects of the viscoelastic parameter and the mixed convection parameter, ξ, on the temperature and velocity fields have also been studied. We notice that with the increase in visco-elastic parameter λ, velocity decreases whereas temperature increases and that velocity gradient is higher than that of temperature. On leave of absence from the Department of Mathematics, University of Dhaka, Bangladesh.     

The optimal dimensions of convective-radiating circular fins

The optimal dimensions of convective-radiating circular fins with variable profile, heat-transfer coefficient and thermal conductivity, as well as internal heat generation are obtained. A profile of the form y=(w/2) [1+(r _o/r) ⁿ ] is studied, while variation of thermal conductivity is of the form k=k _o[1+ɛ((T−T _∞)/ (T _b−T _∞)) ^m ]. The heat-transfer coefficient is assumed to vary according to a power law with distance from the bore, expressed as h=K[(r−r _o)/(r _e−r _o)]^λ. The results for λ=0 to λ=1.9, and −0.4≤ɛ≤0.4, have been expressed by suitable dimensionless parameters. A correlation for the optimal dimensions of a constant and variable profile fins is presented in terms of reduced heat-transfer rate. It is found that a (quadratic) hyperbolic circular fin with n=2 gives an optimum performance. The effect of radiation on the fin performance is found to be considerable for fins operating at higher base temperatures, whereas the effect of variable thermal conductivity on the optimal dimensions is negligible for the variable profile fin. It is also observed, in general, that the optimal fin length and the optimal fin base thickness are greater when compared to constant fin thickness. Received on 22 February 1999     

Combined forced and natural convection in laminar film flow

A mathematical model for the flow and heat transfer in a gravity-driven liquid film is presented, in which the strict Boussinesq approximation is adopted to account for buoyancy. A similarity transformation reduces the governing equations to a coupled set of ordinary differential equations. The resulting two-parameter problem is solved numerically for Prandtl numbers ranging from 1 to 1000. Favourable buoyancy arises when the temperatureT _w of the isothermal surface is lower than the temperatureT ₀ of the incoming fluid, and the principal effects of the aiding buoyancy are to increase the wall shear and heat transfer rate. For unfavourable buoyancy (T _w>T ₀), the buoyancy force and gravity act in opposite directions and the flow in the film boundary layer decelerates, whereas the friction and heat transfer are reduced. The observed effects of buoyancy diminish appreciably for higher Prandtl numbers.     

Exact analytic solutions for free convection boundary layers on a heated vertical plate with lateral mass flux embedded in a saturated porous medium

 The problem of the self-similar boundary flow of a “Darcy-Boussinesq fluid” on a vertical plate with temperature distribution T _w(x) = T _∞+A·x ^λ and lateral mass flux v _w(x) = a·x ^(λ−1)/2, embedded in a saturated porous medium is revisited. For the parameter values λ = 1,−1/3 and −1/2 exact analytic solutions are written down and the characteristics of the corresponding boundary layers are discussed as functions of the suction/ injection parameter in detail. The results are compared with the numerical findings of previous authors. Received on 8 March 1999     

The effect of suction or injection on the boundary layer flows induced by continuous surfaces stretched with prescribed skin friction

The boundary layer flow and heat transfer on a stretched surface moving with prescribed skin friction is studied for permeable surface. Three major cases are studied for isothermal surface (n=0) stretched corresponding to different dimensional skin friction boundary conditions namely; skin friction at the surface scales as (x ^?1/2) at m=0, constant skin friction at m=1/3 and skin friction scales as (x) at m=1. The constants m and n are the indices of the power law velocity and temperature exponent respectively. Similarity solutions are obtained for the boundary layer equations subject to power law temperature and velocity variation. The effect of various governing parameters, such as Prandtl number Pr, suction/injection parameter f _w , m and n are studied. The results show that for isothermal surface increasing m enhances the dimensionless heat transfer coefficient for fixed f _w at the suction case and the reverse is true at the injection case. Furthermore, for fixed m, as f _w increases the dimensionless heat transfer coefficient increases. Large enhancements are observed in the heat transfer coefficient as the temperature boundary condition along the surface changes from uniform to linear where the dimensional skin friction is of order (x) at m=1. This enhancement decreases as the suction increases.     

Nonsimilarity solutions for non-Darcy mixed convection from horizontal surfaces in a porous medium

Non-Darcy mixed convection in a porous medium from horizontal surfaces with variable surface heat flux of the power-law distribution is analyzed. The entire mixed convection regime is divided into two regions. The first region covers the forced convection dominated regime where the dimensionless parameter ζ _f =Ra^* _x /Pe² _x is found to characterize the effect of buoyancy forces on the forced convection with K ^′ U _∞/ν characterizing the effect of inertia resistance. The second region covers the natural convection dominated regime where the dimensionless parameter ζ _n =Pe _x /Ra^*1/2 _x is found to characterize the effect of the forced flow on the natural convection, with (K ^′ U _∞/ν)Ra^*1/2 _x /Pe _x characterizing the effect of inertia resistance. To obtain the solution that covers the entire mixed convection regime the solution of the first regime is carried out for ζ _f =0, the pure forced convection limit, to ζ _f =1 and the solution of the second is carried out for ζ _n =0, the pure natural convection limit, to ζ _n =1. The two solutions meet and match at ζ _f =ζ _n =1, and R ^* _h =G ^* _h . Also a non-Darcy model was used to analyze mixed convection in a porous medium from horizontal surfaces with variable wall temperature of the power-law form. The entire mixed convection regime is divided into two regions. The first region covers the forced convection dominated regime where the dimensionless parameter ξ _f =Ra _x /Pe _x ^3/2 is found to measure the buoyancy effects on mixed convection with Da _x Pe _x /ɛ as the wall effects. The second region covers the natural convection dominated region where ξ _n =Pe _x /Ra _x ^2/3 is found to measure the force effects on mixed convection with Da _x Ra _x ^2/3/ɛ as the wall effects. Numerical results for different inertia, wall, variable surface heat flux and variable wall temperature exponents are presented. Received on 8 July 1996     

The buoyancy effects on the boundary layers induced by continuous surfaces stretched with rapidly decreasing velocities

Heat and mass transfer characteristics of the self-similar boundary layer flows induced by continuous surfaces stretched with rapidly decreasing power law velocities U_w x^m, m < –1 are considered for mixed convection flow. The effect of various governing parameters, such as Prandtl number Pr, temperature exponent n, dimensionless injection/suction velocity f_w, and the mixed convection parameter = s Gr/Re² are studied. These parameters have great effects on velocity and temperature profiles, heat transfer coefficient, and skin friction coefficient at the moving surface. Results show that similarity solutions exist only when the condition n = 2m – 1 is satisfied. Critical values of , Nu/Re^0.5 and C_f Re^0.5 are obtained for predominate natural convection for different Prandtl numbers at m = –2, –6 and n = –5, and –13 respectively. Results also show that the effect of buoyancy is more significant for weak than for strong suction. Furthermore, critical Prandtl numbers where f_w profiles have minimums are obtained for m = –2 and –6. Finally, critical values of , C_f Re^0.5 are also obtained for predominate natural convection for both m = –2 and –6.     

Non-Darcy mixed convection from a vertical plate in saturated porous media-variable surface heat flux

Nonsimilarity solutions for non-Darcy mixed convection from a vertical impermeable surface embedded in a saturated porous medium are presented for variable surface heat flux (VHF) of the power-law form. The entire mixed convection region is divided into two regimes. One region covers the forced convection dominated regime and the other one covers the natural convection dominated regime. The governing equations are first transformed into a dimensionless form by the nonsimilar transformation and then solved by a finite-difference scheme. Computations are based on Keller Box method and a tolerance of iteration of 10⁻⁵ as a criterion for convergence. Three physical aspects are introduced. One measures the strength of mixed convection where the dimensionless parameter Ra^* _x /Pe^3/2 _x characterizes the effect of buoyancy forces on the forced convection; while the parameter Pe _x /Ra^*2/3 _x characterizes the effect of forced flow on the natural convection. The second aspect represents the effect of the inertial resistance where the parameter K′U _∞/ν is found to characterize the effect of inertial force in the forced convection dominated regime, while the parameter (K′U _∞/ν)(Ra^*2/3 _x /Pe _x ) characterizes the effect of inertial force in the natural convection dominated regime. The third aspect is the effect of the heating condition at the wall on the mixed convection, which is presented by m, the power index of the power-law form heating condition. Numerical results for both heating conditions are carried out. Distributions of dimensionless temperature and velocity profiles for both Darcy and non-Darcy models are presented. Received on 26 May 1997     

Steady mixed convection boundary layer flow over a vertical flat plate in a porous medium filled with water at 4°C: case of variable wall temperature

The problem of steady mixed convection boundary layer flow over a vertical impermeable flat plate in a porous medium saturated with water at 4°C (maximum density) when the temperature of the plate varies as x ^m and the velocity outside boundary layer varies as x ^{2 m} , where x measures the distance from the leading edge of the plate and m is a constant is studied. Both cases of the assisting and the opposing flows are considered. The plate is aligned parallel to a free stream velocity U _∞ oriented in the upward or downward direction, while the ambient temperature is T _∞ = T _m (temperature at maximum density). The mathematical models for this problem are formulated, analyzed and simplified, and further transformed into non-dimensional form using non-dimensional variables. Next, the system of governing partial differential equations is transformed into a system of ordinary differential equations using the similarity variables. The resulting system of ordinary differential equations is solved numerically using a finite-difference method known as the Keller-box scheme. Numerical results for the non-dimensional skin friction or shear stress, wall heat transfer, as well as the temperature profiles are obtained and discussed for different values of the mixed convection parameter λ and the power index m. All the numerical solutions are presented in the form of tables and figures. The results show that solutions are possible for large values of λ and m for the case of assisting flow. Dual solutions occurred for the case of opposing flow with limited admissible values of λ and m. In addition, separation of boundary layers occurred for opposing flow, and separation is delayed for the case of water at 4°C (maximum density) compared to water at normal temperature.     

Thermal dispersion effects on non-Darcy natural convection over horizontal plate with surface mass flux

Summary The effect of surface mass flux on the non-Darcy natural convection over a horizontal flat plate in a saturated porous medium is studied using similarity solution technique. Forchheimer extension is considered in the flow equations. The suction/injection velocity distribution has been assumed to have power function form Bx ^l , similar to that of the wall temperature distribution Ax ⁿ , where x is the distance from the leading edge. The thermal diffusivity coefficient has been assumed to be the sum of the molecular diffusivity and the dynamic diffusivity due to mechanical dispersion. The dynamic diffusivity is assumed to vary linearly with the velocity component in the x direction, i.e. along the hot wall. For the problem of constant heat flux from the surface (n=1/2), similarity solution is possible when the exponent l takes the value −1/2. Results indicate that the boundary layer thickness decreases whereas the heat transfer rate increases as the mass flux parameter passes from the injection domain to the suction domain. The increase in the thermal dispersion parameter is observed to favor the heat transfer by reducing the boundary layer thickness. The combined effect of thermal dispersion and fluid suction/injection on the heat transfer rate is discussed. Received 7 December 1995; accepted for publication 7 January 1997