Experimental study on buoyant flow stratification induced by a fire in a horizontal channel

Experiments were carried out in a reduced-scale horizontal channel to investigate the fire-induced buoyant flow stratification behavior, with the effect of the velocity shear between the hot buoyant flow and the cool air flow considered. This shear intensity was controlled and varied by changing the exhaust rate at the ceiling with one of the end of the channel opened. The flow pattern was visualized by the aid of a laser sheet. The horizontal traveling velocity, vertical temperature profile and stratification interface height of the buoyant flow were measured. The stratification pattern was found to fall into three regimes. Buoyancy force and inertia force, as the two factors that dominate the buoyant flow stratification, were correlated through the Froude number and the Richardson number. At Region I (Ri > 0.9 or Fr < 1.2), the buoyant flow stratification was stable, where a distinct interface existed between the upper smoke layer and the lower air layer. At Region II (0.3 < Ri < 0.9 or 1.2 < Fr < 2.4), the buoyant flow stratification was stable but with interfacial instability. At Range III (Ri < 0.3 or Fr > 2.4), the buoyant flow stratification becomes unstable, with a strong mixing between the buoyant flow and the air flow and then a thickened smoke layer.     

Time evolution of double-diffusive convection in a vertical cylinder with radial temperature and axial solutal gradients

The characteristics of transient double-diffusive convection in a vertical cylinder are numerically simulated using a finite element method. Initially the fluid in the cavity is at uniform temperature and solute concentration, then constant temperature and solute concentration, which are lower than their initial values, are imposed along the sidewall and bottom wall, respectively. The time evolution of the double-diffusive convection is investigated for specific parameters, which are the Prandtl number, Pr = 7, the Lewis number, Le = 5, the thermal Grashof number, Gr_T = 10⁷, and the aspect ratio, A = 2, of the enclosure. The objective of the work is to identify the effect of the buoyancy ratio (the ratio of solutal Grashof to thermal Grashof numbers: N = Gr_S/Gr_T) on the evolution of the flow field, temperature and solute field in the cavity. It is found that initially the fluid near the bottom wall is squeezed by the cold flow from the sidewall, a crest of the solute field forms and then pushed to the symmetry line. In the case of N > 0, a domain with higher temperature and weak flow (dead region) forms on the bottom wall near the symmetry line, and the area of dead region increases when N varies from 0.5 to 1.5. More crests of the solute field are formed and the flow near the bottom wall fluctuates continuously for N < 0. The frequency of the fluctuation increases when N varies from −0.5 to −1.5. Corresponding to the variety of the thermal and solutal boundary layers, the average rates of heat transfer (Nu) at the sidewall remain almost unchanged while the average rates of mass transfer (Sh) at the bottom wall change much in the cases of N = 1, 0, −1.     

Prediction of evaporation heat transfer coefficient based on gas–liquid two-phase annular flow regime in horizontal microfin tubes

A physical model of gas–liquid two-phase annular flow regime is presented for predicting the enhanced evaporation heat transfer characteristics in horizontal microfin tubes. The model is based on the equivalence of a periodical distortion of the disturbance wave in the substrate layer. Corresponding to the stratified flow model proposed previously by authors, the dimensionless quantity Fr₀ = G/[gd_eρ_v(ρ_l ? ρ_v)]^0.5 may be used as a measure for determining the applicability of the present theoretical model, which was used to restrict the transition boundary between the stratified-wavy flow and the annular/intermittent flows. Comparison of the prediction of the circumferential average heat transfer coefficient with available experimental data for four tubes and three refrigerants reveals that a good agreement is obtained or the trend is better than that of the previously developed stratified flow model for Fr₀ > 4.0 as long as the partial dry out of tube does not occur. Obviously, the developed annular model is applicable and reliable for evaporation in horizontal microfin tubes under the case of high heat flux and high mass flux.     

Numerical investigation for natural convection of a nanofluid in an inclined L-shaped cavity in the presence of an inclined magnetic field

The problem of natural convection in an inclined L-shaped enclosure filled with Cu/water nanofluid that operates under differentially heated walls in the presence of an inclined magnetic field is presented in this paper. The fully implicit finite difference method is used to solve the governing equations. A comparison with previously published results in special case of the present study is performed and a very good agreement is found. Heat transfer and fluid flow are examined for parameters of the Hartmann number (0 ≤ Ha ≤ 100), the nanoparticles volume fraction (0% ≤ ϕ ≤ 20%), the cavity inclination angle (0° ≤ ϑ ≤ 300°), the magnetic field inclination angle (0° ≤ γ ≤ 270°), the cavity aspect ratio (0.25 ≤ AR ≤ 0.6) and the Rayleigh number (10³ ≤ Ra ≤ 10⁶). It is found that, the presence of the magnetic field in the fluid region causes a significant reduction in the fluid flow and heat transfer characteristics. Also, a good enhancement in the heat transfer rate can be obtained by adding the copper nanoparticles to the base fluid.     

