Mahesh Bandi

Dissolution-driven convection in a Hele-Shaw cell

Motivated by convection in the context of geological carbon-dioxide (CO2) storage, we present an experimental study of dissolution-driven convection in a Hele–Shaw cell for Rayleigh numbers R in the range 100 < R < 1700. We use potassium permanganate (KMnO4) in water as an analog for CO2 in brine and infer concentration profiles at high spatial and temporal resolution and accuracy from transmitted light intensity. We describe behavior from first contact up to 65% average saturation and measure several global quantities including dissolution flux, average concentration, amplitude of perturbations away from pure one-dimensional diffusion, and horizontally averaged concentration profiles. We show that the flow evolves successively through distinct regimes starting with a simple one-dimensional diffusional profile. This is followed by linear growth in which fingers are initiated and grow quasiexponentially, independently of one-another. Once the fingers are well-established, a flux-growth regime begins as fresh fluid is brought to the interface and contaminated fluid removed, with the flux growing to a local maximum. During this regime, fingers still propagate independently. However, beyond the flux maximum, fingers begin to interact and zip together from the root down in a merging regime. Several generations of merging occur before only persistent primary fingers remain. Beyond this, the reinitiation regime begins with new fingers created between primary existing ones before merging into them. Through appropriate scaling, we show that the regimes are universal and independent of layer thickness (equivalently R) until the fingers hit the bottom. At this time, progression through these regimes is interrupted and the flow transitions to a saturating regime. In this final regime, the flux gradually decays in a manner well described by a Howard-style phenomenological model.

Shock-driven jamming and periodic fracture of particulate rafts

A tenuous monolayer of hydrophobic particles at the air-water interface often forms a scum or raft. When such a monolayer is disturbed by the localized introduction of a surfactant droplet, a radially divergent surfactant shock front emanates from the surfactant origin and packs the particles into a jammed, compact, annular band with a packing fraction that saturates at a peak packing fraction *. As the resulting two-dimensional, disordered elastic band grows with time and is driven radially outwards by the surfactant, it fractures to form periodic triangular cracks with robust geometrical features. We find that the number of cracks N and the compaction band radius R* at fracture onset vary monotonically with the initial packing fraction (init). However, the compaction bands width W* is constant for all init. A simple geometric theory that treats the compaction band as an elastic annulus, and accounts for mass conservation allows us to deduce that N2πR*/W*4πRCP/init, a result we verify both experimentally and numerically. We show that the essential ingredients for this phenomenon are an initially low enough particulate packing fraction that allows surfactant-driven advection to cause passive jamming and eventual fracture of the hydrophobic particulate interface.

Time evolution of a fractal distribution: Particle concentrations in free-surface turbulence

Abstract Steady-state turbulence is generated in a tank of water and the trajectories of particles forming a compressible system on the surface are tracked in time. The initial uniformly distributed floating particles coagulate and form a fractal structure, a rare manifestation of a strange attractor observable in real space. The surface pattern reaches a steady state in approximately 1 s. Measurements are made of the fractal dimensions D q ( t ) ( q = 1 to 6) of the floating particles starting with the uniform distribution D q ( 0 ) = 2 for Taylor Microscale Reynolds number R e λ ≃ 160 . Focus is on the time evolution of the correlation dimension D 2 ( t ) as the steady state is approached. This steady state is reached in several large eddy turnover times and does so at an exponential rate.

Sensitivity of Granular Force Chain Orientation to Disorder-Induced Metastable Relaxation.

A two-dimensional system of photoelastic disks subject to vertical tapping against gravity was experimentally monitored from ordered to disordered configurations by varying bidispersity. The packing fraction ϕ, coordination number Z, and an appropriately defined force-chain orientational order parameter S all exhibit as similar sharp transition with a small increase in disorder. A measurable change in S, but not ϕ and Z, was detected under tapping. We find disorder-induced metastability does not show configurational relaxation, but can be detected via force-chain reorientations.

Spectrum of Wind Power Fluctuations

