Low-GWP refrigerants flow boiling heat transfer in a 5 PPI copper foam

This paper reports an experimental investigation of the heat transfer performance of the new low-GWP refrigerants, R1234yf and R1234ze(E), during flow boiling heat transfer inside a horizontal high porosity copper foam with 5 Pores Per Inch (PPI). Metal foams are a class of cellular structured materials consisting of a stochastic distribution of interconnected pores; these materials have been proposed as effective solutions for heat transfer enhancement during both single and two-phase heat transfer. R1234yf and R1234ze(E) refrigerants are appealing alternatives of the more traditional R134a by virtue of their negligible values of GWP and normal boiling temperatures close to that of R134a, which make them suitable solution in several different applications, such as: refrigeration and air conditioning and electronic thermal management. This work compares the two-phase heat transfer behaviour of these new HFO refrigerants, studying the boiling process inside a porous medium and permitting to understand their effective heat transfer capabilities. The experimental measurements were carried out by imposing three different heat fluxes: 50, 75, and 100 kW m⁻², at a constant saturation temperature of 30 °C; the refrigerant mass velocity was varied between 50 and 200 kg m⁻² s⁻¹, whilst the mean vapour quality varied from 0.2 to 0.95. The two-phase heat transfer and pressure drop performance of the two new HFO refrigerants is compared against that of the more traditional R134a.     

3D Numerical study of the effect of aspect ratio on mixed convection air flow in upward solar air heater

A 3D Numerical study of mixed convection air flow in upward solar air heater with large spanwise aspect ratio (A = 10 to 40) was performed using CFD commercial code Fluent 14.5 (ANSYS). The main objective of this study is to investigate the channel height's effect (aspect ratio) on flow pattern and heat transfer in upward solar air heater in the particular case of low Re and high aspect ratio. The bottom plate (absorber) was submitted to Constant Heat Flux (CHF) in the range of 200 to 1000 W/m² and Reynolds number was varied from 50 to 1000. Our results are in concordance with most of authors conclusions about Poiseuille–Rayleigh–Benard flows. In mixed convection, increasing heat flux enhances heat transfer unlike forced convection flows. Simulation results of flow visualizations and Nusselt number calculations have shown that depending on Ri*, the velocity and temperature distributions in SAH vary greatly with the channel's height. The obtained results were different from previous studies. Indeed, our investigation of channel's height was achieved for the same heat flux but different Grashof numbers. For low channel's heights (high aspect ratio), increasing heat flux has not a significant effect but for higher channel's heights, an augmentation of heat flux enhances buoyancy effects in the flow and causes high turbulence. Also, increasing Reynolds number in low channel's heights (high A), can enhance substantially heat transfer. For higher channel's heights (low A), increasing Reynolds number decreases Ri* and thus buoyancy forces. Heat transfer is reduced and so Nusselt number. The obtained results may be very useful for engineers in designing and testing solar collectors.     

Numerical characterization of micro heat exchangers using experimentally tested porous aluminum layers

A microporous heat exchanger device is being developed for cooling high-power electronics. The device uses a mechanically compressed aluminum porous layer to improve the heat transfer at the coolant/solid interface and to provide more uniform cooling of the electronics. The hydraulic characteristics (porosity, permeability, and Forchheimer coefficient) of nine distinct compressed layers are obtained experimentally. These layers have porosity from 0.3 to 0.7 and permeability from 1.8 × 10⁻¹⁰ m2 to 1.2 × 10⁻⁹ m2. The inertia coefficient varies from 0.3 to 0.9. These hydraulic characteristics are used in the numerical simulations of a real microporous heat exchanger for cooling phased-array radars in development. Thermal and hydraulic performances are illustrated in terms of total pressure drop across the heat exchanger, maximum temperature difference in the direction transverse to the electronic modules, and maximum temperature within the coolant passage. Results indicate that the proposed design is capable of achieving a maximum transverse temperature difference of 2°C using polyalphaolephin as coolant.     

Turbulence characteristics over k-type rib roughened porous walls

