Newtonian flow in a triangular duct with slip at the wall

We consider the Newtonian Poiseuille flow in a tube whose cross-section is an equilateral triangle. It is assumed that boundary slip occurs only above a critical value of the wall shear stress, namely the slip yield stress. It turns out that there are three flow regimes defined by two critical values of the pressure gradient. Below the first critical value, the fluid sticks everywhere and the classical no-slip solution is recovered. In an intermediate regime the fluid slips only around the middle of each boundary side and the flow problem is not amenable to analytical solution. Above the second critical pressure gradient non-uniform slip occurs everywhere at the wall. An analytical solution is derived for this case and the results are discussed.     

Slip yield stress effects in start-up Newtonian Poiseuille flows

Analytical solutions are derived for various start-up Newtonian Poiseuille flows assuming that slip at the wall occurs when the wall shear stress exceeds a critical value, known as the slip yield stress. Two distinct regimes characterise the steady axisymmetric and planar flows, which are defined by a critical value of the pressure gradient. If the imposed pressure gradient is below this critical value, the classical no-slip, start-up solution holds. Otherwise, no-slip flow occurs only initially, for a finite time interval determined by a critical time, after which slip does occur. For the annular case, there is an additional intermediate (steady) flow regime where slip occurs only at the inner wall, and hence, there exist two critical values of the pressure gradient. If the applied pressure gradient exceeds both critical values, the velocity evolves initially with no-slip at both walls up to the first critical time, then with slip only along the inner wall up to the second critical time and finally with slip at both walls.     

Slip at the interface between immiscible polymer melts I: method to measure slip

Slip at the interface between immiscible polymer melts remains poorly understood. A method that relies solely on rheological measurements to obtain the interfacial slip velocity uses the slip-induced deviation in the flow variables. To use the method, accurate estimates of the flow variables under the assumption of no-slip are necessary. Although such estimates can be easily derived under some cases, in general, this is not straightforward. Therefore, methods to determine the interfacial slip velocity without using estimates for the flow variables under no-slip conditions are desirable. In this work, we focus on investigations of slip at the interface between two immiscible polymer melts undergoing two-phase coaxial flow. To enable such investigations, we have adapted the Mooney method, usually used to investigate wall slip, to investigate polymer/polymer interfacial slip. Using this method, we have measured the slip velocity at the interface between polypropylene and polystyrene as a function of the interfacial stress. To determine the validity of the modified Mooney method, we also determine the slip velocity using the slip-induced deviation in the flow variables. To enable this determination, we use polypropylene and polystyrene with almost identical shear rate-dependent viscosities over a range of shear rates. The slip velocity obtained from the modified Mooney method displayed excellent agreement with that determined using the deviation from no-slip. In agreement with prior work, the dependence of the slip velocity on the interfacial stress is a power-law. Our investigation spans a sufficiently wide range of interfacial stress to enable the direct observation of two power-law regimes and also the transition between the two regimes. We also find that the power-law exponent of approximately 3 at low stresses decreases to approximately 2 at high stresses.     

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.     

Two-phase flow instabilities in a silicon microchannels heat sink

Two-phase flow instabilities are highly undesirable in microchannels-based heat sinks as they can lead to temperature oscillations with high amplitudes, premature critical heat flux and mechanical vibrations. This work is an experimental study of boiling instabilities in a microchannel silicon heat sink with 40 parallel rectangular microchannels, having a length of 15 mm and a hydraulic diameter of 194 μm. A series of experiments have been carried out to investigate pressure and temperature oscillations during the flow boiling instabilities under uniform heating, using water as a cooling liquid. Thin nickel film thermometers, integrated on the back side of a heat sink with microchannels, were used in order to obtain a better insight related to temperature fluctuations caused by two-phase flow instabilities. Flow regime maps are presented for two inlet water temperatures, showing stable and unstable flow regimes. It was observed that boiling leads to asymmetrical flow distribution within microchannels that result in high temperature non-uniformity and the simultaneously existence of different flow regimes along the transverse direction. Two types of two-phase flow instabilities with appreciable pressure and temperature fluctuations were observed, that depended on the heat to mass flux ratio and inlet water temperature. These were high amplitude/low frequency and low amplitude/high frequency instabilities. High speed camera imaging, performed simultaneously with pressure and temperature measurements, showed that inlet/outlet pressure and the temperature fluctuations existed due to alternation between liquid/two-phase/vapour flows. It was also determined that the inlet water subcooling condition affects the magnitudes of the temperature oscillations in two-phase flow instabilities and flow distribution within the microchannels.     

二阶滑移边界对微型气浮轴承稳态性能的影响

