Polymer solutions in co-rotating Taylor-Couette flow without vorticity

We present experimental results of the flow of dilute and semi-dilute polymer solutions in co-rotating Taylor-Couette cylinders. The experimental set-up consists of a modified Mars II rheometer (Thermo Scientific) with two drive units that are mounted opposite each other. The rotational velocities of the inner and outer cylinders are chosen in a way such that the angular velocity has a 1/r profile and the flow is free of vorticity, but the direction of elongation is not constant, but rotates with the flow. Our particle image velocimetry (PIV) measurements show that for polymer solutions without shear thinning the flow is indeed free of vorticity and is equal to a stagnation point flow at a given position and a given instant in time. In contrast, torque measurements reveal that the stresses are identical to the stresses that are present in a plane shear flow. Thus, we find that for polymer solutions a flow with vorticity and a constant direction of elongation is equal to a flow without vorticity in which the direction of elongation is rotating. Finally, we show that for shear thinning solutions the flow velocity becomes non-monotonic through the gap and resembles a pluglike profile which is known from the Poiseuille flow.     

A novel ultrasound instrument for investigation of arterial mechanics

The study of arterial mechanics concerns functional characteristics depending on wall elasticity and flow profile. Wall elasticity can be investigated through the estimation of parameters like the arterial distensibility, which is of high clinical interest because of its known correlation not only with the advanced atherosclerotic disease, but also with aging and major risk factors for cardiovascular disease. The flow velocity profile is also clinically relevant, because it modulates endothelial function and can be responsible for the development and distribution of atherosclerotic plaques. A clinically relevant variable extracted from the blood velocity profile is the wall shear rate (WSR), which represents the spatial velocity gradient near the vessel wall. This paper describes an integrated ultrasound system, capable of detecting both the velocity profile and the wall movements in human arteries. It basically consists of a PC add-on board including a single high-speed digital signal processor. This is dedicated to the analysis of echo-signals backscattered from 128 range cells located along the axis of the interrogating ultrasound (US) beam. Echoes generated from the walls (characterized by high amplitudes and low Doppler frequencies) and from red blood cells (characterized by low amplitudes and relatively high Doppler frequencies) are independently processed in real-time. Wall velocity is detected through the autocorrelation algorithm, while blood velocity is investigated through a complete spectral analysis of all signals backscattered by erythrocytes and WSR is extracted from the estimated velocity profile. Preliminary applications of the new system, including the simultaneous analysis of blood flow and arterial wall movement in healthy volunteers and in a diseased patient, are discussed, and first results are presented.     

Acoustic power measurements of a damped aeroacoustically driven resonator

Strong self-sustained acoustic oscillations may occur in a gas pipe network under certain gas flow velocities within the network. The pipe network under consideration consists of a main pipe, with a variable mean airflow, with two closed coaxial side branches of variable but equal length joined to the main pipe. Coupling between resonant acoustic standing waves and instabilities of the shear layers separating the flow in the main pipe from the stagnant gas in the closed side branches leads to strong acoustic oscillations at a frequency corresponding to the half-wavelength acoustic mode defined by the total side-branch length. An acoustic damper consisting of a variable acoustic resistance and compliance is used to dissipate power from the resonating mode. The response of the aeroacoustically driven resonator to variable damping will be examined for different fluid flow regimes as well as side-branch geometries.     

Onset of shear thickening in a simple fluid

We report on nonequilibrium molecular-dynamics simulations of the shear-thickening transition in a simple fluid under shear. We relate the shear-thickening transition to the onset of instabilities in the flow profile and to that of dramatic variations in normal stress differences. The dependence of the critical shear rate, which indicates the onset of shear thickening, on density and temperature is rationalized by introducing a ratio between two characteristic times, quantifying the short-time mobility of a particle and the deformation imposed by the applied shear rate, respectively. The shear-thickening transition is shown to occur at a constant value for this ratio for all state points studied. From a structural point of view, this transition is accompanied by the formation of clusters as recently observed in experiments on complex fluids.Received: 26 July 2004, Published online: 21 September 2004PACS: 83.60.Rs Shear rate-dependent structure (shear thinning and shear thickening) - 47.50. + d Non-Newtonian fluid flows - 83.10.Mj Molecular dynamics, Brownian dynamics     

