Stochastically driven instability in rotating shear flows

Origin of hydrodynamic turbulence in rotating shear flows is investigated. The particular emphasis is on flows whose angular velocities decrease but specific angular momenta increase with increasing radial coordinate. Such flows are Rayleigh stable, but must be turbulent in order to explain observed data. Such a mismatch between the linear theory and observations/experiments is more severe when any hydromagnetic/magnetohydrodynamic instability and the corresponding turbulence therein is ruled out. The present work explores the effect of stochastic noise on such hydrodynamic flows. We focus on a small section of such a flow which is essentially a plane shear flow supplemented by the Coriolis effect. This also mimics a small section of an astrophysical accretion disk. It is found that such stochastically driven flows exhibit large temporal and spatial correlations of perturbation velocities, and hence large energy dissipations, that presumably generate instability. A range of angular velocity profiles (for the steady flow), starting with the constant angular momentum to that of the constant circular velocity are explored. It is shown that the growth and roughness exponents calculated from the contour (envelope) of the perturbed flows are all identical, revealing a unique universality class for the stochastically forced hydrodynamics of rotating shear flows. This work, to the best of our knowledge, is the first attempt to understand origin of instability and turbulence in the three-dimensional Rayleigh stable rotating shear flows by introducing additive stochastic noise to the underlying linearized governing equations. This has important implications in resolving the turbulence problem in astrophysical hydrodynamic flows such as accretion disks.

Shapes and Dynamics of Blood Cells in Poiseuille and Shear Flows

The dynamics, shape, deformation, and orientation of red blood cells in microcirculation affect the rheology, flow resistance and transport properties of whole blood. This leads to important correlations of cellular and continuum scales. Furthermore, the dynamics of RBCs subject to different flow conditions and vessel geometries is relevant for both fundamental research and biomedical applications (e.g drug delivery). In this thesis, the behaviour of RBCs is investigated for different flow conditions via computer simulations. We use a combination of two mesoscopic particle-based simulation techniques, dissipative particle dynamics and smoothed dissipative particle dynamics. We focus on the microcapillary scale of several μm. At this scale, blood cannot be considered at the continuum but has to be studied at the cellular level. The connection between cellular motion and overall blood rheology will be investigated. Red blood cells are modelled as viscoelastic objects interacting hydrodynamically with a viscous fluid environment. The properties of the membrane, such as resistance against bending or shearing, are set to correspond to experimental values. Furthermore, thermal fluctuations are considered via random forces. Analyses corresponding to light scattering measurements are performed in order to compare to experiments and suggest for which situations this method is suitable. Static light scattering by red blood cells characterises their shape and allows comparison to objects such as spheres or cylinders, whose scattering signals have analytical solutions, in contrast to those of red blood cells. Dynamic light scattering by red blood cells is studied concerning its suitability to detect and analyse motion, deformation and membrane fluctuations. Dynamic light scattering analysis is performed for both diffusing and flowing cells. We find that scattering signals depend on various cell properties, thus allowing to distinguish different cells. The scattering of diffusing cells allows to draw conclusions on their bending rigidity via the effective diffusion coefficient. The scattering of flowing cells allows to draw conclusions on the shear rate via the scattering amplitude correlation. In flow, a RBC shows different shapes and dynamic states, depending on conditions such as confinement, physiological/pathological state and cell age. Here, two essential flow conditions are studied: simple shear flow and tube flow. Simple shear flow as a basic flow condition is part of any more complex flow. The velocity profile is linear and shear stress is homogeneous. In simple shear flow, we find a sequence of different cell shapes by increasing the shear rate. With increasing shear rate, we find rolling cells with cup shapes, trilobe shapes and quadrulobe shapes. This agrees with recent experiments. Furthermore, the impact of the initial orientation on the dynamics is studied. To study crowding and collective effects, systems with higher haematocrit are set up. Tube flow is an idealised model for the flow through cylindric microvessels. Without cell, a parabolic flow profile prevails. A single red blood cell is placed into the tube and subject to a Poiseuille profile. In tube flow, we find different cell shapes and dynamics depending on confinement, shear rate and cell properties. For strong confinements and high shear rates, we find parachute-like shapes. Although not perfectly symmetric, they are adjusted to the flow profile and maintain a stationary shape and orientation. For weak confinements and low shear rates, we find tumbling slippers that rotate and moderately change their shape. For weak confinements and high shear rates, we find tank-treading slippers that oscillate in a limited range of inclination angles and strongly change their shape. For the lowest shear rates, we find cells performing a snaking motion. Due to cell properties and resultant deformations, all shapes differ from hitherto descriptions, such as steady tank-treading or symmetric parachutes. We introduce phase diagrams to identify flow regimes for the different shapes and dynamics. Changing cell properties, the regime borders in the phase diagrams change. In both flow types, both the viscosity contrast and the choice of stress-free shape are important. For in vitro experiments, the solvent viscosity has often been higher than the cytosol viscosity, leading to a different pattern of dynamics, such as steady tank-treading. The stress-free state of a RBC, which is the state at zero shear stress, is still controversial, and computer simulations enable direct comparisons of possible candidates in equivalent flow conditions.

