The goal of this study is to analyse the steady flow of a Newtonian fluid mixed with spherical particles in a channel for the purpose of modelling proppant transport with gravitational settling in hydraulic fractures. The developments are based on a continuum constitutive model for a slurry, which is approximated by an empirical formula. It is shown that the problem under consideration features a two-dimensional flow and a boundary layer, which effectively introduces slip at the boundary and allows us to describe a transition from Poiseuille flow to Darcy’s law for high proppant concentrations. The expressions for both the outer (i.e. outside the boundary layer) and inner (i.e. within the boundary layer) solutions are obtained in terms of the particle concentration, particle velocity and fluid velocity. Unfortunately, these solutions require the numerical solution of an integral equation, and, as a result, the development of a proppant transport model for hydraulic fracturing based on these results is not practicable. To reduce the complexity of the problem, an approximate solution is introduced. To validate the use of this approximation, the error is estimated for different regimes of flow. The approximate solution is then used to calculate the expressions for the slurry flux and the proppant flux, which are the basis for a model that can be used to account for proppant transport with gravitational settling in a fully coupled hydraulic fracturing simulator.

We perform a local stability analysis of rotational flows in the presence of a constant vertical magnetic field and an azimuthal magnetic field with a general radial dependence. Employing the short-wavelength approximation we develop a unified framework for the investigation of the standard, helical and azimuthal version of the magnetorotational instability (MRI), as well as of current-driven kink-type instabilities. Considering the viscous and resistive setup, our main focus is on the case of small magnetic Prandtl numbers which applies e.g. to liquid-metal experiments but also to the colder parts of accretion disks. We show that the inductionless versions of MRI that were previously thought to be restricted to comparatively steep rotation profiles extend well to the Keplerian case if only the azimuthal field slightly deviates from its current-free (in the fluid) profile. We find an explicit criterion separating the pure azimuthal inductionless MRI from the regime where this instability is mixed with the Tayler instability. We further demonstrate that for particular parameter configurations the azimuthal MRI originates as a result of a dissipation-induced instability of Chandrasekhar’s equipartition solution of ideal magnetohydrodynamics.

We analyse the stability of a mushy layer during the directional solidification of a ternary alloy. Our model includes diffusive and convective transport of heat and solutes, coupled by an equilibrium thermodynamic constraint of the ternary phase diagram. The model contains phase change effects due to latent-heat release, solute rejection and background solidification. We identify novel convective instabilities, both direct and oscillatory, which are present under statically stable conditions. We examine these instabilities asymptotically by simplifying to a thin mushy layer with small growth rates. We also discuss numerical results for the full problem, confirming the asymptotic predictions and providing the stability characteristics outside the small-growth-rate approximation. A physical explanation for these instabilities in terms of parcel arguments is proposed, indicating that the instability mechanisms generally involve different rates of solute diffusion, different rates of solute rejection and different background solute distributions induced by the initial alloy composition.

The physics of the roughness sublayer are studied by direct numerical simulations (DNS) of an open-channel flow with sandgrain roughness. A double-averaging (DA) approach is used to separate the spatial variations of the time-averaged quantities and the turbulent fluctuations. The spatial inhomogeneity of velocity and Reynolds stresses results in an additional production term for the turbulent kinetic energy (TKE) – the ‘wake production’; it is the excess wake kinetic energy (WKE), generated from the work of mean flow against the form drag, that is not directly dissipated into heat, but instead converted into turbulence. The wake production promotes wall-normal turbulent fluctuations and increases the pressure work, which ultimately leads to more homogeneous turbulence in the roughness sublayer, and to the increase of Reynolds shear stress and the drag on the rough wall. In the fully rough regime, roughness directly affects the generation of the wall-normal fluctuations, while in the transitionally rough regime, the region affected by roughness is separated from the region of intense generation of these fluctuations. The budget of the WKE and the connection between the wake and the turbulence suggest strong interactions between the roughness sublayer and the outer layer that are insensitive to the variation of the outer-layer conditions. Furthermore, the present results may have implications for the relationship between the roughness geometry and the flow dynamics in the region directly affected by roughness.

Experimental observations of the flow of a suspension of solid fraction ${\it\phi}\approx 0.084$ over a circular cylindrical post in a shallow microchannel (depth smaller than the cylinder radius) find that the recirculating wake behind the obstacle at moderate Reynolds numbers is depleted or devoid of particles. Particles injected into the wake exit to regain the depleted state. By numerical simulation of the discrete particle motion, the basis for the depletion behind the cylinder is studied; rather than a shallow channel, the numerical simulations consider a periodic domain, mimicking the flow past an infinite cylinder. The Reynolds number is defined, using the average axial velocity ${\bar{U}}$, diameter of the obstacle $D$ and the kinematic viscosity of the suspension ${\it\nu}$, as $Re={\bar{U}}D/{\it\nu}$, and is studied for $Re<30$ in the simulation – conditions for which the pure fluid exhibits an extended steady closed-streamline (recirculating) wake behind the cylinder; unsteadiness is found to be suppressed by the channel walls in the experiments, allowing steady flow at a larger $Re$ than expected for an infinite cylinder (up to at least $Re=300$). The simulations use the lattice-Boltzmann method to determine the motion of the fluid and neutrally buoyant particles. The trajectory of a single particle (small relative to the cylinder) shows migration to a limit cycle inside the wake. With an increase of the number of particles in the wake alone (no particles in the free stream), particles can escape the wake due to velocity fluctuations. Simulation of the flow of suspensions of ${\it\phi}=0.04,0.06$ and 0.08 demonstrates that there is particle exchange between the wake and the free stream; the net flux of particles out of the wake leads to a particle-depleted wake, qualitatively very similar to the experimental observation.