A new model for settling velocity of non-spherical particles

The settlement of non-spherical particles, such as propagules of plants and natural sediments, are commonly observed in riverine ecosystems. The settling process is influenced by both particle properties (size, density and shape) and fluid properties (density and viscosity). Therefore, the drag law of non-spherical particles is a function of both particle Reynolds number and particle shape. Herein, a total of 828 settling data are collected from the literatures, which cover a wide range of particle Reynolds number (0.008–10000). To characterize the influence of particle shapes, sphericity is adopted as the general shape factor, which varies from 0.421 to 1.0. By comparing the measured drag with the standard drag curve of spheres, we modify the spherical drag law with three shape-dependent functions to develop a new drag law for non-spherical particles. Combined with an iterative procedure, a new model is thus obtained to predict the settling velocity of non-spherical particles of various shapes and materials. Further applications in hydrochorous propagule dispersal and sediment transport are projected based on deeper understanding of the settling process.

Parameterizing non-propagating form drag over rough bathymetry

Slowly-evolving stratified flow over rough topography is subject to substantial drag due to internal motions, but often numerical simulations are carried out at resolutions where this “wave” drag must be parameterized. Here we highlight the importance of internal drag from topography with scales that cannot radiate internal waves, but may be highly non-linear, and we propose a simple parameterization of this drag that has a minimum of fit parameters compared to existing schemes. The parameterization smoothly transitions from a quadratic drag law () for low- (linear wave dynamics) to a linear drag law () for high- flows (non-linear blocking and hydraulic dynamics), where N is the stratification, h is the height of the topography, and u0 is the near-bottom velocity; the parameterization does not have a dependence on Coriolis frequency. Simulations carried out in a channel with synthetic bathymetry and steady body forcing indicate that this parameterization accurately predicts drag across a broad range of forcing parameters when the effect of reduced near-bottom mixing is taken into account by reducing the effective height of the topography. The parameterization is also tested in simulations of wind-driven channel flows that generate mesoscale eddy fields, a setup where the downstream transport is sensitive to the bottom drag parameterization and its effect on the eddies. In these simulations, the parameterization replicates the effect of rough bathymetry on the eddies. If extrapolated globally, the sub-inertial topographic scales can account for 2.7 TW of work done on the low-frequency circulation, an important sink that is redistributed to mixing in the open ocean.

Drag of a cylindrical object in a two-dimensional granular environment

The drag of a cylindrical object in a two-dimensional granular environment is numerically studied. It is found that the drag law is fitted by the sum of the yield force and the dynamic force, the latter of which is reproduced by a simple collision model. The angular dependence of the radial stress on the surface of the object is given by the Gaussian below the yield force. The probability of the velocity drops of the object is investigated above the yield force, where this probability is independent of the packing fraction and the drag force.

Direct Measurement of Unsteady Microscale Stokes Flow Using Optically Driven Microspheres

A growing body of work on the dynamics of eukaryotic flagella has noted that their oscillation frequencies are sufficiently high that the viscous penetration depth of unsteady Stokes flow is comparable to the scales over which flagella synchronize. Incorporating these effects into theories of synchronization requires an understanding of the global unsteady flows around oscillating bodies. Yet, there has been no precise experimental test on the microscale of the most basic aspects of such unsteady Stokes flow: the orbits of passive tracers and the position-dependent phase lag between the oscillating response of the fluid at a distant point and that of the driving particle. Here, we report the first such direct Lagrangian measurement of this unsteady flow. The method uses an array of 30 submicron tracer particles positioned by a time-shared optical trap at a range of distances and angular positions with respect to a larger, central particle, which is then driven by an oscillating optical trap at frequencies up to 400 Hz. In this microscale regime, the tracer dynamics is considerably simplified by the smallness of both inertial effects on particle motion and finite-frequency corrections to the Stokes drag law. The tracers are found to display elliptical Lissajous figures whose orientation and geometry are in agreement with a low-frequency expansion of the underlying dynamics, and the experimental phase shift between motion parallel and orthogonal to the oscillation axis exhibits a predicted scaling form in distance and angle. Possible implications of these results for synchronization dynamics are discussed.

A Swing of Beauty: Pendulums, Fluids, Forces, and Computers

While pendulums have been around for millennia and have even managed to swing their way into undergraduate curricula, they still offer a breadth of complex dynamics to which some has still yet to have been untapped. To probe into the dynamics, we developed a computational fluid dynamics (CFD) model of a pendulum using the open-source fluid-structure interaction (FSI) software, IB2d. Beyond analyzing the angular displacements, speeds, and forces attained from the FSI model alone, we compared its dynamics to the canonical damped pendulum ordinary differential equation (ODE) model that is familiar to students. We only observed qualitative agreement after a few oscillation cycles, suggesting that there is enhanced fluid drag during our setup’s initial swing, not captured by the ODE’s linearly-proportional-velocity damping term, which arises from the Stokes Drag Law. Moreover, we were also able to investigate what otherwise could not have been explored using the ODE model, that is, the fluid’s response to a swinging pendulum—the system’s underlying fluid dynamics.

A 2D Multiphase Model of Drop Behavior during Electroslag Remelting

Electroslag remelting is a process extensively used to produce metallic ingots with high quality standards. During the remelting operation, liquid metal droplets fall from the electrode through the liquid slag before entering the liquid pool of the secondary ingot. To better understand the process and help to optimize the operating condition choice, a 2D axisymmetric multiphase model of the slag domain has been developed using a two fluid Eulerian approach. During their fall, droplets hydrodynamic interactions are calculated thanks to an appropriate drag law. Influence of droplets on the electromagnetic field and on the slag hydrodynamics is discussed, as well as their heat exchange with the slag. Even with a small volume fraction, the droplets influence is noticeable. The present investigation shows that small droplets have a large influence on the slag hydrodynamics, due to a great momentum exchange. However heat transfer is more influenced by large drops, which are found to be relatively far from the thermal equilibrium with the slag phase.

The Control of Surface Friction on the Scales of Baroclinic Eddies in a Homogeneous Quasigeostrophic Two-Layer Model

In idealized models of the extratropical troposphere, both β and surface friction can control the equilibrated scales of baroclinic eddies by stopping the inverse cascade. A scaling theory on how surface friction alone sets these scales was proposed by Held in 1999 in the case of a quadratic drag law. However, the theory breaks down when friction is modeled by linear damping, and there are other reasons to suspect that it is oversimplified. An ideal system to test the theory is the homogeneous two-layer quasigeostrophic model in the β = 0 limit with quadratic damping. This study investigates some numerical simulations of the model to analyze two causes of the theory’s breakdown. They are 1) the asymmetry between two layers due to confinement of friction to the lower layer and 2) deviation from a spectrally local inverse energy cascade due to the spread of wavenumbers over which energy is input into the barotropic mode. The former is studied by comparing the simulations with drag appearing asymmetrically or symmetrically between the two layers. The latter is addressed with a heuristic modification of the theory. A regime where eddies equilibrate without an inverse cascade is also examined. A comparison is then made between quadratic and linear drag simulations. The connection to a competing theory based on the dynamics of equivalent barotropic vortices with thermal signatures is further discussed. Finally, we present an example of an inhomogeneous statistically steady state to argue that the diffusivity obtained from the homogeneous model has relevance to more realistic configurations.