To understand the permeability effects on turbulent rib-roughened porous channel flows, particle image velocimetry measurements are performed at the bulk Reynolds number of 5000–20000. Solid impermeable and porous ribs are considered for the rib-roughness whose geometry is categorised in the k-type roughness whose pitch/rib-height is 10. Three isotropic porous media with nearly the same porosity: 0.8, and different permeabilities (0.004, 0.020, 0.033 mm²) are applied. It is observed that the recirculation between the ribs becomes weak and the recirculation vortex submerges into the porous wall as the wall permeability and Reynolds number increase for both solid and porous rib cases while the recirculation vanishes in high permeable cases. These phenomena result in characteristic difference in turbulence quantities. By fitting the mean velocity profiles to the log-law form, the permeability effects of both rib and bottom wall on the log-law parameters and the equivalent sand-grain roughness are discussed. It is concluded that the zero-plane displacement increases while the von Kármán constant and the equivalent sand-grain roughness decrease as the wall and rib permeability increases.     

Experimental study of air natural convection on metallic foam-sintered plate

The natural convection on metallic foam-sintered plate at different inclination angles was experimentally studied. Seven copper foam samples with different pore densities (10–40 pore per inch), porosities (0.90–0.95), and aspect ratios (the ratio of foam thickness to sample length, 0.1–0.5) were measured at inclination angles of 0° (vertical orientation), 15°, 30°, 45°, 60°, 75°, 90° (horizontal orientation). The heat conduction and natural convection inside the foam both contributed to the total heat transfer. Although, the form and viscous drag, which are influenced by permeability and viscous friction in the thermal boundary layer respectively, tend to suppress the natural convection, the heat transfer was finally enhanced by the foam sintered surface due to large surface area extension. Optimum inclination range 60–75° corresponding to maximum average Nu number was found in the heat flux range of 600–1800 W/m². The sintered foam surface with lower porosity and pore density was recommended for heat transfer enhancement. Particularly, the sample with porosity 0.9, pore density of 10 PPI, aspect ratio of 0.5 offered the highest average Nu number among the studied samples. An empirical correlation for modified Nusselt number at isoflux boundary condition considering the foam morphology parameter and inclination angle was proposed within deviation ±15% between the correlation and the experimental data.     

Numerical study of mixed convection and conduction in a 2-D square ventilated cavity with an inlet at the vertical glazing wall and outlet at the top surface

A steady state numerical study of combined laminar mixed convection and conduction heat transfer in a ventilated square cavity is presented. The air inlet gap is located at the bottom of a vertical glazing wall and air exits the cavity via a gap located at the top surface. Three locations for the opening at the top surface: left (case a), center (case b) and right side (case c) are considered. All the remaining surfaces are considered adiabatic. The mass, momentum and energy conservation equations were solved using the finite volume method for different Rayleigh numbers in the interval of 10⁴ < Ra < 10⁶ and Reynolds number in the interval of 100 < Re < 700. Temperature, flow field, and heat transfer rates are analyzed. The effect of the interaction between ambient conditions outside the glazing and the air inlet gap at the bottom for different air outlet gap positions at the top surface modifies the flow structure and temperature distribution of the air inside the cavity. The Nusselt number as a function of the Reynolds number was determined for the three cases. It was found that configuration for case (a) removes a higher amount of heat entering the cavity compared to cases (b) and (c). This is due to the short distance between the main stream and the glass wall surface. Thus, the forced airflow entering the cavity is assisted by the buoyancy forces, and most of the cavity remains at the inlet flow temperature, which should be appropriate for warm climates. These results may provide useful information about the heat transfer and fluid flow for future studies.     

Effect of heater size on ultrasonic enhancement of boiling in water and surfactant solutions

Experiments were performed to study enhancement of heat transfer from the wire of d = 50 µm and the tube of d = 1.5 mm in subcooled pool boiling by ultrasonic waves. The working fluids are clean water and Alkyl (8-16) Glucoside surfactant solutions of different concentrations and bulk temperature 30 °C. The wire resistance was translated to the temperature, using the calibration data, the temperature of the tube was measured by thermocouple. The differences between effect of ultrasonic field on boiling in water for heaters of d = 50 µm and d = 1.5 mm may be summarized as follows: for boiling on the wire of d = 50 µm in subcooled water, T_b = 30 °C, enhancement of heat transfer coefficient due to applied ultrasonic field is about 70% and 20% at heat flux q = 620 kW/m² and q = 1350 kW/m², respectively. For boiling in surfactant solutions at the same boiling conditions enhancement of heat transfer coefficient is in the range of 5–10% at heat flux q = 620 kW/m² and 10–16% at heat flux q = 1350 kW/m² depending on solution concentration. For boiling on the tube of d= 1.5 mm in subcooled water, T_b= 30 ℃, enhancement of heat transfer coefficient due to applied ultrasonic field is about 50% and 45% at heat flux q = 500 kW/m² and q = 2500 kW/m², respectively. The same values are obtained for boiling in surfactant solution of concentration C = 250 ppm. For the wire of d = 50 µm the heat transfer enhancement due to acoustic vibrations in surfactant solutions is not as strong as in water. This fact may be considered as evidence of significant role of relationship between jet flow and ultrasonic field.     

