Wavelet decomposition method decoupled boiling/evaporation oscillation mechanisms over two to three timescales: A study for a microchannel with pin fin structure

Boiling/evaporation heat transfer in a microchannel with pin fin structure was performed with water as the working fluid. Simultaneous measurements of various parameters were performed. The chip wall temperatures were measured by a high spatial-time resolution IR image system, having a sensitivity of 0.02 °C. The flow pattern variations synchronously changed wall temperatures due to ultra-small Bi number. The wavelet decomposition method successfully identified the noise signal and decoupled various temperature oscillations with different amplitudes and frequencies. Three types of temperature oscillations were identified according to heat flux q and mass flux G. The first type of oscillation occurred at q/G < 0.62 kJ/kg. The approximation coefficient of wavelet decomposition decided the dominant cycle period which was ∼3 times of the fluid residence time in the microchannel, behaving the density wave oscillation characteristic. The detail coefficients of wavelet decomposition decided the dominant cycle period, which matched the flow pattern transition determined value well, representing the flow pattern transition induced oscillation. For the second type of oscillation, the wavelet decomposition decoupled the three oscillation mechanisms. The pressure drop oscillation caused the temperature oscillation amplitudes of 5–10 °C and cycle periods of 10–15 s. The density wave oscillation and flow pattern transition induced oscillation are embedded with both the pressure rise and decrease stages of the pressure drop oscillation. The third type of oscillation happened at q/G > 1.13 kJ/kg, having the density wave oscillation coupled with the varied liquid film evaporation induced oscillation. The liquid island, retention bubble induced nucleation sites and cone-shape two-phase developing region are unique features of microchannel boiling with pin fin structure. This study illustrated that pressure drop oscillation and density wave oscillation, usually happened in large size channels, also take place in microchannels. The flow pattern transition and varied liquid film evaporation induced oscillations are specific to microchannel boiling/evaporation flow.     

“Magnetic-ribs” in fully developed laminar liquid–metal channel flow

This paper documents the numerical investigation of the effects of non-uniform magnetic fields, i.e. magnetic-ribs, on a liquid–metal flowing through a two-dimensional channel. The magnetic ribs are physically represented by electric currents flowing underneath the channel walls. The Lorentz forces generated by the magnetic ribs alter the flow field and, as consequence, the convective heat transfer and wall shear stress. The dimensionless numbers characterizing a liquid–metal flow through a magnetic field are the Reynolds (Re) and the Stuart (N) numbers. The latter provides the ratio of the Lorentz forces and the inertial forces. A liquid–metal flow in a laminar regime has been simulated in the absence of a magnetic field (Re_H = 1000, N = 0), and in two different magnetic ribs configurations for increasing values of the Stuart number (Re_H = 1000, N equal to 0.5, 2 and 5). The analysis of the resulting velocity, temperature and force fields has revealed the heat transport phenomena governing these magneto-hydro-dynamic flows. Moreover, it has been noticed that, by increasing the strength of the magnetic field, the convective heat transfer increases with local Nusselt numbers that are as much 27.0% larger if compared to those evaluated in the absence of the magnetic field. Such a convective heat transfer enhancement has been obtained at expenses of the pressure drop, which increases more than twice with respect to the non-magnetic case.     

Solvent-less synthesis of zinc oxide nanostructures from Zn(salen) as precursor and their optical properties

ZnO nanoparticles, 10–20 nm in size, were synthesized by heat treatment in air at 500 °C for 5 h., using [N,N′-bis(salicylaldehydo) ethylene diamine]zinc(II), i.e., Zn(salen), as precursor, which was obtained by a solvent-free solid–solid reaction. Heat-treated products were characterized by X-ray diffraction, Fourier transform infrared spectroscopy, scanning electron microscopy, and transmission electron microscopy. Room temperature photoluminescence spectra of ZnO nanostructures are dominated by green emission attributed to oxygen vacancy related donor–acceptor transition.     