Magnetic and Acoustic Investigations of Turbulent Spherical Couette Flow

Title of dissertation: MAGNETIC AND ACOUSTIC INVESTIGATIONS OF TURBULENT SPHERICAL COUETTE FLOW Matthew M. Adams, Doctor of Philosophy, 2016 Dissertation directed by: Professor Daniel Lathrop Department of Physics This dissertation describes experiments in spherical Couette devices, using both gas and liquid sodium. The experimental geometry is motivated by the Earth's outer core, the seat of the geodynamo, and consists of an outer spherical shell and an inner sphere, both of which can be rotated independently to drive a shear flow in the fluid lying between them. In the case of experiments with liquid sodium, we apply DC axial magnetic fields, with a dominant dipole or quadrupole component, to the system. We measure the magnetic field induced by the flow of liquid sodium using an external array of Hall effect magnetic field probes, as well as two probes inserted into the fluid volume. This gives information about possible velocity patterns present, and we extend previous work categorizing flow states, noting further information that can be extracted from the induced field measurements. The limitations due to a lack of direct velocity measurements prompted us to work on developing the technique of using acoustic modes to measure zonal flows. Using gas as the working fluid in our 60~cm diameter spherical Couette experiment, we identified acoustic modes of the container, and obtained excellent agreement with theoretical predictions. For the case of uniform rotation of the system, we compared the acoustic mode frequency splittings with theoretical predictions for solid body flow, and obtained excellent agreement. This gave us confidence in extending this work to the case of differential rotation, with a turbulent flow state. Using the measured splittings for this case, our colleagues performed an inversion to infer the pattern of zonal velocities within the flow, the first such inversion in a rotating laboratory experiment. This technique holds promise for use in liquid sodium experiments, for which zonal flow measurements have historically been challenging.

Functional heterogeneity of the UpaH autotransporter protein from Uropathogenic Escherichia coli

Uropathogenic Escherichia coli (UPEC) is responsible for the majority of urinary tract infections (UTI). To cause a UTI, UPEC must adhere to the epithelial cells of the urinary tract and overcome the shear flow forces of urine. This function is mediated primarily by fimbrial adhesins, which mediate specific attachment to host cell receptors. Another group of adhesins that contributes to UPEC-mediated UTI is autotransporter (AT) proteins. AT proteins possess a range of virulence properties, such as adherence, aggregation, invasion, and biofilm formation. One recently characterized AT protein of UPEC is UpaH, a large adhesin-involved-in-diffuse-adherence (AIDA-I)-type AT protein that contributes to biofilm formation and bladder colonization. In this study we characterized a series of naturally occurring variants of UpaH. We demonstrate that extensive sequence variation exists within the passenger-encoding domain of UpaH variants from different UPEC strains. This sequence variation is associated with functional heterogeneity with respect to the ability of UpaH to mediate biofilm formation. In contrast, all of the UpaH variants examined retained a conserved ability to mediate binding to extracellular matrix (ECM) proteins. Bioinformatic analysis of the UpaH passenger domain identified a conserved region (UpaHCR) and a hydrophobic region (UpaHHR). Deletion of these domains reduced biofilm formation but not the binding to ECM proteins. Despite variation in the upaH sequence, the transcription of upaH was repressed by a conserved mechanism involving the global regulator H-NS, and mutation of the hns gene relieved this repression. Overall, our findings shed new light on the regulation and functions of the UpaH AT protein.

