Double pulsed particle image velocimeter with directional resolution for complex flows

An automated particle image velocimeter has been constructed to determine fluid velocity fields by directly measuring particle image separations. A new image shifting technique eliminates directional ambiguity from velocity measurements. Theoretical and practical considerations in implementing the PIV system are presented. Using a simple, known flow field to investigate system performance, RMS errors of less than 0.7% of full scale are achieved, and at low seeding density fewer than 10% of velocity measurements produce bad data. Applications of system measurements to the instantaneous velocity field of an axisymmetric jet are presented.A version of this paper was presented at the 10th Symposium on Turbulence, University of Missouri-Rolla, September 22–24, 1986.     

Confined viscoplastic flows with heterogeneous wall slip

The steady, pressure-driven flow of a Herschel-Bulkley fluid in a microchannel is considered, assuming that different power-law slip equations apply at the two walls due to slip heterogeneities, allowing the velocity profile to be asymmetric. Three different flow regimes are observed as the pressure gradient is increased. Below a first critical pressure gradient G ₁, the fluid moves unyielded with a uniform velocity, and thus, the two slip velocities are equal. In an intermediate regime between G ₁ and a second critical pressure gradient G ₂, the fluid yields in a zone near the weak-slip wall and flows with uniform velocity near the stronger-slip wall. Beyond this regime, the fluid yields near both walls and the velocity are uniform only in the central unyielded core. It is demonstrated that the central unyielded region tends towards the midplane only if the power-law exponent is less than unity; otherwise, this region rends towards the weak-slip wall and asymmetry is enhanced. The extension of the different flow regimes depends on the channel gap; in particular, the intermediate asymmetric flow regime dominates when the gap becomes smaller than a characteristic length which incorporates the wall slip coefficients and the fluid properties. The theoretical results compare well with available experimental data on soft glassy suspensions. These results open new routes in manipulating the flow of viscoplastic materials in applications where the flow behavior depends not only on the bulk rheology of the material but also on the wall properties.     

MODELLING TWO-PHASE FLOW-EXCITED DAMPING AND FLUIDELASTIC INSTABILITY IN TUBE ARRAYS

This paper reports the results of an experimental study of the flow-induced vibration of a heat exchanger tube array subjected to two-phase cross-flow of refrigerant 11. The primary concern of the research was to develop a methodology for predicting the critical flow velocities for fluidelastic instability which better characterize the physics of two-phase flows. A new method is proposed for calculating the average fluid density and equivalent flow velocity of the two-phase fluid, using a newly developed void fraction model to account for the difference in velocity between the gas and liquid phases. Additionally, damping measurements in two-phase flow were made and compared with the data of other researchers who used a variety of modelling fluids. The results show that the two-phase damping follows a similar trend with respect to homogeneous void fraction, and when normalized, agree well with the data in the literature. The fluidelastic threshold data of several researchers who used a variety of fluids, is re-examined using the proposed void fraction model, and the results show a remarkable change in trend with flow regime. The data corresponding to the bubbly flow regime shows no significant deviation from the trend established by Connors' theory. However, the data corresponding to the intermittent flow regime show a significant decrease in stability which is nearly independent of the mass-damping parameter. It is believed that the velocity fluctuations that are inherent in the intermittent flow regime are responsible for tripping the instability, causing lower than expected stability of the bundle.     

Flow patterns and interfacial velocities near a moving contact line

Experimental data on velocity fields and flow patterns near a moving contact line is shown to be at variance with existing hydrodynamic theories. The discrepancy points to a new hydrodynamic paradox and suggests that the hydrodynamic approach may be incomplete and further parameters or forces affecting the surfaces may have to be included. A contact line is the line of intersection of three phases: (1) a solid, (2) a liquid, and (3) a fluid (liquid or gas) phase. A moving contact line develops when the contact line moves along the solid surface. A flat plate moved up and down, inside and out of a liquid pool defines a simple, reliable experimental model to characterize dynamic contact lines. Highlighted are three important conclusions from the experimental results that should be prominent in the development of new theoretical models for this flow. First, the velocity along the streamline configuring the liquid–fluid interface is remarkably constant within a distance of a couple of millimeters from the contact line. Second, the relative velocity of the liquid–fluid interface, defined as the ratio of the velocity along the interface to the velocity of the solid surface, is independent of the solid surface velocity. Third, the relative interface velocity is a function of the dynamic contact angle.     

Dynamics of particles floating on liquid flowing in a semicircular open channel

