Influence of the history force on inertial particle advection: gravitational effects and horizontal diffusion

Lagrangian flow networks for passive dispersal: Tracers versus finite-size particles

The transport and distribution of organisms such as larvae, seeds, or litter in the ocean as well as particles in industrial flows is often approximated by a transport of tracer particles. We present a theoretical investigation to check the accuracy of this approximation by studying the transport of inertial particles between different islands embedded in an open hydrodynamic flow aiming at the construction of a Lagrangian flow network reflecting the connectivity between the islands. To this end, we formulate a two-dimensional kinematic flow field which allows the placement of an arbitrary number of islands at arbitrary locations in a flow of prescribed direction. To account for the mixing in the flow, we include a von Kármán vortex street in the wake of each island. We demonstrate that the transport probabilities of inertial particles making up the links of the Lagrangian flow network essentially depend on the properties of the particles, i.e., their Stokes number, the properties of the flow, and the geometry of the setup of the islands. We find a strong segregation between aerosols and bubbles. Upon comparing the mobility of inertial particles to that of tracers or neutrally buoyant particles, it becomes apparent that the tracer approximation may not always accurately predict the probability of movement. This can lead to inconsistent forecasts regarding the fate of marine organisms, seeds, litter, or particles in industrial flows.

On the long-time persistence of hydrodynamic memory

The Basset-Boussinesq-Oseen (BBO) equation correctly describes the nonuniform motion of a spherical particle at a low Reynolds number. It contains an integral term with a singular kernel which accounts for the diffusion of vorticity around the particle throughout its entire history. However, if there are any departures in either rigidity or shape from a solid sphere, besides the integral force with a singular kernel, the Basset history force, we should add a second history force with a non-singular kernel, related to the shape or composition of the particle. In this work, we introduce a fractional generalized Basset-Boussinesq-Oseen equation which includes both history terms as fractional derivatives. Using the Laplace transform, an integral representation of the solution is obtained. For a driven single particle, the solution shows that memory effects persist indefinitely under rather general driving conditions.

Oxygen-Selective Diffusion-Bubbling Membranes with Core-Shell Structure: Bubble Dynamics and Unsteady Effects

Oxygen is the second-largest-volume industrial gas that is mainly produced using cryogenic air separation. However, the state-of-the-art cryogenic technology thermodynamic efficiency has approached a theoretical limit as near as is practicable. Therefore, there is stimulus to develop an alternative technology for efficient oxygen separation from air. Mixed ionic electronic-conducting (MIEC) ceramic membrane-based oxygen separation technology could become this alternative, but commercialization aspects, including cost, have revealed inadequacies in ceramic membrane materials. Currently, diffusion-bubbling molten oxide membrane-based oxygen separation technology is being developed. It is a potentially disruptive technology that would propose an improvement in oxygen purity and a reduction in capital costs. Bubbles play an important role in ensuring the oxygen mass transfer in diffusion-bubbling membranes. However, there is not sufficient understanding of the bubble dynamics. This understanding is important to be able to control transport properties of these membranes and assess their potential for technological application. The aim of this feature article is to highlight the progress made in developing this understanding and specify the directions for future research.

Hydrodynamic memory can boost enormously driven nonlinear diffusion and transport

Hydrodynamic memory force or Basset force has been known since the 19th century. Its influence on Brownian motion remains, however, mostly unexplored. Here we investigate its role in nonlinear transport and diffusion within a paradigmatic model of tilted washboard potential. In this model, a giant enhancement of driven diffusion over its potential-free limit [Phys. Rev. Lett. 87, 010602 (2001)PRLTAO0031-900710.1103/PhysRevLett.87.010602] presents a well-established paradoxical phenomenon. In the overdamped limit, it occurs at a critical tilt of vanishing potential barriers. However, for weak damping, it takes place surprisingly at another critical tilt, where the potential barriers are clearly expressed. Recently we showed [Phys. Rev. Lett. 123, 180603 (2019)PRLTAO0031-900710.1103/PhysRevLett.123.180603] that Basset force could make such a diffusion enhancement enormously large. In this paper, we discover that even for moderately strong damping, where the overdamped theory works very well when the memory effects are negligible, substantial hydrodynamic memory unexpectedly makes a strong impact. First, the diffusion boost occurs at nonvanishing potential barriers and can be orders of magnitude larger. Second, transient anomalous diffusion regimes emerge over many time decades and potential periods. Third, particles' mobility can also be dramatically enhanced, and a long transient supertransport regime emerges.

