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.

Stability analysis of core-annular flow and neutral stability wave number

A general transition criterion is proposed in order to locate the core-annular flow pattern in horizontal and vertical oil-water flows. It is based on a rigorous one-dimensional two-fluid model of liquid-liquid two-phase flow and considers the existence of critical interfacial wave numbers related to a non-negligible interfacial tension term to which the linear stability theory still applies. The viscous laminar-laminar flow problem is fully resolved and turbulence effects on the stability are analyzed through experimentally obtained shape factors. The proposed general transition criterion includes in its formulation the inviscid Kelvin-Helmholtz`s discriminator. If a theoretical maximum wavelength is considered as a necessary condition for stability, a stability criterion in terms of the Eotvos number is achieved. Effects of interfacial tension, viscosity ratio, density difference, and shape factors on the stability of core-annular flow are analyzed in detail. The more complete modeling allowed for the analysis of the neutral-stability wave number and the results strongly suggest that the interfacial tension term plays an indispensable role in the correct prediction of the stable region of core-annular flow pattern. The incorporation of a theoretical minimum wavelength into the transition model produced significantly better results. The criterion predictions were compared with recent data from the literature and the agreement is encouraging. (C) 2007 American Institute of Chemical Engineers.

CONCENTRATION FIELDS NEAR AIR-WATER INTERFACES DURING INTERFACIAL MASS TRANSPORT: OXYGEN TRANSPORT AND RANDOM SQUARE WAVE ANALYSIS

Mass transfer across a gas-liquid interface was studied theoretically and experimentally, using transfer of oxygen into water as the gas-liquid system. The experimental results support the conclusions of a theoretical description of the concentration field that uses random square waves approximations. The effect of diffusion over the concentration records was quantified. It is shown that the peak of the normalized rills concentration fluctuation profiles must be lower than 0.5, and that the position of the peak of the rms value is an adequate measure of the thickness of the diffusive layer. The position of the peak is the boundary between the regions more subject to molecular diffusion or to turbulent transport of dissolved mass.

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.

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.

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.

Computational Method to Determine Reflected Ultrasonic Signals From Arbitrary-Geometry Targets

A computational method based on the impulse response and on the discrete representation computational concept is proposed for the determination of the echo responses from arbitrary-geometry targets. It is supposed that each point of the transducer aperture can be considered as a source radiating hemispherical waves to the reflector. The local interaction with each of the hemispherical waves at the reflector surface can be modeled as a plane wave impinging on a planar surface, using the respective reflection coefficient. The method is valid for all field regions and can be performed for any excitation waveform radiated from an arbitrary acoustic aperture. The effects of target geometry, position, and material on both the amplitude and the shape of the echo response are studied. The model is compared with experimental results obtained using broadband transducers together with plane and cylindrical concave rectangular reflectors (aluminum, brass, and acrylic), as well as a circular cavity placed on a plane surface, in a water medium. The method can predict the measured echoes accurately. This paper shows an improved approach of the method, considering the reflection coefficient for all incident hemispherical waves arriving at each point of the target surface.

Viscosity measurement of Newtonian liquids using the complex reflection coefficient

This work presents the implementation of the ultrasonic shear reflectance method for viscosity measurement of Newtonian liquids using wave mode conversion from longitudinal to shear waves and vice versa. The method is based on the measurement of the complex reflection coefficient (magnitude and phase) at a solid-liquid interface. The implemented measurement cell is composed of an ultrasonic transducer, a water buffer, an aluminum prism, a PMMA buffer rod, and a sample chamber. Viscosity measurements were made in the range from 1 to 3.5 MHz for olive oil and for automotive oils (SAE 40, 90, and 250) at 15 and 22.5 degrees C, respectively. Moreover, olive oil and corn oil measurements were conducted in the range from 15 to 30 degrees C at 3.5 and 2.25 MHz, respectively. The ultrasonic measurements, in the case of the less viscous liquids, agree with the results provided by a rotational viscometer, showing Newtonian behavior. In the case of the more viscous liquids, a significant difference was obtained, showing a clear non-Newtonian behavior that cannot be described by the Kelvin-Voigt model.

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.

Animated solid model of the lung constructed from unsynchronized MR sequential images

This work discusses a 4D lung reconstruction method from unsynchronized MR sequential images. The lung, differently from the heart, does not have its own muscles, turning impossible to see its real movements. The visualization of the lung in motion is an actual topic of research in medicine. CT (Computerized Tomography) can obtain spatio-temporal images of the heart by synchronizing with electrocardiographic waves. The FOV of the heart is small when compared to the lung`s FOV. The lung`s movement is not periodic and is susceptible to variations in the degree of respiration. Compared to CT, MR (Magnetic Resonance) imaging involves longer acquisition times and it is not possible to obtain instantaneous 3D images of the lung. For each slice, only one temporal sequence of 2D images can be obtained. However, methods using MR are preferable because they do not involve radiation. In this paper, based on unsynchronized MR images of the lung an animated B-Repsolid model of the lung is created. The 3D animation represents the lung`s motion associated to one selected sequence of MR images. The proposed method can be divided in two parts. First, the lung`s silhouettes moving in time are extracted by detecting the presence of a respiratory pattern on 2D spatio-temporal MR images. This approach enables us to determine the lung`s silhouette for every frame, even on frames with obscure edges. The sequence of extracted lung`s silhouettes are unsynchronized sagittal and coronal silhouettes. Using our algorithm it is possible to reconstruct a 3D lung starting from a silhouette of any type (coronal or sagittal) selected from any instant in time. A wire-frame model of the lung is created by composing coronal and sagittal planar silhouettes representing cross-sections. The silhouette composition is severely underconstrained. Many wire-frame models can be created from the observed sequences of silhouettes in time. Finally, a B-Rep solid model is created using a meshing algorithm. Using the B-Rep solid model the volume in time for the right and left lungs were calculated. It was possible to recognize several characteristics of the 3D real right and left lungs in the shaded model. (C) 2007 Elsevier Ltd. All rights reserved.

Reproductive cycles in Bos indicus cattle

Several studies using transrectal ovarian ultrasonic scanning in Bos taurus (B. taurus) cattle and more recently in Bos indicus (B. Indicus) females evaluated the reproductive cycles of heifers and cows under different conditions. In general, B. indicus cattle have more follicles and more follicular waves during the estrous cycle and ovulate from smaller follicles than B. taurus. Consequently B. indicus females have smaller corpora lutea and it is assumed circulating concentrations of estradiol and progesterone are also less. However, these findings may vary depending on the nutritional status and regimen in which the animals are managed. Moreover, there are significant differences between B. taurus and B. indicus regarding follicle size at the time of deviation of the dominant follicle. These differences in ovarian function between B. indicus and B. taurus, e.g. greater antral follicle population are, probably, the main reasons for the great success of in vitro embryo production programs in Zebu cattle, especially in Brazil. (C) 2011 Elsevier B.V. All rights reserved.

Synthesis and characterization of the [Ru(3)O(CH(3)COO)(6)(py)(2)(BPE)Ru(bpy)(2)Cl](PF(6))(2) dimer

The polymetallic [Ru(3)O(CH(3)COO)(6)(py)(2)(BPE)Ru( bpy)(2)Cl](PF(6))(2) complex (bpy = 2,2`-bipyridine, BPE = trans- 1,2-bis(4-pyridil) ethylene and py = pyridine) was assembled by the combination of an electroactive [Ru(3)O] moiety with a [ Ru( bpy) 2( BPE) Cl] photoactive centre, and its structure was determined using positive ion electrospray (ESI-MS) and tandem mass (ESI-MS/MS) spectrometry. The [Ru(3)O(CH(3)COO)(6)(py)(2)(BPE)Ru(bpy)(2)Cl] (2+) doubly charged ion of m/z 732 was mass-selected and subject to 15 eV collision-induced dissociation, leading to a specific dissociation pattern, diagnostic of the complex structure. The electronic spectra display broad bands at 409, 491 and 692 nm ascribed to the [Ru(bpy)(2)(BPE)] charge-transfer bands and to the [Ru(3)O] internal cluster transitions. The cyclic voltammetry shows five reversible waves at - 1.07 V, 0.13 V, 1.17 V, 2.91 V and - 1.29 V (vs SHE) assigned to the [Ru(3)O](-1/0/+ 1/+ 2/+3) and to the bpy (0/-1) redox processes; also a wave is observed at 0.96 V, assigned to the Ru (+2/+ 3) pair. Despite the conjugated BPE bridge, the electrochemical and spectroelectrochemical results indicate only a weak coupling through the pi-system, and preliminary photophysical essays showed the compound decomposes under visible light irradiation.

