An experimental investigation of flow boiling in an asymmetrically heated rectangular microchannel

By using unique experimental techniques and carefully constructed experimental apparatus, the characteristics of flow boiling of water in microscale were investigated using a single horizontal rectangular microchannel. A polydimethylsiloxane rectangular microchannel (D_h = 103.5 and 133 μm) was fabricated by using the replica molding technique, a kind of soft lithography. A piecewise serpentine platinum microheater array on a Pyrex substrate was fabricated with the surface micromachining MEMS technique. Real time flow visualization of the phase change phenomena inside the microchannel was performed using a high speed CCD camera with microscope. The experimental local boiling heat transfer coefficients were studied, and single bubble inception, growth, and departure, as well as elongated bubble behavior were analyzed to elucidate the microscale heat transfer mechanisms. Tests were performed for mass fluxes of 77.5, 154.9, and 309.8 kg/m² s and heat fluxes of 180–500 kW/m². The effects of mass flux, heat flux, and vapor qualities on flow boiling heat transfer in a microchannel were studied.     

Rheological properties of suspensions of TiO2 particles in polymer solutions. 1. Shear viscosity

This paper presents results on the rheological behaviour of suspensions of two kinds of TiO₂ particles in two different polymer solutions. The particles differ in their hydrophilic or hydrophobic properties. The dispersing media are a solution of high molecular weight polyisobutylene in decalin and a solution of a low molecular weight polybutene in decalin. The concentrations of polymer are adjusted in order to get the same zero shear viscosity. The shear viscosity measurements display an apparent yield stress in some cases. The existence and the values of the yield stress depend on the volume fraction of solid particles and on the type of particles. The evolution of the intrinsic viscosity and of the maximum packing fraction vs the shear rate is interpreted in terms of evolution of the size and of the shape of aggregates of particles under shear. The effect of temperature on the development of the yield stress is also discussed. The results are completed by microscopic observations.     

Burnett simulations of gas flow in microchannels

The Burnett equations with slip boundary conditions are used to model the gas flow in microchannels in transition flow regime. As the Navier-Stokes equations are not appropriate to model the gas flow in this regime, the higher-order Burnett equations are adopted in the present study. In earlier studies, convergent solutions of the Burnett equations of microPoiseuille flow could only be obtained when Knudsen number is less than 0.2. By using a relaxation method on the boundary values, convergent solutions of the Burnett equations can be obtained even when Knudsen number reaches 0.4. The solutions of Burnett equations agree very well with experimental data and direct simulation Monte Carlo (DSMC) results. The pressure distributions and velocity profiles are then discussed in detail.     

Experimental studies of adiabatic flow boiling in fractal-like branching microchannels

Experimental results of adiabatic boiling of water flowing through a fractal-like branching microchannel network are presented and compared to numerical model simulations. The goal is to assess the ability of current pressure loss models applied to a bifurcating flow geometry. The fractal-like branching channel network is based on channel length and width ratios between adjacent branching levels of 2^−1/2. There are four branching sections for a total flow length of 18 mm, a channel height of 150 μm and a terminal channel width of 100 μm. The channels were Deep Reactive Ion Etched (DRIE) into a silicon disk. A Pyrex disk was anodically bonded to the silicon to form the channel top to allow visualization of the flow within the channels. The flow rates ranged from 100 to 225 g/min and the inlet subcooling levels varied from 0.5 to 6 °C. Pressure drop along the flow network and time averaged void fraction in each branching level were measured for each of the test conditions. The measured pressure drop ranged from 20 to 90 kPa, and the measured void fraction ranged from 0.3 to 0.9. The measured pressure drop results agree well with separated flow model predictions accounting for the varying flow geometry. The measured void fraction results followed the same trends as the model; however, the scatter in the experimental results is rather large.     

Lubricated optical rheometer for the study of two-dimensional complex flows of polymer melts

We describe a novel optical cross-slot channel rheometer generating two-dimensional and isothermal complex flows of polymer melts. This is made possible by lubricating the channel front and back viewing windows. Flow-induced birefringence and particle tracking velocimetry are reviewed and used to investigate the cross-slot flow of a low density polyethylene melt involving mixed shear and planar extensional deformations. This new device solves the issue of end effects in flow birefringence experiments where no variations of the optical properties along the light path are expected. It greatly facilitates the interpretation of stress field data by providing reliable measurements of the polymer melt extinction angle χ$\chi$ and retardation δ$\delta$, with a spatial resolution of one tenth of a millimeter. At the same time, it offers an enhanced temperature control and an increased optical accuracy due to an improved laser beam shaping. The capabilities and performances of this unique type of lubricated rheometer are discussed in detail and compared with previous approaches.     

