Heat transfer in a power-law fluid film over a unsteady stretching sheet

In this paper the problem of momentum and heat transfer in a thin liquid film of power-law fluid on an unsteady stretching surface has been studied. Numerical solutions are obtained for some representative values of the unsteadiness parameter S and the power-law index n for a wide range of the generalized Prandtl number, 0.001 ≤ Pr ≤ 1000. Typical temperature and velocity profiles, the dimensionless film thickness, free-surface temperature, and the surface heat fluxes are presented at selected controlling parameters. The results show that increasing the value of n tends to increase the boundary-layer thickness and broadens the temperature distributions. The free-surface temperature of a shear thinning fluid is larger than that of a Newtonian fluid, but the opposite trend is true for a shear thickening fluid. For small generalized Prandtl numbers, the surface heat flux increases with a decrease in n, but the impacts of n on the heat transfer diminish for Pr greater than a moderate value (approximately 1 ≤ Pr ≤ 10, depending on the magnitude of S).     

Self-similar solution of the unsteady mixed convection flow in the stagnation point region of a rotating sphere

An analysis is performed to present a new self-similar solution of unsteady mixed convection boundary layer flow in the forward stagnation point region of a rotating sphere where the free stream velocity and the angular velocity of the rotating sphere vary continuously with time. It is shown that a self-similar solution is possible when the free stream velocity varies inversely with time. Both constant wall temperature and constant heat flux conditions have been considered in the present study. The system of ordinary differential equations governing the flow have been solved numerically using an implicit finite difference scheme in combination with a quasilinearization technique. It is observed that the surface shear stresses and the surface heat transfer parameters increase with the acceleration and rotation parameters. For a certain value of the acceleration parameter, the surface shear stress in x-direction vanishes and due to further reduction in the value of the acceleration parameter, reverse flow occurs in the x–component of the velocity profiles. The effect of buoyancy parameter is to increase the surface heat transfer rate for buoyancy assisting flow and to decrease it for buoyancy opposing flow. For a fixed buoyancy force, heating by constant heat flux yields a higher value of surface heat transfer rate than heating by constant wall temperature.     

Transient MHD stagnation flow of a non-Newtonian fluid due to impulsive motion from rest

The transient boundary layer flow and heat transfer of a viscous incompressible electrically conducting non-Newtonian power-law fluid in a stagnation region of a two-dimensional body in the presence of an applied magnetic field have been studied when the motion is induced impulsively from rest. The non-linear partial differential equations governing the flow and heat transfer have been solved by the homotopy analysis method and by an implicit finite-difference scheme. For some cases, analytical or approximate solutions have also been obtained. The special interest are the effects of the power-law index, magnetic parameter and the generalized Prandtl number on the surface shear stress and heat transfer rate. In all cases, there is a smooth transition from the transient state to steady state. The shear stress and heat transfer rate at the surface are found to be significantly influenced by the power-law index N except for large time and they show opposite behaviour for steady and unsteady flows. The magnetic field strongly affects the surface shear stress, but its effect on the surface heat transfer rate is comparatively weak except for large time. On the other hand, the generalized Prandtl number exerts strong influence on the surface heat transfer. The skin friction coefficient and the Nusselt number decrease rapidly in a small interval 0<t^*<1 and reach the steady-state values for t^*≥4.     

Investigation of the different base fluid effects on the nanofluids heat transfer and pressure drop

A numerical study of laminar forced convective flows of three different nanofluids through a horizontal circular tube with a constant heat flux condition has been performed. The effect of Al₂O₃ volume concentration 0 ≤ φ ≤ 0.09 in the pure water, water-ethylene glycol mixture and pure ethylene glycol as base fluids, and Reynolds number of 100 ≤ Re ≤ 2,000 for different power inputs in the range of 10 ≤ Q(W) ≤ 400 have been investigated. In this study, all of the nanofluid properties are temperature and nanoparticle volume concentration dependent. The governing equations have been solved using finite volume approach with the SIMPLER algorithm. The results indicate an increase in the averaged heat transfer coefficient with increasing the mass of ethylene glycol in the water base fluid, solid concentration and Reynolds number. From the investigations it can be inferred that, the pressure drop and pumping power in the nanofluids at low solid volumetric concentration (φ < 3%) is approximately the same as in the pure base fluid in the various Reynolds numbers, but the higher solid nanoparticle volume concentration causes a penalty drop in the pressure. Moreover, this study shows it is possible to achieve a higher heat transfer rate with lower wall shear stress with the use of proper nanofluids.     

