Direct numerical simulation of the turbulent flows of power-law fluids in a circular pipe

Fully developed turbulent pipe flows of power-law fluids are studied by means of direct numerical simulation. Two series of calculations at generalised Reynolds numbers of approximately 10000 and 20000 were carried out. Five different power law indexes n from 0.4 to 1 were considered. The distributions of components of Reynolds stress tensor, averaged viscosity, viscosity fluctuations, and measures of turbulent anisotropy are presented. The friction coefficient predicted by the simulations is in a good agreement with the correlation obtained from experiment. Flows of power-law fluids exhibit stronger anisotropy of the Reynolds stress tensor compared with the flow of Newtonian fluid. The turbulence anisotropy becomes more significant with the decreasing flow index n. An increase in apparent viscosity away from the wall leads to the damping of the wall-normal velocity pulsations. The suppression of the turbulent energy redistribution between the Reynolds stress tensor components observed in the simulations leads to a strong domination of the axial velocity pulsations. The damping of wall-normal velocity pulsations leads to a reduction of the fluctuating transport of momentum from the core toward the wall, which explains the effect of drag reduction.     

Mathematical model for a Herschel-Bulkley fluid flow in an elastic tube

The constitution of blood demands a yield stress fluid model, and among the available yield stress fluid models for blood flow, the Herschel-Bulkley model is preferred (because Bingham, Power-law and Newtonian models are its special cases). The Herschel-Bulkley fluid model has two parameters, namely the yield stress and the power law index. The expressions for velocity, plug flow velocity, wall shear stress, and the flux flow rate are derived. The flux is determined as a function of inlet, outlet and external pressures, yield stress, and the elastic property of the tube. Further when the power-law index n = 1 and the yield stress τ ₀ → 0, our results agree well with those of Rubinow and Keller [J. Theor. Biol. 35, 299 (1972)]. Furthermore, it is observed that, the yield stress and the elastic parameters (t ₁ and t ₂) have strong effects on the flux of the non-Newtonian fluid flow in the elastic tube. The results obtained for the flow characteristics reveal many interesting behaviors that warrant further study on the non-Newtonian fluid flow phenomena, especially the shear-thinning phenomena. Shear thinning reduces the wall shear stress.     

分形油藏非牛顿幂律流体低速非达西不稳定渗流的组合数学模型

结合分形理论与渗流理论,对分形油藏非牛顿幂律流体低速非达西不稳定渗流的试井分析问题的数学模型进行了推导.该分形油藏模型由内域为非牛顿幂律流体低速非达西渗流,外域为非牛顿幂律流体达西渗流的同心圆域组成.在考虑井筒储存、表皮效应影响下,建立了该油藏的不稳定渗流有效井径组合数学模型,在3种外边界条件下求出了两个区域内压力在Laplace空间的解析解,应用Stehfest数值反演方法求得井底的无因次压力,分析了井底压力动态特征和参数影响.非牛顿幂律流体的幂律指数、分形参数均对典型曲线产生较大的影响,呈现出与牛顿流体和均质油藏明显不同的特征.这对非均质油藏非牛顿流体的不稳定试井分析及研究其非线性渗流特征均十分重要.     

Bulk viscous L.R.S. Bianchi type V tilted stiff fluid cosmological model in general relativity

Locally rotationally symmetric (L.R.S.) Bianchi type V bulk viscous tilted stiff fluid cosmological model is investigated. To get the deterministic model of the universe, we have also assumed a condition A=Bⁿ between metric potentials A, B where n is the constant. The behaviour of the model in presence and absence of bulk viscosity and singularities in the model are also discussed. In general, the models represent accelerating, shearing, tilted and non-rotating universe. The models have point type singularity in presence and absence of bulk viscosity both.     

The inertial dynamics of thin film flow of non-Newtonian fluids

Consider the flow of a thin layer of non-Newtonian fluid over a solid surface. I model the case where the viscosity depends nonlinearly on the shear-rate; power law fluids are an important example, but the analysis here is for general nonlinear dependence. The modelling allows for large changes in film thickness provided the changes occur over a relatively large enough lateral length scale. Modifying the surface boundary condition for tangential stress forms an accessible foundation for the analysis where flow with constant shear is a neutral critical mode, in addition to a mode representing conservation of fluid. Perturbatively removing the modification then constructs a model for the coupled dynamics of the fluid depth and the lateral momentum. For example, the results model the dynamics of gravity currents of non-Newtonian fluids when the flow is not creeping.     

