Wall slip and multi-tier yielding in capillary suspensions

Impact of wall slip on the yield stress measurement is examined for capillary suspensions consisting of cocoa powder as the dispersed phase, vegetable oil as the continuous primary fluid, and water as the secondary fluid using smooth and serrated parallel plates. Using dynamic oscillatory measurements, we investigated the yielding behavior of this ternary solid-fluid-fluid system with varying particle volume fraction, ?, from 0.45 to 0.65 and varying water volume fraction, ?_w, from 0.02 to 0.08. Yield stress is defined as the maximum in the elastic stress (G^′γ), which is obtained by plotting the product of elastic modulus (G^′) and strain amplitude (γ) as a function of applied strain amplitude. With serrated plates, which offer minimal slippage, capillary suspensions with ? ≥?0.45 and a fixed ?_w =?0.06 showed a two-step yielding behavior as indicated by two peaks in the plots of elastic stress as a function of strain amplitude. On the other hand with smooth plates, the capillary suspensions showed strong evidence of wall slip as evident by the presence of three distinct peaks and lowered first yield stresses for all ? and ?_w. These results can be interpreted based on the fact that a particle-depleted layer, which is known to be responsible for slip, is present in the vicinity of the smooth surfaces. The slip layer presents itself as an additional “pseudo-microstructure” (characteristic length scale) besides the two microstructures, aqueous bridges and solid particle agglomerates, that may occur in the system. With serrated plates, both the yield stresses (σ₁σ₂) and storage moduli plateau at lower strain (before the first yield point) and at higher strain (before the second yield point) (G\(^{\prime }_{p1}\), G\(^{\prime }_{p2}\)) were found to increase with ? (at a fixed ?_w =?0.06) following power-law dependences. Similarly with increasing ?_w (0.02 – 0.08) at a fixed ? =?0.62, the system behaved as a solid-like material in a jammed state with particles strongly held together as manifested by rapidly increasing σ₁ and σ₂. The usage of smooth surfaces primarily affected σ₁ which was reflected by an approximately 70–90% decrement in the measured σ₁ for all values of ?. By contrast, σ₂ and G\(^{\prime }_{p2}\) were found to be unaffected as shown by close agreement of values obtained using serrated geometry due to vanishing slip layers at higher strain amplitudes.     

Direct Numerical Simulation and Theory of a Wall-Bounded Flow with Zero Skin Friction

We study turbulent plane Couette-Poiseuille (CP) flows in which the conditions (relative wall velocity ΔU _w ≡ 2U _w , pressure gradient dP/dx and viscosity ν) are adjusted to produce zero mean skin friction on one of the walls, denoted by APG for adverse pressure gradient. The other wall, FPG for favorable pressure gradient, provides the friction velocity u _τ , and h is the half-height of the channel. This leads to a one-parameter family of one-dimensional flows of varying Reynolds number Re ≡ U _w h/ν. We apply three codes, and cover three Reynolds numbers stepping by a factor of two each time. The agreement between codes is very good, and the Reynolds-number range is sizable. The theoretical questions revolve around Reynolds-number independence in both the core region (free of local viscous effects) and the two wall regions. The core region follows Townsend’s hypothesis of universal behavior for the velocity and shear stress, when they are normalized with u _τ and h; on the other hand universality is not observed for all the Reynolds stresses, any more than it is in Poiseuille flow or boundary layers. The FPG wall region obeys the classical law of the wall, again for velocity and shear stress. For the APG wall region, Stratford conjectured universal behavior when normalized with the pressure gradient, leading to a square-root law for the velocity. The literature, also covering other flows with zero skin friction, is ambiguous. Our results are very consistent with both of Stratford’s conjectures, suggesting that at least in this idealized flow turbulence theory is successful like it was for the classical logarithmic law of the wall. We appear to know the constants of the law within a 10% bracket. On the other hand, that again does not extend to Reynolds stresses other than the shear stress, but these stresses are passive in the momentum equation.     

Dual solutions of Casson fluid flows over a power-law stretching sheet

A steady boundary layer flow of a non-Newtonian Casson fluid over a power-law stretching sheet is investigated. A self-similar form of the governing equation is obtained, and numerical solutions are found for various values of the governing parameters. The solutions depend on the fluid material parameter. Dual solutions are obtained for some particular range of these parameters. The fluid velocity is found to decrease as the power-law stretching parameter β in the rheological Casson equation increases. At large values of β, the skin friction coefficient and the velocity profile across the boundary layer for the Casson fluid tend to those for the Newtonian fluid.     

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.     

Passive Boundary Layer Flow Control Using Porous Lamination

