Evaluation of flow patterns and elongated bubble characteristics during the flow boiling of halocarbon refrigerants in a micro-scale channel

In the present study, quasi-diabatic two-phase flow pattern visualizations and measurements of elongated bubble velocity, frequency and length were performed. The tests were run for R134a and R245fa evaporating in a stainless steel tube with diameter of 2.32 mm, mass velocities ranging from 50 to 600 kg/m(2) s and saturation temperatures of 22 degrees C, 31 degrees C and 41 degrees C. The tube was heated by applying a direct DC current to its surface. Images from a high-speed video-camera (8000 frames/s) obtained through a transparent tube just downstream the heated sections were used to identify the following flow patterns: bubbly, elongated bubbles, churn and annular flows. The visualized flow patterns were compared against the predictions provided by Barnea et al. (1983) [1], Felcar et al. (2007) [10], Revellin and Thome (2007) [3] and Ong and Thome (2009) [11]. From this comparison, it was found that the methods proposed by Felcar et al. (2007) [10] and Ong and Thome (2009) [1] predicted relatively well the present database. Additionally, elongated bubble velocities, frequencies and lengths were determined based on the analysis of high-speed videos. Results suggested that the elongated bubble velocity depends on mass velocity, vapor quality and saturation temperature. The bubble velocity increases with increasing mass velocity and vapor quality and decreases with increasing saturation temperature. Additionally, bubble velocity was correlated as linear functions of the two-phase superficial velocity. (C) 2010 Elsevier Inc. All rights reserved.

Film thickness measurement techniques applied to micro-scale two-phase flow systems

Recently semi-empirical models to estimate flow boiling heat transfer coefficient, saturated CHF and pressure drop in micro-scale channels have been proposed. Most of the models were developed based on elongated bubbles and annular flows in the view of the fact that these flow patterns are predominant in smaller channels. In these models, the liquid film thickness plays an important role and such a fact emphasizes that the accurate measurement of the liquid film thickness is a key point to validate them. On the other hand, several techniques have been successfully applied to measure liquid film thicknesses during condensation and evaporation under macro-scale conditions. However, although this subject has been targeted by several leading laboratories around the world, it seems that there is no conclusive result describing a successful technique capable of measuring dynamic liquid film thickness during evaporation inside micro-scale round channels. This work presents a comprehensive literature review of the methods used to measure liquid film thickness in macro- and micro-scale systems. The methods are described and the main difficulties related to their use in micro-scale systems are identified. Based on this discussion, the most promising methods to measure dynamic liquid film thickness in micro-scale channels are identified. (C) 2009 Elsevier Inc. All rights reserved.

Hybrid modeling of convective laminar flow in a permeable tube associated with the cross-flow process

The confined flows in tubes with permeable surfaces arc associated to tangential filtration processes (microfiltration or ultrafiltration). The complexity of the phenomena do not allow for the development of exact analytical solutions, however, approximate solutions are of great interest for the calculation of the transmembrane outflow and estimate of the concentration, polarization phenomenon. In the present work, the generalized integral transform technique (GITT) was employed in solving the laminar and permanent flow in permeable tubes of Newtonian and incompressible fluid. The mathematical formulation employed the parabolic differential equation of chemical species conservation (convective-diffusive equation). The velocity profiles for the entrance region flow, which are found in the connective terms of the equation, were assessed by solutions obtained from literature. The velocity at the permeable wall was considered uniform, with the concentration at the tube wall regarded as variable with an axial position. A computational methodology using global error control was applied to determine the concentration in the wall and concentration boundary layer thickness. The results obtained for the local transmembrane flux and the concentration boundary layer thickness were compared against others in literature. (C) 2007 Elsevier B.V. All rights reserved.

Verification of a mixed high-order accurate DNS code for laminar turbulent transition by the method of manufactured solutions