Simulation of heat and chemical reactions on Reiner Rivlin fluid model for blood flow through a tapered artery with a stenosis

In the present article, we have analyzed the effects of heat and mass transfer on Reiner Rivlin fluid model for blood flow through a tapered artery with a stenosis. The constitutive equations for a Reiner Rivlin fluid have been modelled in cylindrical coordinates. A perturbation series in dimensionless Reiner Rivlin fluid parameter (λ ₁ ≪ 1) have been used to obtain explicit forms for the velocity, temperature, concentration, resistance impedance, wall shear stress and shearing stress at the stenosis throat. The graphical results of different type of tapered arteries i.e. converging tapering, diverging tapering, non-tapered artery have been examined for different parameters of interest.     

Numerical simulation of fluid flow and heat transfer characteristics in channel with V corrugated plates

A detailed numerical study is carried out to investigate fluid flow and heat transfer characteristics in a channel with heated V corrugated upper and lower plates. The parameters studied include the Reynolds number (Re = 2,000–5,500), angles of V corrugated plates (θ = 20°, 40°, 60°), and constant heat fluxs (q″ = 580, 830, 1,090 W/m²). Numerical results have been validated using the experimented data reported by Naphon, and a good agreement has been found. The angles of V corrugated plates (θ) and the Reynolds number are demonstrated to significantly affect the fluid flow and the heat transfer rate. Increasing the angles of V corrugated plates can make the heat transfer performance become better. The increasing Reynolds number leads to a more complex fluid flow and heat transfer rate. The numerical calculations with a non-equilibrium wall function have a better accuracy than with a standard wall function for solving high Reynolds numbers or complex flow problems.     

Natural Convection in a Cavity Filled with Porous Layers on the Top and Bottom Walls

Natural convection in a partially filled porous square cavity is numerically investigated using SIMPLEC method. The Brinkman-Forchheimer extended model was used to govern the flow in the porous medium region. At the porous-fluid interface, the flow boundary condition imposed is a shear stress jump, which includes both the viscous and inertial effects, together with a continuity of normal stress. The thermal boundary condition is continuity of temperature and heat flux. The results are presented with flow configurations and isotherms, local and average Nusselt number along the cold wall for different Darcy numbers from 10⁻¹ to 10⁻⁶, porosity values from 0.2 to 0.8, Rayleigh numbers from 10³ to 10⁷, and the ratio of porous layer thickness to cavity height from 0 to 0.50. The flow pattern inside the cavity is affected with these parameters and hence the local and global heat transfer. A modified Darcy–Rayleigh number is proposed for the heat convection intensity in porous/fluid filled domains. When its value is less than unit, global heat transfer keeps unchanged. The interfacial stress jump coefficients β ₁ and β ₂ were varied from  −1 to +1, and their effects on the local and average Nusselt numbers, velocity and temperature profiles in the mid-width of the cavity are investigated.     

On the Shear Number Effect in Stratified Shear Flow

The influence of the shear number on the turbulence evolution in a stably stratified fluid is investigated using direct numerical simulations on grids with up to 512 × 256 × 256 points. The shear number SK/ε is the ratio of a turbulence time scale K/ε to the shear time scale 1/S. Simulations are performed at two initial values of the Reynolds number Re _Λ= 44.72 and Re _Λ= 89.44. When the shear number is increased from small to moderate values, the nondimensional growth rate γ= (1/SK)dK/dt of the turbulent kinetic energy K increases since the shear forcing and its associated turbulence production is larger. However, a further increase of the shear number from moderate to large values results in a reduction of the growth rate γ and the turbulent kinetic energy K shows long-time decay for sufficiently large values of the shear number. The inhibition of turbulence growth at large shear numbers occurs for both initial values of the Reynolds number and can be explained with the predominance of linear effects over nonlinear effects when the shear number is sufficiently high. It is found that, at the higher initial value of the Reynolds number, the reduction of the growth rate occurs at a higher value of the shear number. The shear number is found to affect spectral space dynamics. Turbulent transport coefficients decrease with increasing shear number. Received 23 June 1998 and accepted 25 February 1999     

