Search

Laboratory experiments and particle-resolved simulations are employed to investigate the settling dynamics of a pair of rigidly connected spherical particles of unequal size. They yield a detailed picture of the transient evolution and the terminal values of the aggregate's orientation angle and its settling and drift velocities as functions of the aspect ratio and the Galileo number $Ga$ , which denotes the ratio of buoyancy and viscous forces acting on the aggregate. At low to moderate values of $Ga$ , the aggregate's orientation and velocity converge to their terminal values monotonically, whereas for higher $Ga$ -values the aggregate tends to undergo a more complex motion. If the aggregate assumes an asymmetric terminal orientation, it displays a non-zero terminal drift velocity. For diameter ratios much larger than one and small $Ga$ , the terminal orientation of the aggregate becomes approximately vertical, whereas when $Ga$ is sufficiently large for flow separation to occur, the aggregate orients itself such that the smaller sphere is located at the separation line. Empirical scaling laws are obtained for the terminal settling velocity and orientation angle as functions of the aspect ratio and $Ga$ for diameter ratios from 1 to 4 and particle-to-fluid density ratios from 1.3 to 5. An analysis of the accompanying flow field shows the formation of vortical structures exhibiting complex topologies in the aggregate's wake, and indicates the formation of a horizontal pressure gradient across the larger sphere, which represents the main reason for the emergence of the drift velocity.

Abstract:

Abstract:

Abstract:

We develop a physical and computational model for performing fully coupled, grain-resolved direct numerical simulations of cohesive sediment, based on the immersed boundary method. The model distributes the cohesive forces over a thin shell surrounding each particle, thereby allowing for the spatial and temporal resolution of the cohesive forces during particle–particle interactions. The influence of the cohesive forces is captured by a single dimensionless parameter in the form of a cohesion number, which represents the ratio of cohesive and gravitational forces acting on a particle. We test and validate the cohesive force model for binary particle interactions in the drafting–kissing–tumbling (DKT) configuration. Cohesive sediment grains can remain attached to each other during the tumbling phase following the initial collision, thereby giving rise to the formation of flocs. The DKT simulations demonstrate that cohesive particle pairs settle in a preferred orientation, with particles of very different sizes preferentially aligning themselves in the vertical direction, so that the smaller particle is drafted in the wake of the larger one. This preferred orientation of cohesive particle pairs is found to remain influential for systems of higher complexity. To this end, we perform large simulations of 1261 polydisperse settling particles starting from rest. These simulations reproduce several earlier experimental observations by other authors, such as the accelerated settling of sand and silt particles due to particle bonding, the stratification of cohesive sediment deposits, and the consolidation process of the deposit. They identify three characteristic phases of the polydisperse settling process, viz. (i) initial stir-up phase with limited flocculation, (ii) enhanced settling phase characterized by increased flocculation, and (iii) consolidation phase. The simulations demonstrate that cohesive forces accelerate the overall settling process primarily because smaller grains attach to larger ones and settle in their wakes. For the present cohesive number values, we observe that settling can be accelerated by up to 29 %. We propose physically based parametrization of classical hindered settling functions introduced by earlier authors, in order to account for cohesive forces. An investigation of the energy budget shows that, even though the work of the collision forces is much smaller than that of the hydrodynamic drag forces, it can substantially modify the relevant energy conversion processes.

Abstract: Fingering is a hydrodynamic instability that occurs when a more mobile fluid displaces a fluid of lower mobility. When the primary source of the mobility difference is viscosity, the instability is termed viscous fingering. Viscous fingering is often, though not always, undesirable in industrial processes, particularly secondary petroleum recovery. Linear stability analysis by Hejazi et al. has indicated that the production of a non-monotonic viscosity profile can stabilize the interface. Herein, we use step-growth polymerization at the interface between two miscible monomers as a model system. In particular, a dithiol monomer displaced a diacrylate that reacted to form a linear polymer that behaves as a Newtonian fluid. Viscous fingering was imaged in a horizontal Hele-Shaw cell via Schlieren, which is sensitive to changes in index of refraction, and therefore polymer conversion. By varying reaction rate via initiator concentration along with flow rate via a syringe pump, we were able to demonstrate increasing stabilization of the flow with increasing Damköhler number. Graphical abstract: [Figure not available: see fulltext.].