Fractional Hydrodynamic Memory and Superdiffusion in Tilted Washboard Potentials

Diffusion in tilted washboard potentials can paradoxically exceed free normal diffusion. The effect becomes much stronger in the underdamped case due to inertial effects. What happens upon inclusion of usually neglected fractional hydrodynamics memory effects (Basset-Boussinesq frictional force), which result in a heavy algebraic tail of the velocity autocorrelation function of the potential-free diffusion making it transiently superdiffusive? Will a giant enhancement of diffusion become even stronger, and the transient superdiffusion last even longer? These are the questions that we answer in this Letter based on an accurate numerical investigation. We show that a resonancelike enhancement of normal diffusion becomes indeed much stronger and sharper. Moreover, a long-lasting transient regime of superdiffusion, including Richardson-like diffusion, ⟨δx^{2}(t)⟩∝t^{3} and ballistic supertransport, ⟨δx(t)⟩∝t^{2}, is revealed.

Inhomogeneities and caustics in the sedimentation of noninertial particles in incompressible flows

In an incompressible flow, fluid density remains invariant along fluid element trajectories. This implies that the spatial distribution of non-interacting noninertial particles in such flows cannot develop density inhomogeneities beyond those that are already introduced in the initial condition. However, in certain practical situations, density is measured or accumulated on (hyper-) surfaces of dimensionality lower than the full dimensionality of the flow in which the particles move. An example is the observation of particle distributions sedimented on the floor of the ocean. In such cases, even if the initial distribution of noninertial particles is uniform but its support is finite, advection in an incompressible flow will give rise to inhomogeneities in the observed density. In this paper, we analytically derive, in the framework of an initially homogeneous particle sheet sedimenting toward a bottom surface, the relationship between the geometry of the flow and the emerging distribution. From a physical point of view, we identify the two processes that generate inhomogeneities to be the stretching within the sheet and the projection of the deformed sheet onto the target surface. We point out that an extreme form of inhomogeneity, caustics, can develop for sheets. We exemplify our geometrical results with simulations of particle advection in a simple kinematic flow, study the dependence on various parameters involved, and illustrate that the basic mechanisms work similarly if the initial (homogeneous) distribution occupies a more general region of finite extension rather than a sheet.

Turbulence induces metabolically costly behaviors and inhibits food capture in oyster larvae, causing net energy loss

Planktotrophic invertebrate larvae require energy to develop, disperse and settle successfully, and it is unknown how their energetics are impacted by turbulence. Ciliated larvae gain metabolic energy from their phytoplankton food to offset the energetic costs of growth, development and ciliary activity for swimming and feeding. Turbulence may affect the energetic balance by inducing behaviors that alter the metabolic costs and efficiency of swimming, by raising the encounter rate with food particles and by inhibiting food capture. We used experiments and an empirical model to quantify the net rate of energy gain, swimming efficiency and food capture efficiency for eyed oyster larvae (Crassostrea virginica) in turbulence. At dissipation rates representative of coastal waters, larvae lost energy even when food concentrations were very high. Both feeding activity and turbulence-induced behaviors incurred high metabolic costs. Swimming efficiency was concave up versus dissipation rate, suggesting that ciliary activity for food handling became more costly while swimming became more efficient with turbulence intensity. Though counter-intuitive, swimming may have become more efficient in turbulence because vorticity-induced rotation caused larvae to swim more horizontally, which requires less effort than swimming vertically against the pull of gravity. Overall, however, larvae failed to offset high activity costs with food energy gains because turbulence reduced food capture efficiency more than it enhanced food encounter rates. Younger, smaller larvae may have some energetic advantages, but competent larvae would lose energy at turbulence intensities they experience frequently, suggesting that turbulence-induced starvation may account for much of oysters' high larval mortality.