This paper presents results on a verification test of a Direct Numerical Simulation code of mixed high-order of accuracy using the method of manufactured solutions (MMS). This test is based on the formulation of an analytical solution for the Navier-Stokes equations modified by the addition of a source term. The present numerical code was aimed at simulating the temporal evolution of instability waves in a plane Poiseuille flow. The governing equations were solved in a vorticity-velocity formulation for a two-dimensional incompressible flow. The code employed two different numerical schemes. One used mixed high-order compact and non-compact finite-differences from fourth-order to sixth-order of accuracy. The other scheme used spectral methods instead of finite-difference methods for the streamwise direction, which was periodic. In the present test, particular attention was paid to the boundary conditions of the physical problem of interest. Indeed, the verification procedure using MMS can be more demanding than the often used comparison with Linear Stability Theory. That is particularly because in the latter test no attention is paid to the nonlinear terms. For the present verification test, it was possible to manufacture an analytical solution that reproduced some aspects of an instability wave in a nonlinear stage. Although the results of the verification by MMS for this mixed-order numerical scheme had to be interpreted with care, the test was very useful as it gave confidence that the code was free of programming errors. Copyright (C) 2009 John Wiley & Sons, Ltd.

Slip ratio in dispersed viscous oil-water pipe flow

In this article, dispersed flow of viscous oil and water is investigated. The experimental work was performed in a 26.2-mm-i.d. 12-m-long horizontal glass pipe using water and oil (viscosity of 100 mPa s and density of 860 kg/m(3)) as test fluids. High-speed video recording and a new wire-mesh sensor based on capacitance (permittivity) measurements were used to characterize the flow. Furthermore, holdup data were obtained using quick-closing-valves technique (QCV). An interesting finding was the oil-water slip ratio greater than one for dispersed flow at high Reynolds number. Chordal phase fraction distribution diagrams and images of the holdup distribution over the pipe cross-section obtained via wire-mesh sensor indicated a significant amount of water near to the pipe wall for the three different dispersed flow patterns identified in this study: oil-in-water homogeneous dispersion (o/w H), oil-in-water non-homogeneous dispersion (o/w NH) and Dual continuous (Do/w & Dw/o). The phase slip might be explained by the existence of a water film surrounding the homogeneous mixture of oil-in-water in a hidrofilic-oilfobic pipe. (C) 2010 Elsevier Inc. All rights reserved.

Roughness and surface material effects on nucleate boiling heat transfer from cylindrical surfaces to refrigerants R-134a and R-123

This paper presents results of an experimental investigation carried out to determine the effects of the surface roughness of different materials on nucleate boiling heat transfer of refrigerants R-134a and R-123. Experiments have been performed over cylindrical surfaces of copper, brass and stainless steel. Surfaces have been treated by different methods in order to obtain an average roughness, Ra, varying from 0.03 mu m to 10.5 mu m. Boiling curves at different reduced pressures have been raised as part of the investigation. The obtained results have shown significant effects of the surface material, with brass being the best performing and stainless steel the worst. Polished surfaces seem to present slightly better performance than the sand paper roughened. Boiling on very rough surfaces presents a peculiar behavior characterized by good thermal performance at low heat fluxes, the performance deteriorating at high heat fluxes with respect to smoother surfaces. (C) 2008 Elsevier Inc. All rights reserved.

Modeling and experimental study of nucleate boiling on a vertical array of horizontal plain tubes

An investigation of nucleate boiling on a vertical array of horizontal plain tubes is presented in this paper. Experiments were performed with refrigerant RI 23 at reduced pressures varying from 0.022 to 0.64, tube pitch to diameter ratios of 1.32, 1.53 and 2.00, and heat fluxes from 0.5 to 40 kW/m(2). Brass tubes with external diameters of 19.05 mm and average roughness of 0.12 mu m were used in the experiments. The effect of the tube spacing on the local heat transfer coefficient along the tube array was negligible within the present range of experimental conditions. For partial nucleate boiling, characterized by low heat fluxes, and low reduced pressures, the tube positioning shows a remarkable effect on the heat transfer coefficient. Based on these data, a general correlation for the prediction of the nucleate boiling heat transfer coefficient on a vertical array of horizontal tubes under flooded conditions was proposed. According to this correlation, the ratio between the heat transfer coefficients of a given tube and the lowest tube in the array depends only on the tube row number, the reduced pressure and the heat flux. By using the proposed correlation, most of the experimental heat transfer coefficients obtained in the present study were predicted within +/- 15%. The new correlation compares reasonably well with independent data from the literature. (C) 2008 Elsevier Inc. All rights reserved.