The influence of pipe length on thermal statistics computed from DNS of turbulent heat transfer

We present results from direct numerical simulation of turbulent heat transfer in pipe flow at a bulk flow Reynolds number of 5000 and Prandtl numbers ranging from 0.025 to 2.0 in order to examine the effect of streamwise pipe length (πδ ≡ πD/2 ⩽ L ⩽ 12πδ) on the convergence of thermal turbulence statistics. Various lower and higher order thermal statistics such as mean temperature, rms of fluctuating temperature, turbulent heat fluxes, two-point auto and cross-correlations, skewness and flatness were computed and it is found that the value of L required for convergence of the statistics depends on the Prandtl number: larger Prandtl numbers requires comparatively shorter pipe length for convergence of most of the thermal statistics.     

The influence of pipe length on thermal statistics computed from DNS of turbulent heat transfer

We present results from direct numerical simulation of turbulent heat transfer in pipe flow at a bulk flow Reynolds number of 5000 and Prandtl numbers ranging from 0.025 to 2.0 in order to examine the effect of streamwise pipe length (πδ ≡ πD/2 ? L ? 12πδ) on the convergence of thermal turbulence statistics. Various lower and higher order thermal statistics such as mean temperature, rms of fluctuating temperature, turbulent heat fluxes, two-point auto and cross-correlations, skewness and flatness were computed and it is found that the value of L required for convergence of the statistics depends on the Prandtl number: larger Prandtl numbers requires comparatively shorter pipe length for convergence of most of the thermal statistics.     

Periodic boiling in parallel micro-channels at low vapor quality

Based on experimental investigations the present study evaluates instability and heat transfer phenomenon under condition of periodic flow boiling of water and ethanol in parallel triangular micro-channels. Tests were performed in the range of hydraulic diameter 100–220 μm, mass flux 32–200 kg/m² s, heat flux 120–270 kW/m², vapor quality x = 0.01–0.08. The period between successive events depends on the boiling number and decreases with an increase in the boiling number. The initial film thickness decreases with increasing heat flux. When the liquid film reached the minimum initial film thickness CHF regime occurred. Temporal variations of pressure drop, fluid and heater temperatures were periodic. Oscillation frequency is the same for the pressure drop, for the fluid temperature at the outlet manifold, and for the mean and maximum heater temperature fluctuations. All these fluctuations are in phase. The CHF phenomenon is different from that observed in a single channel of conventional size. A key difference between micro-channel heat sink and single conventional channel is amplification of parallel-channel instability prior to CHF. The dimensionless experimental values of the heat transfer coefficient are presented as the Nusselt number dependence on the Eotvos number and the boiling number.     

Developing structure of two-phase flow in a large diameter pipe at low liquid flow rate

In order to develop the interfacial area transport equation for the interfacial transfer terms in the two-fluid model, accurate data sets on axial development of local parameters such as void fraction, interfacial area concentration, interfacial gas velocity and Sauter mean diameter are indispensable to verify the modeled source and sink terms in the interfacial area transport equation. From this point of view, local measurements of both group 1 spherical/distorted bubbles and group 2 cap/slug bubbles in vertical upward air–water two-phase flow in a large diameter pipe with 200 mm in inner diameter and 26 m in height were performed at three axial locations of z/D = 41.5, 82.8 and 113 as well as 11 radial locations from r/R = 0–0.95 by using four-sensor probe method. Here, z, r, D and R are the axial distance from the inlet, radial distance from the pipe center, pipe diameter and pipe radius, respectively. The liquid flow rate and the void fraction ranged from 0.0505 m/s to 0.312 m/s and from 1.98% to 32.6%, respectively in the present experiment. The flow condition covered extensive region of bubbly flow, cap turbulent flow as well as their transition. The extensive analysis on the radial profiles of local flow parameters and their axial developments demonstrate the development of interfacial structures along the flow direction due to the bubble coalescence and breakup and the gas expansion. The significant decrease in void faction and interfacial area concentration and the increase in Sauter mean diameter and interfacial velocity were observed when the gradual flow regime transition occurred. Finally, the net change in the interfacial area concentration due to the bubble coalescence and breakup was quantitatively investigated in the present paper to reflect the true transfer mechanisms in observed two-phase flows.     