考虑微型气浮轴承的尺寸特征，内部气流不再满足连续流的假设，根据Knudsen数可确定内 部气流为滑移流. 分别利用一阶速度滑移模型和二阶速度滑移模型对连续流的状态方程进行 修改，得到一阶滑移流和二阶滑移流机制下修正的雷诺方程. 利用有限差分法对连续流、一 阶滑移流和二阶滑移流的雷诺方程分别求解，得出相应的承载力和偏位角. 经过对比分析， 发现采用滑移流模型得到的轴承的稳态力学性能与连续流机制的结果存在较大差异，一阶滑 移流与二阶滑移流间的差异随偏心率增加而增加. 说明在MEMS环境下必须考虑滑移流效应 对微型气浮轴承稳态力学性能的影响. 在大偏心率工作状态下，二阶滑移流模型能够得到最 好的结果.     

Cessation of Newtonian circular and plane Couette flows with wall slip and non-zero slip yield stress

We solve analytically the cessation flows of a Newtonian fluid in circular and plane Couette geometries assuming that wall slip occurs provided that the wall shear stress exceeds a critical threshold, the slip yield stress. In steady-state, slip occurs only beyond a critical value of the angular velocity of the rotating inner cylinder in circular Couette flow or of the speed of the moving upper plate in plane Couette flow. Hence, in cessation, the classical no-slip solution holds if the corresponding wall speed is below the critical value. Otherwise, slip occurs only initially along both walls. Beyond a first critical time, slip along the fixed wall ceases, and beyond a second critical time slip ceases also along the initially moving wall. Beyond this second critical time no slip is observed and the decay of the velocity is faster. The velocity decays exponentially in all regimes and the decay is reduced with slip. The effects of slip and the slip yield stress are discussed.     

The Effect of the Mach Number on a Turbulent Backward-Facing Step Flow

The flow around a backward-facing step in the sub-, trans- and supersonic regimes was investigated at the Trisonic Wind Tunnel Munich with particle image velocimetry and dynamic pressure measurements. These two techniques were combined to simultaneously measure and correlate the velocity fluctuations in a streamwise vertical plane with the pressure fluctuations on the reattachment surface. The results show that the dynamic loads on the reattachment surface increase from subsonic up to the transonic regime while the mean reattachment location moves downstream. As soon as the flow becomes locally supersonic aft of the backward-facing step, the mean reattachment location suddenly moves upstream while the normalized dynamic loads drastically decrease. By correlating the velocity and the dynamic pressure data, it was shown that a clear separation between outer flow and the flow close to the surface aft of the step is responsible for the drastic load reduction. Due to the large difference in pressure/density, the disturbances from the locally supersonic flow do not have an effect on the flow close to the surface. This is also reflected in the power spectral densities of the pressure fluctuations on the surface, showing that at supersonic free-stream Mach numbers a low-frequency pumping motion of the locally subsonic flow is the dominant mode, while in sub-/transonic flow Kelvin-Helmholtz instabilities and a cross-pumping motion of the shear layer dominate the dynamic loads.     

Particle dynamics calculations of wall stresses and slip velocities for couette flow of smooth inelastic spheres

The objective of this work is to characterize the effects of boundary geometry on the flow of dry granular materials composed of smooth, inelastic spheres between parallel, bumpy walls in the absence of gravity. Particle dynamic simulations are done in which wall stresses and slip velocities are computed over a wide range of parameters, including shear gap height, geometry of the wall particles, wall to free particle diameter ratio, and normal restitution coefficient. Calculated wall stresses and slip velocities are found to be highly sensitive to boundary roughness, which is characterized in terms of the mean spacing value for the wall particles. Most noticeable is a pronounced stress drop for dense flows with an associated large slip velocity in a system having a narrow shear gap height.     

On the stability of the shear flow of a viscoelastic fluid with slip along the fixed wall

In this paper we solve the time-dependent shear flow of an Oldroyd-B fluid with slip along the fixed wall. We use a non-linear slip model relating the shear stress to the velocity at the wall and exhibiting a maximum and a minimum. We assume that the material parameters in the slip equation are such that multiple steady-state solutions do not exist. The stability of the steady-state solutions is investigated by means of a one-dimensional linear stability analysis and by numerical calculations. The instability regimes are always within or coincide with the negative-slope regime of the slip equation. As expected, the numerical results show that the instability regimes are much broader than those predicted by the linear stability analysis. Under our assumptions for the slip equation, the Newtonian solutions are stable everywhere. The interval of instability grows as one moves from the Newtonian to the upper-convected Maxwell model. Perturbing an unstable steady-state solution leads to periodic solutions. The amplitude and the period of the oscillations increase with elasticity.     