Deformation of a single red blood cell in a microvessel

Red blood cells (RBCs) are the most common type of cells in human blood and they exhibit different types of motions and deformed shapes in capillary flows. The behaviour of the RBCs should be studied in order to explain the RBC motion and deformation mechanism. This article presents a numerical simulation method for RBC deformation in microvessels. A two dimensional spring network model is used to represent the RBC membrane, where the elastic stretch/compression energy and the bending energy are considered with the constraint of constant RBC surface area. The forces acting on the RBC membrane are obtained from the principle of virtual work. The whole fluid domain is discretized into a finite number of particles using smoothed particle hydrodynamics concepts and the motions of all the particles are solved using Navier--Stokes equations. Minimum energy concepts are used to simulate the deformed shape of the RBC model. To verify the model, the motion of a single RBC is simulated in a Poiseuille flow and the characteristic parachute shape of the RBC is observed. Further simulations reveal that the RBC shows a tank treading motion when it flows in a linear shear flow.

Decompression melting driving intraplate volcanism in Australia: Evidence from magnetotelluric sounding

A long-period magnetotelluric (MT) survey, with 39 sites covering an area of 270 by 150 km, has identified melt within the thinned lithosphere of Pleistocene-Holocene Newer Volcanics Province (NVP) in southeast Australia, which has been variously attributed to mantle plume activity or edge-driven mantle convection. Two-dimensional inversions from the MT array imaged a low-resistivity anomaly (10-30Ωm) beneath the NVP at ∼40-80 km depth, which is consistent with the presence of ∼1.5-4% partial melt in the lithosphere, but inconsistent with elevated iron content, metasomatism products or a hot spot. The conductive zone is located within thin juvenile oceanic mantle lithosphere, which was accreted onto thicker Proterozoic continental mantle lithosphere. We propose that the NVP owes its origin to decompression melting within the asthenosphere, promoted by lithospheric thickness variations in conjunction with rapid shear, where asthenospheric material is drawn by shear flow at a "step" at the base of the lithosphere.

Behaviour and strength of hollow flange channel sections under torsion and bending

Hollow flange channel section is a cold-formed high-strength and thin-walled steel section with a unique shape including two rectangular hollow flanges and a slender web. Due to its mono-symmetric characteristics, it will also be subjected to torsion when subjected to transverse loads in practical applications. Past research on steel beams subject to torsion has concentrated on open sections while very few steel design standards give suitable design rules for torsion design. Since the hollow flange channel section is different from conventional open sections, its torsional behaviour remains unknown to researchers. Therefore the elastic behaviour of hollow flange channel sections subject to uniform and non-uniform torsion, and combined torsion and bending was investigated using the solutions of appropriate differential equilibrium equations. The section torsion shear flow, warping normal stress distribution, and section constants including torsion constant and warping constant were obtained. The results were compared with those from finite element analyses that verified the accuracy of analytical solutions. Parametric studies were undertaken for simply supported beams subject to a uniformly distributed torque and a uniformly distributed transverse load applied away from the shear centre. This paper presents the details of this research into the elastic behaviour and strength of hollow flange channel sections subject to torsion and bending and the results.

Dynamics of a dilute sheared inelastic fluid. I. Hydrodynamic modes and velocity correlation functions