Numerical prediction of flow and heat transfer of power-law fluids in a plane channel with a built-in heated square cylinder

Structure of unsteady laminar flow and heat transfer of power-law fluids in two-dimensional horizontal plane channel with a built-in heated square cylinder is studied numerically. The governing equations are solved using a control volume finite element method (CVFEM) adapted to the staggered grid. Computations are performed over a range of Reynolds and Richardson numbers from Re = 20 to 200 and from Ri = 0 to 8, respectively at fixed Prandtl number Pr = 50 and blockage ratio value β′ = 1/8. Three different values of the power-law index (n = 0.5, 1 and 1.4) are considered in this study to show its effect on the value of the critical Reynolds number defining the transition between two different flow regimes (symmetrical and periodic flows), the variations of Strouhal number, drag and lift coefficients and the heat transfer from the square cylinder as function of Reynolds number. Heat transfer correlations are obtained through forced convection. A discussion about the buoyancy effect on the flow pattern and the heat transfer for different power-law index is also presented.     

Energy separation in the wake of a cylinder: Effect of Reynolds number and acoustic resonance

Energy separation is a spontaneous redistribution of total energy in a flowing fluid without external work or heat transfer. The energy separation mechanism in the vortex field behind an adiabatic circular cylinder in a cross flow of air is investigated. Time-averaged velocity and temperature measurements taken one diameter downstream of the cylinder (Re ～ 10⁵, M_∞ ～ 0.25) indicate flow reversal. The measured recovery temperature, expressed as distribution of energy separation factor indicates that energy separation is caused by the vortex flow in the wake, enhanced by acoustic excitation, and is insensitive to Reynolds number in the sub critical range studied.     

Onset of entrainment in transitional round fountains

It is of fundamental interest to understand the behavior of transitional fountains with intermediate Froude and Reynolds numbers, together with the associated entrainment and turbulence. In this work, the transient behavior of axisymmetric fountains with 1 ? Fr ? 8 and 200 ? Re ? 800 is studied by direct numerical simulation. It is found that at Re ? 200, there is little entrainment present at the upflow–downflow interface and at the downflow–ambient interface, even for a value of Fr as high as 8; however, at Re > 200, entrainment is present at these interfaces and the extent increases with Re, which clearly demonstrates that entrainment is strongly dependent on Re whereas the contribution from the Fr effect is relatively much smaller. The DNS results also show that z_m, which is the maximum fountain penetration height, fluctuates, even when the flow reaches full development, due to the entrainment at the upflow–downflow and the downflow–ambient interfaces, and the averaged z_m scales with ${\mathit{Fr}}^{\frac{3}{2}}{\mathit{Re}}^{\frac{1}{4}}$ for 1 < Fr ? 8 and 100 ? Re ? 800.     

Mixed convection flow and heat transfer across a square cylinder under the influence of aiding buoyancy at low Reynolds numbers

In this paper, mixed convection flow and heat transfer around a long cylinder of square cross-section under the influence of aiding buoyancy are investigated in the vertical unconfined configuration (Reynolds number, Re = 1–40 and Richardson number, Ri = 0–1). The semi-explicit finite volume method implemented on the collocated grid arrangement is used to solve the governing equations along with the appropriate boundary conditions. The onset of flow separation occurs between Re = 1–2, between Re = 2–3 and between Re = 3–4 for Ri = 0, 0.5 and 1, respectively. The flow is found to be steady for the range of conditions studied here. The friction, pressure and total drag coefficients are found to increase with Richardson number, i.e., as the influence of aiding buoyancy increases drag coefficients increase at the constant value of the Reynolds number. The temperature field around the obstacle is presented by isotherm contours at the Prandtl number of 0.7 (air). The local and average Nusselt numbers are calculated to give a detailed study of heat transfer over each surface of the square cylinder and an overall heat transfer rate and it is found that heat transfer increases with increase in Reynolds number and/or Richardson number. The simple expressions for the wake length and average cylinder Nusselt number are obtained for the range of conditions covered in this work.     

Heat transfer enhancement by flow bifurcations in asymmetric wavy wall channels