Aeroacoustics of the swinging corrugated tube: voice of the Dragon

When one swings a short corrugated pipe segment around one's head, it produces a musically interesting whistling sound. As a musical toy it is called a "Hummer" and as a musical instrument, the "Voice of the Dragon." The fluid dynamics aspects of the instrument are addressed, corresponding to the sound generation mechanism. Velocity profile measurements reveal that the turbulent velocity profile developed in a corrugated pipe differs notably from the one of a smooth pipe. This velocity profile appears to have a crucial effect both on the non-dimensional whistling frequency (Strouhal number) and on the amplitude of the pressure fluctuations. Using a numerical model based on incompressible flow simulations and vortex sound theory, excellent predictions of the whistling Strouhal numbers are achieved. The model does not provide an accurate prediction of the amplitude. In the second part of the paper the sound radiation from a Hummer is discussed. The acoustic measurements obtained in a semi-anechoic chamber are compared with a theoretical radiation model. Globally the instrument behaves as a rotating (Leslie) horn. The effects of Doppler shift, wall reflections, bending of the tube, non-constant rotational speed on the observed frequency, and amplitude are discussed.     

应用光纤布喇格光栅传感器监测地下管道腐蚀的新方法研究

地下管线工作环境恶劣,腐蚀是引起它破坏和失效的主要原因之一.本文提出了一种应用光纤光栅传感器监测管道腐蚀的新方法.输送石油和天然气等一些物质的管道内部存在的压力会使管道发生膨胀,进而导致管道外壁产生环向变形.假定管道为无限长,管道内部压力导致的管道外壁轴向应变可忽略不计,当管线发生腐蚀时,会导致管壁变薄,管道内部的压力又可以看作基本不变,这样就会增大管道外壁的环向应变.当管道内部压力保持恒定时,管道壁厚与应变之间存在一定的反比例关系,本文方法就是将封装好的光纤光栅传感器粘贴在管道外壁,通过监测外壁环向变形来判断管道腐蚀的情况和程度.通过理论值、模拟值和试验值三者的对比分析来说明了这种新方法是可行的,并且值得被广泛应用,为埋地管道的防灾减灾提供了一种对策.     

Intraluminal fiber-optic Doppler imaging catheter for structural and functional optical coherence tomography

We describe a miniature fiber-optic Doppler imaging catheter for integrated functional and structural optical coherence tomography (OCT) imaging. The Doppler catheter can map blood flow within a vessel as well as image vessel wall structures. A prototype Doppler catheter has been developed and demonstrated for measuring the intraluminal velocity profile in a vessel phantom (conduit). A simple mathematical model is demonstrated to estimate the total flow rate. This estimation technique also enables the spatial range of flow measurements to be extended by approximately two times the normal OCT image-penetration depth. The Doppler OCT catheter could be a powerful device for cardiovascular imaging.     

Flow Rate Profiler: an instrument to measure blood velocity profiles

In this paper we present Flow Rate Profiler (FRP), an instrument for measuring the blood velocity by means of ultrasound-based techniques. The velocity is directly related to the shear rate, which is in turn proportional to the shear stress, a parameter expressing the pressure exerted by the blood on the vessel walls. The knowledge of this value is important in medicine to establish the state of the vessels, directly related to vascular diseases. FRP provides a non-invasive measure of the blood velocity by exploiting the red corpuscles property of diffusing ultrasound waves: in practice blood velocity is determined by a cross-correlation technique, which analyses the time shift between correlated subsequent echo waves, instead of frequency shift characteristic of the Doppler technique. The acquired data are then processed on a personal computer by means of mathematical techniques based on the evaluation of the correlation function, giving a reconstructed velocity profile and showing a good adherence with experimental data, since the average error is nearly the 10%. The reconstructed profile is displayed to the operator, who can follow the vessel status in real time. A few comparisons between the reconstructed and the experimental profiles are also presented, together with a study on a small set of patients suffering from artery hypertension.     