Finite Element Modelling of Flow Through a Porous Medium Between Two Parallel Plates Using The Brinkman Equation

Despite the widespread use of the Darcy equation to model porous flow, it is well known that this equation is inconsistent with commonly prescribed no slip conditions at flow domain walls or interfaces between different sections. Therefore, in cases where the wall effects on the flow regime are expected to be significant, the Darcy equation which is only consistent with perfect slip at solid boundaries, cannot predict velocity and pressure profiles properly and alternative models such as the Brinkman equation need to be considered. This paper is devoted to the study of the flow of a Newtonian fluid in a porous medium between two impermeable parallel walls at different Darcy parameters (Da). The flow regime is considered to be isothermal and steady. Three different flow regimes can be considered using the Brinkman equation: free flow (Da > 1), porous flow (high permeability, 1 > Da > 10⁻⁶) and porous flow (low permeability Da < 10⁻⁶). In the present work the described bench mark problem is used to study the effects of solid walls for a range of low to high Darcy parameters. Both no-slip and slip conditions are considered and the results of these two cases are compared. The range of the applicability of the Brinkman equation and simulated results for different cases are shown.     

Influence of wall slip on the hydrodynamic behavior of a two-dimensional slider bearing

In the present paper, a multi-linearity method is used to address the nonlinear slip control equation for the hydrodynamic analysis of a two-dimensional (2-D) slip gap flow. Numerical analysis of a finite length slider bearing with wall slip shows that the surface limiting shear stress exerts complicated influences on the hydrodynamic behavior of the gap flow. If the slip occurs at either the stationary surface or the moving surface (especially at the stationary surface), there is a transition point in the initial limiting shear stress for the proportional coefficient to affect the hydrodynamic load support in two opposite ways: it increases the hydrodynamic load support at higher initial limiting shear stresses, but decreases the hydrodynamic load support at lower initial limiting shear stresses. If the slip occurs at the moving surface only, no fluid pressure is generated in the case of null initial limiting shear stress. If the slip occurs at both the surfaces with the same slip property, the hydrodynamic load support goes off after a critical sliding speed is reached. A small initial limiting shear stress and a small proportionality coefficient always give rise to a low friction drag. The project supported by the National Natural Science Foundation of China (10421002, 10332010), the National Basic Research Program of China (2006CB601205), and the Science Research Foundation of Liaoning Province (20052178). The English text was polished by Yunming Chen.     

Local characteristics of upward gas-liquid flows

The results of two-phase flow structure measurements in an upward gas-liquid flow in a 86.4 mm i.d. tube by the electrochemical and conductivity techniques are presented. Measurements were made in bubble and slug flow regimes at liquid flow rates ranging from 0.2 to 2 m/s.The flow instability and ambiguity in a bubble regime at low velocities is shown to exist. Great discrepancy between measured wall shear stress values and those predicted by the Lockhart-Martinelli model are due to the nonuniform distribution of gas over the tube cross section. Measurements of intensity of wall shear stress and liquid velocity fluctuations in a two-phase flow are presented.     

疏液壁面上亲液杂质对滑移特性的影响

纳米流动系统具有高效、经济等优势,在众多领域具有广泛的应用前景.因该类系统具有极高的表面积体积比,致使界面滑移效应对流动具有显著影响.论文采用非平衡分子动力学模拟方法,研究了疏液壁面表层混入少量亲液原子时纳米通道内液体的滑移特性,并基于分子动理论解释了其影响机制.研究发现,亲液杂质(均布或集中)对液体法向密度振荡程度影响较弱,但会显著改变壁面附近类固体层的分布和液体滑移规律;随亲液杂质占比增加,液体类固体现象更趋明显,壁面处液体接触密度也线性增大,但通道内液体平均速度逐渐降低,滑移长度也迅速减小;相比于集中的亲液杂质,均匀分布亲液杂质对滑移的弱化效应更强.如当亲液杂质占比为28%时,其滑移长度比单纯疏液表面的降低率从50%扩大至56%.基于分子动理论的分析发现,亲液杂质会导致杂质原子附近第一液体层内的原子发生跃迁的能垒加大,即弱化了液体原子的流向跃迁行为,从而降低了滑移量;相比于集中杂质,均匀分布的杂质还会降低固液间的非公度性,致使滑移特性破坏更严重.     

Molecular dynamics simulation of liquid argon flow at platinum surfaces

