Active Viscoelastic Matter: From Bacterial Drag Reduction to Turbulent Solids

A paradigm for internally driven matter is the active nematic liquid crystal, whereby the equations of a conventional nematic are supplemented by a minimal active stress that violates time-reversal symmetry. In practice, active fluids may have not only liquid-crystalline but also viscoelastic polymer degrees of freedom. Here we explore the resulting interplay by coupling an active nematic to a minimal model of polymer rheology. We find that adding a polymer can greatly increase the complexity of spontaneous flow, but can also have calming effects, thereby increasing the net throughput of spontaneous flow along a pipe (a ``drag-reduction'' effect). Remarkably, active turbulence can also arise after switching on activity in a sufficiently soft elastomeric solid.

Acoustic response of vortex breakdown modes in a coaxial isothermal unconfined swirling jet

The present experimental work is concerned with the study of amplitude dependent acoustic response of an isothermal coaxial swirling jet. The excitation amplitude is increased in five distinct steps at the burner's Helmholtz resonator mode (i.e., 100 Hz). Two flow states are compared, namely, sub-critical and super-critical vortex breakdown (VB) that occur before and after the critical conical sheet breakdown, respectively. The geometric swirl number is varied in the range 2.14-4.03. Under the influence of external pulsing, global response characteristics are studied based on the topological changes observed in time-averaged 2D flow field. These are obtained from high resolution 2D PIV (particle image velocimetry) in the longitudinal-mid plane. PIV results also illustrate the changes in the normalized vortex core coordinates (r(vcc)/(r(vcc))(0) (Hz), y(vcc)/(y(vcc))(0) (Hz)) of internal recirculation zone (IRZ). A strong forced response is observed at 100 Hz (excitation frequency) in the convectively unstable region which get amplified based on the magnitude of external forcing. The radial extent of this forced response region at a given excitation amplitude is represented by the acoustic response region (b). The topological placement of the responsive convectively unstable region is a function of both the intensity of imparted swirl (characterized by swirl number) and forcing amplitude. It is observed that for sub-critical VB mode, an increase in the excitation amplitude till a critical value shifts the vortex core centre (particularly, the vortex core moves downstream and radially outwards) leading to drastic fanning-out/widening of the IRZ. This is accompanied by similar to 30% reduction in the recirculation velocity of the IRZ. It is also observed that b < R (R: radial distance from central axis to outer shear layer-OSL). At super-critical amplitudes, the sub-critical IRZ topology transits back (the vortex core retracts upstream and radially inwards) and finally undergoes a transverse shrinkage ((r(vcc))/(r(vcc))(0 Hz) decreases by similar to 20%) when b >= R. In contrast, the vortex core of super-critical breakdown mode consistently spreads radially outwards and is displaced further downstream. Finally, the IRZ fans-out at the threshold excitation amplitude. However, the acoustic response region b is still less than R. This is explained based on the characteristic geometric swirl number (S-G) of the flow regimes. The super-critical flow mode with higher S-G (hence, higher radial pressure drop due to rotational effect which scales as Delta P similar to rho u theta(2) and acts inwards towards the center line) compared to sub-critical state imposes a greater resistance to the radial outward spread of b. As a result, the acoustic energy supplied to the super-critical flow mode increases the degree of acoustic response at the pulsing frequency and energizes its harmonics (evident from power spectra). As a disturbance amplifier, the stronger convective instability mode within the flow structure of super-critical VB causes the topology to widen/fan-out severely at threshold excitation amplitude. (C) 2015 AIP Publishing LLC.

Jet substructure and probes of CP violation in Vh production

We analyse the hVV (V = W, Z) vertex in a model independent way using Vh production. To that end, we consider possible corrections to the Standard Model Higgs Lagrangian, in the form of higher dimensional operators which parametrise the effects of new physics. In our analysis, we pay special attention to linear observables that can be used to probe CP violation in the same. By considering the associated production of a Higgs boson with a vector boson (W or Z), we use jet substructure methods to define angular observables which are sensitive to new physics effects, including an asymmetry which is linearly sensitive to the presence of CP odd effects. We demonstrate how to use these observables to place bounds on the presence of higher dimensional operators, and quantify these statements using a log likelihood analysis. Our approach allows one to probe separately the hZZ and hWW vertices, involving arbitrary combinations of BSM operators, at the Large Hadron Collider.

Associated jet and subjet rates in light-quark and gluon jet discrimination