Transport, diffusion, and energy studies in the Arnold-Beltrami-Childress map

We study the transport and diffusion properties of passive inertial particles described by a six-dimensional dissipative bailout embedding map. The base map chosen for the study is the three-dimensional incompressible Arnold-Beltrami-Childress (ABC) map chosen as a representation of volume preserving flows. There are two distinct cases: the two-action and the one-action cases, depending on whether two or one of the parameters (A,B,C) exceed 1. The embedded map dynamics is governed by two parameters (α,γ), which quantify the mass density ratio and dissipation, respectively. There are important differences between the aerosol (α<1) and the bubble (α>1) regimes. We have studied the diffusive behavior of the system and constructed the phase diagram in the parameter space by computing the diffusion exponents η. Three classes have been broadly classified-subdiffusive transport (η<1), normal diffusion (η≈1), and superdiffusion (η>1) with η≈2 referred to as the ballistic regime. Correlating the diffusive phase diagram with the phase diagram for dynamical regimes seen earlier, we find that the hyperchaotic bubble regime is largely correlated with normal and superdiffusive behavior. In contrast, in the aerosol regime, ballistic superdiffusion is seen in regions that largely show periodic dynamical behaviors, whereas subdiffusive behavior is seen in both periodic and chaotic regimes. The probability distributions of the diffusion exponents show power-law scaling for both aerosol and bubbles in the superdiffusive regimes. We further study the Poincáre recurrence times statistics of the system. Here, we find that recurrence time distributions show power law regimes due to the existence of partial barriers to transport in the phase space. Moreover, the plot of average particle kinetic energies versus the mass density ratio for the two-action case exhibits a devil's staircase-like structure for higher dissipation values. We explain these results and discuss their implications for realistic systems.

Experimental study on dynamics of coherent structures formed by inertial solid particles in three-dimensional periodic flows

We present experimental results obtained under normal gravity on the dynamics of solid particles in periodic oscillatory thermocapillary-driven flows in a non-isothermal liquid bridge made of decane. Inertial particles of different densities and in the size range approximately 0.75-75 μm are able to form stable coherent structures (particle accumulation structures, or PASs). Two image processing techniques were developed and successfully applied to compute time required for an ensemble of particles to form a structure. It is shown that the formation time grows with the decrease of the Stokes number. The observations indicate the probable irrelevance of the memory term for these experiments. Two types of PAS were observed-single (SL-I) and double-loop (SL-II)-which sometimes co-existed. Only large or very dense particles may form an SL-II type structure. A number of novel features of the system were perceived. In some cases, intermittently stable structures emerged (their dynamics is characterized by alternating time intervals during which a structure exists and is destroyed). Whereas in most experiments we observed a conventional symmetric and centered PAS, there were cases when a long-term stable asymmetric structure appeared. Experiments wherein two different types of PAS-forming particles were used simultaneously revealed the destructive role of collisions between the particles on formation of structures.

Critical Stokes number for the capture of inertial particles by recirculation cells in two-dimensional quasisteady flows