A numerical study of the steady forced convection heat transfer from an unconfined circular cylinder

Forced convection heat transfer from an unconfined circular cylinder in the steady cross-flow regime has been studied using a finite volume method (FVM) implemented on a Cartesian grid system in the range as 10 ≤ Re ≤ 45 and 0.7 ≤ Pr ≤ 400. The numerical results are used to develop simple correlations for Nusselt number as a function of the pertinent dimensionless variables. In addition to average Nusselt number, the effects of Re, Pr and thermal boundary conditions on the temperature field near the cylinder and on the local Nusselt number distributions have also been presented to provide further physical insights into the nature of the flow. The rate of heat transfer increases with an increase in the Reynolds and/or Prandtl numbers. The uniform heat flux condition always shows higher value of heat transfer coefficient than the constant wall temperature at the surface of the cylinder for the same Reynolds and Prandtl numbers. The maximum difference between the two values is around 15–20%.     

Experimental investigation of natural convection heat transfer in horizontal and inclined annular fluid layers

In the present study, an experimental investigation of heat transfer and fluid flow characteristics of buoyancy-driven flow in horizontal and inclined annuli bounded by concentric tubes has been carried out. The annulus inner surface is maintained at high temperature by applying heat flux to the inner tube while the annulus outer surface is maintained at low temperature by circulating cooling water at high mass flow rate around the outer tube. The experiments were carried out at a wide range of Rayleigh number (5 × 10⁴ < Ra < 5 × 10⁵) for different annulus gap widths (L/D _o = 0.23, 0.3, and 0.37) and different inclination of the annulus (α = 0°, 30° and 60°). The results showed that: (1) increasing the annulus gap width strongly increases the heat transfer rate, (2) the heat transfer rate slightly decreases with increasing the inclination of the annulus from the horizontal, and (3) increasing Ra increases the heat transfer rate for any L/D _o and at any inclination. Correlations of the heat transfer enhancement due to buoyancy driven flow in an annulus has been developed in terms of Ra, L/D _o and α. The prediction of the correlation has been compared with the present and previous data and fair agreement was found.     

The role of the mathematical mean value theorem (MVT) in rheometry: an easy way to convert apparent flow curves into correct ones

The mean value theorem of integral calculus guarantees that the apparent viscosity η _a can easily be converted into the correct viscosity η. For ordinary liquids there is a direct identity between η _a and η but the apparent shear rate (or apparent shear stress) has to be shifted to the representative shear rate γ˙^ (or representative shear stress τ^). A model free approximation scheme is introduced which implies a constant shift factor. The corresponding approximation for η is acceptable for liquids most commonly encountered. For plastic fluids the relation between η and η _a is more complex since it involves a function depending upon α; the yield stress relative to the maximum stress within the viscometer. Using the same approximation scheme as before the shift factor will involve α as well. The corresponding approximation of η is shown to be acceptable for the whole range of α. Received: 7 February 2000/Accepted: 15 February 2000     

Unsteady compressible boundary layer flow over a circular cone near a plane of symmetry