We show that in studies of light quark- and gluon-initiated jet discrimination, it is important to include the information on softer reconstructed jets (associated jets) around a primary hard jet. This is particularly relevant while adopting a small radius parameter for reconstructing hadronic jets. The probability of having an associated jet as a function of the primary jet transverse momentum (PT) and radius, the minimum associated jet pi, and the association radius is computed up to next-to-double logarithmic accuracy (NDLA), and the predictions are compared with results from Herwig++, Pythia6 and Pythia8 Monte Carlos (MC). We demonstrate the improvement in quark-gluon discrimination on using the associated jet rate variable with the help of a multivariate analysis. The associated jet rates are found to be only mildly sensitive to the choice of parton shower and hadronization algorithms, as well as to the effects of initial state radiation and underlying event. In addition, the number of k(t) subjets of an anti-k(t) jet is found to be an observable that leads to a rather uniform prediction across different MC's, broadly being in agreement with predictions in NDLA, as compared to the often used number of charged tracks observable.

Behaviour of turbulent Prandtl/Schmidt number in compressible mixing layer

The behaviour of turbulent Prandtl/Schmidt number is explored through the model-free simulation results. It has been observed that compressibility affects the Reynolds scalar flux vectors. Reduced peak values are also observed for compressible convective Mach number mixing layer as compared with the incompressible convective Mach number counterpart, indicating a reduction in the mixing of enthalpy and species. Pr-t and Sc-t variations also indicate a reduction in mixing. It is observed that unlike the incompressible case, it is difficult to assign a constant value to these numbers due to their continuous variation in space. Modelling of Pr-t and Sc-t would be necessary to cater for this continuous spatial variation. However, the turbulent Lewis number is evaluated to be near unity for the compressible case, making it necessary to model only one of the Pr-t and Sc-t..

NUMERICAL STUDIES OF DYNAMO ACTION IN A TURBULENT SHEAR FLOW. I.

We perform numerical experiments to study the shear dynamo problem where we look for the growth of a large-scale magnetic field due to non-helical stirring at small scales in a background linear shear flow in previously unexplored parameter regimes. We demonstrate the large-scale dynamo action in the limit where the fluid Reynolds number (Re) is below unity while the magnetic Reynolds number (Rm) is above unity; the exponential growth rate scales linearly with shear, which is consistent with earlier numerical works. The limit of low Re is particularly interesting, as seeing the dynamo action in this limit would provide enough motivation for further theoretical investigations, which may focus attention on this analytically more tractable limit of Re < 1 compared to the more formidable limit of Re > 1. We also perform simulations in the regimes where (i) both (Re, Rm) < 1, and (ii) Re > 1 and Rm < 1, and compute all of the components of the turbulent transport coefficients (alpha(ij) and alpha(ij)) using the test-field method. A reasonably good agreement is observed between our results and the results of earlier analytical works in similar parameter regimes.

A model supersonic buried-nozzle jet: instability and acoustic wave scattering and the far-field sound

We consider sound source mechanisms involving the acoustic and instability modes of dual-stream isothermal supersonic jets with the inner nozzle buried within an outer shroud-like nozzle. A particular focus is scattering into radiating sound waves at the shroud lip. For such jets, several families of acoustically coupled instability waves exist, beyond the regular vortical Kelvin-Helmholtz mode, with different shapes and propagation characteristics, which can therefore affect the character of the radiated sound. In our model, the coaxial shear layers are vortex sheets while the incident acoustic disturbances are the propagating shroud modes. The Wiener-Hopf method is used to compute their scattering at the sharp shroud edge to obtain the far-field radiation. The resulting far-field directivity quantifies the acoustic efficiency of different mechanisms, which is particularly important in the upstream direction, where the results show that the scattered sound is more intense than that radiated directly by the shear-layer modes.

Effect of H-2/CO Composition on Extinction Strain Rates of Counterflow Syngas Flames

This work presents a detailed experimental and numerical investigation of the effect of H-2/CO composition on extinction characteristics of premixed and nonpremixed syngas flames. Experimental measurements of local and global extinction strain rates in counterflow diffusion flames have been reported at atmospheric pressure for six different compositions of syngas fuel. The concentration of H-2 was varied from 5 to 20% with a 3% increment, and correspondingly, CO was decreased from 35 to 20% in steps of 3%. Particle imaging velocimetry has been used to determine the local extinction strain rates. Local extinction strain rates increased with an increase in the H-2/CO ratio in both nonpremixed and premixed flames. The predicted extinction strain rates for both nonpremixed and premixed counterflow flames using five different mechanisms available in the literature were compared with measurements. The Davis H-2/CO and Ranzi H-2/CO mechanisms predicted extinction strain rates within 10% of experimental values irrespective of the H-2/CO ratio. In the nonpremixed case, the Cl mechanism by Li et al., GRI 3.0, and the Ranzi H-2/CO mechanism predicted extinction strain rates well for low H-2/CO ratios (from 5:35 to 14:26) but deviated from experiments for higher H-2/CO values (17:23 and 20:20). In addition to kinetics, preferential diffusion effects were found to affect the reaction zone significantly and create distinct localized reaction zone structures in nonpremixed flames, which could contribute to discrepancies in extinction predictions.