Adhesion and fingering in the lifting Hele-Shaw cell: Role of the substrate

The lifting Hele-Shaw cell (LHSC) is used to study adhesion as well as viscous fingering. In the present paper we report a series of observations of development of the interface for different viscous fluids, both Newtonian and non-Newtonian, in a LHSC operated at a constant lifting force. Glass and perspex are used as the plates in two different sets of experiments. The objectives are 1) to measure the time required to separate the plates as a function of the lifting force and 2) to note the force above which viscous fingering appears. We find that for the Newtonian fluids, the plate separation time follows a universal power law with the lifting force, irrespective of fluid and substrate. The non-Newtonian fluids too, with proper scaling obey the same power law. The appearance of fingering, however, depends on the properties of the fluid as well as the substrate. We suggest a modified form of the capillary number which controls the onset of fingering; this new quantity, termed the “fingering parameter” involves the dielectric constants of the substrate and fluid in addition to the viscosity and surface tension.     

Marangoni flow and mass transfer of power-law non-Newtonian fluids over a disk with suction and injection

We scrutinize the approximate analytical solutions by the optimal homotopy analysis method (OHAM) for the flow and mass transfer within the Marangoni boundary layer of power-law fluids over a disk with suction and injection in the present paper. Concentration distribution on the surface of a disk varies in a power-law form. The non-Newtonian fluid flow is due to the surface concentration gradient without considering gravity and buoyancy. According to the conservation of mass, momentum and concentration, the governing partial differential equations are established, and the appropriate generalized Kármán transformation is found to reduce them to the nonlinear ordinary differential equations. OHAM is used to access the approximate analytical solution. The influences of Marangoni the number, suction/injection parameters and power-law exponent on the flow and mass transfer are examined.     

Rapid communication: Transverse spin with coupled plasmons

The steady two-dimensional flow and heat transfer of a non-Newtonian power-law nanofluid over a stretching surface under convective boundary conditions and temperature-dependent fluid viscosity has been numerically investigated. The power-law rheology is adopted to describe non-Newtonian characteristics of the flow. Four different types of nanoparticles, namely copper (Cu), silver (Ag), alumina (Al ₂ O ₃) and titanium oxide (TiO ₂) are considered by using sodium alginate (SA) as the base non-Newtonian fluid. Lie symmetry group transformations are used to convert the boundary layer equations into non-linear ordinary differential equations. The transformed equations are solved numerically by using a shooting method with fourth-order Runge–Kutta integration scheme. The results show that the effect of viscosity on the heat transfer rate is remarkable only for relatively strong convective heating. Moreover, the skin friction coefficient and the rate of heat transfer increase with an increase in Biot number.     

微扩张管道内幂律流体非定常电渗流动

研究了电渗驱动下幂律流体在有限长微扩张管道内非稳态流动特性.基于Ostwald-de Wael幂律模型,采用高精度紧致差分离散二维Poisson-Nernst-Planck方程及修正的Cauchy动量方程,数值模拟了初始及稳态时刻微扩张管道内幂律流体电渗流流场分布情况,研究了管道截面改变对幂律流体无量纲剪切应变率及无量纲表观黏度的影响,以及无量纲表观黏度对拟塑性流体与胀流型流体流速分布的影响.数值模拟结果显示,当扩张角和无量纲电动宽度一定时,电场驱动下的幂律流体在近壁区域速度响应都很快;初始时刻,近壁处表观黏度的变化受到剪切应变率变化的影响,从而影响了三种幂律流体速度峰值的分布,出现拟塑性流体流速在扩张段上游及扩张段近壁处速度峰值均为幂律流体中最大、而在扩张段下游三种幂律流体速度峰值相近的现象;稳态时刻,幂律流体速度剖面呈现塞型分布,且满足连续性条件下,幂律流体流速随扩张管半径增大而减小,牛顿流体流动规律与宏观尺度下流动规律相同;初始时刻,在相同电动宽度、不同壁面电势作用下,幂律流体在扩张管近壁处剪切应变率分布的差异导致表观黏度分布的差异,并最终导致拟塑性流体与胀流型流体流速分布的差异.     