Effects of inclination angle on natural convection in enclosures filled with Cu–water nanofluid

Effects of inclination angle on natural convection heat transfer and fluid flow in a two-dimensional enclosure filled with Cu-nanofluid has been analyzed numerically. The performance of nanofluids is tested inside an enclosure by taking into account the solid particle dispersion. The angle of inclination is used as a control parameter for flow and heat transfer. It was varied from  = 0° to  = 120°. The governing equations are solved with finite-volume technique for the range of Rayleigh numbers as 10³  Ra  10⁵. It is found that the effect of nanoparticles concentration on Nusselt number is more pronounced at low volume fraction than at high volume fraction. Inclination angle can be a control parameter for nanofluid filled enclosure. Percentage of heat transfer enhancement using nanoparticles decreases for higher Rayleigh numbers.     

Conjugate heat transfer in the unbounded flow of a viscoelastic fluid past a sphere

This work addresses the conjugate heat transfer of a simplified PTT fluid flowing past an unbounded sphere in the Stokes regime (Re = 0.01). The problem is numerically solved with the finite-volume method assuming axisymmetry, absence of natural convection and constant physical properties. The sphere generates heat at a constant and uniform rate, and the analysis is conducted for a range of Deborah (0 ≤ De ≤ 100), Prandtl (10⁰ ≤ Pr ≤ 10⁵) and Brinkman (0 ≤ Br ≤ 100) numbers, in the presence or absence of thermal contact resistance at the solid–fluid interface and for different conductivity ratios (0.1 ≤ κ ≤ 10). The drag coefficient shows a monotonic decrease with De, whereas the normalized stresses on the sphere surface and in the wake first increase and then decrease with De. A negative wake was observed for the two solvent viscosity ratios tested (β = 0.1 and 0.5), being more intense for the more elastic fluid. In the absence of viscous dissipation, the average Nusselt number starts to decrease with De after an initial increase. Heat transfer enhancement relative to an equivalent Newtonian fluid was observed for the whole range of conditions tested. The dimensionless temperature of the sphere decreases and becomes more homogeneous when its thermal conductivity increases in relation to the conductivity of the fluid, although small changes are observed in the Nusselt number. The thermal contact resistance at the interface increases the average temperature of the sphere, without affecting significantly the shape of the temperature profiles inside the sphere. When viscous dissipation is considered, significant changes are observed in the heat transfer process as Br increases. Overall, a simplified PTT fluid can moderately enhance heat transfer compared to a Newtonian fluid, but increasing De does not necessarily improve heat exchange.     

Experimental investigations of flow field and heat transfer characteristics due to periodically pulsating impinging air jets

The paper concentrates on increasing convective heat transfer due to periodically pulsating impinging air jets. A maximum enhancement rate of cooling effectiveness up to 20% could be detected at an excitation Strouhal number of Sr = 0.82 when using a high pulsation magnitude. Reductions up to 5% occured at low Strouhal numbers with coincident high pulsation magnitudes as well. The thermal results were completed with phase-locked flow field investigations by means of PIV and surface visualizations using the oil film method.     

Effect of force fields on pool boiling flow patterns in normal and reduced gravity

This paper reports the observations of boiling flow patterns in FC-72, performed during a microgravity experiment, recently flown aboard of Foton-M2 satellite, in some instances with the additional aid of an electrostatic field to replace the buoyancy force. The heater consisted of a flat plate, 20 × 20 mm², directly heated by direct current. Several levels of liquid subcooling (from 20 to 6 K) and heat fluxes up to 200 kW/m² were tested. A complete counterpart test, carried out on ground before the mission, allowed direct comparison with terrestrial data. The void fraction in microgravity revealed much larger than in normal gravity condition: this may be attributed to increased bubble coalescence that hinders vapor condensation in the bulk of the subcooled fluid. In several cases, an oscillatory boiling behavior was detected, leading to periodical variation of average wall overheating of some degrees. The electric field confirmed to be very effective, even at low values of applied voltage, in reducing bubble size, thus improving their condensation rate in the bulk fluid, and in enhancing heat transfer performance, suppressing the boiling oscillations and preventing surface dryout.     