Experimental study of a Tollmien-Schlichting wave interacting with a shallow 3D roughness element

The paper is devoted to an experimental study of the effect of a shallow 3D roughness element on the evolution of a 2D Tollmien-Schlichting wave in a Blasius boundary layer. The experiments were carried out under controlled disturbance conditions on an airfoil section which could provide a long run with zero pressure gradient flow. A pneumatically driven slit source was used to introduce the Tollmien-Schilichting wave upstream of the lower branch of the neutral stability curve. A few wavelengths downstream, the T-S wave interacts with a cylindrical roughness element. The height of the roughness was slowly oscillating in time, which allows a continuous measurement of the T-S wave response downstream the roughness. The oscillation frequency was approximately 1500 times lower than the frequency of the studied Tollmien-Schlichting wave and therefore, behaved as a steady roughness with respect to the T-S wave. Hot wire anemometry was used to measure wall normal profiles and spanwise scans close to the maximum of the eigenfunction of the T-S wave. The oscillation of the roughness and the synchronization of all-equipments permitted the use of ensemble average techniques. Two different amplitudes of T-S waves with a non-dimensional frequency of F120E-06 were studied. They show a strong amplification of the disturbances in a small spanwise wave number range. The analysis of the wall normal T-S profiles suggests the growth of oblique modes.

On the wake-induced vibration of tandem circular cylinders: the vortex interaction excitation mechanism

The mechanism of wake-induced vibrations (WIV) of a pair of cylinders in a tandem arrangement is investigated by experiments. A typical WIV response is characterized by a build-up of amplitude persisting to high reduced velocities; this is different from a typical vortex-induced vibration (VIV) response, which occurs in a limited resonance range. We suggest that WIV of the downstream cylinder is excited by the unsteady vortex-structure interactions between the body and the upstream wake. Coherent vortices interfering with the downstream cylinder induce fluctuations in the fluid force that are not synchronized with the motion. A favourable phase lag between the displacement and the fluid force guarantees that a positive energy transfer from the flow to the structure sustains the oscillations. If the unsteady vortices are removed from the wake of the upstream body then WIV will not be excited. An experiment performed in a steady shear flow turned out to be central to the understanding of the origin of the fluid forces acting on the downstream cylinder.

Numerical solution of the two-phase expansion of a metastable flashing liquid jet using the dispersion-controlled dissipative scheme

This paper presents a study of the stationary phenomenon of superheated or metastable liquid jets, flashing into a two-dimensional axisymmetric domain, while in the two-phase region. In general, the phenomenon starts off when a high-pressure, high-temperature liquid jet emerges from a small nozzle or orifice expanding into a low-pressure chamber, below its saturation pressure taken at the injection temperature. As the process evolves, crossing the saturation curve, one observes that the fluid remains in the liquid phase reaching a superheated condition. Then, the liquid undergoes an abrupt phase change by means of an oblique evaporation wave. Across this phase change the superheated liquid becomes a two-phase high-speed mixture in various directions, expanding to supersonic velocities. In order to reach the downstream pressure, the supersonic fluid continues to expand, crossing a complex bow shock wave. The balance equations that govern the phenomenon are mass conservation, momentum conservation, and energy conservation, plus an equation-of-state for the substance. A false-transient model is implemented using the shock capturing scheme: dispersion-controlled dissipative (DCD), which was used to calculate the flow conditions as the steady-state condition is reached. Numerical results with computational code DCD-2D vI have been analyzed. Copyright (C) 2009 John Wiley & Sons, Ltd.

Secondary instabilities in the flow around two circular cylinders in tandem