The flow over a porous laminated flat plate is investigated from a flow control perspective through experiments and computations. A square array of circular cylinders is used to model the porous lamination. We determine the velocities at the fluid–porous interface by solving the two-dimensional Navier–Stokes and the continuity equations using a staggered flow solver and using LDV in experiments. The control parameters for the porous region are porosity, \(\phi \) and Reynolds number, Re, based on the diameter of the circular cylinders used to model the porous lamination. Computations are conducted for \(0.4< \phi < 0.9\) and \(25< Re < 1000\), and the experiments are conducted for \(\phi = 0.65\) and 0.8 at \(Re \approx 391,\ 497\) and 803. The permeability of the porous lamination is observed to induce a slip velocity at the interface, effectively making it a slip wall. The slip velocity is seen to be increasing functions of \(\phi \) and Re. For higher porosities at higher Re, the slip velocity shows non-uniform and unsteady behavior and a breakdown Reynolds number is defined based on this characteristic. A map demarcating the two regimes of flow is drawn from the computational and experimental data. We observe that the boundary layer over the porous lamination is thinner than the Blasius boundary layer and the shear stress is higher at locations over the porous lamination. We note that the porous lamination helps maintain a favorable pressure gradient at the interface which delays separation. The suitable range of porosities for effective passive separation control is deduced from the results.     

The “Twin capillary” a simple device to separate shear- and slip-flow of fluids

When the flow behaviour of fluids is investigated with capillary-or rotational rheometers, adhesion of the fluid to the wall is normally one of the boundary conditions. For many fluids, especially for suspensions, this assumption is not valid. These fluids tend to slip at the wall. Therefore the normal evaluation of rheometer measurements leads to apparent but not compatible flow functions. The flow behaviour of these fluids can be characterized with two material functions which describe separately slipping in the boundary layer and shearing within the fluid. Only if both functions are known, correct predictions of flow processes are possible. A simple equipment to separate the shear function and the slip function is described.List of symbols Y^* apparent shear rate - Y _w ^* apparent wall shear rate - Y_w wall shear rate corrected with Rabinowitsch and Weissenberg correction - Y_s reduced shear rate (slip corrected) - Y_ws reduced wall shear rate (slip corrected) - * (r) velocity distribution in a capillary - _G slip velocity (at the wall) - * (r) velocity distribution in a capillary (without slip) - shear stress - _w wall shear stress - V_S total volume rate - V_G shear volume rate - V_G slip volume rate - p ₁ pressure in the reservoir channel of the capillary rheometer - p ₀ athmospheric pressure - L capillary length - R capillary radius     

Steady flow of a power-law non-Newtonian fluid across an unconfined square cylinder

A two-dimensional flow of a non-Newtonian power-law fluid directed normally to a horizontal cylinder with a square cross section is considered in the present paper. The problem is investigated numerically with a finite volume method by using the commercial code Ansys Fluent with a very large computational domain so that the flow could be considered unbounded. The investigation covers the power-law index from 0.1 to 2.0 and the Reynolds number range from 0.001 to 45.000. It is found that the drag coefficient for low Reynolds numbers and low power-law index (n ≤ 0.5) obeys the relationship C_D = A/Re. An equation for the quantity A as a function of the power-law index is derived. The drag coefficient becomes almost independent of the power-law index at high Reynolds numbers and the wake length changes nonlinearly with the Reynolds number and power-law index.     

Hydrodynamic properties of squirmer swimming in power-law fluid near a wall

Hydrodynamic properties of squirmer swimming in power-law fluid near a wall considering the interaction between squirmer and wall are numerically studied with an immersed boundary-lattice Boltzmann method. The power-law index, Reynolds number, initial orientation angle of squirmer, and initial distance of squirmer from the wall are all taken into account to investigate the swimming characteristics for pusher (β?<?0), neutral squirmer (β?=?0), and puller (β?>?0) (three kinds of swimmer types) near the no-slip boundary. Four new kinds of swimming modes are found. Results show that, for the pushers and pullers, the wall displays an increasing attraction with increasing power-law index n, which differs from the neutral squirmer who always departs from the wall after the first collision with the wall. Both the initial orientation angle and initial distance from the wall only affect the moving situations rather than the moving modes of the squirmers. However, the squirmers depart from the wall as the Reynolds number increases and chaotic orbits appear for some squirmers at Re?=?5. Several typical flow fields are analyzed and the power consumption and torque for different kinds of flows are also studied. It is found that, as the absolute value of β increases, the power consumption generally increases in shear-thinning (n?=?0.4), Newtonian (n?=?1), and shear-thickening (n?=?1.6) fluids. Moreover, the pushers (β?<?0) and the pullers (β?>?0) expend almost the same power if the absolute value of β remains the same. In addition, the power consumption of the squirmers is highly dependent on the power-law index n.     