The hydrodynamic modes and the velocity autocorrelation functions for a dilute sheared inelastic fluid are analyzed using an expansion in the parameter epsilon=(1-e)(1/2), where e is the coefficient of restitution. It is shown that the hydrodynamic modes for a sheared inelastic fluid are very different from those for an elastic fluid in the long-wave limit, since energy is not a conserved variable when the wavelength of perturbations is larger than the ``conduction length.'' In an inelastic fluid under shear, there are three coupled modes, the mass and the momenta in the plane of shear, which have a decay rate proportional to k(2/3) in the limit k -> 0, if the wave vector has a component along the flow direction. When the wave vector is aligned along the gradient-vorticity plane, we find that the scaling of the growth rate is similar to that for an elastic fluid. The Fourier transforms of the velocity autocorrelation functions are calculated for a steady shear flow correct to leading order in an expansion in epsilon. The time dependence of the autocorrelation function in the long-time limit is obtained by estimating the integral of the Fourier transform over wave number space. It is found that the autocorrelation functions for the velocity in the flow and gradient directions decay proportional to t(-5/2) in two dimensions and t(-15/4) in three dimensions. In the vorticity direction, the decay of the autocorrelation function is proportional to t(-3) in two dimensions and t(-7/2) in three dimensions.

Structure and dynamics of two-dimensional sheared granular flows

The structure and dynamics of the two-dimensional linear shear flow of inelastic disks at high area fractions are analyzed. The event-driven simulation technique is used in the hard-particle limit, where the particles interact through instantaneous collisions. The structure (relative arrangement of particles) is analyzed using the bond-orientational order parameter. It is found that the shear flow reduces the order in the system, and the order parameter in a shear flow is lower than that in a collection of elastic hard disks at equilibrium. The distribution of relative velocities between colliding particles is analyzed. The relative velocity distribution undergoes a transition from a Gaussian distribution for nearly elastic particles, to an exponential distribution at low coefficients of restitution. However, the single-particle distribution function is close to a Gaussian in the dense limit, indicating that correlations between colliding particles have a strong influence on the relative velocity distribution. This results in a much lower dissipation rate than that predicted using the molecular chaos assumption, where the velocities of colliding particles are considered to be uncorrelated.

The machinability of a cast aluminium alloy-graphite particle composite

The microstructure of a cast Al---Si alloy-graphite particle composite is examined using optical and analytical scanning electron microscopy. Specimens containing different percentages of graphite were machined by orthogonal planning with 25° and 45° rake angle tools at both 6.5 and 13.2 m min−1. The machining forces are reported and the chip-rake-face friction coefficients and shear flow stresses are calculated. It is shown that the reduction in machining forces with increasing graphite content is due mostly to a decrease in the shear flow stress rather than to lower chip-rake-face friction. Both the polished and the machined surfaces of the composite are rougher than those of the simple alloy, apparently owing to the greater porosity, the tearing out of graphite particles, or the opening of cracks at the graphite particles in the wake of the tool.

Dynamics of dense sheared granular flows. Part 1. Structure and diffusion

Shear flows of inelastic spheres in three dimensions in the Volume fraction range 0.4-0.64 are analysed using event-driven simulations.Particle interactions are considered to be due to instantaneous binary collisions, and the collision model has a normal coefficient of restitution e(n) (negative of the ratio of the post- and pre-collisional relative velocities of the particles along the line joining the centres) and a tangential coefficient of restitution e(t) (negative of the ratio of post- and pre-collisional velocities perpendicular to the line Joining the centres). Here, we have considered both e(t) = +1 and e(t) = e(n) (rough particles) and e(t) =-1 (smooth particles), and the normal coefficient of restitution e(n) was varied in the range 0.6-0.98. Care was taken to avoid inelastic collapse and ensure there are no particle overlaps during the simulation. First, we studied the ordering in the system by examining the icosahedral order parameter Q(6) in three dimensions and the planar order parameter q(6) in the plane perpendicular to the gradient direction. It was found that for shear flows of sufficiently large size, the system Continues to be in the random state, with Q(6) and q(6) close to 0, even for volume fractions between phi = 0.5 and phi = 0.6; in contrast, for a system of elastic particles in the absence of shear, the system orders (crystallizes) at phi = 0.49. This indicates that the shear flow prevents ordering in a system of sufficiently large size. In a shear flow of inelastic particles, the strain rate and the temperature are related through the energy balance equation, and all time scales can be non-dimensionalized by the inverse of the strain rate. Therefore, the dynamics of the system are determined only by the volume fraction and the coefficients of restitution. The variation of the collision frequency with volume fraction and coefficient of estitution was examined. It was found, by plotting the inverse of the collision frequency as a function of volume fraction, that the collision frequency at constant strain rate diverges at a volume fraction phi(ad) (volume fraction for arrested dynamics) which is lower than the random close-packing Volume fraction 0.64 in the absence of shear. The volume fraction phi(ad) decreases as the coefficient of restitution is decreased from e(n) = 1; phi(ad) has a minimum of about 0.585 for coefficient of restitution e(n) in the range 0.6-0.8 for rough particles and is slightly larger for smooth particles. It is found that the dissipation rate and all components of the stress diverge proportional to the collision frequency in the close-packing limit. The qualitative behaviour of the increase in the stress and dissipation rate are well Captured by results derived from kinetic theory, but the quantitative agreement is lacking even if the collision frequency obtained from simulations is used to calculate the pair correlation function used In the theory.