Turbulent states and their transitions in mathematical models for ventricular tissue: The effects of random interstitial fibroblasts

We study the dynamical behaviors of two types of spiral-and scroll-wave turbulence states, respectively, in two-dimensional (2D) and three-dimensional (3D) mathematical models, of human, ventricular, myocyte cells that are attached to randomly distributed interstitial fibroblasts; these turbulence states are promoted by (a) the steep slope of the action-potential-duration-restitution (APDR) plot or (b) early afterdepolarizations (EADs). Our single-cell study shows that (1) the myocyte-fibroblast (MF) coupling G(j) and (2) the number N-f of fibroblasts in an MF unit lower the steepness of the APDR slope and eliminate the EAD behaviors of myocytes; we explore the pacing dependence of such EAD suppression. In our 2D simulations, we observe that a spiral-turbulence (ST) state evolves into a state with a single, rotating spiral (RS) if either (a) G(j) is large or (b) the maximum possible number of fibroblasts per myocyte N-f(max) is large. We also observe that the minimum value of G(j), for the transition from the ST to the RS state, decreases as N-f(max) increases. We find that, for the steep-APDR-induced ST state, once the MF coupling suppresses ST, the rotation period of a spiral in the RS state increases as (1) G(j) increases, with fixed N-f(max), and (2) N-f(max) increases, with fixed G(j). We obtain the boundary between ST and RS stability regions in the N-f(max)-G(j) plane. In particular, for low values of N-f(max), the value of G(j), at the ST-RS boundary, depends on the realization of the randomly distributed fibroblasts; this dependence decreases as N-f(max) increases. Our 3D studies show a similar transition from scroll-wave turbulence to a single, rotating, scroll-wave state because of the MF coupling. We examine the experimental implications of our study and propose that the suppression (a) of the steep slope of the APDR or (b) EADs can eliminate spiral-and scroll-wave turbulence in heterogeneous cardiac tissue, which has randomly distributed fibroblasts.

COOL CORE CYCLES: COLD GAS AND AGN JET FEEDBACK IN CLUSTER CORES

Using high-resolution 3D and 2D (axisymmetric) hydrodynamic simulations in spherical geometry, we study the evolution of cool cluster cores heated by feedback-driven bipolar active galactic nuclei (AGNs) jets. Condensation of cold gas, and the consequent enhanced accretion, is required for AGN feedback to balance radiative cooling with reasonable efficiencies, and to match the observed cool core properties. A feedback efficiency (mechanical luminosity approximate to epsilon(M) over dot(acc)c(2); where (M) over dot(acc). is the mass accretion rate at 1 kpc) as small as 6 x 10(-5) is sufficient to reduce the cooling/accretion rate by similar to 10 compared to a pure cooling flow in clusters (with M-200 less than or similar to 7 x 10(14) M-circle dot). This value is much smaller compared to the ones considered earlier, and is consistent with the jet efficiency and the fact that only a small fraction of gas at 1 kpc is accreted onto the supermassive black hole (SMBH). The feedback efficiency in earlier works was so high that the cluster core reached equilibrium in a hot state without much precipitation, unlike what is observed in cool-core clusters. We find hysteresis cycles in all our simulations with cold mode feedback: condensation of cold gas when the ratio of the cooling-time to the free-fall time (t(cool)/t(ff)) is less than or similar to 10 leads to a sudden enhancement in the accretion rate; a large accretion rate causes strong jets and overheating of the hot intracluster medium such that t(cool)/t(ff) > 10; further condensation of cold gas is suppressed and the accretion rate falls, leading to slow cooling of the core and condensation of cold gas, restarting the cycle. Therefore, there is a spread in core properties, such as the jet power, accretion rate, for the same value of core entropy t(cool)/t(ff). A smaller number of cycles is observed for higher efficiencies and for lower mass halos because the core is overheated to a longer cooling time. The 3D simulations show the formation of a few-kpc scale, rotationally supported, massive (similar to 10(11) M-circle dot) cold gas torus. Since the torus gas is not accreted onto the SMBH, it is largely decoupled from the feedback cycle. The radially dominant cold gas (T < 5 x 10(4) K; vertical bar v(r)vertical bar >vertical bar v(phi vertical bar)) consists of fast cold gas uplifted by AGN jets and freely infalling cold gas condensing out of the core. The radially dominant cold gas extends out to 25 kpc for the fiducial run (halo mass 7 x 10(14) M-circle dot and feedback efficiency 6 x 10(-5)), with the average mass inflow rate dominating the outflow rate by a factor of approximate to 2. We compare our simulation results with recent observations.

