Transition to magnetorotational turbulence in Taylor-Couette flow with imposed azimuthal magnetic field

The magnetorotational instability (MRI) is thought to be a powerful source of turbulence and momentum transport in astrophysical accretion discs, but obtaining observational evidence of its operation is challenging. Recently, laboratory experi-ments of Taylor–Couette flow with externally imposed axial and azimuthal magnetic fields have revealed the kinematic and dynamic properties of the MRI close to the instability onset. While good agreement was found with linear stability analyses, little is known about the transition to turbulence and transport properties of the MRI. We here report on a numerical investigation of the MRI with an imposed azimuthal magnetic field. We show that the laminar Taylor–Couette flow becomes unstable to a wave rotating in the azimuthal direction and standing in the axial direction via a supercritical Hopf bifurcation. Subsequently, the flow features a catastrophic transition to spatio-temporal defects which is mediated by a subcritical subharmonic Hopf bifurcation. Our results are in qualitative agreement with the PROMISE ex-periment and dramatically extend their realizable parameter range. We find that as the Reynolds number increases defects accumulate and grow into turbulence, yet the momentum transport scales weakly

Dynamo action in a quasi-Keplerian Taylor-Couette flow

We numerically compute the flow of an electrically conducting fluid in a Taylor-Couette geometry where the rotation rates of the inner and outer cylinders satisfy \Omega_o/\Omega_i=(r_o/r_i)-3/2. In this quasi-Keplerian regime a non-magnetic system would be Rayleigh-stable for all Reynolds numbers Re, and the resulting purely azimuthal flow incapable of kinematic dynamo action for all magnetic Reynolds numbers Rm. For Re=10^4 and Rm=10^5 we demonstrate the existence of a finite-amplitude dynamo, whereby a suitable initial condition yields mutually sustaining turbulence and magnetic fields, even though neither could exist without the other. This dynamo solution results in significantly increased outward angular momentum transport, with the bulk of the transport being by Maxwell rather than Reynolds stresses

The rise of fully turbulent flow

Over a century of research into the origin of turbulence in wallbounded shear flows has resulted in a puzzling picture in which turbulence appears in a variety of different states competing with laminar background flow. At slightly higher speeds the situation changes distinctly and the entire flow is turbulent. Neither the origin of the different states encountered during transition, nor their front dynamics, let alone the transformation to full turbulence could be explained to date. Combining experiments, theory and computer simulations here we uncover the bifurcation scenario organising the route to fully turbulent pipe flow and explain the front dynamics of the different states encountered in the process. Key to resolving this problem is the interpretation of the flow as a bistable system with nonlinear propagation (advection) of turbulent fronts. These findings bridge the gap between our understanding of the onset of turbulence and fully turbulent flows.Comment: 31 pages, 9 figure

Stabilisation and drag reduction of pipe flows by flattening the base profile

Recent experimental observations (Kuehnen et al., 2018) have shown that flattening a turbulent streamwise velocity profile in pipe flow destabilises the turbulence so that the flow relaminarises. We show that a similar phenomenon exists for laminar pipe flow profiles in the sense that the nonlinear stability of the laminar state is enhanced as the profile becomes more flattened. Significant drag reduction is also observed for the turbulent flow when triggered by sufficiently large disturbances. The flattening is produced by an artificial body force designed to mimick a baffle used in the experiments of Kuehnen et al. (2018) and the nonlinear stability measured by the size of the energy of the initial perturbations needed to trigger transition. In order to make the latter computation more efficient, we examine how indicative the minimal seed for transition is in measuring transition thresholds. We first show that the minimal seed is relatively robust to base profile changes and spectral filtering. We then compare the (unforced) transition behaviour of the minimal seed with several forms of randomised initial conditions in the range of Reynolds numbers Re=2400 to 10000 and find that the energy of the minimal seed after the Orr and oblique phases of its evolution is close to that of a localised random disturbance. In this sense, the minimal seed at the end of the oblique phase can be regarded as a good proxy for typical disturbances (here taken to be the localised random ones) and is thus used as initial condition in the simulations with the body force. The enhanced nonlinear stability and drag reduction predicted in the present study are an encouraging first step in modelling the experiments of Kuehnen et al. and should motivate future developments to fully exploit the benefits of this promising direction for flow control

Designing a minimal baffle to destabilise turbulence in pipe flows