The bimolecular sensitization of nitric oxide release from weak interacting ruthenium units

This work reports on the bimolecular sensitization of nitric oxide release from cis-[Ru(bpy)(2)(iso)-NO](PF(6))(3) (1) (iso = isoquinoline and bpy = 2,2`- bipyridine) by irradiating the MLCT transition of the chloro analog cis-[Ru(bpy) 2(iso) Cl] PF6 (2). The compounds displayed peaks in the ESI-MS spectra at m/z 749.1 and m/z 578.1 ascribed, respectively, to ([1(NO(o))-2PF(6)center dot CH(3)OH](2+)) and ([2-PF(6)](+)). In the cyclic voltammograms, the nitrosyl complex presented two redox waves related to the NO ligand at 0.48 and -0.37 V (versus Ag/AgCl, NO(+/0/-1) processes), while the sensitizer showed two reversible waves at 0.79 and -1.46 V (versus Ag/AgCl, Ru(2+/3+) and bpy(0/-1), respectively). The most important feature of this system is that the nitrosyl compound does not have significant absorption in the visible region, while the sensitizer has an intense band centered at 496 nm. The irradiation of an equimolar mixture of the two compounds in an ethanol: water solution (v: v) with light of lambda > 500 nm leads to NO release, as probed by amperometric measurements. The variational method was applied, showing that the two compounds self-assembly in solution with a 1: 1 stoichiometry. Fluorescence spectra acquired at 77 K provided the E(0-0) for the system and, from the thermodynamic cycle it was estimated that the photoinduced electron transfer between the species has a Delta G value of -1.59 eV. (C) 2011 Elsevier B. V. All rights reserved.

Excitability changes in human forearm corticospinal projections and spinal reflex pathways during rhythmic voluntary movement of the opposite limb

Rhythmic movements brought about by the contraction of muscles on one side of the body give rise to phase-locked changes in the excitability of the homologous motor pathways of the opposite limb. Such crossed facilitation should favour patterns of bimanual coordination in which homologous muscles are engaged simultaneously, and disrupt those in which the muscles are activated in an alternating fashion. In order to examine these issues, we obtained responses to transcranial magnetic stimulation (TMS), to stimulation of the cervicomedullary junction (cervicomedullary-evoked potentials, CMEPs), to peripheral nerve stimulation (H-reflexes and f-waves), and elicited stretch reflexes in the relaxed right flexor carpi radialis (FCR) muscle during rhythmic (2 Hz) flexion and extension movements of the opposite (left) wrist. The potentials evoked by TMS in right FCR were potentiated during the phases of movement in which the left FCR was most strongly engaged. In contrast, CMEPs were unaffected by the movements of the opposite limb. These results suggest that there was systematic variation of the excitability of the motor cortex ipsilateral to the moving limb. H-reflexes and stretch reflexes recorded in right FCR were modulated in phase with the activation of left FCR. As the f-waves did not vary in corresponding fashion, it appears that the phasic modulation of the H-reflex was mediated by presynaptic inhibition of Ia afferents. The observation that both H-reflexes and f-waves were depressed markedly during movements of the opposite indicates that there may also have been postsynaptic inhibition or disfacilitation of the largest motor units. Our findings indicate that the patterned modulation of excitability in motor pathways that occurs during rhythmic movements of the opposite limb is mediated primarily by interhemispheric interactions between cortical motor areas.

Modulational Instability in Periodic Quadratic Nonlinear Materials

We investigate the modulational instability of plane waves in quadratic nonlinear materials with linear and nonlinear quasi-phase-matching gratings. Exact Floquet calculations, confirmed by numerical simulations, show that the periodicity can drastically alter the gain spectrum but never completely removes the instability. The low-frequency part of the gain spectrum is accurately predicted by an averaged theory and disappears for certain gratings. The high-frequency part is related to the inherent gain of the homogeneous non-phase-matched material and is a consistent spectral feature.