Wind power fluctuations for an individual turbine and plant have been widely reported to follow the Kolmogorov spectrum of atmospheric turbulence; both vary with a fluctuation time scale τ as τ^{2/3}. Yet, this scaling has not been explained through turbulence theory. Using turbines as probes of turbulence, we show the τ^{2/3} scaling results from a large scale influence of atmospheric turbulence. Owing to this long-range influence spanning 100s of kilometers, when power from geographically distributed wind plants is summed into aggregate power at the grid, fluctuations average (geographic smoothing) and their scaling steepens from τ^{2/3}→τ^{4/3}, beyond which further smoothing is not possible. Our analysis demonstrates grids have already reached this τ^{4/3} spectral limit to geographic smoothing.

Craig's XY distribution and the statistics of Lagrangian power in two-dimensional turbulence.

We examine the probability distribution function (PDF) of the energy injection rate (power) in numerical simulations of stationary two-dimensional (2D) turbulence in the Lagrangian frame. The simulation is designed to mimic an electromagnetically driven fluid layer, a well-documented system for generating 2D turbulence in the laboratory. In our simulations, the forcing and velocity fields are close to Gaussian. On the other hand, the measured PDF of injected power is very sharply peaked at zero, suggestive of a singularity there, with tails which are exponential but asymmetric. Large positive fluctuations are more probable than large negative fluctuations. It is this asymmetry of the tails which leads to a net positive mean value for the energy input despite the most probable value being zero. The main features of the power distribution are well described by Craigs XY distribution for the PDF of the product of two correlated normal variables. We show that the power distribution should exhibit a logarithmic singularity at zero and decay exponentially for large absolute values of the power. We calculate the asymptotic behavior and express the asymmetry of the tails in terms of the correlation coefficient of the force and velocity. We compare the measured PDFs with the theoretical calculations and briefly discuss how the power PDF might change with other forcing mechanisms.

Linear stability analysis for monami in a submerged seagrass bed

Amala Mahadevan Woods Hole Oceanographic Institution, Woods Hole MA 02543 USA Abstract The onset of monami – the synchronous waving of sea grass beds driven by a steady flow – is modeled as a linear instability of the flow. Our model treats the drag exerted by the grass in establishing the steady flow profile, and in damping out perturbations to it. This damping leads to a finite threshold flow for the instability, which agrees with experimental observations. This role of vegetation drag differentiates our mechanism from the previous hypothesis that the Kelvin-Helmholtz instability underlies monami.

A pendulum in a flowing soap film

We consider the dynamics of a pendulum made of a rigid ring attached to an elastic filament immersed in a flowing soap film. The system shows an oscillatory instability whose onset is a function of the flow speed, length of the supporting string, the ring mass, and ring radius. We characterize this system and show that there are different regimes where the frequency is dependent or independent of the pendulum length depending on the relative magnitude of the added-mass. Although the system is an infinite-dimensional, we can explain many of our results in terms of a one degree-of-freedom system corresponding to a forced pendulum. Indeed, using the vorticity measured via particle imaging velocimetry allows us to make the model quantitative, and a comparison with our experimental results shows we can capture the basic phenomenology of this system.

Probability distribution of power fluctuations in turbulence

We study local power fluctuations in numerical simulations of stationary, homogeneous, isotropic turbulence in two and three dimensions with Gaussian forcing. Due to the near-Gaussianity of the one-point velocity distribution, the probability distribution function (pdf) of the local power is well modeled by the pdf of the product of two joint normally distributed variables. In appropriate units, this distribution is parametrized only by the mean dissipation rate, epsilon. The large deviation function for this distribution is calculated exactly and shown to satisfy a fluctuation relation (FR) with a coefficient which depends on epsilon. This FR is entirely statistical in origin. The deviations from the model pdf are most pronounced for positive fluctuations of the power and can be traced to a slightly faster than Gaussian decay of the tails of the one-point velocity pdf. The resulting deviations from the FR are consistent with several recent experimental studies.

A minimal description of morphological hierarchy in two-dimensional aggregates

A dimensionless parameter Λ is proposed to describe a hierarchy of morphologies in two-dimensional (2D) aggregates formed due to varying competition between short-range attraction and long-range repulsion. Structural transitions from finite non-compact to compact to percolated structures are observed in the configurations simulated by molecular dynamics at a constant temperature and density. Configurational randomness across the transition, measured by the two-body excess entropy S2, exhibits data collapse with the average potential energy [small epsilon, Greek, macron] of the systems. Independent master curves are presented among S2, the reduced second virial coefficient B2* and Λ, justifying this minimal description. This work lays out a coherent basis for the study of 2D aggregate morphologies relevant to diverse nano- and bio-processes.