A parameter investigation of shear-induced coalescence in semidilute PIB–PDMS polymer blends: effects of shear rate, shear stress volume fraction, and viscosity

In this work, drop coalescence of polymer blends under shear flow in a parallel flow apparatus was investigated by optical sectioning microscopy. In each experiment, shear rate was set at values low enough to avoid any break-up phenomena. The time evolution of the drop size distribution was determined by motorized sample scanning and iterative acquisition of stacks of images along sample depth. Drop size and location in the acquired images was found by automated image analysis techniques. A systematic experimental campaign to investigate the effects of shear rate (in the range 0.1–0.5 s⁻¹), volume fraction (2.5–10%), and viscosity of the two phases (3–63 Pa s) at different viscosity ratio (0.1–2.3) was carried out. By comparing data from different experiments, it was found that at any strain value, the average drop size decreases monotonically with the shear stress, calculated as the product of shear rate and matrix viscosity. Furthermore, the coalescence rate slowed down with increasing viscosity ratio. Overall, these results provide an extensive set of data, which can be used as a benchmark for modeling shear-induced coalescence in polymer blends.Paper presented at the Annual Meeting of the European Society of Rheology, Grenoble, April 2005.     

A procedure to determine the viscosity function from experimental data of capillary flow

The accurate calculation of the viscosity η as function of the shear rate &γdot; from capillary viscometry is still a matter of debate in the literature. In fact, this problem involves the inversion of an integral equation, which leads to multiple solutions due to the unavoidable noise present in the experimental data. The purpose of this work is to develop an efficient procedure to determine the viscosity function from experimental data of capillary flow without presenting the difficulties inherent in other methods discussed previously in the literature. The system identification procedure is used here to estimate the parameters of a viscosity model, which is appropriately selected for the fluid under study through preliminary calculations involving the apparent shear rate – shear stress data. Once the model is chosen by satisfying criteria for the fit goodness and its parameters are evaluated, a smooth and continuous function η(γdot;) is obtained in the range of experimental shear rates. The procedure proposed is also applicable to fluids in shear flow that present two Newtonian plateaus, as it is typically found in macromolecular dilute solutions. The mean value theorem of continuous functions is used to reduce significantly the computational time. Received: 15 November 1999 Accepted: 7 November 2000     

An experimental investigation of outflow of liquids from single-outlet vessels

This paper considers the outflow of a liquid from a single outlet vessel, i.e. a vessel in which the outflowing liquid is displaced by another fluid which enters the vessel through the same opening. The simplest possible arrangement is investigated: a sealed axisymmetric cylindrical vessel with an outlet in its base, in which water is displaced by air.

It is shown experimentally that the average liquid discharge velocity is independent of the liquid level in the vessel and the shape of the outlet for the range of outlets employed; it increases weakly with both the diameter of the vessel and the diameter of the outlet.     

Brownian dynamics simulations of shear-thickening in dilute polymer solutions