Inertial particles are often observed to be trapped, temporarily or permanently, by recirculation cells which are ubiquitous in natural or industrial flows. In the limit of small particle inertia, determining the conditions of trapping is a challenging task, as it requires a large number of numerical simulations or experiments to test various particle sizes or densities. Here, we investigate this phenomenon analytically and numerically in the case of heavy particles (e.g., aerosols) at low Reynolds number, to derive a trapping criterion that can be used both in analytical and numerical velocity fields. The resulting criterion allows one to predict the characteristics of trapped particles as soon as single-phase simulations of the flow are performed. Our analysis is valid for two-dimensional particle-laden flows in the vertical plane, in the limit where the particle inertia, the free-fall terminal velocity, and the flow unsteadiness can be treated as perturbations. The weak unsteadiness of the flow generally induces a chaotic tangle near heteroclinic or homoclinic cycles if any, leading to the apparent diffusion of fluid elements through the boundary of the cell. The critical particle Stokes number St_{c} below which aerosols also enter and exit the cell in a complex manner has been derived analytically, in terms of the flow characteristics. It involves the nondimensional curvature-weighted integral of the squared velocity of the steady fluid flow along the dividing streamline of the recirculation cell. When the flow is unsteady and St>St_{c}, a regular motion takes place due to gravity and centrifugal effects, like in the steady case. Particles driven towards the interior of the cell are trapped permanently. In contrast, when the flow is unsteady and St<St_{c}, particles wander in a chaotic manner in the vicinity of the border of the cell, and can escape the cell.

Effect of the vertical component of diffusion on passive scalar transport in an isolated vortex model

On the basis of the ellipsoidal vortex model and a Monte-Carlo-type diffusion simulation, we examine the flux and ensuing distribution of passive fluid particles through the boundary of an idealized geophysical vortex. Our focus is on features that the horizontal and vertical diffusion components introduce into the fluid particle transport. We examine the concurrent effect of both components, and we compare it with the only horizontal diffusion impact. We analyze the ellipsoid vortex model in two cases: (i) the steady state when the ellipsoid is motionless, i.e., there is no variation in its axes' lengths, and consequently the exterior fluid is not being stirred; (ii) the perturbed case when the ellipsoid rotates periodically, varying it axes' lengths, which results in the appearance of stirred fluid outside the ellipsoid. Influenced by diffusion, a fluid particle is now permitted to move to another vertical horizon, thus there is an increased possibility that the particle will eventually be located in the exterior stirred region rather than in the ellipsoid vortex with the regular dynamics. This is because the area of the horizontal section of the ellipsoid vortex decreases with depth, but the region of stirred exterior fluid extends significantly deeper. Numerical calculations show that factoring in the vertical component of diffusion significantly affects scalar spreading in the horizontal plane in the perturbed case, while in the steady state the vertical component of diffusion only induces dispersion linear growth according to a Gaussian distribution.

Dynamics of impurities in a three-dimensional volume-preserving map

We study the dynamics of inertial particles in three-dimensional incompressible maps, as representations of volume-preserving flows. The impurity dynamics has been modeled, in the Lagrangian framework, by a six-dimensional dissipative bailout embedding map. The fluid-parcel dynamics of the base map is embedded in the particle dynamics governed by the map. The base map considered for the present study is the Arnold-Beltrami-Childress (ABC) map. We consider the behavior of the system both in the aerosol regime, where the density of the particle is larger than that of the base flow, as well as the bubble regime, where the particle density is less than that of the base flow. The phase spaces in both the regimes show rich and complex dynamics with three types of dynamical behaviors--chaotic structures, regular orbits, and hyperchaotic regions. In the one-action case, the aerosol regime is found to have periodic attractors for certain values of the dissipation and inertia parameters. For the aerosol regime of the two-action ABC map, an attractor merging and widening crisis is identified using the bifurcation diagram and the spectrum of Lyapunov exponents. After the crisis an attractor with two parts is seen, and trajectories hop between these parts with period 2. The bubble regime of the embedded map shows strong hyperchaotic regions as well as crisis induced intermittency with characteristic times between bursts that scale as a power law behavior as a function of the dissipation parameter. Furthermore, we observe a riddled basin of attraction and unstable dimension variability in the phase space in the bubble regime. The bubble regime in the one-action case shows similar behavior. This study of a simple model of impurity dynamics may shed light upon the transport properties of passive scalars in three-dimensional flows. We also compare our results with those seen earlier in two-dimensional flows.