Analysis of the flow of non-Newtonian visco-elastic fluids in fractal reservoir with the fractional derivative

In both the oil reservoir engineering and seepage flow mechanics, heavy oil with relaxation property shows non-Newtonian rheological characteristics. The relationship between shear rate g& and shear stress t is nonlinear. Because of the relaxation phenomena of heavy oil flow in porous media, the equation of motion can be written as[1] 2,rrvpqkppqtrrtll秏骣+=-+琪抖桫 (1) where lv and lp are velocity relaxation and pressure retardation times. For most porous media, the above motion equation (1)…     

Scaling of hysteresis in a multi-dimensional all-optical bistable system

We analyse the hysteresis enlargements of an optical bistable system involving three dynamical variables. We investigate, both experimentally and numerically, the local dynamics of the up- and down-switching process versus the sweeping frequency of the control parameter. In particular, we delineate the domain of validity of the scaling law predicted for one-dimensional systems. At high sweeping frequency, we show the appearance of another asymptotic scaling low in . Thereafter, we analyse the global evolution of the hysteresis loop induced by these processes. At low frequency, a scaling law is retrieved, whereas at high frequency, the dynamical behaviour is shown to strongly depend on the particular shape of the bistability curve. Received: 14 September 1998 / Received in final form: 15 February 1999     

Dynamical mean-field theory of a simplified double-exchange model

Simplified double-exchange model including transfer of the itinerant electrons with spin parallel to the localized spin in the same site and the indirect interaction J of kinetic type between localized spins is comprihensively investigated. The model is exactly solved in infinite dimensions. The exact equations describing the main ordered phases (ferromagnetic and antiferromagnetic) are obtained for the Bethe lattice with (z is the coordination number) in analytical form. The exact expression for the generalized paramagnetic susceptibility of the localized-spin subsystem is also obtained in analytical form. It is shown that temperature dependence of the uniform and the staggered susceptibilities has deviation from Curie-Weiss law. Dependence of Curie and Néel temperatures on itinerant-electron concentration is discussed to study instability conditions of the paramagnetic phase. Anomalous temperature behaviour of the chemical potential, the thermopower and the specific heat is investigated near the Curie point. It is found for J=0 that the system is unstable towards temperature phase separation between ferromagnetic and paramagnetic states. A phase separation connected with antiferromagnetic and the paramagnetic phases can occur only at . Zero-temperature phase diagram including the phase separation between ferromagnetic and antiferromagnetic states is given. Received 28 May 1999 and Received in final form 14 July 1999     

Modelling capillary filling dynamics using lattice Boltzmann simulations

We investigate the dynamics of capillary filling using two lattice Boltzmann schemes: a liquid-gas model and a binary model. The simulation results are compared to the well-known Washburn's law, which predicts that the filled length of the capillary scales with time as l ∝ t ^1/2. We find that the liquid-gas model does not reproduce Washburn's law due to condensation of the gas phase at the interface, which causes the asymptotic behaviour of the capillary penetration to be faster than t ^1/2. The binary model, on the other hand, captures the correct scaling behaviour when the viscosity ratio between the two phases is sufficiently high.     

Gravity driven instabilities in miscible non-Newtonian fluid displacements in porous media

Gravity driven instabilities in model porous packings of 1 mm diameter spheres are studied by comparing the broadening of the displacement front between fluids of slightly different densities in stable and unstable configurations. Water, water–glycerol and water–polymer solutions are used to vary independently viscosity and molecular diffusion and study the influence of shear-thinning properties. Both injected and displaced solutions are identical but for a different concentration of NaNO₃ salt used as an ionic tracer and to introduce the density contrast. Dispersivity in stable configuration increases with polymer concentration – as already reported for double porosity packings of porous grains. Gravity-induced instabilities are shown to develop below a same threshold Péclet number Pe for water and water–glycerol solutions of different viscosities and result in considerable increases of the dispersivity. Measured threshold Pe values decrease markedly on the contrary with polymer concentration. The quantitative analysis demonstrates that the development of the instabilities is controlled by viscosity through a characteristic gravity number G (ratio between hydrostatic and viscous pressure gradients). A single threshold value of G accounts for results obtained on Newtonian and non-Newtonian solutions.     

Generalized second law of thermodynamics for FRW cosmology with power-law entropy correction