Effect of side ratio on fluid flow and heat transfer from rectangular cylinders using the PANS method

A computational study of heat transfer from rectangular cylinders is carried out. Rectangular cylinders are distinguished based on the ratio of the length of streamwise face to the height of the cross-stream face (side ratio, R). The simulations were performed to understand the heat transfer in a flow field comprising separation, reattachment, vortex shedding and stagnation. The Partially-Averaged Navier–Stokes (PANS) modeling approach is used to solve the turbulent flow physics associated and the wall resolve approach is used for the near wall treatment because of the flow separation involved. The simulations were performed using a finite volume based opensource software, OpenFOAM, at Reynolds number (Re) = 22,000 for rectangular cylinder at constant temperature kept in an air stream. Two critical side ratios were obtained, R = 0.62 and 3.0. At R = 0.62, the maximum value of the drag coefficient (C_d) = 2.681 was observed which gradually reduced by 54% at R = 4.0. The base pressure coefficient and global Nusselt number also attained the maximum value at R = 0.62 and from R = 2.5 to 3.0 a sharp discontinuous increase by 140% in the Strouhal number was observed. At R = 0.62, it was observed that the separated flow reattaches at the trailing edge after rolling over the side face and therefore increases the overall Nusselt number. The phase averaging was also performed to analyze the unsteady behavior of heat transfer.     

Study of material agglomeration during nozzle injection of multiviscosity liquid into a fluidized bed reactor by the electric conductance method

Fluidized bed agglomeration is an important and challenging problem for thermal cracking in fluid cokers. A low coker temperature can be problematic because the bitumen is injected into the fluidized bed with a different viscosity, resulting in formation of agglomerates of varying sizes, which slows the cracking reactions. In the present study, the bed material agglomeration process during nozzle injection of multiviscosity liquid was investigated in a fluidized bed operated at different mass ratios of the atomization gas to the liquid jets (GLR = 1%–3.5%) and gas velocities (3.9U_mf and 5.9U_mf) based on a conductance method using a water–sand system to simulate the hot bitumen–coke system at room temperature. During the tests of liquid-jet dispersion throughout the bed, different agglomeration stages are observed at both gas velocities. The critical amount of tert-butanol in the liquid jets that could lead to severe agglomeration of the bed materials (poor fluidization) at GLR = 1% is about 10 wt% at the low fluidizing gas velocity (3.9U_mf) and 18 wt% at the high gas velocity (5.9U_mf). This study provides a new approach for on-line monitoring of bed agglomeration during liquid injection to guarantee perfect contact between the atomized liquid and the bed particles.     

Effect of magnetic field on 3D flow and heat transfer during solidification from a melt

We present the effect of a magnetic field on three-dimensional fluid flow and heat transfer during solidification from a melt in a cubic enclosure. The walls of the enclosure are considered perfectly electrically conducting and the magnetic field is applied separately in three directions. The finite-volume method with enthalpy formulation is used to solve the mathematical model in the solid and liquid phases. The results obtained by our computer code are compared with the numerical and experimental data found in the literature. For Gr = 5 × 10⁵ and Ha = 0, 25, 50, 75, and 100 (where Gr and Ha are the Grashof and Hartmann numbers, respectively), the effects of magnetic field on flow and thermal fields, and on solid/liquid interface shape are presented and discussed. The interface is localized with and without magnetic field. The results show a strong dependence between the interface shape and the intensity and orientation of magnetic field. When the magnetic field is applied along the X-direction, the magnetic stability diagrams (V_max–Ha) and (Nu_avg–Ha) show the strongest stabilization of the flow field and heat transfer.     