The dynamic characteristics of surface-floating particles in liquids flowing in a two-dimensional, semicircular open channel is studied experimentally. For high visibility in the experiments, relatively large particles are employed whose particle-liquid density ratio is either equal to or less than unity. Particles of different size and geometry are tested in a water-glycerin mixture. A video camera traces the pathline of each particle from which the velocity and direction of particle motion are evaluated. Liquid velocity distribution is determined by hot-film anemometry. A modified dynamics (Basset-Boussinesq-Oseen) equation is derived and numerically solved by means of a finite-difference technique to determine fluid velocity. A new dimensionless parameter is disclosed which is pertinent to both particle geometry and fluid flow conditions. It correlates particle trajectory and velocity, trajectory dispersion and fluid-particle velocity ratio.Visiting Scholar on leave from Department of Mechanical Engineering, Fukuyama University, Fukujama, Japan     

A fractal analysis of subcooled flow boiling heat transfer

A fractal model for the subcooled flow boiling heat transfer is proposed in this paper. The analytical expressions for the subcooled flow boiling heat transfer are derived based on the fractal distribution of nucleation sites on boiling surfaces. The proposed fractal model for the subcooled flow boiling heat transfer is found to be a function of wall superheat, liquid subcooling, bulk velocity of fluid (or Reynolds number), fractal dimension, the minimum and maximum active cavity size, the contact angle and physical properties of fluid. No additional/new empirical constant is introduced, and the proposed model contains less empirical constants than the conventional models. The proposed model takes into account all the possible mechanisms for subcooled flow boiling heat transfer. The model predictions are compared with the existing experimental data, and fair agreement between the model predictions and experimental data is found for different bulk flow rates.     

Photoviscous analysis of two-dimensional laminar flow in an expanding jet

This paper presents experimental measurements of the velocity profiles in an expanding fluid jet. The data were obtained from analysis of birefringence measurements over the flow field. A new method of photoviscous analysis to determine the pertinent flow information from isochromate data alone was developed as part of the investigation. A numerical method for obtaining the flow data from the optical measurements is presented.     

Reconstruction of a quasi-instantaneous image of coherent structures from hotwire signals obtained by a multi-point simultaneous measurement system

The turbulent velocity components (u, v) at 11 points in a reciprocating oscillatory turbulent flow have been measured simultaneously by a set of eleven X-type hotwire probes located in a plane perpendicular to the mean flow. Using a conditional sampling technique and a new method of data analysis for the inverse estimation of flow fields called the “virtual plate/load and MASCON model”, a quasi-instantaneous three-dimensional image of coherent structures of turbulence was first reconstructed directly from the experimental velocity data. The quasi-instantaneous image was expressed in terms of the velocity components u, v, w and the vorticity components ω_x ω_y, ω_z and we found that the large-scale coherent structure was composed of a pair of counter-rotating fluid motions with asymmetry which was quite different from that of the ensemble-averaged one. Flow patterns induced by the large-scale structure have been clarified by perspective representations visualized by computer simulations that produce timelines and streaklines of fluid particle traces. Results showed that the new experimental method was applicable for investigating the three-dimensional feature of coherent structures including asymmetry.     

On translational and rotational relative velocities of fibers and fluid in a turbulent channel flow with a backward-facing step

A dilute fiber suspension in a turbulent channel with a backward-facing step is investigated by means of Feature Tracking. Its combination with a phase-discrimination methodology, which is described in detail, allows simultaneous and separate measurement of carrier and dispersed phases velocity fields, the orientation and rotation rate of fibers as well as the fiber–fluid translational and rotational slip velocities. The patterns of fibers concentration, angular velocity and the probability distribution of fibers velocity appear to be dominated by the mechanical interactions with the wall and the local high shear rather than by near-wall turbulent structures. The translational slip velocity obtained from instantaneous data shows that fibers move faster than the surrounding fluid inside the buffer layer, the velocity gap reducing gradually when approaching the channel centerline. On the other hand, the rotational slip profile suggests a gradual decoupling of the translational and rotational dynamics. Downstream of the step, the excess of streamwise velocity displayed by fibers is still observed and extends in the free-shear region, whereas the rotation rate slip decreases at a relatively short distance from the step, as the effect of the wall presence fades away.     

