Physics

Bubbles in complex fluids are often desirable, and sometimes simply inevitable, in the processing of formulated products. Bubbles can rise by buoyancy, grow or dissolve by mass transfer, and readily respond to changes in pressure, thereby applying a deformation to the surrounding complex fluid. The deformation field around a stationary, spherical bubble undergoing a change in radius is simple and localised, thus making it suitable for rheological measurements. This article reviews emerging approaches to extract information on the rheology of complex fluids by analysing bubble dynamics. The focus is on three phenomena: changes in radius by mass transfer, harmonic oscillations driven by an acoustic wave, and bubble collapse. These phenomena cover a broad range of deformation frequencies, from 10 −4 to 10 6 Hz, thus paving the way to broadband microrheology using bubbles as active probes. The outstanding challenges that need to be overcome to achieve a robust technique are also discussed.

In acoustofluidics, acoustic resonance modes for fluid and microparticle handling are traditionally excited by bulk piezoelectric transducers. In this work, we demonstrate by numerical simulation in three dimensions (3D) that integrated piezoelectric thin-film transducers constituting less than 0.1% of the device work equally well. The simulations are done using a well-tested and experimentally validated numerical model. Our proof-of-concept example is a water-filled straight channel embedded in a mm-sized glass chip with a 1-um thick thin-film transducer made of (Al,Sc)N. We compute the acoustic energy, streaming, and radiation force, and show that it is comparable to that of a conventional silicon-glass device actuated by a bulk PZT transducer. The ability of the thin-film transducer to create the desired acoustofluidic effects in bulk acoustofluidic devices rely on three physical aspects: The in-plane-expansion of the thin-film transducer under the orthogonal applied electric field, the acoustic whole-system resonance of the device, and the high Q-factor of the elastic solid constituting the bulk part of the device. Consequently, the thin-film device is surprisingly insensitive to the Q-factor and resonance properties of the thin-film transducer.

When a rising bubble in a Newtonian liquid reaches the liquid-air interface, it can burst, leading to the formation of capillary waves and a jet on the surface. Here, we numerically study this phenomenon in a yield stress fluid. We show how viscoplasticity controls the fate of these capillary waves and their interaction at the bottom of the cavity. Unlike Newtonian liquids, the free surface converges to a non-flat final equilibrium shape once the driving stresses inside the pool fall below the yield stress. Details of the dynamics, including the flow's energy budgets, are discussed. The work culminates in a regime map with four main regimes with different characteristic behaviours.

This paper reports on the investigation of combustion instabilities in a Methane air nonpremixed half dump combustor for different flow Reynolds number using Computational aero acoustic simulation. In order to simulate the flow physics under turbulent flow conditions, Detached Eddy Simulation in conjunction with generalized eddy dissipation model is adopted, while a CAA formulation based on Lighthills acoustic analogy is used for computing the acoustic field. It is observed that the unsteady pressure signals either predominantly arise from the natural acoustic modes of the duct or the local flow fluctuations in the vortex shedding process downstream of the dump plane, giving rise to different dominant frequencies at different spatial locations at lower Re and a single dominant lockedon frequency at higher Re. The nondimensional numbers (Helmholtz and Strouhal numbers) are used to characterize the duct acoustic modes from the vortex shedding modes. Under reacting flow conditions, unsteady heat release and pressure oscillations are monitored to compute Rayleigh index in order to verify if the instability is driving or damping at different frequency levels. Moreover, the predicted Helmholtz and Strouhal numbers are found to be in excellent agreement with the experimental data available in open literature for wide range of Re.

In a previous work (Li et al. Science 364, 170) [1], we proposed a heat transfer system that preserves the anti-parity-time (APT) symmetry, and observe the rest-to-motion phase transition during the symmetry breaking. Recently, it was suggested (Zhao et al. arXiv:1906.08431) [2] that the behaviours of the system can be understood using scaling analysis based on the Péclet and Nusselt numbers (Pe and Nu). It was further proposed that there exists a third regime in the phase diagram in addition to the symmetric and symmetry broken phases. Although we appreciate the proposal to characterize the contributions of coupling, diffusion, and advection with dimensionless numbers, here we show that they do not help to predict or interpret the behaviours of the APT system. The dimensionless numbers do not provide enough details about the system to conclude that there is a motionless phase, a phase transition, to find the critical point, or to give the correct phase diagram with only two regimes.