The enhancement characteristics of heat transfer, through a transition scenario of flow bifurcations, in asymmetric wavy wall channels, are investigated by direct numerical simulations of the mass, momentum and energy equations, using the spectral element method. The heat transfer characteristics, flow bifurcation and transition scenarios are determined by increasing the Reynolds numbers for three geometrical aspect ratios r = 0.25, 0.375, and 0.5, and Prandtl numbers 1.0 and 9.4. The transition scenarios to transitional flow regimes depend on the aspect ratio. For the aspect ratios r = 0.25 and 0.5, the transition scenario is characterized by one Hopf flow bifurcation. For the aspect ratio r = 0.375, the transition scenario is characterized by a first Hopf flow bifurcation from a laminar to a periodic flow, and a second Hopf flow bifurcation from a periodic to quasi-periodic flow. The periodic and quasi-periodic flows are characterized by fundamental frequencies ω₁, and ω₁ and ω₂, respectively. For all the aspect ratios and Prandtl numbers, the time-average mean Nusselt number and heat transfer enhancement increases with the Reynolds number as the flow evolves from a laminar to a transitional regime. For both Prandtl numbers, the highest increase in the Nusselt number occurs for the aspect ratio r = 0.5; whereas, the lowest increases happen to r = 0.25. The increase of the Nusselt number occurs at the expense of a higher pumping power, which, for both Prandtl numbers, grows as the aspect ratio increases from r = 0.25 to r = 0.5 for reaching a specific Nusselt number. This enhancement is obtained without the necessity of high volumetric flow rates associated with turbulent flow regimes, which demand much higher pumping powers. Significant heat transfer enhancements are obtained when the asymmetric wavy channel is operated in the appropriate transitional Reynolds number range.     

Frequency-dependent heat transfer enhancement from rectangular heated block array in a pulsating channel flow

The effect of pulsating flow on convective heat transfer from periodically spaced blocks in tandem on a channel wall is experimentally investigated. The spacing l between repeated blocks varied from l/L = 0.3 to 0.6 where L is the block pitch. The experiments are carried out in the range of 10 Hz < f_F < 100 Hz and 0.2 < A < 0.3. A pulsating flow is imposed by an acoustic woofer at the channel inlet and a constant heat is generated at each protruding block. The impact of the important governing parameters such as the Reynolds number, the Strouhal number and the inter-block spacing on heat transfer rate from heated blocks is examined in detail. The experimental results show that thermal transport from the blocks is greatly affected by the frequency, the amplitude of the flow pulsation, the inter-block spacing and the Reynolds number. Thermal resonance frequency which shows a maximum heat transfer coincides well with the inverse of traveling time of a fluid parcel that can be determined from the block periodicity and the Reynolds number.     

Effects of moving lid direction on MHD mixed convection in a linearly heated cavity

Effects of moving lid-direction on MHD mixed convection in a cavity with the bottom wall being linearly heated are analyzed using a numerical technique. Vertical walls of the enclosure are adiabatic and the sliding wall at the top has constant temperature. The lid moves in the negative and positive x-direction. Finite volume method has been used to solve the governing equations. Results are presented for different values of Hartmann number (0 ? Ha ? 30), Reynolds number (100 ? Re ? 1000) and Grashof number (10⁴ ? Gr ? 10⁶). It is found that direction of lid is more effective on heat transfer and fluid flow in the case of mixed convection than it is the case in forced convection. Heat transfer is also decreased with increasing of magnetic field for all studied parameters.     

Experiments on heat transfer enhancement from a heated square cylinder in a pulsating channel flow

Experiments have been performed to investigate heat transfer enhancement from a heated square cylinder in a channel by pulsating flow. For all the experiments, the amplitude of the pulsating flow is fixed at A = 0.05. The effects of the Reynolds number based on the mean flow velocity (Re = 350 and 540), the pulsating frequency (0 Hz < f_p < 60 Hz) and the blockage ratio of the square cylinder (β = 1/10, 1/8, and 1/6) on convective heat transfer are examined. The measured Strouhal numbers of shedding vortices for non-pulsating (A = 0) steady inlet flow are compared with the previously published data, and good agreement is found. The “lock-on” phenomenon is clearly observed for a square cylinder in the present flow pulsation. When the pulsating frequency is within the lock-on regime, heat transfer from the square cylinder is substantially enhanced. In addition, the influence of the Reynolds number and the blockage ratio on the lock-on occurrence is discussed in detail.     

Laminar free convection from a horizontal semi-circular cylinder to power-law fluids

