Mixed convection flow investigation in a rectangular horizontal tube stenosis via liquid crystal thermography and planar PIV

The temperature and velocity field in a horizontal convergent-divergent rectangular channel heated from below is studied experimentally for a Reynolds range 8-120, Grashof numbers from 0.44 × 10⁵ to 2.56 × 10⁵ and Richardson numbers from 3 to 4000, using water as working fluid. The duct aspect ratio (width/height) varies from 1 at its inlet to 2.28 at the stenosis neck, and both the upper and bottom walls are tilted with an angle of ±15.7° with respect to the horizontal. The temperature of the bottom wall is kept constant above that of the issuing fluid. The temperature field is recorded by liquid crystals in the vertical mid plane whereas the velocity field is measured in this plane as well as in four cross sections of the divergent passage by planar PIV, revealing the characteristics of the secondary velocity field. For all the examined cases the flow in the convergent passage is free of thermal plumes, and the thermal boundary layer is thin. In contrast, the divergent passage is characterized by a thermal plume which is shifted upstream with increasing Gr or reducing Re. Both transversal and longitudinal rolls emerge in this diffuser the strength of which depend on Re and Gr influencing the streamwise distribution of Nusselt which for low Re presents a minimum.     

Natural convection in inclined two dimensional rectangular cavities

Steady two-dimensional natural convection in fluid filled cavities is numerically investigated. The channel is heated from below and cooled from the top with insulated side walls and the inclination angle is varied. The field equations for a Newtonian Boussinesq fluid are solved numerically for three cavity height based Rayleigh numbers, Ra = 10⁴, 10⁵ and 10⁶, and several aspect ratios. The calculations are in excellent agreement with previously published benchmark results. The effect of the inclination of the cavity to the horizontal with the angle varying from 0° to 180° and the effect of the startup conditions on the flow pattern, temperature distribution and the heat transfer rates have been investigated. Flow admits different configurations at different angles as the angle of inclination is increased depending on the initial conditions. Regardless of the initial conditions Nusselt number Nu exhibits discontinuities triggered by gradual transition from multiple cell to a single cell configuration. The critical angle of inclination at which the discontinuity occurs is strongly influenced by the assumed startup field. The hysteresis effect previously reported is not always present when the calculations are reversed from 90° to 0°. A comprehensive study of the flow structure, the Nu variation with varying angle of inclination, the effect of the initial conditions and the hysteresis effect are presented.     

Numerical study of simultaneous natural convection heat transfer from both surfaces of a uniformly heated thin plate with arbitrary inclination

 Numerical studies were conducted to investigate the natural convection heat transfer around a uniformly heated thin plate with arbitrary inclination in an infinite space. The numerical approach was based on the finite volume technique with a nonstaggered grid arrangement. For handling the pressure–velocity coupling the SIMPLE-algorithm was used. QUICK scheme and first order upwind scheme were employed for discretization of the momentum and energy convective terms respectively. Plate width and heating rate were used to vary the modified Rayleigh number over the range of 4.8×10⁶ to 1.87×10⁸. Local and average heat transfer characteristics were compared with regarding to the inclination angle. The empirical expressions for local and average Nusselt number were correlated. It has been found that for inclination angle less than 10^∘, the flow and heat transfer characteristics are complicated and the average Nusselt number can not be correlated by one equation while for inclination angle larger than 10^∘, the average Nusselt number can be correlated into an elegant correlation. Received on 18 April 2001 / Published online: 29 November 2001     

Numerical investigation of natural convection inside an inclined parallel‐walled channel

Steady two‐dimensional natural convection in an inclined parallel‐walled channel was investigated numerically. The full elliptic forms of conservation equations were solved together and the velocity vectors, temperature contours and local and average Nusselt number distribution were obtained. The comparisons of local and average Nusselt number with published experimental and numerical results indicate very good agreement. Results are presented for a single aspect ratio, L/b=24, over the range of Rayleigh number of 3–1000 and angle of inclination 0–90°. The results indicate that the overall channel average Nusselt number is reduced as the inclination angle is increased. Significant reductions in the overall Nusselt number are exhibited at high angle of channel inclination. Copyright © 2005 John Wiley & Sons, Ltd.     