Echo tracer dispersion in flows of polymer solutions through porous media: A tool for detecting weak permeability heterogeneities?

Tracer dispersion in Newtonian and shear-thinning fluids (scleroglucan-water polymer solutions) flowing through single and double porosity grain packings has been studied experimentally using both classical transmission dispersion and echo dispersion (in the latter, the concentration variation front is pumped back through a detector at the inlet after penetrating for a chosen distance into the sample). Transmission dispersion increases markedly in both types of samples with the shear thinning index of the fluid at all Péclet numbers (except when molecular diffusion is dominant). Echo and transmission experiments give nearly identical dispersivity values for Newtonian fluids while echo dispersivity is lower than transmission for shear thinning ones. The normalized dispersivity difference has same order of magnitude for single and double porosity samples and increases with the shear thinning exponent α (by a factor of 2 between α = 0.35 and α = 0.60). This difference may be due to heterogeneities inducing permeability variations of small amplitude over distances of the order of the sample section : their influence on tracer dispersion is partly reversible with respect to a change of the flow direction and is only detectable if it is amplified by the shear-thinning properties of the fluid. Received 19 September 2001     

Properties of a soft-sphere liquid from non-Newtonian molecular dynamics

A soft-sphere, inverse-12 liquid is simulated in both the isokinetic-isochoric and the isokinetic-isobaric ensemble using nonequilibrium molecular dynamics. The simulation for the isobaric ensemble is discussed in detail. The non-Newtonian characteristics of the liquid are clearly demonstrated; namely, the shear-rate-dependent pressure and density (shear dilatancy), the shear-rate-dependent shear viscosity (shear thinning), and evidence of normal pressure differences. For the first time, it is clearly shown that a significant component of isobaric shear thinning is due to shear dilatancy. The isochoric and isobaric results are checked for consistency. Simple empirical relations for the equation of state and transport properties of the fluid are presented.This is a publication in part of the National Institute of Standards and Technology (formerly NBS) and is not subject to copyright.     

Simplification of formulas for compression wave and oscillating flow in circular pipe

We had already obtained the analytical solutions for the compression wave and steady-state oscillating flow in a pipe with a circular cross section [Sato Y, Kanki H. Formulas for compression wave and oscillating flow in circular pipe. To be published in Appl Acoust [accepted 11 Sept. 2006]]. This work contains three key-components. The first key is to simplify the formulas using the unique mathematic technique without losing the accuracy. Simplifying the formulas is the one of the most important factors for formulas used in engineering use. The results will enable us to greatly reduce the work and computation costs. The second is to verify the flow distribution calculated by our formulas. The third is to study the possibility of application of our method to the analysis in turbulent flow region. Kawata et al. have represented the validity of one-dimensional quasi-analysis in turbulent flow region, using the method of D’Souza et al. and the shear viscosity. Therefore, in this paper the validity of the analysis in a turbulence region was verified by proving theoretically that the methods of D’Souza et al. and ours are intrinsically equivalent. The proof of equivalence was accomplished using the formulas simplified in this paper.     

Doppler ultrasound signals simulation from vessels with various stenosis degrees

As a non-invasive method, the Doppler ultrasound technique is used to detect the vessel stenosis. To search for characteristics of Doppler ultrasound signals sensitive to the stenosis, a computer simulation approach is proposed in this paper to generate Doppler ultrasound signals from vessels with various stenosis degrees. The blood flow velocity distribution in a stenosed vessel is firstly calculated using the transient finite element method (FEM). Then the power spectral density of Doppler signals is estimated using the overall-distribution nonparametric estimation method. Finally Doppler signals are generated using the cosine-superposed method. The proposed approach is proved to be useful for simulating Doppler ultrasound signals from vessels with various stenosis degrees. It is also shown that characteristics of Doppler ultrasound signals may be used to estimate the vessel's stenosis degree.     