Numerical analysis of turbulent convection heat transfer from an array of perforated fins

Numerical investigation is made for three-dimensional fluid flow and convective heat transfer from an array of solid and perforated fins that are mounted on a flat plate. Incompressible air as working fluid is modeled using Navier–Stokes equations and RNG based k ? ? turbulent model is used to predict turbulent flow parameters. Temperature field inside the fins is obtained by solving Fourier’s conduction equation. The conjugate differential equations for both solid and gas phase are solved simultaneously by finite volume procedure using SIMPLE algorithm. Perforations such as small channels with square cross section are arranged streamwise along the fin’s length and their numbers varied from 1 to 3. Flow and heat transfer characteristics are presented for Reynolds numbers from 2 × 10⁴ to 4 × 10⁴ based on the fin length and Prandtl number is taken Pr = 0.71. Numerical computations are validated with experimental studies of the previous investigators and good agreements were observed. Results show that fins with longitudinal pores, have remarkable heat transfer enhancement in addition to the considerable reduction in weight by comparison with solid fins.     

The influence of oil on nucleate pool boiling heat transfer

The influence of various oil contents in R134a is investigated for nucleate pool boiling on copper tubes either sandblasted or with enhanced heating surfaces (GEWA-B tube). Polyolester oils (POE) (Reniso Triton) with medium viscosity 55 cSt (SE55) and high viscosity 170 cSt (SE170) were used. Heat transfer coefficients were obtained for boiling temperatures between −28.6 and +20.1°C. The oil content varied from 0 to 5% mass fraction. For the sandblasted tube and the SE55 oil the heat transfer coefficients for the refrigerant/oil-mixture can be higher or lower than those for the pure refrigerant, depending on oil mass fraction, boiling temperature and heat flux. In some cases the highest heat transfer coefficients were obtained at a mass fraction of 3%. For the 170 cSt oil there is a clear decrease in heat transfer for all variations except for a heat flux 4,000 W/m² and −10.1°C at 0.5% oil content. The heat transfer coefficients are compared to those in the literature for a smooth stainless steel tube and a platinum wire. For the enhanced tube and 55 cSt oil the heat transfer coefficients are clearly below those for pure refrigerant in all cases. The experimental results for the sandblasted tube are compared with the correlation by Jensen and Jackman. The calculated values are within +20 and −40% for the medium viscosity oil and between +50% and −40% for the high viscosity oil. A correlation for predicting oil-degradation effects on enhanced surfaces does not exist.     

Augmentation of heat transfer in a bubble agitated vertical rectangular channel

This paper presents the results of an experimental study of convective heat transfer between three parallel vertical plates symmetrically spaced with and without bubble agitation to ascertain the degree of augmentation of the heat transfer coefficients due to agitation. The centre plate was electrically heated, while the other side plates were water-cooled forming two successive parallel vertical rectangular channels of dimensions 20 cm × 3.5 cm × 35 cm (length W, gap L, height H) each. At the bottom of the hot and cold plates air spargers were fitted. Water/ethylene glycol (100%) was used to fill the channels. The superficial gas velocity ranged from 0.0016 to 0.01 m/s. Top, bottom and sides of the channels were open to the water/ethylene glycol in the chamber which is the novel aspect of this study. Experimental data have been correlated as under: Natural convective heat transfer: Nu = 0.60 Gr ^0.29, r = 0.96, σ = 0.186, 1.17 E6 < Gr < 1.48 E7; Bubble agitated heat transfer: St = 0.11(ReFrPr ²)^−0.23, r = 0.82, σ = 0.002, 1.20 E−2 < (ReFrPr ²) < 1.36 E2.     

Comparison of experimental data with results of some drying models for regularly shaped products