Quantitative NMR velocity imaging of a main-chain liquid crystalline polymer flowing through an abrupt contraction

 The flow of isotropic and liquid crystalline (LC) hydroxypropylcellulose (HPC) aqueous solutions into an abrupt axisymmetric contraction has been quantitatively measured by pulsed-field-gradient NMR techniques. Steady-state axial velocity profiles, acquired upstream of the contraction, reveal a large contraction entry length for the LC solution. This entry flow field exists over an order of magnitude change in flow rate and is attributed to elasticity that is associated with polydomain liquid crystallinity. Pronounced, off-centerline velocity maxima (in an axisymmetric flow field) were present upstream of the contraction, in the entry flow region. Apparently, a more viscous and elastic core of fluid was present along the centerline; this fluid resisted elongational strain more than the fluid closer to the walls. Quantitative velocity profiles were extracted from displacement distributions and corrected for elongational dispersion. The isotropic solution velocity profiles matched those obtained from viscoelastic simulations using an approximate Doi-Edwards model, parameterized with independent rheological data. Received: 29 April 1999/Accepted: 30 August 1999     

Rotationally symmetric flow over a rotating disk

The rotationally symmetric flow over a rotating disk in an incompressible viscous fluid is analyzed by a new method when the fluid at infinity is in a state of rigid rotation (in the same or in the opposite sense) about the same axis as that of the disk. Asymptotic expansions for the velocity field over the entire flow field are obtained for the general class of one-parameter rotationally symmetric flows. This method is further extended to the case when a uniform suction or injection is assumed at the rotating disk. Fluid motion induced by oscillatory suction of small amplitude at the rotating disk is also discussed.An initial-value analysis reveals that resonance is possible only when the angular velocity of the rotating fluid is greater than that of the rotating disk.     

An analytical poroelastic model for laboratorial mechanical testing of the articular cartilage (AC)

The articular cartilage (AC) can be seen as a biphasic poroelastic material. The cartilage deformation under compression mainly leads to an interstitial fluid flow in the porous solid phase. In this paper, an analytical poroelastic model for the AC under laboratorial mechanical testing is developed. The solutions of interstitial fluid pressure and velocity are obtained. The results show the following facts. (i) Both the pressure and fluid velocity amplitudes are proportional to the strain loading amplitude. (ii) Both the amplitudes of pore fluid pressure and velocity in the AC depend more on the loading amplitude than on the frequency. Thus, in order to obtain the considerable fluid stimulus for the AC cell responses, the most effective way is to increase the loading amplitude rather than the frequency. (iii) Both the interstitial fluid pressure and velocity are strongly affected by permeability variations. This model can be used in experimental tests of the parameters of AC or other poroelastic materials, and in research of mechanotransduction and injury mechanism involved interstitial fluid flow.     

输液曲管平面内振动的波动方法研究

采用Flügge曲梁模拟弯曲管道,推导了管内流体的加速度,在总体轴线不可伸长假定的基础上建立了曲管平面内振动的动力学方程;采用波动方法,获得了曲管内振动波的传播和反射矩阵,提出了计算曲管平面内振动固有频率的数值方法。算例分析中,通过计算两端固定半圆形曲管的临界流速并与已有文献结果对比,验证了本文方法的正确性。最后,计算了两端固定半圆形曲管在四种不同流速下的前四阶固有频率,结果表明,管内流速的增大会降低管道的固有频率,当流速增大到某一特定值时,管道的一阶固有频率消失。     

Direct numerical simulation of particle behaviour in the wall region of turbulent flows in horizontal channels

The motion of small particles in the wall region of turbulent channel flows has been investigated using direct numerical simulation. It is assumed that the particle concentration is low enough to allow the use of one-way coupling in the calculations, i.e. the fluid moves the particles but there is no feedback from the particles on the fluid motion. The velocity of the fluid is calculated by using a pseudospectral, direct solution of the Navier-Stokes equations. The calculations indicate that particles tend to segregate into the low-speed regions of the fluid motion near the wall. The segregation tendency depends on the time constant of the particle made non-dimensional with the wall shear velocity and kinematic viscosity. For very small and very large time constants, the particles are distributed more uniformly. For intermediate time constants (of the order 3), the segregation into the low-speed fluid regions is the highest. The finding that segregation occurs for a range of particle time constants is supported by experimental results. The findings regarding the more uniform distributions, however, still remain to be verified against experimental data which is not yet available. For horizontal channel flows, it is also found that particles are resuspended by ejections (of portions of the low-speed streaks) from the wall and are, therefore, primarily associated with low-speed fluid. The smaller particles are flung further upwards and, as they fall back towards the wall, they tend to be accelerated close to the fluid velocity. The larger particles have greater inertia and, consequently, accelerate to lower velocities giving higher relative velocities. This velocity difference, as a function of wall-normal distance, follows the same trend as in experiments but is always somewhat smaller in the calculations. This appears to be due to the Reynolds number for the numerical simulation being smaller than that in the experiment. It is concluded that the average particle velocity depends not only on the wall variables for scaling, but also on outer variables associated with the mean fluid velocity and fluid depth in the channel. This is because fluid depth in combination with the wall shear velocity determines how much time a particle, of a given size and density, spends in the outer flow and, hence, how close it gets to the local fluid velocity.     