A versatile model describing the shear thickening behaviour of dilute polymer solutions in high shear flows is presented. The polymer macromolecules are modelled as Hookean elastic dumbbells which deform affinely during flow. In addition, the dumbbells feel a retractive anisotropic hydrodynamic drag and an isotropic Brownian force. Furthermore, it is assumed that high shear rate increases the probability of molecules forming associations and this is described through expressions for the frequencies of association and dissociation, without explicitly accounting for finite extensibility, hydrodynamic interaction or excluded volume effects. Thus, a reversible kinetic process is incorporated into the model, which results in two diffusion equations for the associated and dissociated dumbbells. Numerical simulations predict shear thickening for specific range of parameters, which are physically meaningful and related to molecular characteristics of the polymer. A comparison against experimental data reported in the literature revealed very promising results, thus confirming the ability of this model to predict shear thickening under a wide range of conditions, for various polymer models.Nomenclature A A factor in the frequency of association - B Frequency of dissociation - B ₀ Reference frequency of dissociation - c Concentration of polymer solution - c _i Concentration of singlets (i = 1) and doublets (i = 2) in the solution - c ^* The overlap concentration - D _t Translation coefficient of molecule - F _i (Q) Spring force for a singlet (i = 1) and for a doublet (i = 2) - F Frequency of association - F ₀ Reference frequency of association - H _i Dumbbell spring constant for a singlet (i = 1) and for a doublet (i = 2) - k Boltzman's constant - k _H Huggins constant - MW Molecular weight - MW _c Critical molecular weight for formation of entanglements - n Number density of molecules in the polymer solution - n ₀ Number density of dumbbells at equilibrium - n _i Number density of singlets (i = 1) and doublets (i = 2) - Q Vector defining the size and orientation of a dumbbell - t Time - T Absolute temperature - x Degree of multimerization - W Interaction energy between the two components of a doublet Greek letters a Dimensionless anisotropy parameter - Shear rate - _i Friction coefficient of singlets (i = 1) and doublets (i = 2) - _i Intrinsic viscosity of singlets (i = 1) and doublets (i = 2) - _red Reduced viscosity of solution - _sp Specific viscosity - Viscosity of the polymer solution of concentration c - _s Viscosity of the solvent - (t) White noise - K ^T Velocity gradient tensor - _Hi Time constant of a singlet (i = 1) and a doublet (i = 2) - ₁ Length scale of singlets (standard deviation of singlet lengths at equilibrium) - ₂ Length scale of doublets - T _p Stress tensor - T _xy Shear Stress (xy element of T _p ) - T _pi Contributions to the stress tensor of singlets (i = 1) and doublets (i = 2) - ₀ Equilibrium configuration distribution function of Q - _i Configuration distribution function of singlets (i = 1) and doublets (i = 2)     

铁电陶瓷PZT53复杂力电耦合行为的实验研究

通过实验研究了平行和垂直于极化方向的正应力对铁电陶瓷锆钛酸铅(PZT53)的电滞回线(E3-P3)和电致应变曲线(E-ε)的影响. 实验发现平行于极化方向的压应力对PZT53陶瓷的电滞回线、电致应变曲线形状以及矫顽场大小都有明显的影响,但是垂直于极化方向的拉、压应力只对PZT53陶瓷的电致应变曲线形状有明显的影响,但对电滞回线形状和矫顽场大小都没有显著影响. 采用畴翻转的模型详细解释了观察到的实验现象,所得结果为建立铁电陶瓷的多轴力、电耦合本构模型,提供了物理基础.     

Coarse-grained model of entangled polymer melts in non-equilibrium

A coarse-grained model developed for entangled polymeric systems and calibrated to represent melts in equilibrium (Rakshit, Picu, J Chem Phys 125:164907(1)–(10), 2006) is used to model shear flows. The model is a hybrid between multimode and mean-field representations: chain inner blobs are constrained to move along the chain backbone and the end blobs are free to move in 3D and continuously redefine the diffusion path for the inner blobs. Therefore, contour length fluctuations and reptation are captured. Constraint release is implemented by tracing the position of chain ends and performing a local relaxation of the chain backbones once end retraction is detected. This algorithm takes advantage of the multi-body nature of the model and requires no phenomenological parameters other than the length of an entanglement segment. The model is used to study start-up and step strain shear flows and reproduces features observed experimentally such as the overshoot during start-up shear flow, the Lodge–Meissner law, the monotonicity of the steady state shear stress with the strain rate, and shear thinning at large . These simulations are performed in conditions in which using a fully refined model of the same system would have been extremely computationally demanding or simply impossible with the current methods.     

The role of elongational viscosity in the mechanism of drag reduction by polymer additives

The role of elongational viscosity in the mechanism of drag reduction by polymer additives is investigated qualitatively by means of direct numerical simulations of a turbulent pipe flow. For the polymer solution, a generalised Newtonian constitutive model is utilised in which the viscosity depends on the second and third invariant of the rate-of-strain tensor via an elongation parameter. This elongation parameter is capable of identifying elongational type of regions within the flow. The simulations show that complementary to stretching of the polymers, also compression must be incorporated to have drag reduction, contrary to many suggestions done in the literature on the mechanism which assume that stretching of the polymers is most important.     

Direct comparison of equibiaxial elongational viscosity measurements from lubricated squeezing flow and the MultiAxiales Dehnrheometer