Natural convection in inclined partitioned enclosures

The problem of steady, laminar, natural convective flow of a viscous fluid in an inclined enclosure with partitions is considered. Transverse gradient of temperature is applied on the two opposing regular walls of the inclined enclosure while the other walls are maintained adiabatic. The problem is formulated in terms of the vorticity-stream function procedure. A numerical solution based on the finite volume method is obtained. Representative results illustrating the effects of the enclosure inclination angle and the degree of irregularity on the contour maps of the streamlines and temperature are reported and discussed. In addition, results for the average Nusselt number at the heated wall of the enclosure and the difference of extreme stream-function values are presented and discussed for various Rayleigh numbers, inclination angles and dimensionless partition heights.     

Natural convection due to heat and mass transfer in a composite system

A numerical study is performed to analyse heat and mass transfer phenomena due to natural convection in a composite cavity containing a fluid layer overlying a porous layer saturated with the same fluid. The flow in the porous region is modelled using Brinkman–Forchheimer-extended Darcy model that includes both the effect of macroscopic shear (Brinkman effect) and flow inertia (Forchheimer effect). The vertical walls of the two-dimensional enclosure are isothermal whilst the horizontal walls are adiabatic. The two regions are coupled by equating the velocity and stress components at the interface. The resulting coupled equations in non-dimensional form are solved by an alternating direction implicit method by transforming them into parabolic form by the addition of false transient terms. The numerical results show that the amount of fluid penetration into the porous layer depends strongly upon the Darcy, thermal and solutal Rayleigh numbers. Average Nusselt number decreases while average Sherwood number increases with an increase of the Lewis number. The transfer of heat and mass on the heated wall near the interface depends strongly on the Darcy number. Received on 11 May 1998     

Fully Developed Mixed Convection Flow in a Vertical Channel Containing Porous and Fluid Layer with Isothermal or Isoflux Boundaries

An analysis of fully developed combined free and forced convective flow in a fluid saturated porous medium channel bounded by two vertical parallel plates is presented. The flow is modeled using Brinkman equation model. The viscous and Darcy dissipation terms are also included in the energy equation. Three types of thermal boundary conditions such as isothermal–isothermal, isoflux–isothermal, and isothermal–isoflux for the left–right walls of the channel are considered. Analytical solutions for the governing ordinary differential equations are obtained by perturbation series method. In addition, closed form expressions for the Nusselt number at both the left and right channel walls are derived. Results have been presented for a wide range of governing parameters such as porous parameter, ratio of Grashof number and Reynolds number, viscosity ratio, width ratio, and conductivity ratio on velocity, and temperature fields. It is found that the presence of porous matrix in one of the region reduces the velocity and temperature.     

Magnetohydrodynamic mixed convection in a vertical channel

The problem of combined free and forced convective magnetohydrodynamic flow in a vertical channel is analysed by taking into account the effect of viscous and ohmic dissipations. The channel walls are maintained at equal or at different constant temperatures. The velocity field and the temperature field are obtained analytically by perturbation series method and numerically by finite difference technique. The results are presented for various values of the Brinkman number and the ratio of Grashof number to the Reynolds number for both equal and different wall temperatures. Nusselt number at the walls is determined. It is found that the viscous dissipation enhances the flow reversal in the case of downward flow while it counters the flow in the case of upward flow. It is also found that the analytical and numerical solutions agree very well for small values of ε.     

Turbulent natural convection flow in a vertical channel with anti-symmetric heating

This paper presents the fluid flow characteristics of natural convection flow in an anti-symmetrically heated parallel-plate vertical channel using the Particle Image Velocimetry (PIV) system. The channel walls were subjected to uniform temperature conditions, with one wall heated above ambient and the opposing wall maintained below ambient temperature. Velocity measurements were conducted at three different sections on the horizontal plane to validate two-dimensionality of the flow and at three different heights in the vertical plane to establish vertical mean velocity profiles. The results indicate that the flow recirculation in the channel exhibits a similar character to that of a closed cavity and the induced flow rate is much lower than the one for a channel with both walls heated. A correlation for the dimensionless flow rate is developed.     