Direct stability analysis and numerical simulations have been employed to identify and characterize secondary instabilities in the wake of the flow around two identical circular cylinders in tandem arrangements. The centre-to-centre separation was varied from 1.2 to 10 cylinder diameters. Four distinct regimes were identified and salient cases chosen to represent the different scenarios observed, and for each configuration detailed results are presented and compared to those obtained for a flow around an isolated cylinder. It was observed that the early stages of the wake transition changes significantly if the separation is smaller than the drag inversion spacing. The onset of the three-dimensional instabilities were calculated and the unstable modes are fully described. In addition, we assessed the nonlinear character of the bifurcations and physical mechanisms are proposed to explain the instabilities. The dependence of the critical Reynolds number on the centre-to-centre separation is also discussed.

Phase-Locked Loop design applied to frequency-modulated atomic force microscope

The atomic force microscope (AFM) introduced the surface investigation with true atomic resolution. In the frequency modulation technique (FM-AFM) both the amplitude and the frequency of oscillation of the micro-cantilever must be kept constant even in the presence of tip-surface interaction forces. For that reason, the proper design of the Phase-Locked Loop (PLL) used in FM-AFM is vital to system performance. Here, the mathematical model of the FM-AFM control system is derived considering high order PLL In addition a method to design stable third-order Phase-Locked Loops is presented. (C) 2010 Elsevier B.V. All rights reserved.

Design constraints for third-order PLL nodes in master-slave clock distribution networks

Clock signal distribution in telecommunication commercial systems usually adopts a master-slave architecture, with a precise time basis generator as a master and phase-locked loops (PLLs) as slaves. In the majority of the networks, second-order PLLs are adopted due to their simplicity and stability. Nevertheless, in some applications better transient responses are necessary and, consequently, greater order PLLs need to be used, in spite of the possibility of bifurcations and chaotic attractors. Here a master-slave network with third-order PLLs is analyzed and conditions for the stability of the synchronous state are derived, providing design constraints for the node parameters, in order to guarantee stability and reachability of the synchronous state for the whole network. Numerical simulations are carried out in order to confirm the analytical results. (C) 2009 Elsevier B.V. All rights reserved.

Network of phase-locking oscillators and a possible model for neural synchronization

In order to model the synchronization of brain signals, a three-node fully-connected network is presented. The nodes are considered to be voltage control oscillator neurons (VCON) allowing to conjecture about how the whole process depends on synaptic gains, free-running frequencies and delays. The VCON, represented by phase-locked loops (PLL), are fully-connected and, as a consequence, an asymptotically stable synchronous state appears. Here, an expression for the synchronous state frequency is derived and the parameter dependence of its stability is discussed. Numerical simulations are performed providing conditions for the use of the derived formulae. Model differential equations are hard to be analytically treated, but some simplifying assumptions combined with simulations provide an alternative formulation for the long-term behavior of the fully-connected VCON network. Regarding this kind of network as models for brain frequency signal processing, with each PLL representing a neuron (VCON), conditions for their synchronization are proposed, considering the different bands of brain activity signals and relating them to synaptic gains, delays and free-running frequencies. For the delta waves, the synchronous state depends strongly on the delays. However, for alpha, beta and theta waves, the free-running individual frequencies determine the synchronous state. (C) 2011 Elsevier B.V. All rights reserved.

Optimized network structure for full-synchronization

A network of Kuramoto oscillators with different natural frequencies is optimized for enhanced synchronizability. All node inputs are normalized by the node connectivity and some important properties of the network Structure are determined in this case: (i) optimized networks present a strong anti-correlation between natural frequencies of adjacent nodes: (ii) this anti-correlation should be as high as possible since the average path length between nodes is maintained as small as in random networks: and (iii) high anti-correlation is obtained without any relation between nodes natural frequencies and the degree of connectivity. We also propose a network construction model with which it is shown that high anti-correlation and small average paths may be achieved by randomly rewiring a fraction of the links of a totally anti-correlated network, and that these networks present optimal synchronization properties. (C) 2008 Elsevier B.V. All rights reserved.