This paper presents an experimental and theoretical investigation of drying of moist slab, cylinder and spherical products to study dimensionless moisture content distributions and their comparisons. Experimental study includes the measurement of the moisture content distributions of slab and cylindrical carrot, slab and cylindrical pumpkin and spherical blueberry during drying at various temperatures (e.g., 30, 40, 50 and 60°C) at specific constant velocity (U = 1 m/s) and the relative humidity φ = 30%. In theoretical analysis, two moisture transfer models are used to determine drying process parameters (e.g., drying coefficient and lag factor) and moisture transfer parameters (e.g., moisture diffusivity and moisture transfer coefficient), and to calculate the dimensionless moisture content distributions. The calculated results are then compared with the experimental moisture data. A considerably high agreement is obtained between the calculations and experimental measurements for the cases considered. The effective diffusivity values were evaluated between 0.741 × 10⁻⁵ and 5.981 × 10⁻⁵ m²/h for slab products, 0.818 × 10⁻⁵ and 6.287 × 10⁻⁵ m²/h for cylindrical products and 1.213 × 10⁻⁷ and 7.589 × 10⁻⁷ m²/h spherical products using the Model-I and 0.316 × 10⁻⁵–5.072 × 10⁻⁵ m²/h for slab products, 0.580 × 10⁻⁵–9.587 × 10⁻⁵ m²/h for cylindrical products and 1.408 × 10⁻⁷–13.913 × 10⁻⁷ m²/h spherical products using the Model-II.     

Heat transfer and fluid dynamics of air–water two-phase flow in micro-channels

Heat transfer, pressure drop, and void fraction were simultaneously measured for upward heated air–water non-boiling two-phase flow in 0.51 mm ID tube to investigate thermo–hydro dynamic characteristics of two-phase flow in micro-channels. At low liquid superficial velocity j_l frictional pressure drop agreed with Mishima–Hibiki’s correlation, whereas agreed with Chisholm–Laird’s correlation at relatively high j_l. Void fraction was lower than the homogeneous model and conventional empirical correlations. To interpret the decrease of void fraction with decrease of tube diameter, a relation among the void fraction, pressure gradient and tube diameter was derived. Heat transfer coefficient fairly agreed with the data for 1.03 and 2.01 mm ID tubes when j_l was relatively high. But it became lower than that for larger diameter tubes when j_l was low. Analogy between heat transfer and frictional pressure drop was proved to hold roughly for the two-phase flow in micro-channel. But satisfactory relation was not obtained under the condition of low liquid superficial velocity.     

Effect of porous disc receiver configurations on performance of solar parabolic trough concentrator

In this article, heat transfer enhancement of line focus solar collector with porous disc receiver is studied with water and therminol oil. A three dimensional (3-D) numerical simulation of porous disc enhanced receiver is carried out using commercial CFD software Fluent 6.3 to evolve the optimum configuration. The 3-D numerical model is solved by renormalization-group based k-ε turbulent model associated with standard wall function. The effect of porous disc receiver configurations (solid disc at bottom; porous disc at bottom; porous disc at top; and alternative porous disc) on performance of the trough concentrator is investigated. The effect of porous disc geometric parameters (φ, θ, W, H and t) and fluid parameters (Pr and m) on heat transfer enhancement of the receiver is also studied. The numerical simulation results show that the flow pattern around the solid and porous discs are entirely different and it significantly influences the local heat transfer coefficient. The porous disc receiver experiences low pressure drop as compared to that of solid disc receiver due to less obstruction. The optimum configuration of porous disc receiver enhances the heat transfer rate of 221 W m⁻¹ and 13.5% with pumping penalty of 0.014 W m⁻¹ for water and for therminol oil-55, heat transfer rate enhances of 575 W m⁻¹ and 31.4% with pumping penalty of 0.074 W m⁻¹ as compared to that of tubular receiver at the mass flow rate of 0.5 kg s⁻¹. The Nusselt number and friction factor correlations are proposed for porous disc receiver to calculate heat transfer characteristics. The porous disc receiver can be used to increase the performance of solar parabolic trough concentrator.     

Air-side performance of louver-finned flat aluminum heat exchangers at a low velocity region

The heat transfer and pressure drop characteristics of heat exchangers having louver fins were experimentally investigated. The samples had small fin pitches (1.0–1.4 mm), and experiments were conducted up to a very low frontal air velocity (as low as 0.3 m/s). Below a certain Reynolds number (critical Reynolds number), the fall-off of the heat transfer coefficient curve was observed. The critical Reynolds number was insensitive to the louver angle, and decreased as the louver pitch to fin pitch ratio (L _p /F _p ) decreased. Existing correlations on the critical Reynolds number did not adequately predict the data. The heat transfer coefficient curves crossed over as the Reynolds number decreased. Possible explanation is provided considering the louver pattern between neighboring rows. Different from the heat transfer coefficient, the friction factor did not show the fall-off characteristic. The reason was attributed to the form drag by louvers, which offsets the decreased skin friction at low Reynolds numbers. The friction factor increased as the fin pitch decreased and the louver angle increased. A new correlation predicted 92% of the heat transfer coefficient and 94% of the friction factor within ±10%.