Extensive numerical results on the flow and thermal fields are presented for free convection from a semi-circular cylinder (flat base upward) immersed in quiescent power-law fluids for the following ranges of conditions: Grashof number, 10 ? Gr ? 10⁵, Prandtl number, 0.72 ? Pr ? 100, and power-law index, 0.2 ? n ? 1.8. The heat transfer characteristics are analyzed in terms of the isotherm patterns, local and average Nusselt number as functions of the pertinent dimensionless parameters. The flow field is visualized in terms of the streamline patterns adjacent to the surface of the cylinder for a range of values of the Grashof number, Prandtl number and power-law index. A separated flow region forms at as low values of the Prandtl number as Pr = 0.72 for n ? 1 (Newtonian and shear-thickening fluids); whereas for shear-thinning fluids (n < 1), the flow remains attached to the cylinder surface over the range of conditions encompassed here. The bubble size grows with Grashof number and it shrinks with Prandtl number. In order to quantify the deviation from the Newtonian behaviour, the normalized values of average Nusselt number are analyzed as a function of the power-law index. In addition, a correlation is proposed for average Nusselt number as a function of the Grashof number, Prandtl number and power-law index. In general terms, shear-thinning fluid behaviour enhances heat transfer whereas shear-thickening has adverse influence on it.     

Experimental study of thermal behaviors in a rectangular channel with baffle of pores

This experimental study attempts to explore the local heat transfer in rectangular channel with baffles, and analyzes the experimental results of baffles with different heights and pores in the event of five Reynolds numbers and three heating quantities. Apart from increasing the perturbation of flow field, the channel's flow field with baffles, which is similar to a backward-facing step flow field, is very helpful to heat transfer. To obtain an optimized baffle and increase the perturbation of flow field, this experiment employed baffles with five heights (H = 10–50 mm) and different numbers of pores (N = 1–3), as well as heat flux: Q = 40–100 l/min, Reynolds number: 702–1752, and heating quantity: q_in = 90–750 W/m². In addition to measurement of overall temperature distribution, emphasis is also placed on analysis of local heat transfer coefficient. Furthermore, heat transfer distribution of channel can be applied to explain how the baffles of pores have an influence upon backward-facing step flow field, shear layer, recirculation region, reattachment region and redeveloped boundary layer. Finally, some empirical formulas derived form experimental results may provide a reference for future design.     

Analysis of combined thermal and magnetic convection ferrofluid flow in a cavity

A numerical analysis of the magnetic gradient and thermal buoyancy induced cavity ferrofluid flow is conducted by a semi-implicit finite element method. The physical model for a square cavity containing two different temperature side walls and a magnet near bottom wall is described by mass, momentum and energy equations. Conditions for the fixed Prandtl number, Rayleigh number and different ferro-hydrodynamic interaction parameter are studied for 5 × 10⁸ ≤ β ≤ 1.6 × 10¹⁰. Results show the flow strength increases with the strengthening magnetic field. However, the side-wall heat transfer rate presents a decrease trend with the increase in magnetic field strength, for the magnet located near the bottom central area evokes the circulation to move toward the central portion. In summary, a proper choice of magnet strength and location can adjust the flow field and local heat transfer rate to fit the practical application.     

Numerical study of swirling flows in a cylindrical container with co-/counter-rotating end disks under stable temperature difference

A numerical study was conducted to investigate swirling flows of a Boussinesq fluid confined in a cylindrical container with co-/counter-rotating end disks. A vertically stable temperature gradient is imposed, with the stationary sidewall assumed as adiabatic. Flows are studied for a range of parameters: the Reynolds number, Re, 100 ⩽ Re ⩽ 2000; the Richardson number, Ri, 0 ⩽ Ri ⩽ 1.0; and the Prandtl number, Pr, Pr = 1.0. The ratio of the angular velocity of the top disk to the bottom disk, s, −1.0 ⩽ s ⩽ 1.0. The cylinder aspect ratio: h = 2.0. For the case of negligibly small temperature difference (Ri ∼ 0) and high Re, interior fluid rotates with an intermediate angular velocity of both end disks when they are co-rotating (s > 0). When end disks are counter-rotating (s < 0), shearing flow with meridional recirculation is created. For the case of large temperature difference (Pr · Ri ∼ O(1)), the Ekman suction is suppressed and the sidewall boundary layer disappears at mid-height of the cylinder. For all the values of s considered in the present study, the bulk of the fluid is brought close to rest with the fluid in the vicinity of both end disks rotating in each direction. The secondary flow in the meridional section of the cylindrical container exhibits various types of vortices as the governing parameters are varied. These flow patterns are presented in the form of diagrams on the (s, Re) plane and (s, Ri) plane. The average Nusselt number is computed and presented as functions of Ri, Re and s.     