Bulk flow parameters in turbulent wall boundary layers near separation

Experiments were conducted in a turbulent boundary layer near separation along a flat plate. The pressure gradient in flow direction was varied such that three significant boundary layer configurations could be maintained. The flow in the test section thus had simultaneously a region of favourable pressure gradient, a region of strong adverse pressure gradient with boundary layer separation and a region of reattached boundary layer. Specially designed fine probes facilitated the measurements of skin friction and velocity distribution very close to the wall. Bulk flow parameters such as skin friction coefficient C _f, Reynold's number Re_δ2 and shape factors H and G, which are significant characteristics of wall boundary layers were evaluated. The dependence of these parameters on the Reynolds number and along the test section was explored and the values were compared with other empirical and analytical formulae known in the literature.     

The Dean instability in power-law and Bingham fluids in a curved rectangular duct

The laminar flow of power-law and yield-stress fluids in 180° curved channels of rectangular cross section was studied experimentally and numerically in order to understand the effect of rheological fluid behavior on the Dean instability that appears beyond a critical condition in the flow. This leads to the apparition of Dean vortices that differ from the two corner vortices created by the channel wall curvature.Flow visualizations showed that the Dean vortices develop first in the near-wall zone on the concave (outer) wall, where the shear rate is higher and the viscosity weaker; then they penetrate into the centre of the channel cross section where power-law fluids have high viscosity and Bingham fluids are unyielded in laminar flow. Based on the complete formation on the concave wall of the new pairs of counter-rotating vortices (Dean vortices), the critical value of the Dean number decreases as the power-law index increases for the power-law fluids, and the Bingham number decreases for the Bingham fluids. For power-law fluids, a diagram of critical Dean numbers, based on the number of Dean vortices formed, was established for different axial positions. For the same flow conditions, the critical Dean number obtained using the axial velocity gradient criterion was smaller then that obtained with the visualization technique.     

Stability of Traveling Waves of Nonlinear Schrödinger Equation with Nonzero Condition at Infinity

We study effective elastic behavior of the incompatibly prestrained thin plates, where the prestrain is independent of thickness and uniform through the plate’s thickness h. We model such plates as three-dimensional elastic bodies with a prescribed pointwise stress-free state characterized by a Riemannian metric G, and seek the limiting behavior as \({h \to 0}\). We first establish that when the energy per volume scales as the second power of h, the resulting \({\Gamma}\) -limit is a Kirchhoff-type bending theory. We then show the somewhat surprising result that there exist non-immersible metrics G for whom the infimum energy (per volume) scales smaller than h². This implies that the minimizing sequence of deformations carries nontrivial residual three-dimensional energy but it has zero bending energy as seen from the limit Kirchhoff theory perspective. Another implication is that other asymptotic scenarios are valid in appropriate smaller scaling regimes of energy. We characterize the metrics G with the above property, showing that the zero bending energy in the Kirchhoff limit occurs if and only if the Riemann curvatures R₁₂₁₃, R₁₂₂₃ and R₁₂₁₂ of G vanish identically. We illustrate our findings with examples; of particular interest is an example where \({G_{2 \times 2}}\), the two-dimensional restriction of G, is flat but the plate still exhibits the energy scaling of the Föppl–von Kármán type. Finally, we apply these results to a model of nematic glass, including a characterization of the condition when the metric is immersible, for \({G = Id_{3} + \gamma n \otimes n}\) given in terms of the inhomogeneous unit director field distribution \({ n \in \mathbb{R}^3}\).     

Time-dependent drag reduction and ageing in aqueous solutions of a cationic surfactant

The turbulent pipe flow of a highly dilute aqueous cationic surfactant solution is investigated by means of a pulsed ultrasound Doppler method with special emphasis on the wall boundary layer. The velocity profiles are recorded for several Reynolds numbers at varying ages of the solution. The wall shear stress velocities u _τ used for the normalization of the velocity profiles are determined by fitting the measured profiles to the universal linear velocity profile in the viscous sublayer. The theoretical pressure loss is then calculated from the numerical values of u _τ and compared to the experimental values. Two different scaling methods are discussed for the velocity fluctuations concerning the correlation of the root-mean square values with the effect and the amount of drag reduction. It is shown that outer scaling with the mean velocity is appropriate for the detection of drag reduction in surfactant solutions, rather than inner scaling with the wall shear stress velocity, which is common practice in investigations of 'usual' turbulent flows.     

Spontaneous rotation in the exact solution of magnetohydrodynamic equations for flow between two stationary impermeable disks