We present the first direct comparisons of rheological data from the lubricated squeezing flow (LSF) technique and the MultiAxiales Dehnrheometer (MAD) instrument developed by Meissner and coworkers (J Rheol 47:989–1010, 2003). Comparisons of transient equibiaxial elongational viscosity are carried out at strain rates well into the nonlinear regime on low-density polyethylene and polystyrene melts. We find data obtained using LSF deviate from the MAD data when the Hencky strain reaches a value of approximately 1, which we interpret as a failure of the LSF technique. The strain at which the LSF technique fails is relatively insensitive to experimental parameters including strain rate. For Hencky strains larger than 1, LSF data display behavior that could easily be mistaken for the phenomenon of strain hardening.     

An experimental investigation of pulsatile flow through a smooth constriction

Measurements of flow disturbances in the downstream region of modeled stenoses in a rigid tube, with upstream pulsatile flow are reported. Experiments were conducted over physiologically relevant mean Reynolds numbers of 600; based on the tube diameter and the time-averaged value of upstream centerline velocity. Contoured constrictions with 25%, 50% and 75% area reductions were investigated and velocity data were obtained from ensemble averaging techniques (phase-locked waveform). Experimental data over extensive spatial regions of poststenotic fields were taken, employing a two-component laser Doppler velocimeter LDV. Constant time sampling techniques for performing data or frequency analyses were used to avoid velocity bias and to study the evolution of poststenotic flow disturbances. It is found that different types of flow disturbances exist downstream of the constriction. Data analysis methods with the aid of flow visualization allow accurate classification of the disturbances which are sensitive indicators of mild to moderate constrictions. Although the present study was motivated by a biological situation, sufficient data were reported in detail that they may also be used by investigators working in computational fluid dynamics.     

Determination of the elongational viscosity of polymer melts by melt spinning experiments. A comparison with different experimental techniques

In this work, melt spinning experiments were tentatively used for the determination of the elongational viscosity of polymer melts at different levels of tensile strain and strain rate. The materials examined were two high-density polyethylene grades for blow moulding with similar number-average molecular mass but different polydispersity index. The data from melt spinning tests were compared with transient extensional viscosity data obtained by uniform isothermal tensile tests, performed by means of an extensional rheometer, as well as with those produced by converging flow tests (Cogswell model). The results showed that for high strain and strain rate levels, the melt spinning experiments provide elongational viscosity data quite close to the transient extensional viscosity values obtained from the tensile tests.     

An experimental investigation of the shear behaviour of polyethylenes with different structures

Rheological properties in shear flow are presented for four different polyethylene samples: a high density, a linear low density and two low density polyethylenes manufactured using different techniques. Tests have been performed with the aid of capillary types of instrument equipped with capillaries of various lengths at three different temperatures. End correction factors have been determined and true flow curves obtained. Swelling ratios for both unannealed and annealed samples have been determined as well as the shear rate and shear stress at which irregularities begin. In some cases generalized plots have been prepared and in all cases the rheological response is discussed in terms of molecular characteristics, in particular the average molecular weight, molecular weight distribution and degree of branching.     

An experimental investigation of heat transfer to pulsating pipe air flow with different amplitudes

Heat transfer characteristics to both laminar and turbulent pulsating pipe flows under different conditions of Reynolds number, pulsation frequency, pulsator location and tube diameter were experimentally investigated. The tube wall of uniform heat flux condition was considered for both cases. Reynolds number varied from 750 to 12,320 while the frequency of pulsation ranged from 1 to 10 Hz. With locating the pulsator upstream of the inlet of the test section tube, results showed an increase in heat transfer rate due to pulsation by as much as 30% with flow Reynolds number of 1,643 and pulsation frequency of 1 Hz, depending on the upstream location of the pulsator valve. Closer the valve to the tested section inlet, the better improvement in the heat transfer coefficient is achieved. Upon comparing the heat transfer results of the upstream and the downstream pulsation, at Reynolds number of 1,366 and 1,643, low values of the relative mean Nusselt number were obtained with the upstream pulsation. Comparing the heat transfer results of the two studied test sections tubes for Reynolds number range from 8,000 to 12,000 and pulsation frequency range from 1.0 to 10 Hz showed that more improvement in heat transfer rate was observed with a larger tube diameter. For Reynolds number ranging from 8,000 to 12,000 and pulsation frequency of 10 Hz, an improvement in the relative mean Nusselt number of about 50% was obtained at Reynolds number of 8,000 for the large test section diameter of 50 mm. While, for the small test section diameter of 15 mm, at same conditions of Reynolds number and frequency, a reduction in the relative mean Nusselt number of up to 10% was obtained.