A theoretical analysis is presented for turbulent flows, applicable for canonical (channel, boundary-layer and free jet) geometries. Momentum and energy balance for a control volume moving at the local mean velocity decouples the fluctuation from the mean velocities, resulting in a symmetric set of transport equations for the Reynolds normal and shear stresses. In this formalism, gradients of the fluctuating velocities represent flux vectors, easily verifiable using the available DNS data. A derivative of this transport concept is the scaling for the Reynolds stresses in the dissipation space. Combining with the statistical energy distribution function, a full prescription of turbulent flows is enabled in the basic canonical geometries. Based on this theoretical foundation, more complex flow configurations may be addressed with far more efficient algorithms.

We present an experimental and theoretical study of capillary surfers, a new fluid-mediated active system that bridges the gap between dissipation- and inertia-dominated regimes. Surfers are wave-driven particles that self-propel and interact on a fluid interface via an extended field of surface waves. A surfer's speed and interaction with its environment can be tuned broadly through the particle, fluid, and vibration parameters. The wave nature of interactions among surfers allows for multistability of interaction modes and promises a number of novel collective behaviors.

In laboratory studies and numerical simulations, we observe clear signatures of unstable time-periodic solutions in a moderately turbulent quasi-two-dimensional flow. We validate the dynamical relevance of such solutions by demonstrating that turbulent flows in both experiment and numerics transiently display time-periodic dynamics when they shadow unstable periodic orbits (UPOs). We show that UPOs we computed are also statistically significant, with turbulent flows spending a sizable fraction of the total time near these solutions. As a result, the average rates of energy input and dissipation for the turbulent flow and frequently visited UPOs differ only by a few percent.

At high Reynolds number, the interaction between two vortex tubes leads to intense velocity gradients, which are at the heart of fluid turbulence. This vorticity amplification comes about through two different instability mechanisms of the initial vortex tubes, assumed anti-parallel and with a mirror plane of symmetry. At moderate Reynolds number, the tubes destabilize via a Crow instability, with the nonlinear development leading to strong flattening of the cores into thin sheets. These sheets then break down into filaments which can repeat the process. At higher Reynolds number, the instability proceeds via the elliptical instability, producing vortex tubes that are perpendicular to the original tube directions. In this work, we demonstrate that these same transition between Crow and Elliptical instability occurs at moderate Reynolds number when we vary the initial angle β between two straight vortex tubes. We demonstrate that when the angle between the two tubes is close to ?/2 , the interaction between tubes leads to the formation of thin vortex sheets. The subsequent breakdown of these sheets involves a twisting of the paired sheets, followed by the appearance of a localized cloud of small scale vortex structures. At smaller values of the angle β between the two tubes, the breakdown mechanism changes to an elliptic cascade-like mechanism. Whereas the interaction of two vortices depends on the initial condition, the rapid formation of fine-scales vortex structures appears to be a robust feature, possibly universal at very high Reynolds numbers.

Compressible vortex rings, usually formed at the open end of a shock tube, often show interesting phenomena during their formation, evolution, and propagation depending on the shock Mach number (Ms) and exit flow conditions. The Mach number of the translating compressible vortex rings (Mv) investigated so far in the literature is subsonic as, the shock tube pressure ratio (PR) considered is relatively low. In this numerical study we focus on low to high vortex ring Mach numbers (0.31 < Mv < 1.08) cases with a particular focus on very high Mv cases that are not been reported in experiments as, it is difficult to obtain in laboratory. Using hydrogen as a driver section gas inside the shock tube, a supersonic compressible vortex ring (Mv > 1) is obtained for first time. It is established that the SST k-{\omega} based DES turbulent model replicates the experimental observation better than the previously published results at different stages of development of the vortex ring. DES, which is an inbuilt hybrid of LES and RANS approaches is evoked that can automatically switch to the sub-grid scale (SGS) model in the LES regions (i.e. with different scale vortical structures) and to a RANS model in the rest of the region (i.e. where the grid spacing is greater than the turbulent length scale). The DES model can predict characteristics of the shear layer vortices as well as counter-rotating vortex rings (CRVRs) as reported in the experimental measurements. Formation of multiple triple points and the corresponding slip-stream shear layers and thus multiple CRVRs behind the primary vortex ring at different radial locations, in addition to the usual CRVRs, appears to be a unique characteristic for high Mach number vortex rings. For high PR, H2, case during formation stage, a vortex layer of reverse circulation (that of primary vortex ring) is formed (contd...)