Magnetohydrodynamic (MHD) flow of a viscous electrically conducting incompressible fluid between two stationary impermeable disks is considered. A homogeneous electric current density vector normal to the surface is specified on the upper disk, and the lower disk is nonconducting. The exact von Karman solution of the complete system of MHD equations is studied in which the axial velocity and the magnetic field depend only on the axial coordinate. The problem contains two dimensionless parameters: the electric current density on the upper plate Y and the Batchelor number (magnetic Prandtl number). It is assumed that there is no external source that produces an axial magnetic field. The problem is solved for a Batchelor number of 0–2. Fluid flow is caused by the electric current. It is shown that for small values of Y, the fluid velocity vector has only axial and radial components. The velocity of motion increases with increasing Y, and at a critical value of Y, there is a bifurcation of the new steady flow regime with fluid rotation, while the flow without rotation becomes unstable. A feature of the obtained new exact solution is the absence of an axial magnetic field necessary for the occurrence of an azimuthal component of the ponderomotive force, as is the case in the MHD dynamo. A new mechanism for the bifurcation of rotation in MHD flow is found.     

Effects of Knudsen number and geometry on gaseous flow and heat transfer in a constricted microchannel

A flow and heat transfer numerical simulation is performed for a 2D laminar incompressible gas flow through a constricted microchannel in the slip regime with constant wall temperature. The effects of rarefaction, creeping flow, first order slip boundary conditions and hydrodynamically/thermally developing flow are assumed. The effects of Knudsen number and geometry on thermal and hydrodynamic characteristics of flow in a constricted microchannel are explored. SIMPLE algorithm in curvilinear coordinate is used to solve the governing equations including continuity, energy and momentum with the temperature jump and velocity slip conditions at the solid walls in discretized form. The resulting velocity and temperature profiles are then utilized to obtain the microchannel C _f Re and Nusselt number as a function of Knudsen number and geometry. The results show that Knudsen number has declining effect on the C _f Re and Nusselt number in the constricted microchannel. In addition, the temperature jump on wall and slip velocity increase with increasing Knudsen number. Moreover, by decreasing the throttle area, the fluid flow characteristics experience more intense variations in the constricted region. To verify the code a comparison is carried out with available results and good agreement is achieved.     

Experimental Control of Turbulent Boundary Layers with In-plane Travelling Waves

The experimental control of turbulent boundary layers using streamwise travelling waves of spanwise wall velocity, produced using a novel active surface, is outlined in this paper. The innovative surface comprises a pneumatically actuated compliant structure based on the kagome lattice geometry, supporting a pre-tensioned membrane skin. Careful design of the structure enables waves of variable length and speed to be produced in the flat surface in a robust and repeatable way, at frequencies and amplitudes known to have a favourable influence on the boundary layer. Two surfaces were developed, a preliminary module extending 152 mm in the streamwise direction, and a longer one with a fetch of 2.9 m so that the boundary layer can adjust to the new surface condition imposed by the forcing. With a shorter, 1.5 m portion of the surface actuated, generating an upstream-travelling wave, a drag reduction of 21.5% was recorded in the boundary layer with Re _τ =?1125. At the same flow conditions, a downstream-travelling produced a much smaller drag reduction of 2.6%, agreeing with the observed trends in current simulations. The drag reduction was determined with constant temperature hot-wire measurements of the mean velocity gradient in the viscous sublayer, while simultaneous laser Doppler vibrometer measurements of the surface recorded the wall motion. Despite the mechanics of the dynamic surface resulting in some out-of-plane motion (which is small in comparison to the in-plane streamwise movement), the positive drag reduction results are encouraging for future investigations at higher Reynolds numbers.     

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.     

Further remarks on numerical investigation on gas displacement of a shear-thinning liquid and a visco-plastic material in capillary tubes

In this note, we present some results concerning the gas displacement of power-law liquids and visco-plastic Papanastasiou’s materials improving the understanding of the problem considered in Sousa et al. (2007) [1]. Specifically, we present: the fraction of mass attached to the wall for a viscosity-thickening power-law fluid, different transition patterns between by pass and full-recirculating flow regimes, expressions for the critical (in the sense proposed by Soares and Thompson (2009) [2]) fraction of residual mass as a function of the rheological parameter of interest, and fields of yielded and unyielded zones for the visco-plastic material.     

Effects of slip velocity and surface acoustic wave on the laminar flow stability

The stability of slip flows when a surface acoustic wave (SAW) propagating along the walls of a microchannel in the laminar flow regime is investigated. The governing equation which was derived by considering the weakly nonlinear coupling between the deformable wall and streaming slip flow is linearized and then the eigenvalue problem is solved by a numerical code together with the associated interface and slip velocity boundary conditions. The value of the critical Reynolds number was found to be near 1,441 for a Knudsen number being 0.001 (associated with a physical parameter K ₀ characterizing the SAW effect) which is much smaller than the static-wall case for conventional pressure-driven flows.