Importance of depletion interactions for structure and dynamics of ferrofluids

The influence of attractive depletion forces on the structure and dynamics of ferrofluids is studied by computer simulations. In the presence of a magnetic field, we find that sufficiently strong depletion forces lead to an assembly of particle chains into columnar structures with hexagonal ordering inside the columns. In a planar shear flow, this ordering is destroyed, leading to strong shear thinning behavior. A pronounced anisotropy of the shear viscosity is observed. The viscosity is found to be largest when the magnetic field is oriented in the gradient direction of the flow.     

Molecular origin of shear thickening in transient polymer networks: A molecular dynamics study

This paper proposes a simple model of transient networks of telechelic associating polymers for molecular simulations and reports the main results obtained by molecular dynamics on the rheological properties of the transient networks. The steady shear viscosity obtained by the non-equilibrium molecular dynamics simulation exhibits shear thickening at moderate shear rates and shear thinning at larger shear rates. The behavior is similar to that observed in experiments of telechelic associating polymers. By analyzing the distribution function of the end-to-end vector of bridge chains as a function of the shear rate, we find that shear thickening is mainly caused by the stress from the bridge chains highly stretched by shear flow. We also find that fracture of the transient network occurs in the shear-thinning regime at high shear rates.     

On the acoustic modes in a duct containing a parabolic shear flow

The exact solution of the acoustic wave equation in an unidirectional shear flow with a parabolic velocity profile is obtained, representing sound propagation in a plane, parallel walled duct, with two boundary layers over rigid or impedance walls. It is shown that there are four cases, depending on the critical level(s) where the Doppler shifted frequency vanishes: (i) for propagation upstream the critical levels are outside the duct (case II); (ii) for propagation downstream there may be two (case IV), one (case I) or no (case III) critical level inside the duct. The acoustic wave equation is transformed in each of the four cases to particular forms of the extended hypergeometric equation, which has power series solutions, some involving logarithmic singularities. In the cases where critical levels occur, at real or ‘imaginary’ distance, matching of two or three pairs of solutions, valid over regions each overlapping the next, is needed. The particular case of the parabolic velocity profile is used to address general properties of sound in unidirectional shear flows. For example, it is shown that for ducted shear flows, there exist a pair of even and odd eigenfunctions, in the absence of critical levels. It is also proved, in more than one instance, that there is no single set of eigenvalues and eigenfunctions valid across one or two shear layers. This leads to the general conjecture, considering the acoustics of shear flows in ducts, that critical levels separate regions with distinct sets of eigenvalues and eigenfunctions.     

Existence of a Hartmann layer in the peristalsis of Sisko fluid

Analytical solutions for the peristaltic flow of a magneto hydrodynamic(MHD) Sisko fluid in a channel, under the effects of strong and weak magnetic fields, are presented. The governing nonlinear problem, for the strong magnetic field,is solved using the matched asymptotic expansion. The solution for the weak magnetic field is obtained using a regular perturbation method. The main observation is the existence of a Hartman boundary layer for the strong magnetic field at the location of the two plates of the channel. The thickness of the Hartmann boundary layer is determined analytically. The effects of a strong magnetic field and the shear thinning parameter of the Sisko fluid on the velocity profile are presented graphically.     

Shear thinning and orientational ordering of wormlike micelles

Shear thinning and orientation of cylindrical surfactant and block copolymer micelles was investigated by rheo-SANS (small-angle neutron scattering) experiments. Shear thinning and orientation occur for shear rates (.)gamma tau(dis)>1, where tau(dis) is the disentanglement time of the micelles. Micelles align in the flow direction with an orientational distribution that can be well described by an Onsager-type distribution function. Over nearly the whole range of concentrations and for all cylindrical micelles investigated, the shear viscosity eta follows a simple eta approximately e(-aS) behavior as a function of the orientational order parameter S with the same prefactor a.     