In this work, we have considered the power-law correction of entropy on the horizon. If the flat FRW Universe is filled with the n components fluid with interactions, the GSL of thermodynamics for apparent and event horizons have been investigated for equilibrium and non-equilibrium cases. If we consider a small perturbation around the de Sitter spacetime, the general conditions of the validity of GSL have been found. Also if a phantom dominated Universe has a pole-like type scale factor, the validity of GSL has also been analyzed. Further we have obtained constraints on the power-law parameter α in the phantom and quintessence dominated regimes. Finally we obtain conditions under which GSL breaks down in a cosmological background.     

Stability of fluid flow past a membrane

The stability of the flow of a fluid past a solid membrane of infinitesimal thickness is investigated using a linear stability analysis. The system consists of two fluids of thicknesses R and H R and bounded by rigid walls moving with velocities and , and separated by a membrane of infinitesimal thickness which is flat in the unperturbed state. The fluids are described by the Navier-Stokes equations, while the constitutive equation for the membrane incorporates the surface tension, and the effect of curvature elasticity is also examined for a membrane with no surface tension. The stability of the system depends on the dimensionless strain rates and in the two fluids, which are defined as and for a membrane with surface tension , and and for a membrane with zero surface tension and curvature elasticity K. In the absence of fluid inertia, the perturbations are always stable. In the limit , the decay rate of the perturbations is O(k ³ ) smaller than the frequency of the fluctuations. The effect of fluid inertia in this limit is incorporated using a small wave number asymptotic analysis, and it is found that there is a correction of smaller than the leading order frequency due to inertial effects. This correction causes long wave fluctuations to be unstable for certain values of the ratio of strain rates and ratio of thicknesses H. The stability of the system at finite Reynolds number was calculated using numerical techniques for the case where the strain rate in one of the fluids is zero. The stability depends on the Reynolds number for the fluid with the non-zero strain rate, and the parameter , where is the surface tension of the membrane. It is found that the Reynolds number for the transition from stable to unstable modes, , first increases with , undergoes a turning point and a further increase in the results in a decrease in . This indicates that there are unstable perturbations only in a finite domain in the plane, and perturbations are always stable outside this domain. Received: 29 May 1997 / Revised: 9 October 1997 / Accepted: 26 November 1997     

Asymptotic analysis of wall modes in a flexible tube

The stability of wall modes in a flexible tube of radius R surrounded by a viscoelastic material in the region R < r < H R in the high Reynolds number limit is studied using asymptotic techniques. The fluid is a Newtonian fluid, while the wall material is modeled as an incompressible visco-elastic solid. In the limit of high Reynolds number, the vorticity of the wall modes is confined to a region of thickness in the fluid near the wall of the tube, where the small parameter , and the Reynolds number is , and are the fluid density and viscosity, and V is the maximum fluid velocity. The regime is considered in the asymptotic analysis, where G is the shear modulus of the wall material. In this limit, the ratio of the normal stress and normal displacement in the wall, , is only a function of H and scaled wave number . There are multiple solutions for the growth rate which depend on the parameter .In the limit , which is equivalent to using a zero normal stress boundary condition for the fluid, all the roots have negative real parts, indicating that the wall modes are stable. In the limit , which corresponds to the flow in a rigid tube, the stable roots of previous studies on the flow in a rigid tube are recovered. In addition, there is one root in the limit which does not reduce to any of the rigid tube solutions determined previously. The decay rate of this solution decreases proportional to in the limit , and the frequency increases proportional to . Received: 5 November 1997 / Revised: 10 March 1998 / Accepted: 29 April 1998     

Faraday instability with a polymer solution

We present an experimental study of the Faraday instability in which we compare the behavior of a Newtonian fluid (water-glycerine mixture) with that of a semi-dilute non-Newtonian solution of high molecular weight polymer. We show that although the dispersion relation of surface waves, derived for a layer of inviscid fluid, remains valid in that particular non-Newtonian case, the behavior of the instability threshold with frequency strongly differs from the Newtonian case. We explain this effect as a result of a frequency-dependent viscosity. The linear stability analysis of the non-Newtonian case shows a perfect agreement with the experimental results both for the dispersion relation and for the reduction of the instability threshold. We discuss the use of the characteristics of the Faraday experiment as a measurement tool to determine frequency dependent properties of non-Newtonian fluids. Received 5 January 1999