An analysis has been performed to study the unsteady laminar compressible boundary layer governing the hypersonic flow over a circular cone at an angle of attack near a plane of symmetry with either inflow or outflow in the presence of suction. The flow is assumed to be steady at time t=0 and at t>0 it becomes unsteady due to the time-dependent free stream velocity which varies arbitrarily with time. The nonlinear coupled parabolic partial differential equations under boundary layer approximations have been solved by using an implicit finite-difference method. It is found that suction plays an important role in stabilising the fluid motion and in obtaining unique solution of the problem. The effect of the cross flow parameter is found to be more pronounced on the cross flow surface shear stress than on the streamwise surface shear stress and surface heat transfer. Beyond a certain value of the cross flow parameter overshoot in the cross flow velocity occurs and the magnitude of this overshoot increases with the cross flow parameter. The time variation of the streamwise surface shear stress is more significant than that of the cross flow surface shear stress and surface heat transfer. The suction and the total enthalpy at the wall exert strong influence on the streamwise and cross flow surface shear stresses and the surface heat transfer except that the effect of suction on the cross flow surface shear stress is small.     

Laminar forced convection in two-dimensional impacting tee junctions

This study presents numerical predictions of the laminar fluid-flow and heat-transfer characteristics in planar (two-dimensional) impacting tee junctions. The applicable Navier-Stokes equations and the energy equation were solved for airflow at two inlet Reynolds numbers (Re₁) and a wide range of the mass split ratio (β). The results include wall shear stress distributions, streamlines showing the number, location, and size of the re-circulation zones, the pressure loss coefficient, wall heat flux distributions, isotherms, and the overall rate of heat transfer. These results indicate that two re-circulation zones always form on the inside-bend wall of the tee at all values of β and Re₁. Two more re-circulation zones may form on the impacting wall of the tee depending on the values of β and Re₁. It was also found that the pressure loss coefficient reaches a minimum and the overall rate of heat transfer reaches a maximum at even mass split (β = 0.5).     

Heat transfer due to a round jet impinging normal to a circular cylinder

An experimental investigations of heat transfer for a stationary isothermal circular cylinder exposed normal to an impinging round air-jet has been reported. The circumferential heat transfer distributions as well as axial Nusselt number is measured. The measurements are taken as a function of the Reynolds number ranging from 3.8 × 10³ to 4 × 10⁴, the cylinder separation distance to the nozzle diameter (z/d) varying from 7 to 30, and the nozzle to cylinder diameter ratio (d/D) changing from 0.06 to 0.14. The output results indicated that the axial and radial distributions of the local heat transfer peaked at the impingement point. The heat transfer rate increases as the values of z decreases, for the same d and Re. The drop-off of the Nusselt number with increasing axial distance or radial angle from the impingement point was more pronounced for smaller z and d. The peripheral and surface average Nusselt numbers were determined by integration. The experimental data was used to produce correlations for both average and stagnation point heat transfer. Received on 4 January 1999     

Hydromagnetic flow between two porous disks rotating about non-coincident axes

Hydromagnetic flow between two porous disks rotating with same angular velocity Ω about two noncoincident axes has been studied in the presence of a uniform transverse magnetic field. An exact solution of the governing equations has been obtained in a closed form. It is found that the primary velocity f/Ωl increases and the secondary velocity g/Ωl decreases with increase in either Reynolds number Re or the Hartman number M. It is also found that the torque at the disk η= 0 increases with increase in either M^2 or K^2. On the other hand there is no torque at the disk η= 1 for large M^2 and K^2. The heat transfer characteristic has also been studied on taking viscous and Joule dissipation into account. It is seen that the temperature increases with increase in either M^2 or K^2. It is found that the rate of heat transfer at the disk η= 0 increases with increase in either M or K. On the other hand the rate of heat transfer at the disk η= 1 increases with increase in K but decreases with increase in M.     

Average and extremal properties of heat transfer and shear stress on a wall surface in Rayleigh–Bénard convection

In this study surface-averaged and extremal properties of heat transfer and shear stress on the upper wall surface of Rayleigh–Bénard convection are numerically examined. The Prandtl number was raised up to 10³, and the Rayleigh number was changed between 10⁴ and 10⁷. As a result, average Nusselt number Nu and shear rate τ/Pr depends on Pr, Ra, and the entire numerical results are distributed between two correlation equations corresponding to small and large Pr. The small and large Pr equations are closely related to steady and unsteady flow regimes, respectively. Nevertheless, a single relation τ/Pr ~ Nu ^3.0 exists to explain the entire results. Similarly the change of local maximal properties Nu _max and τ _max/Pr depends on Pr, Ra, and these values are also distributed between two correlation equations corresponding to small and large Pr cases. Despite such complicated dependence we can obtain a correlation equation as a form of τ _max/Pr ~ Nu _max^2.6, which has not been obtained theoretically.     