Motivated by the results of recent experiments (K\"uhnen et al., Flow Turb. Combust., vol. 100, 2018, pp. 919-943), we consider the problem of designing a baffle (an obstacle to the flow) to relaminarise turbulence in pipe flows. Modelling the baffle as a spatial distribution of linear drag $\mathbf{f}(\mathbf{x},t)=-\chi(\mathbf{x})\mathbf{u}_{tot}(\mathbf{x},t)$ within the flow ($\mathbf{u}_{tot}$ is the total velocity field and $\chi \ge 0$ a scalar field), two different optimisation problems are considered to design $\chi$ at a Reynolds number $Re=3000$. In the first, the smallest baffle defined in terms of a $L_1$ norm of $\chi$ is sought which minimises the viscous dissipation rate of the flow. In the second, a baffle which minimises the total energy consumption of the flow is treated. Both problems indicate that the baffle should be axisymmetric and radially localised near the pipe wall, but struggle to predict the optimal streamwise extent. A manual search finds an optimal baffle one radius long which is then used to study how the amplitude for relaminarisation varies with $Re$ up to $15\,000$. Large stress reduction is found at the pipe wall, but at the expense of an increased pressure drop across the baffle. Estimates are then made of the break even point downstream of the baffle where the stress reduction at the wall due to the relaminarised flow compensates for the extra drag produced by the baffle.EPSR

Order-of-magnitude speedup for steady states and traveling waves via Stokes preconditioning in Channelflow and Openpipeflow

Steady states and traveling waves play a fundamental role in understanding hydrodynamic problems. Even when unstable, these states provide the bifurcation-theoretic explanation for the origin of the observed states. In turbulent wall-bounded shear flows, these states have been hypothesized to be saddle points organizing the trajectories within a chaotic attractor. These states must be computed with Newton's method or one of its generalizations, since time-integration cannot converge to unstable equilibria. The bottleneck is the solution of linear systems involving the Jacobian of the Navier-Stokes or Boussinesq equations. Originally such computations were carried out by constructing and directly inverting the Jacobian, but this is unfeasible for the matrices arising from three-dimensional hydrodynamic configurations in large domains. A popular method is to seek states that are invariant under numerical time integration. Surprisingly, equilibria may also be found by seeking flows that are invariant under a single very large Backwards-Euler Forwards-Euler timestep. We show that this method, called Stokes preconditioning, is 10 to 50 times faster at computing steady states in plane Couette flow and traveling waves in pipe flow. Moreover, it can be carried out using Channelflow (by Gibson) and Openpipeflow (by Willis) without any changes to these popular spectral codes. We explain the convergence rate as a function of the integration period and Reynolds number by computing the full spectra of the operators corresponding to the Jacobians of both methods.Comment: in Computational Modelling of Bifurcations and Instabilities in Fluid Dynamics, ed. Alexander Gelfgat (Springer, 2018

Search For Heavy Pointlike Dirac Monopoles

We have searched for central production of a pair of photons with high transverse energies in $p\bar p$ collisions at $\sqrt{s} = 1.8$ TeV using $70 pb^{-1}$ of data collected with the D\O detector at the Fermilab Tevatron in 1994--1996. If they exist, virtual heavy pointlike Dirac monopoles could rescatter pairs of nearly real photons into this final state via a box diagram. We observe no excess of events above background, and set lower 95% C.L. limits of $610, 870, or 1580 GeV/c^2$ on the mass of a spin 0, 1/2, or 1 Dirac monopole.Comment: 12 pages, 4 figure

Search for High Mass Photon Pairs in p-pbar --> gamma-gamma-jet-jet Events at sqrt(s)=1.8 TeV

A search has been carried out for events in the channel p-barp --> gamma gamma jet jet. Such a signature can characterize the production of a non-standard Higgs boson together with a W or Z boson. We refer to this non-standard Higgs, having standard model couplings to vector bosons but no coupling to fermions, as a "bosonic Higgs." With the requirement of two high transverse energy photons and two jets, the diphoton mass (m(gamma gamma)) distribution is consistent with expected background. A 90(95)% C.L. upper limit on the cross section as a function of mass is calculated, ranging from 0.60(0.80) pb for m(gamma gamma) = 65 GeV/c^2 to 0.26(0.34) pb for m(gamma gamma) = 150 GeV/c^2, corresponding to a 95% C.L. lower limit on the mass of a bosonic Higgs of 78.5 GeV/c^2.Comment: 9 pages, 3 figures. Replacement has new H->gamma gamma branching ratios and corresponding new mass limit

Ratio of the Isolated Photon Cross Sections at \sqrt{s} = 630 and 1800 GeV

The inclusive cross section for production of isolated photons has been measured in \pbarp collisions at $\sqrt{s} = 630$ GeV with the \D0 detector at the Fermilab Tevatron Collider. The photons span a transverse energy ($E_T$) range from 7-49 GeV and have pseudorapidity $|\eta| < 2.5$. This measurement is combined with to previous \D0 result at $\sqrt{s} = 1800$ GeV to form a ratio of the cross sections. Comparison of next-to-leading order QCD with the measured cross section at 630 GeV and ratio of cross sections show satisfactory agreement in most of the $E_T$ range.Comment: 7 pages. Published in Phys. Rev. Lett. 87, 251805, (2001

Molecular Interactions of Prodiginines with the BH3 Domain of Anti-Apoptotic Bcl-2 Family Members

Prodigiosin and obatoclax, members of the prodiginines family, are small molecules with anti-cancer properties that are currently under preclinical and clinical trials. The molecular target(s) of these agents, however, is an open question. Combining experimental and computational techniques we find that prodigiosin binds to the BH3 domain in some BCL-2 protein families, which play an important role in the apoptotic programmed cell death. In particular, our results indicate a large affinity of prodigiosin for MCL-1, an anti-apoptotic member of the BCL-2 family. In melanoma cells, we demonstrate that prodigiosin activates the mitochondrial apoptotic pathway by disrupting MCL-1/BAK complexes. Computer simulations with the PELE software allow the description of the induced fit process, obtaining a detailed atomic view of the molecular interactions. These results provide new data to understand the mechanism of action of these molecules, and assist in the development of more specific inhibitors of anti-apoptotic BCL-2 proteins.Spanish government and the European Union (FIS-PI10/00338) and from the ERC-2009-Adg 25027-PELE European project