Experimental investigation of mixed convection in a rectangular duct with a heated plate in the middle of cross section

The flow and heat transfer in an inclined and horizontal rectangular duct with a heated plate longitudinally mounted in the middle of cross section was experimentally investigated. The heated plate and rectangular duct were both made of highly conductive materials, and the heated plate was subjected to a uniform heat flux. The heat transfer processes through the test section were under various operating conditions: Pr ≈ 0.7, inclination angle ϕ = −60° to +60°, Reynolds number Re = 334–1,911, Grashof number Gr = 5.26 × 10²–5.78 × 10⁶. The experimental results showed that the average Nusselt number in the entrance region was 1.6–2 times as large as that in the fully developed region. The average Nusselt numbers and pressure drops increased with the Reynolds number. The average Nusselt numbers and pressure drops decreased with an increase in the inclination angle from −60° to +60° when the Reynolds number was less than 1,500. But when the Reynolds number increased to over about 1,800, the heat transfer coefficients and pressure drops were independent of inclination angles.     

A two-phase model for laminar film condensation from steam-air mixtures in vertical parallel-plate channels

Detailed results are presented for laminar film condensation from steam-air mixtures flowing downward in vertical flat-plate channels. The mixture flow is laminar and saturation conditions prevail at the inlet. A fully coupled implicit numerical approach is used that achieves excellent convergence behavior, even for high inlet gas mass fractions. The detailed results include velocity, temperature, and gas mass fraction profiles, as well as axial variations of film thickness, pressure gradient and Nusselt number. The effects of a wide range of changes in the four independent variables (the inlet-to-wall temperature difference and the inlet values of gas concentration, Reynolds number, and pressure) on the film thickness, axial pressure gradient, and the local and average Nusselt numbers are carefully examined. It was found that increases in inlet concentration of noncondensable gas caused significant decreases in the film thickness, local Nusselt number, and axial pressure gradient. An analytical solution for the film thickness and velocity field at the end of condensation path was developed and shown to be the asymptotic value of the numerical results for large distances along the channel.     

Two-point correlation and integral scale of spatially advancing curved channel flow at friction-velocity-based Reynolds number 550

Two-point correlation and integral scale of curved channel turbulent flow were examined for onset and growth of large-scale vortex. Away from the two walls of the straight channel, the two-point correlation of streamwise velocity has large negative values around local minimum at y/δ = 1.6 and those of radial and spanwise components take small negative values at the same area. The large-scale motion in the straight channel is thus indicated to have conspicuous streamwise component and weak secondary motion. In the curved channel, two-point correlation of radial velocity decreased around the minimum away from the wall (y/δ = 1.6) and large negative values extended to the whole channel width. It is thus indicated that the large-scale vortex originated from the large-scale structures with conspicuous streamwise component in the outer-wall side, grew and expanded to the full channel width of the curved part. The integral scale integrating two-point correlation conditionally on the negative values greatly decreased to denote the onset and growth of the large-scale motion quantitatively. The local minima of two-point correlation for radial velocity were used as indicator of the large-scale vortex and near-wall small-scale vortex. The magnitude of it at the channel center area was as low as one thirds of those of near-wall area in the straight channel, and however increased and surpassed those of near-wall area in the curved channel. Therefore, large-scale vortices were weaker than the small-scale motion in the straight part of the channel, and they grew and surpassed the small-scale motion in strength in the curved part of the channel.     

Experimental research on heat transfer of natural convection in vertical rectangular channels with large aspect ratio