Vortex dynamics of a trapezoidal bluff body placed inside a circular pipe

The influence of Reynolds number and blockage ratio on the vortex dynamics of a trapezoidal bluff body placed inside a circular pipe is studied experimentally and numerically. Low aspect ratio, high blockage ratio, curved end conditions (junction of pipe and bluff body), axisymmetric upstream flow with shear and turbulence are some of the intrinsic features of this class of bluff body flows which have been scarcely addressed in the literature. A large range (200:200,000) of Reynolds number (Re_D) is covered in this study, encompassing all the three pipe flow regimes (laminar, transition and turbulent). Four different flow regimes are defined based on the distinct features of Strouhal number (St)–Re_D relation: steady, laminar irregular, transition and turbulent. The wake in the steady regime is stationary with no oscillations in the shear layer. The laminar regime is termed as irregular owing to irregular vortex shedding. The vortex shedding in this regime is observed to be symmetric. The emergence of separation bubble downstream of the bluff body on either side is another interesting feature of this regime, which is further observed to be symmetric. Two pairs of mean streamwise vortices are noticed in the near-wake regime, which are termed as reverse dipole-type wake topology. Beyond the irregular laminar regime, the Strouhal number falls gradually and vortex shedding becomes more periodic. This regime is named transition and occurs close to the Reynolds number at which transition to turbulence takes place in a fully developed pipe. The turbulent regime is characterised by a nearly constant Strouhal number. Typical Karman-type vortex shedding is noticed in this regime. The convection velocity, wake width formation length and irrecoverable pressure loss are quantified to highlight the influence of blockage ratio. These results will be useful to develop basic understanding of vortex dynamics of confined bluff body flow for several practical applications.     

光学相干层析多普勒成像功能拓展研究

光学多普勒成像(Optical Doppler tomography,ODT)是一种结合了光学相干层析成像技术(Opticalcoherence tomography,OCT)和多普勒流速仪的非侵入、非接触的成像技术,能够实现对高散介质组织内部的血管分布和血液流速的探测。阐述了基于数字希尔伯特变换的相位分离多普勒光学相干层析成像技术的工作原理,并且通过对玻璃毛细管和生物芯片微通道管中聚苯乙烯溶液流速的实验测量,准确测量管内微粒缓慢移动时的多普勒频移量,获得了玻璃管内和生物芯片微通道管中流速分布曲线,证实了所提方法的可行性。获取的多普勒图像具有较高的空间分辨力和速度分辨力,在未来的临床应用中有潜在的应用价值。     

Flow spectra from spectral power density calculations for pulsed Doppler

Range gated pulsed Doppler can be used to make localized velocity measurements within a blood vessel. Both the transducer and the sample volume are of finite size, and this prohibits the measurement of velocity at a point. A spectral flow profile can be created by stepping a sufficiently small sample volume across the lumen of a vessel. However no such set of spectra will correspond directly to the true velocity profile. In this study we developed a systematic theoretical treatment which allows Doppler spectral power density (SPD) functions to be calculated under a very wide range of conditions. Simulated flow spectra were created from sets of these spectra. The model is based on the beam intensity weighted volume method and incorporates, through the idea of a spread function, Guidi's individual flow line spectrum. Our method can be applied for different spread functions; with beam profiles which are uniform, Gaussian or arbitrarily narrow (needle beam); with range gated sample volumes which can be maximal (CW-type) or minimal (PW-type); and for beams which intersect the flow tube axis, or are off centre. Under all conditions we find the spread function parameter k, equal to the ratio of the central Doppler shift to half the bandwidth, plays a key role. After formulating the model analytically, we sought simplifications to allow results to be obtained from simple, practical formulae. Spread and unspread SPD functions are in most cases given as single integrals which contain measurable physical parameters and can be easily evaluated numerically. Model results are presented for flow spectra of parabolic flow, illustrating the interplay between different factors in determining the appearance of spectral flow profiles.