Near-wall velocity measurements by particle-shadow tracking

We report a new method to measure the velocity of a fluid in the vicinity of a wall. The method, that we call particle-shadow tracking (PST), simply consists in seeding the fluid with a small number of fine tracer particles of density close to that of the fluid. The position of each particle and of its shadow on the wall are then tracked simultaneously, allowing one to accurately determine the distance separating tracers from the wall and therefore to extract the velocity field. We present an application of the method to the determination of the velocity profile inside a laminar density current flowing along an inclined plane.     

Theory on transverse migration of particles in a turbulent two-phase suspension flow due to turbulent diffusion—I

A new theoretical model has been developed to explain the behavior of transverse particle transport in turbulent flow of a dilute two-phase suspension due to turbulent diffusion. This model is based on the ability of a particle to respond to surrounding fluid motion and depends on particle size and density relative to the carrier fluid, the fractional variation in particle concentration in the transverse direction as well as the existing turbulence structure of the surrounding fluid. The model developed in this investigation has been formulated by dividing the transverse fluid velocity, as seen by a particular particle, into two superimposed components representing, respectively, the transverse turbulent fluid fluctuations and an apparent transverse local fluid drifting velocity due to the effect on the transverse oscillatory component of fluid motion by the transverse concentration distribution of particles. A subsequent paper will show that the theory (together with other new results on the concentration effects on particle drag and lift and fluid turbulence properties) can help to explain the phenomena measured previously.     

The surface marker and micro cell method

A new method is presented for the simulation of two-dimensional, incompressible, free surface fluid flow problems. The surface marker and micro cell (SMMC) method is capable of simulating transient free surface fluid flow problems that include multivalued free surfaces, impact of free surfaces with solid obstacles and converging fluid fronts (including wave breaking). New approaches are presented for the advection of the free surface, the calculation of the tentative velocity, final velocity and pressure fields and the use of multivalued velocities to treat converging fluid fronts. Simulation results are compared with experimental results for water sloshing in a tank to demonstrate the validity of the new method. Convergence of the new method is demonstrated by a grid refinement study. © 1997 John Wiley & Sons, Ltd.     

Experimental investigation of dynamic stability of a cantilever pipe aspirating fluid

The dynamic stability of a submerged cantilever pipe conveying fluid from the free end to the fixed one is considered as one of the unresolved issues in the area of fluid–structure interaction. There is a contradiction between theoretical predictions and experiments. Reported experiments did not show any instability, while theory predicts instability beyond a critical fluid velocity. Recently, several papers appeared, improving the theoretical modelling of pipe dynamics. All theories predict instability, either oscillatory or static, referred to here as flutter and divergence, respectively. A new test set-up was designed to investigate the hypothesis that previous experimental set-ups could not allow observations of pipe instability or the pipe aspirating water is unconditionally stable. In this new test set-up, the fluid velocity could exceed the theoretically predicted critical velocities. A cantilever pipe of about 5 m length was partly submerged in water. The free open end of the pipe was in the water, whereas the fixed end was above the waterline. The experiments clearly showed that the cantilever pipe aspirating water is unstable beyond a critical velocity of water convection through the pipe. Below this velocity the pipe is stable, whereas above it the pipe shows a complex motion that consists of two alternating phases. The first phase is a nearly periodic orbital motion with maximum amplitude of a few pipe diameters, whereas the second one is a noise-like vibration with very small amplitudes. Increasing the internal fluid velocity results in a larger amplitude of the orbital motion, but does not change the pipe motion qualitatively.     

混合炸药爆速预报的新方法

为了提高混合炸药爆速尤其是含铝混合炸药爆速的计算精度,基于体积加权法思路构造一种新的混合炸药爆速计算公式,对近50种常见的混合炸药配方进行了计算,并与实验值进行对比,结果很好吻合,最大误差低于3%,平均误差低于1%,与BKW和Urizar常用方法相比,计算精度有显著提高。