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

Abstract

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. Sedimentation of small particles in complex flows is an outstanding problem in science and technology. In particular, the sinking of biogenic particles from the marine surface to the bottom is a fundamental process of the biological carbon pump, playing a key role in the global carbon cycle. A complete understanding of this problem is still lacking. It has been recently shown that despite fluid incompressibility, sedimenting particles moving as noninertial particles in the ocean on large scales show density inhomogeneities when accumulated on some bottom surface. Here, we analytically derive the relation between the geometry of the flow and the emerging distribution for an initially homogeneous sheet of particles. 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 conditions under which an extreme form of inhomogeneity, caustics, can develop.