When particles settle through a stable temperature or salinity gradient they can drive an instability known as sedimentary fingering convection. This phenomenon is thought to occur beneath sediment-rich river plumes in lakes and oceans, in the context of marine snow where decaying organic materials serve as the suspended particles or in the atmosphere in the presence of aerosols or volcanic ash. Laboratory experiments of Houk & Green (Deep-Sea Res., vol. 20, 1973, pp. 757–761) and Green (Sedimentology, vol. 34(2), 1987, pp. 319–331) have shown sedimentary fingering convection to be similar to the more commonly known thermohaline fingering convection in many ways. Here, we study the phenomenon using three-dimensional direct numerical simulations. We find evidence for layer formation in sedimentary fingering convection in regions of parameter space where it does not occur for non-sedimentary systems. This is due to two complementary effects. Sedimentation affects the turbulent fluxes and broadens the region of parameter space unstable to the $\unicode[STIX]{x1D6FE}$-instability (Radko, J. Fluid Mech., vol. 497, 2003, pp. 365–380) to include systems at larger density ratios. It also gives rise to a new layering instability that exists in $\unicode[STIX]{x1D6FE}$-stable regimes. The former is likely quite ubiquitous in geophysical systems for sufficiently large settling velocities, while the latter probably grows too slowly to be relevant, at least in the context of sediments in water.

We analyse the linear stability of stably stratified fluids whose density depends on two scalar fields where one of the scalar fields is unstably stratified and involves a settling velocity. Such conditions may be found, for example, in flows involving the transport of sediment in addition to heat or salt. A linear stability analysis for constant-gradient base states demonstrates that the settling velocity generates a phase shift between the perturbation fields of the two scalars, which gives rise to a novel, settling-driven instability mode. This instability mechanism favours the growth of waves that are inclined with respect to the horizontal. It is active for all density and diffusivity ratios, including for cases in which the two scalars diffuse at identical rates. If the scalars have unequal diffusivities, it competes with the elevator mode waves of the classical double-diffusive instability. We present detailed linear stability results as a function of the governing dimensionless parameters, including for lateral gradients of the base state density fields that result in predominantly horizontal intrusion instabilities. Highly resolved direct numerical simulation results serve to illustrate the nonlinear competition of the various instabilities for such flows in different parameter regimes.

We extend the vorticity-based modelling approach of Borden & Meiburg (Phys. Fluids, vol. 25 (10), 2013, 101301) to non-Boussinesq gravity currents and derive an analytical expression for the Froude number without the need for an energy closure or any assumptions about the pressure. The Froude-number expression we obtain reduces to the correct form in the Boussinesq limit and agrees closely with simulation data. Via detailed comparisons with simulation results, we furthermore assess the validity of three key assumptions underlying both our as well as earlier models: (i) steady-state flow in the moving reference frame; (ii) inviscid flow; and (iii) horizontal flow sufficiently far in front of and behind the current. The current approach does not require an assumption of zero velocity in the current.

Interfaces between different fluids can be unstable with regard to hydrodynamic instabilities such as viscous fingering or buoyancy-driven convection. To study such instabilities experimentally for transparent fluids, dyes or chemical indicators are most often used to track the dynamics. While the interfacial deformation can easily be tracked by color changes, it is difficult to have access to the internal flow structure for comparison with theoretical predictions. To overcome this problem, a modification of a Schlieren technique is introduced to image 3D flows during viscously driven instabilities in a horizontal Hele-Shaw cell without using any dye or chemical indicator. The method is exquisitely sensitive, readily yielding information about 3D flows in gaps under a millimeter and allowing imaging of the flow structure internal to the fingers, rather than merely imaging the flow boundary. Following a description of the technique, visualization of dynamics for nonreactive water–glycerol and reactive displacements is presented revealing previously unobserved internal flows. These flows are tentatively interpreted in terms of known theoretical predictions.