An experimental study of non-isothermal miscible displacements in a Hele-Shaw cell

Non-isothermal miscible displacements in a radial Hele-Shaw cell were experimentally investigated using a scheme in which room temperature liquids of relatively high viscosity were displaced by high-temperature (80 °C), less-viscous liquids. Fundamental characteristics have been presented regarding how the effect of a non-isothermal field on miscible displacement patterns varies in terms of factors such as the viscosity ratio of the more- and less-viscous liquids at 20 °C, M₂₀, the rate of an increase in the pattern’s area, R, and the gap width of the cell, b. The concept of area density was used to quantitatively evaluate the effect of the non-isothermal fields on the patterns. We have found that the effect of the non-isothermal field on the patterns does not monotonically vary with M₂₀ and b. In contrast, it increases with R in the present experimental condition. The experimental results can be explained by introducing an assumption in which heat is transferred mainly to the plates of the cell, in other words, the temperature of the more-viscous liquid remains constant, whereas that of the less-viscous liquid spatiotemporally decreases and the viscosity of it increases along with the temperature decrease. Visualization of non-isothermal field in the cell has been done by means of a thermo sheet and the results support the assumption mentioned above.     

Transient heating effects in high pressure Diesel injector nozzles

The tendency of today’s fuel injection systems to reach injection pressures up to 3000 bar in order to meet forthcoming emission regulations may significantly increase liquid temperatures due to friction heating; this paper identifies numerically the importance of fuel pressurization, phase-change due to cavitation, wall heat transfer and needle valve motion on the fluid heating induced in high pressure Diesel fuel injectors. These parameters affect the nozzle discharge coefficient (C_d), fuel exit temperature, cavitation volume fraction and temperature distribution within the nozzle. Variable fuel properties, being a function of the local pressure and temperature are found necessary in order to simulate accurately the effects of depressurization and heating induced by friction forces. Comparison of CFD predictions against a 0-D thermodynamic model, indicates that although the mean exit temperature increase relative to the initial fuel temperature is proportional to (1 − C_d²) at fixed needle positions, it can significantly deviate from this value when the motion of the needle valve, controlling the opening and closing of the injection process, is taken into consideration. Increasing the inlet pressure from 2000 bar, which is the pressure utilized in today’s fuel systems to 3000 bar, results to significantly increased fluid temperatures above the boiling point of the Diesel fuel components and therefore regions of potential heterogeneous fuel boiling are identified.     

Study on premixed combustion in cylindrical micro combustors: Transient flame behavior and wall heat flux

The micro combustor is a key component of the micro thermophotovoltaic (TPV) system. Improving the wall temperature of the micro combustor is an effective way to elevate the system efficiency. An experimental study on the wall temperature and radiation heat flux of a series of cylindrical micro combustors (with a backward-facing step) was carried out. For the micro combustors with d = 2 mm, the regime of successful ignition (under the cold wall condition) was identified for different combustor lengths. Acoustic emission was detected for some cases and the emitted sound was recorded and analyzed. Under the steady-state condition, the effects of the combustor diameter (d), combustor length (L), flow velocity (u₀) and fuel–air equivalence ratio (Ф) on the wall temperature distribution were investigated by measuring the detailed wall temperature profiles. In the case that the micro combustor is working as an emitter, the optimum efficiency was found at Ф ≈ 0.8, independent of the combustor dimensions (d and L) and the flow velocity. Under the experimental conditions employed in the present study, the positions of the peak wall temperature were found to be about 8–11 mm and 4–6 mm from the step for the d = 3 mm and d = 2 mm micro combustors, respectively, which are 8–11 and 8–12 times of their respective step heights. This result suggests that the backward-facing step employed in the combustor design is effective in stabilizing the flame position.     

Flow structures and heat transfer on dimples in a staggered arrangement