This work presents the experimental research on the steady laminar natural convection heat transfer of air in three vertical thin rectangular channels with different gap clearance. The much higher ratio of width to gap clearance (60–24) and the ratio of length to gap clearance (800–320) make the rectangular channels similar with the coolant flow passage in plate type fuel reactors. The vertical rectangular channels were composed of two stainless steal plates and were heated by electrical heating rods. The wall temperatures were detected with the K-type thermocouples which were inserted into the blind holes drilled in the steal plates. Also the air temperatures at the inlet and outlet of the channel were detected. The wall heat fluxes added to the air flow were calculated by the Fourier heat conduction law. The heat transfer characteristics were analyzed, and the average Nusselt numbers in all the three channels could be well correlated with the Rayleigh number or the modified Rayleigh number in a uniform correlation. Furthermore, the maximum wall temperatures were investigated, which is a key parameter for the fuel’s integrity during some accidents. It was found that even the wall heat flux was up to 1500 W/m², the maximum wall temperature was lower than 350 °C. All this work is valuable for the plate type reactor’s design and safety analysis.     

Nanoscale Poiseuille flow and effects of modified Lennard–Jones potential function

Numerical simulation of Poiseuille flow of liquid Argon in a nanochannel using the non-equilibrium molecular dynamics simulation (NEMD) is performed. The nanochannel is a three-dimensional rectangular prism geometry where the concerned numbers of Argon atoms are 2,700, 2,550 and 2,400 at 102, 108 and 120 K. Poiseuille flow is simulated by embedding the fluid particles in a uniform force field. An external driving force, ranging from 1 to 11 PN (Pico Newton), is applied along the flow direction to inlet fluid particles during the simulation. To obtain a more uniform temperature distribution across the channel, local thermostating near the wall are used. Also, the effect of other mixing rules (Lorenthz–Berthelot and Waldman–Kugler rules) on the interface structure are examined by comparing the density profiles near the liquid/solid interfaces for wall temperatures 108 and 133 K for an external force of 7 PN. Using Kong and Waldman–Kugler rules, the molecules near the solid walls were more randomly distributed compared to Lorenthz–Berthelot rule. These mean that the attraction between solid–fluid atoms was weakened by using Kong rule and Waldman–Kugler rule rather than the Lorenthz–Berthelot rule. Also, results show that the mean axial velocity has symmetrical distribution near the channel centerline and an increase in external driving force can increase maximum and average velocity values of fluid. Furthermore, the slip length and slip velocity are functions of the driving forces and they show an arising trend with an increase in inlet driving force and no slip boundary condition is satisfied at very low external force (<1 PN).     

Numerical investigation of static flow instability in a low-pressure subcooled boiling channel

A three-dimensional two-fluid model to predict subcooled boiling flow at low pressure is presented. The model is adopted to investigate the two-phase flow and heat transfer characteristics in a heated channel. The presence of bubbles as a consequence of heating flow through a vertical rectangular channel has a significant effect on the overall pressure drop along the channel. Numerical results were compared against a series experimental data performed at various conditions – mass flux, heat flux, inlet temperature and exit pressure. Good agreement on the overall pressure drop was achieved. The onset of flow instability velocity was also accurately determined when compared against measurements. Predicted results of void fraction provided useful information towards a more fundamental understanding of the occurrence of onset of nucleate boiling, onset of significant voiding and onset of flow instability. The phenomenon of boiling onset oscillations was also predicted through the use of the two-fluid model.     

Radiation effect on Darcy free convection flow along an inclined surface placed in porous media

The paper investigates the effect of radiation on Darcy's buoyancy induced flow of an optically dense viscous incompressible fluid along a heated inclined flat surface maintained at uniform temperature placed in a saturated porous medium with Rosseland diffusion approximation employing the implicit finite difference method together with Keller box elimination technique. Both the streamwise and normal components of the buoyancy force are retained in the momentum equations. The numerical results show that as the buoyancy parameter, ξ, increases the local Nusselt number increases. The results for the locally nonsimilar solutions are compared with the locally similar solutions for small angle of inclination and approximate similar solutions along vertical surface. The effect of the conduction-radiation parameter, R _d , and the surface temperature excess ration, θ _w , on the local Nusselt number, the tangential velocity distribution and the temperature distribution are also shown graphically.     

On the modeling of aiding mixed convection in vertical channels