The micro Poiseuille flow for liquid argon flowing in a nanoscale channel formed by two solid walls was studied in the present paper. The solid wall material was selected as platinum, which has well established interaction potential. We consider the intermolecular force not only among the liquid argon molecules, but also between the liquid argon atoms and the solid wall particles, therefore three regions, i.e. the liquid argon computation domain, the top and bottom solid wall regions are included for the force interaction. The present MD (Molecular Dynamics) simulation was performed without any assumptions at the wall surface. The objective of the study is to find how the flow and the slip boundaries at the wall surface are affected by the applied gravity force, or the shear rate. The MD simulations are performed in a nondimensional unit system, with the periodic boundary conditions applied except in the channel height direction. Once the steady state is reached, the macroscopic parameters are evaluated using the statistical mechanics approach. For all the cases tested numerically in the present paper, slip boundaries occur, and such slip velocity at the stationary wall surface increases with increasing the applied gravity force, or the shear rate. The slip length, which is defined as the distance that the liquid particles shall travel beyond the wall surfaces to reach the same velocity as the wall surface, sharply decreases at small shear rate, then slightly decreases with increasing the applied shear rate. We observe that the liquid viscosity remains nearly constant at small shear rates, and the Newtonian flow occurs. However, with increasing the shear rate, the viscosity increases and the non-Newtonian flow appears.     

Unsteadiness of Flow Separation and End-Effects Regime in a Thrust-Optimized Contour Rocket Nozzle

Turbulent flow separation in over-expanded rocket nozzles is investigated experimentally in a sub-scale model nozzle fed with cold air and having a thrust-optimized contour. Depending upon the pressure ratio either a free shock separation (FSS) or a restricted shock separation (RSS) is observed with a significant hysteresis between these two flow regimes. It is shown that the RSS configuration may involve several separated regions. Analysis of wall pressure fluctuations give quantitative information on the fluctuating pressure field directly connected with the occurrence of significant side loads. Direct measurements of the evolution of the side loads with respect to the pressure ratio show the occurrence of three distinct peaks which are explained by the wall pressure fluctuations measurements.     

Stick-slip flow of high density polyethylene in a transparent slit die investigated by laser Doppler velocimetry

The oscillating flow instability of a molten linear high-density polyethylene is carefully studied using a single screw extruder equipped with a transparent slit die. Experiments are performed using laser Doppler velocimetry in order to obtain the local velocities field across the entire die width. At low flow rate, the extrusion is stable and steady state velocity profiles are obtained. During the instability, the velocity oscillates between two steady state limits, suggesting a periodic stick-slip transition mechanism. At high flow rate, the flow is mainly characterized by a pronounced wall slip. We show that wall slip occurs all along the die land. An investigation of the slip flow conditions shows that wall slip is not homogeneous in a cross section of the slit die, and that pure plug flow occurs only for very high flow rates. A numerical computation of the profile assuming wall slip boundary conditions is done to obtain the true local wall slip velocity. It confirms that slip velocities are of the same order of magnitude as those measured with a capillary rheometer.     

Slope effect on core-annular flow

To facilitate the flow of heavy viscous oils in a pipe, a water-lubricated transport is generally used. The water migrates into the regions of high shear at the pipe wall where it lubricates the flow. The pumping pressures are balanced by wall shear stresses in the water, the process therefore requires pressures comparable to pumping water alone, with no dependence on the viscosity of the oil. This means that significant savings in pumping power can be derived from this process, provided that it is well monitored. Indeed, the flow of a water/oil mixture in a pipe has two main characteristics. First, the fluids can adopt different spatial arrangements called flow regimes, and second, the presence of a water layer at the channel wall significantly reduces the global pressure drop. In this paper, an experimental investigation was performed on the effect of pipe slope and fluids flow rates on flow regimes, pressure drop and interfacial instability.     

Numerical study of aerosol particle deposition in simple and converging-diverging micro-channels with a slip boundary condition at the wall

The design of micro-devices involving aerosol transport requires the study of the deposition of aerosols in micro-channels. In this study, the slip and no-slip boundary conditions for the gas flow regime were applied to the Navier-Stokes equations to obtain the particle deposition in simple and converging-diverging micro-channels. The equation of particle motion included inertial, viscous, Brownian, and gravity terms. It was found that the ratio of gravity to inertial effects controls the deposition of particles with diameters of 0.1-1 μm, and the ratio of diffusion to inertial effects controls the deposition of particles with diameters of 0.01-0.001 i~m. Comparison between the no-slip and slip flow regimes showed that the deposition of 0.1- to 1-μm-diameter particles was less and the deposition of 0.01- to 0.001-1μm-diameter particles was greater for the slip flow regime. There was no significant difference between slip and no-slip flow regimes for the deposition of 0.01- to 0.1-μm-diameter particles. Finally, it was shown that the stagnated gas in the corners of the converging-diverging micro-channel produced similar gas velocity profiles under the slip and no-slip flow regimes.