Vortex structures and heat transfer enhancement mechanism of turbulent flow over a staggered array of dimples in a narrow channel have been investigated using Large Eddy Simulation (LES), Laser Doppler Velocimetry (LDV) and pressure measurements for Reynolds numbers Re_H = 6521 and Re_H = 13,042.The flow and temperature fields are calculated by LES using dynamic mixed model applied both for the velocity and temperature. Simulations have been validated with experimental data obtained for smooth and dimpled channels and empiric correlations. The flow structures determined by LES inside the dimple are chaotic and consist of small eddies with a broad range of scales where coherent structures are hardly to detect. Proper Orthogonal Decomposition (POD) method is applied on resolved LES fields of pressure and velocity to identify spatial–temporal structures hidden in the random fluctuations. For both Reynolds numbers it was found that the dimple package with a depth h to diameter D ratio of h/D = 0.26 provides the maximum thermo-hydraulic performance. The heat transfer rate could be enhanced up to 201% compared to a smooth channel.     

Flow boiling heat transfer of R134a,R236fa and R245fa in a horizontal 1.030 mm circular channel

This research focuses on acquiring accurate flow boiling heat transfer data and flow pattern visualization for three refrigerants, R134a, R236fa and R245fa in a 1.030 mm channel. We investigate trends in the data, and their possible mechanisms, for mass fluxes from 200 to 1600 kg/m²s, heat fluxes from 2.3 kW/m² to 250 kW/m² at T_sat = 31 °C and ΔT_sub from 2 to 9 K. The local saturated flow boiling heat transfer coefficients display a heat flux and a mass flux dependency but no residual subcooling influence. The changes in heat transfer trends correspond well with flow regime transitions. These were segregated into the isolated bubble (IB) regime, the coalescing bubble (CB) regime, and the annular (A) regime for the three fluids. The importance of nucleate boiling and forced convection in these small channels is still relatively unclear and requires further research.     

A numerical study of the longitudinal thermoconvective rolls in a mixed convection flow in a horizontal channel with a free surface

This paper presents a numerical study of three-dimensional laminar mixed convection within a liquid flowing on a horizontal channel heated uniformly from below. The upper surface is free and assumed to be flat. The coupled Navier–Stokes and energy equations are solved numerically by the finite volume method taking into account the thermocapillary effects (Marangoni effect). When the strength of the buoyancy, thermocapillary effects and forced convective currents are comparable (Ri ⋍ O(1) and Bd = Ra/Ma ⋍ O(1)), the results show that the development of instabilities in the form of steady longitudinal convective rolls is similar to those encountered in the Poiseuille–Rayleigh–Bénard flow. The number and spatial distribution of these rolls along the channel depend on the flow conditions. The objective of this work is to study the influence of parameters, such as the Reynolds, Rayleigh and Biot numbers, on the flow patterns and heat transfer characteristics. The effects of variations in the surface tension with temperature gradients (Marangoni effect) are also considered.     

CFD based mini- vs. micro-system delineation in elongated bubble flow regime

Delineation of mini- and micro-scale channels with respect to two-phase flow has been the subject of many research papers. There is no consensus on when the small channel can be characterized as a mini-channel or micro-channel. The idea proposed by this paper is to use the normalized bubble nose radius, liquid film thickness top over bottom ratio, and bubble shape contour, which are found under normal gravity conditions in slug flow through a horizontal adiabatic channel, as the delineation criteria. The input parameters are bubble nose radius and bubble nose velocity as the characteristic length scale and characteristic velocity scale respectively. 3D numerical simulation with ANSYS FLUENT was used to obtain the necessary data. Following CFD practice, a mesh independence study and a numerical model validation against published experimental data were both conducted. Analysis of the numerical simulation results showed that channels with D ⩽ 100 μm can be characterized as a micro-system, while channels with D ⩾ 400 μm belong to mini-systems. The region 200 μm ⩽ D ⩽ 300 μm represents a transition from the micro-scale to mini-scale.     

Saturated critical heat flux in a multi-microchannel heat sink fed by a split flow system