Dynamics of dense sheared granular flows. Part II. The relative velocity distributions

The distribution of relative velocities between colliding particles in shear flows of inelastic spheres is analysed in the Volume fraction range 0.4-0.64. Particle interactions are considered to be due to instantaneous binary collisions, and the collision model has a normal coefficient of restitution e(n) (negative of the ratio of the post- and pre-collisional relative velocities of the particles along the line joining the centres) and a tangential coefficient of restitution e(t) (negative of the ratio of post- and pre-collisional velocities perpendicular to line joining the centres). The distribution or pre-collisional normal relative velocities (along the line Joining the centres of the particles) is Found to be an exponential distribution for particles with low normal coefficient of restitution in the range 0.6-0.7. This is in contrast to the Gaussian distribution for the normal relative velocity in all elastic fluid in the absence of shear. A composite distribution function, which consists of an exponential and a Gaussian component, is proposed to span the range of inelasticities considered here. In the case of roughd particles, the relative velocity tangential to the surfaces at contact is also evaluated, and it is found to be close to a Gaussian distribution even for highly inelastic particles.Empirical relations are formulated for the relative velocity distribution. These are used to calculate the collisional contributions to the pressure, shear stress and the energy dissipation rate in a shear flow. The results of the calculation were round to be in quantitative agreement with simulation results, even for low coefficients of restitution for which the predictions obtained using the Enskog approximation are in error by an order of magnitude. The results are also applied to the flow down an inclined plane, to predict the angle of repose and the variation of the volume fraction with angle of inclination. These results are also found to be in quantitative agreement with previous simulations.

A rationalisation of fatigue thresholds in pearlitic steels using a theoretical model

From a detailed re-examination of results in the literature, the effects of microstructure sizes, namely interlamellar spacing, pearlitic colony size and the prior austentitic grain size on the thresholds for fatigue crack growth (ΔKth) and crack closure (Kcl, th) have been illustrated. It is shown that while interlamellar spacing explicitly controls yield strength, a similar effect on ΔKth cannot be expected. On the other hand, the pearlitic colony size is shown to strongly influence ΔKth and Kcl, th through the deflection and retardation of cracks at colony boundaries. Consequently, an increase in ΔKth and Kcl, th with colony size has been found. The development of a theoretical model to illustrate the effects of colony size, shear flow stress in the slip band and macroscopic yield strength on Kcl, th and ΔKth is presented. the model assumes colony boundaries as potential sites for slip band pile-up formation and subsequent crack deflection finally leading to zig-zag crack growth. Using the concepts of roughness induced crack closure, the magnitude of Kcl, th is quantified as a function of colony size. In deriving the model, the flow stress in the slip band has been considered to represent the work hardened state in pearlite. Comparison of the theoretically predicted trend with the experimental data demonstrates very good agreement. Further, the intrinsic or closure free component of the fatigue threshold, ΔKeff, th is found to be insensitive to colony size and interlamellar spacing. Using a criterion for intrinsic fatigue threshold which considers the attainment of a critical fracture stress over a characteristic distance corresponding to interlamellar spacing, ΔKth values at high R values can be estimated with reasonable accuracy. The magnitude of ΔKth as a function of colony size is then obtained by summing up the average value of experimentally obtained ΔKeff, th values and the predicted Kcl, th values as a function of colony size. Again, very good agreement of the theoretically predicted ΔKth values with those experimentally obtained has been demonstrated.