Coupled heat and mass transfer in mixed convection over a vertical flat plate embedded in saturated porous media: PST/PSC or PHF/PMF

A boundary layer analysis is used to investigate the heat and mass transfer characteristics of mixed convection about a vertical flat plate embedded in a saturated porous medium under the coupled effects of thermal and mass diffusion. The plate is maintained at prescribed surface temperature/concentration (PST/PSC) or prescribed heat/mass flux (PHF/PMF). The nonsimilar governing equations are obtained by using a suitable transformation and solved by Keller box method. Numerical results for the local heat transfer rate and the local mass transfer rate are presented for various parameters. The local heat and mass transfer rates increase with increasing n and m and buoyancy parameter ξ. When buoyancy parameter ξ is very small (large) the value of local Nusselt and the local Sherwood number correspond with the pure forced (free) convection, respectively. Increasing buoyancy ratio N (or N ^*) increases the local heat and mass transfer rates. It is apparent that Lewis number has a pronounced effect on the local mass transfer rate than it does on the local heat transfer rate. Furthermore, increasing Lewis number decreases (increases) the local heat (mass) transfer rate. Received on 8 December 1997     

Forced convection flow over a microsphere

Forced convection heat transfer characteristics around a microsphere subjected to uniform heat flux boundary condition is numerically investigated in this study. Moderate to high values of Reynolds number and a wide range of Prandtl number are considered. The analysis assumes that the continuity assumption is valid and hence the Navier–Stokes equations are solved for the range of Knudsen number of 0.001 ≤ Kn ≤ 0.1. The appropriate boundary conditions at the surface of the microsphere; the velocity slip and temperature jump are applied. The effect of the flow parameters: Re, Pr and Kn on the velocity and temperature distribution is presented and hence a better control on the boundary layer thickness can be achieved in the microscale level. Furthermore, the effect of the controlling parameters on the delay of flow separation, reduced shear stress, drag coefficient and on the Nusselt number profiles is also presented in the results.     

Flow and heat transfer in a moving fluid over a moving flat surface

In this paper, a numerical analysis of the momentum and heat transfer of an incompressible fluid past a parallel moving sheet based on composite reference velocity U is carried out. A single set of equations has been formulated for both momentum and thermal boundary layer problems containing the following parameters: r the ratio of the free stream velocity to the composite reference velocity, σ (Prandtl number) the ratio of the momentum diffusivity of the fluid to its thermal diffusivity, and E _c (E _ck ) (Eckert number). The present study has been carried out in the domain 0 ≤ r ≤ 1. It is found that the direction of the wall shear changes in such an interval and an increase of the parameter r yields an increase in temperature.      

Evaluation and comparison of routes to obtain pressure coefficients from high-pressure capillary rheometry data

A capillary rheometer equipped with a pressure chamber is used to measure the pressure-dependent viscosity of polymethylmethacrylate (PMMA), poly-α-methylstyrene-co-acrylonitrile (PαMSAN), and low-density polyethylene (LDPE). Data analysis schemes are discussed to obtain pressure coefficients at constant shear rate and at constant shear stress. It is shown that the constant shear stress pressure coefficients have the advantage of being shear stress independent for the three polymers. The constant shear rate pressure coefficients, on the other hand, turn out to depend on shear rate, which makes them less suitable for use, e.g., in process simulations. In addition to the commonly used superposition method, a direct calculation method for the pressure coefficients is tested. Values obtained from both methods are equivalent. However, the latter requires less experimental and calculational efforts. From the obtained pressure coefficients, it is clear that PMMA and PαMSAN have a very similar pressure dependence, while LDPE is less sensitive to pressure.