This paper aims at showing that to prescribe a flow rate at the inlet section of a vertical channel with heated walls leads to surprising and counterintuitive physical solutions, especially when the problem is modeled as elliptical. Such an approach can give rise to the onset of recirculation cells in the entry region while the heat transfer is slightly increased under the influence of the buoyancy force. We suggest an alternative model based on more realistic boundary conditions based on a prescribed total pressure at the inlet and a fixed pressure at the outlet sections. In this case, the pressure and buoyancy forces act effectively in the same direction and, the concept of buoyancy aiding convection makes sense. The numerical results emphasize the large differences between solutions based on prescribed inlet velocity and those obtained with the present pressure-based boundary conditions.     

Influence of inclination angle on buoyancy-driven convection in triangular enclosure filled with a fluid-saturated porous medium

The present investigation deals with the numerical analysis of steady-state laminar buoyancy-driven convection in an inclined triangular enclosure filled with fluid saturated porous media using the Darcy law equation. One wall of the enclosure is isothermally heated and the other is cooled, while the remaining wall is adiabatic. The effect of inclination angle on natural convection is investigated by varying the angle of inclination (φ) between 0° and 360°. The governing transformed equations are solved numerically using a finite-difference method. Obtained results are shown in the form of streamlines, isotherms, mean Nusselt numbers and dimensionless stream function for different values of the Rayleigh number Ra in the range 100 ≤ Ra ≤ 1,000. It is found that the values of the maximum and minimum mean Nusselt number are reached for φ = 330° and φ = 210° , respectively. However, the lowest flow strength is formed at φ = 240° for all values of Ra.     

Indirect measurements of streamwise solid fraction variations of granular flows accelerating down a smooth rectangular chute

In this study, we detail a method for estimating the flux-averaged solid fraction of a steady granular flows moving down an inclined rectangular chute using velocity measurements from along the perimeter cross section, combined with knowledge of the mass flow rate through the cross section. The chute is 5 cm wide and 150 cm long with an adjustable inclination angle. Four inclination angles, from 27° to 36° at 3° intervals, are tested. This angle range overlaps the internal friction angle of the glass beads, which are 4 mm nominal in diameter. Two slender mirrors are installed at the top and the bottom of the transparent chute to reflect images of the flow down the chute of the two surfaces. This allows photographic recording of the flow with a PIV imaging system and measurement of the flow depth. The mass flow rate can be calibrated simultaneously by collecting the accumulated mass at the chute exit. A linear interpolation scheme is proposed to interpolate the volume flow rate in each section of the chute. Sensitivity analysis suggests that the relative standard deviation of this scheme is about ±6%, i.e., the resultant solid volume fraction is only moderately dependent on the interpolation scheme for the tested cases. This is further confirmed by a direct intercepting method. Compared to the sophisticated magnetic resonance imaging (MRI) or the radioactive positron emission particle tracking (PEPT) methods, the present method is verified as a cost-effective and nonhazardous alternative for ordinary laboratories. Two distinct groups of streamwise dependence of the solid fractions are found. They are separated by the inclination angle of the chute and agreed with the internal friction angle. In the experiments using the two smaller inclination angles, the solid fraction ratios are found to be linear functions of the streamwise distance, while for the two larger inclination angles, the ratios have a nonlinear concave shape. All decrease with growing downstream distance.     

Conjugate natural convection in a square porous cavity

 Steady-state conjugate natural convection in a square cavity filled with a porous medium is studied numerically in this paper. The enclosure consists of two horizontal conductive walls of finite thickness and two vertical walls at different uniform temperatures. The focus is on the role of solid-fluid conductivity ratio, k, on the flow and heat transfer characteristics and the average Nusselt number, , over the vertical hot and cold walls of the cavity for a limited set of particular parameters. It was shown that the interface temperature, θ_w, along the top of the solid wall decreases with the increase in the wall conductivity k. Also, the values of decreases with the increase of the values of the parameter k. Comparison with known results from the open literature when the wall thickness of the horizontal solid walls is neglected (non-conjugate problem) is excellent. Received on 4 April 2000