No Arabic abstract
The seminal Batchelor-Greens (BG) theory on the hydrodynamic interaction of two spherical particles of radii a suspended in a viscous shear flow neglects the effect of the boundaries. In the present paper we study how a plane wall modifies this interaction. Using an integral equation for the surface traction we derive the expression for the particles relative velocity as a sum of the BGs velocity and the term due to the presence of a wall at finite distance, z_0. Our calculation is not the perturbation theory of the BG solution, so the contribution due to the wall is not necessarily small. The distance at which the wall significantly alters the particles interaction scales as z_0^{3/5}. The phase portrait of the particles relative motion is different from the BG theory, where there are two singly-connected regions of open and closed trajectories both of infinite volume. For finite z_0, there is a new domain of closed (dancing) and open (swapping) trajectories. The width of this region behaves as 1/z_0. Along the swapping trajectories, that have been previously observed numerically, the incoming particle is turning back after the encounter with the reference particle, rather than passing it by, as in the BG theory. The region of dancing trajectories has infinite volume. We found a one-parameter family of equilibrium states, overlooked previously, whereas the pair of spheres flows as a whole without changing its configuration. These states are marginally stable and their perturbation yields a two-parameter family of the dancing trajectories, where the particle is orbiting around a fixed point in a frame co-moving with the reference particle. We suggest that the phase portrait obtained at z_0>>a is topologically stable and can be extended down to rather small z_0 of several particle diameters. We confirm this by direct numerical simulations of the Navier-Stokes equations with z_0=5a.
We numerically investigate the effect of entrance condition on the spatial and temporal evolution of multiple three-dimensional vortex pairs and wall shear stress distribution in a curved artery model. We perform this study using a Newtonian blood-analog fluid subjected to a pulsatile flow with two inflow conditions. The first flow condition is fully developed while the second condition is undeveloped (i.e. uniform). We discuss the connection along the axial direction between regions of organized vorticity observed at various cross-sections of the model and compare results between the different entrance conditions. We model a human artery with a simple, rigid $180^circ$ curved pipe with circular cross-section and constant curvature, neglecting effects of taper, torsion and elasticity. Numerical results are computed from a discontinuous high-order spectral element flow solver. The flow rate used in this study is physiological. We observe differences in secondary flow patterns, especially during the deceleration phase of the physiological waveform where multiple vortical structures of both Dean-type and Lyne-type coexist. We highlight the effect of the entrance condition on the formation of these structures and subsequent appearance of abnormal inner wall shear stresses - a potentially significant correlation since cardiovascular disease is known to progress along the inner wall of curved arteries under varying degrees of flow development.
A comprehensive study of the effect of wall heating or cooling on the linear, transient and secondary growth of instability in channel flow is conducted. The effect of viscosity stratification, heat diffusivity and of buoyancy are estimated separately, with some unexpected results. From linear stability results, it has been accepted that heat diffusivity does not affect stability. However, we show that realistic Prandtl numbers cause a transient growth of disturbances that is an order of magnitude higher than at zero Prandtl number. Buoyancy, even at fairly low levels, gives rise to high levels of subcritical energy growth. Unusually for transient growth, both of these are spanwise-independent and not in the form of streamwise vortices. At moderate Grashof numbers, exponential growth dominates, with distinct Rayleigh-Benard and Poiseuille modes for Grashof numbers upto $sim 25000$, which merge thereafter. Wall heating has a converse effect on the secondary instability compared to the primary, destabilising significantly when viscosity decreases towards the wall. It is hoped that the work will motivate experimental and numerical efforts to understand the role of wall heating in the control of channel and pipe flows.
This paper presents a method for calculating the wall shear rate in pipe turbulent flow. It collapses adequately the data measured in laminar flow and turbulent flow into a single flow curve and gives the basis for the design of turbulent flow viscometers. Key words: non-Newtonian, wall shear rate, turbulent, rheometer
In this study, we combine experiments and direct numerical simulations to investigate the effects of the height of transverse ribs at the walls on both global and local flow properties in turbulent Taylor-Couette flow. We create rib roughness by attaching up to 6 axial obstacles to the surfaces of the cylinders over an extensive range of rib heights, up to blockages of 25% of the gap width. In the asymptotic ultimate regime, where the transport is independent of viscosity, we emperically find that the prefactor of the $Nu_{omega} propto Ta^{1/2}$ scaling (corresponding to the drag coefficient $C_f(Re)$ being constant) scales with the number of ribs $N_r$ and by the rib height $h^{1.71}$. The physical mechanism behind this is that the dominant contribution to the torque originates from the pressure forces acting on the rib which scale with rib height. The measured scaling relation of $N_r h^{1.71}$ is slightly smaller than the expected $N_r h^2$ scaling, presumably because the ribs cannot be regarded as completely isolated but interact. In the counter-rotating regime with smooth walls, the momentum transport is increased by turbulent Taylor vortices. We find that also in the presence of transverse ribs these vortices persist. In the counter-rotating regime, even for large roughness heights, the momentum transport is enhanced by these vortices.
Driven or active suspensions can display fascinating collective behavior, where coherent motions or structures arise on a scale much larger than that of the constituent particles. Here, we report experiments and numerical simulations revealing that red blood cells (RBCs) assemble into regular patterns in a confined shear flow. The order is of pure hydrodynamic and inertialess origin, and emerges from a subtle interplay between (i) hydrodynamic repulsion by the bounding walls which drives deformable cells towards the channel mid-plane and (ii) intercellular hydrodynamic interactions which can be attractive or repulsive depending on cell-cell separation. Various crystal-like structures arise depending on RBC concentration and confinement. Hardened RBCs in experiments and rigid particles in simulations remain disordered under the same conditions where deformable RBCs form regular patterns, highlighting the intimate link between particle deformability and the emergence of order. The difference in structuring ability of healthy (deformable) and diseased (stiff) RBCs creates a flow signature potentially exploitable for diagnosis of blood pathologies.