The impact of latent heating on the location and strength of the tropical easterly jet

The tropical easterly jet (TEJ) is a prominent atmospheric circulation feature observed during the Asian summer monsoon. It is generally assumed that sensible heating over the Tibetan Plateau directly influences the location of the TEJ. However, other studies have suggested the importance of latent heating in determining the jet location. In this paper, the relative importance of latent heating on the maintenance of the TEJ is explored through simulations with a general circulation model. The simulation of the TEJ by the Community Atmosphere Model, version 3.1 is discussed in detail. These simulations showed that the location of the TEJ is well correlated with the location of the precipitation. Significant zonal shifts in the location of the precipitation resulted in similar shifts in the zonal location of the TEJ. These zonal shifts had minimal effect on the large-scale structure of the jet. Further, provided that precipitation patterns were relatively unchanged, orography did not directly impact the location of the TEJ. These changes were robust even with changes in the cumulus parameterization. This suggests the potential important role of latent heating in determining the location and structure of the TEJ. These results were used to explain the significant differences in the zonal location of the TEJ in the years 1988 and 2002. To understand the contribution of the latitudinal location of latent heating on the strength of the TEJ, aqua-planet simulations were carried out. It has been shown that for similar amounts of net latent heating, the jet is stronger when heating is in the higher tropical latitudes. This may partly explain the reason for the jet to be very strong during the JJA monsoon season.

Hierarchy of structure functions for passive scalars advected by turbulent flows

A hierarchical model is proposed for the joint moments of the passive scalar dissipation and the velocity dissipation in fluid turbulence. This model predicts that the joint probability density function (PDF) of the dissipations is a bivariate log-Poisson. An analytical calculation of the scaling exponents of structure functions of the passive scalar is carried out for this hierarchical model, showing a good agreement with the results of direct numerical simulations and experiments.

Direct Numerical Simulation of a Spatially Evolving Supersonic Turbulent Boundary Layer at Ma=6

Direct numerical simulation is carried out for a spatially evolving supersonic turbulent boundary layer at free-stream Mach number 6. To overcome numerical instability, the seventh-order WENO scheme is used for the convection terms of Navier-Stokes equations, and fine mesh is adopted to minimize numerical dissipation. Compressibilty effects on the near-wall turbulent kinetic energy budget are studied. The cross-stream extended self-similarity and scaling exponents including the near-wall region are studied. In high Mach number flows, the coherence vortex structures are arranged to be smoother and streamwised, and the hair-pin vortices are less likely to occur.

An experimental study on turbulent coherent structures near a sheared air-water interface

The turbulence structures near a sheared air-water interface were experimentally investigated with the hydrogen bubble visualization technique. Surface shear was imposed by an airflow over the water flow which was kept free from surface waves. Results show that the wind shear has the main influence on coherent structures under air-water interfaces. Low- and high- speed streaks form in the region close to the interface as a result of the imposed shear stress. When a certain airflow velocity is reached, "turbulent spots" appear randomly at low-speed streaks with some characteristics of hairpin vortices. At even higher shear rates, the flow near the interface is dominated primarily by intermittent bursting events. The coherent structures observed neat sheared air-water interfaces show qualitative similarities with those occurring in near-wall turbulence. However, a few distinctive phenomena were also observed, including the fluctuating thickness of the instantaneous boundary layer and vertical vortices in bursting processes, which appear to be associated with the characteristics of air-water interfaces.

Burst Detection in Turbulent Channel Flows Based on Large Eddy Simulation Databases

Reliable turbulent channel flow databases at several Reynolds numbers have been established by large eddy simulation (LES), with two of them validated by comparing with typical direct numerical simulation (DNS) results. Furthermore, the statistics, such as velocity profile, turbulent intensities and shear stress, were obtained as well as the temporal and spatial structure of turbulent bursts. Based on the LES databases available, the conditional sampling methods are used to detect the structures of burst events. A method to deterimine the grouping parameter from the probability distribution function (pdf) curve of the time separation between ejection events is proposed to avoid the errors in detected results. And thus, the dependence of average burst period on thresholds is considerably weakened. Meanwhile, the average burst-to-bed area ratios are detected. It is found that the Reynolds number exhibits little effect on the burst period and burst-to-bed area ratio.