An extensive experimental campaign has been carried out for the measurement of saturated critical heat flux in a multi-microchannel copper heat sink. The heat sink was formed by 29 parallel channels that were 199 μm wide and 756 μm deep. In order to increase the critical heat flux and reduce the two-phase pressure drop, a split flow system was implemented with one central inlet at the middle of the channels and two outlets at either end. The base critical heat flux was measured using three HFC Refrigerants (R134a, R236fa and R245fa) for mass fluxes ranging from 250 to 1500 kg/m² s, inlet subcoolings from ?25 to ?5 K and saturation temperatures from 20 to 50 °C. The parametric effects of mass velocity, saturation temperature and inlet subcooling were investigated. The analysis showed that significantly higher CHF was obtainable with the split flow system (one inlet–two outlets) compared to the single inlet–single outlet system, providing also a much lower pressure drop. Notably several existing predictive methods matched the experimental data quite well and quantitatively predicted the benefit of higher CHF of the split flow.     

Influences of peripherally-cut twisted tape insert on heat transfer and thermal performance characteristics in laminar and turbulent tube flows

Effects of peripherally-cut twisted tape insert on heat transfer, friction loss and thermal performance factor characteristics in a round tube were investigated. Nine different peripherally-cut twisted tapes with constant twist ratio (y/W = 3.0) and different three tape depth ratios (DR = d/W = 0.11, 0.22 and 0.33), each with three different tape width ratios (WR = w/W = 0.11, 0.22 and 0.33) were tested. Besides, one typical twisted tape was also tested for comparison. The measurement of heat transfer rate was conducted under uniform heat flux condition while that of friction factor was performed under isothermal condition. Tests were performed with Reynolds number in a range from 1000 to 20,000, using water as a working fluid. The experimental results revealed that both heat transfer rate and friction factor in the tube equipped with the peripherally-cut twisted tapes were significantly higher than those in the tube fitted with the typical twisted tape and plain tube, especially in the laminar flow regime. The higher turbulence intensity of fluid in the vicinity of the tube wall generated by the peripherally-cut twisted tape compared to that induced by the typical twisted tape is referred as the main reason for achieved results. The obtained results also demonstrated that as the depth ratio increased and width ratio decreased, the heat transfer enhancement increased. Over the range investigated, the peripherally-cut twisted tape enhanced heat transfer rates in term of Nusselt numbers up to 2.6 times (turbulent regime) and 12.8 times (laminar regime) of that in the plain tube. These corresponded to the maximum performance factors of 1.29 (turbulent regime) and 4.88 (laminar regime).     

Conjugate turbulent heat transfer in a square cavity with a solar control coating deposited to a vertical semitransparent wall

This paper presents a numerical study of the conjugate heat transfer (natural convection, surface thermal radiation and conduction) in a square cavity with turbulent flow. The cavity has one vertical isothermal wall, two horizontal adiabatic walls and one vertical semitransparent wall with a selective coating applied to the inner side to control the solar radiation transmission. Later on the semitransparent wall is replaced with another one without the selective coating. The mathematical model for the turbulent flow in the cavity was solved using the finite volume method. The system had the following conditions: the uniform temperature in the isothermal wall was 21 °C, the external ambient temperature was fixed at 35 °C and on the semitransparent wall the direct normal solar irradiation of 750 W/m² was considered constant. The Rayleigh number was varied in the range of 10⁹ ? Ra ? 10¹² by changing the lengths of the cavity from 0.70 m to 6.98 m, respectively. The results show that, even though the air temperature of the cavity with the solar control film coating semitransparent wall (case A) is higher compared with the one without solar film coating (case B), the total amount of heat going through the cavity is lower compared to the one going through the cavity without solar control film. The total amount of energy transferred to the air in cavity for the case A was 41.98% less than for the case B. A set of correlations for the Nusselt number was obtained for both cases considering the conjugate heat transfer.