Linear stability analysis of Poiseuille–Bénard–Marangoni flow in a horizontal infinite liquid film

We present the first linear stability analysis of a Poiseuille–Bénard–Marangoni flow, which refers to a horizontal infinite liquid film flowing in one direction with uniform heating from below. This study concerns the two limiting cases of pure buoyancy effect (Ma = 0) and pure thermocapillary effect (Ra = 0). The stability thresholds of the flow and their variation with the control parameters (Biot, Reynolds and Prandtl numbers) are given and compared with those for a Poiseuille–Rayleigh–Bénard flow. The spatial structures of the flow are presented, and it is shown that the centers of the rolls are shifted upwards compared to the PRB case and that there is a loss of symmetry with respect to the vertical axis for the transverse rolls. These effects are directly linked to thermocapillary convection.     

Numerical simulation of natural convection heat transfer in a porous cavity heated from below using a non-Darcian and thermal non-equilibrium model

The present paper investigates the numerical simulation of steady laminar incompressible natural convection heat transfer in an enclosed cavity that is filled with a fluid-saturated porous medium. The bottom wall is subjected to a relatively higher temperature than the top wall while the vertical walls are considered to be insulated. The flow field is modeled upon incorporating different non-Darcian effects, such as the convective term, Brinkman effect and Forchhiemer quadratic inertial effect. Moreover the two-equation model is used to separately account for the local fluid and solid temperatures. The numerical solution is obtained through the application of the finite volume method. The appraisals of the sought objectives are performed upon identifying key dimensionless groups of parameters. These dimensionless groups along with their operating domains are: Rayleigh number 1 ⩽ Ra ⩽ 400, Darcy number 10⁻⁴ ⩽ Da ⩽ 10⁻³, effective fluid-to-solid thermal conductivity ratio 0.1 ⩽ κ ⩽ 1.0, and the modified Biot number 1 ⩽ χ ⩽ 100. The non-Darcian effects are first examined over a broad range of Rayleigh number. Next, the implications of the group of parameters on the flow circulation intensity, local thermal non-equilibrium (LTNE) and average Nusselt number are highlighted and pertinent observations are documented.     

Mixed convective transport in a lid-driven cavity containing a nanofluid and a rotating circular cylinder at the center

The mixed convective transport of Cu-H₂O nanofluid in a differentially heated and lid-driven square enclosure in the presence of a rotating circular cylinder is investigated numerically. The top wall of the enclosure is sliding from left to right at a uniform speed while all other walls are stationary. A thermally insulated circular cylinder is placed centrally within the enclosure. The cylinder can rotate about its centroidal axis. The top and bottom walls are kept isothermal at different temperatures while the side walls are assumed adiabatic. Simulations are performed for, Richardson number 1 ≤ Ri ≤ 10, dimensionless rotational speed 0 ≤ Ω ≤ 5 and nanoparticle concentration 0 ≤ ϕ ≤ 0.20 keeping the Grashof number fixed as Gr = 10⁴. The flow and thermal fields are analyzed through streamline and isotherm plots for various Ω and Ri. Furthermore, the drag coefficient of the moving lid and Nusselt number of the hot wall are also computed to understand the effects of Ω and Ri on them. It is observed that the heat transfer greatly depends on the rotational speed of the cylinder, mixed convective strength and the nanoparticle concentration.     

Numerical study of transient conjugate heat transfer of a turbulent impinging jet

This study presents the numerical study of transient conjugate heat transfer in a high turbulence air jet impinging over a flat circular disk. The numerical simulation of transient, two-dimensional cylindrical coordinate, turbulent flow and heat transfer is adopted to test the accuracy of the theoretical model. The turbulent governing equations are resolved by the control-volume based finite-difference method with a power-low scheme, and the well-known low-Re κ–ω turbulence model to describe the turbulent structure. The SIMPLE algorithm is adopted to solve the pressure–velocity coupling. The parameters studied include turbulent flow Reynolds number (Re = 16,100–29,600), heated temperature of a circular disk (T_h = 373 K) or heat flux (q″ = 63–189 kW/m²), and orifice to heat-source spacing (H/D = 4–10). The numerical results of the transient impinging process indicate that the jet Reynolds number has a significant effect on the hydrodynamics and heat transfer, particularly in the stagnation region of an impinging jet. High turbulence values lead to greater heat transfer coefficients in the stagnation region and cause a bypass of the laminar-to-turbulent transition region in the wall jet region. Induced turbulence from the environment around the jet also influences the variation of the stagnation heat transfer. The modeling approach used here effectively captures both the stagnation region behavior and the transition to turbulence, thus forming the basis of a reliable turbulence model.