Source author record

Bruno Eckhardt

Bruno Eckhardt appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

36works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

36 published item(s)

preprint2022arXiv

Interpreted machine learning in fluid dynamics: Explaining relaminarisation events in wall-bounded shear flows

Machine Learning (ML) is becoming increasingly popular in fluid dynamics. Powerful ML algorithms such as neural networks or ensemble methods are notoriously difficult to interpret. Here, we introduce the novel Shapley Additive Explanations (SHAP) algorithm (Lundberg & Lee, 2017), a game-theoretic approach that explains the output of a given ML model, in the fluid dynamics context. We give a proof of concept concerning SHAP as an explainable AI method providing useful and human-interpretable insight for fluid dynamics. To show that the feature importance ranking provided by SHAP can be interpreted physically, we first consider data from an established low-dimensional model based on the self-sustaining process (SSP) in wall-bounded shear flows, where each data feature has a clear physical and dynamical interpretation in terms of known representative features of the near-wall dynamics, i.e. streamwise vortices, streaks and linear streak instabilities. SHAP determines consistently that only the laminar profile, the streamwise vortex, and a specific streak instability play a major role in the prediction. We demonstrate that the method can be applied to larger fluid dynamics datasets by a SHAP evaluation on plane Couette flow in a minimal flow unit focussing on the relevance of streaks and their instabilities for the prediction of relaminarisation events. Here, we find that the prediction is based on proxies for streak modulations corresponding to linear streak instabilities within the SSP. That is, the SHAP analysis suggests that the break-up of the self-sustaining cycle is connected with a suppression of streak instabilities.

preprint2021arXiv

Purely elastic linear instabilities in parallel shear flows with free-slip boundary conditions

We perform a linear stability analysis of viscoelastic plane Couette and plane Poiseuille flows with free-slip boundary conditions. The fluid is described by the Oldroyd-B constitutive model, and the flows are driven by a suitable body force. We find that both types of flow become linearly unstable, and we characterise the spatial structure of the unstable modes. By performing a boundary condition homotopy from the free-slip to no-slip boundaries, we demonstrate that the unstable modes are directly related to the least stable modes of the no-slip problem, destabilised under the free-slip situation. We discuss how our observations can be used to study recently discovered purely elastic turbulence in parallel shear flows.

preprint2020arXiv

Ab initio study of time-dependent dynamics in strong-field triple ionization

An ab initio analysis of strong-field three-electron ionization in a restricted-dimensionality model reveals the dynamics of the ionization process and the dominant channels for double (DI) and triple ionization (TI). Simulations using wave functions that respect the Pauli principle show that the most likely channel is a sequence of single ionization (SI) and DI, while direct TI has a much lower probability. The dominant DI process has the highest probability for a singlet of up- and down-spin electrons. The results demonstrate the significance of the Pauli principle for the selection of dominant pathways in ionization and possibly other many-electron processes in strong fields.

preprint2020arXiv

Linear feedback control of invariant solutions in channel flow

Considering channel flow at Reynolds numbers below the linear stability threshold of the laminar profile as a generic example system showing a subcritical transition to turbulence connected with the existence of simple invariant solutions, we here discuss issues that arise in the application of linear feedback control of invariant solutions of the Navier-Stokes equations. We focus on the simplest possible problem, that is, travelling waves with one unstable direction. In view of potential experimental applicability we construct a pressure-based feedback strategy and study its effect on the stable, marginal and unstable directions of these solutions in different periodic cells. Even though the original instability can be removed, new instabilities emerge as the feedback procedure affects not only the unstable but also the stable directions. We quantify these adverse effects and discuss their implications for the design of successful control strategies. In order to highlight the challenges that arise in the application of feedback control methods in principle and concerning potential applications in the search for simple invariant solutions of the Navier-Stokes equations in particular, we consider an explicitly constructed analogue to closed-loop linear optimal control that leaves the stable directions unaffected.

preprint2020arXiv

Using Machine Learning to predict extreme events in the Hénon map

Machine Learning (ML) inspired algorithms provide a flexible set of tools for analyzing and forecasting chaotic dynamical systems. We here analyze the performance of one algorithm for the prediction of extreme events in the two-dimensional Hénon map at the classical parameters. The task is to determine whether a trajectory will exceed a threshold after a set number of time steps into the future. This task has a geometric interpretation within the dynamics of the Hénon map, which we use to gauge the performance of the neural networks that are used in this work. We analyze the dependence of the success rate of the ML models on the prediction time $T$ , the number of training samples $N_T$ and the size of the network $N_p$. We observe that in order to maintain a certain accuracy, $N_T \propto exp(2 h T)$ and $N_p \propto exp(hT)$, where $h$ is the topological entropy. Similar relations between the intrinsic chaotic properties of the dynamics and ML parameters might be observable in other systems as well.

preprint2019arXiv

Condensate formation and multiscale dynamics in two-dimensional active suspensions

The collective effects of microswimmers in active suspensions result in active turbulence, a spatiotemporally chaotic dynamics at mesoscale, which is characterized by the presence of vortices and jets at scales much larger than the characteristic size of the individual active constituents. To describe this dynamics, Navier-Stokes-based one-fluid models driven by small-scale forces have been proposed. Here, we provide a justification of such models for the case of dense suspensions in two dimensions (2d). We subsequently carry out an in-depth numerical study of the properties of one-fluid models as a function of the active driving in view of possible transition scenarios from active turbulence to large-scale pattern, referred to as condensate, formation induced by the classical inverse energy cascade in Newtonian 2d turbulence. Using a one-fluid model it was recently shown (Linkmann et al., Phys. Rev. Lett. (in press)) that two-dimensional active suspensions support two non-equilibrium steady states, one with a condensate and one without, which are separated by a subcritical transition. Here, we report further details on this transition such as hysteresis and discuss a low-dimensional model that describes the main features of the transition through nonlocal-in-scale coupling between the small-scale driving and the condensate.

preprint2019arXiv

Dynamic feedback control through wall suction in shear flows

Flow control is of interest in many open and wall-bounded shear flows in order to reduce drag or to avoid sudden large fluctuations that may lead to material failure. An established means of control is the application of suction through a porous wall. Here, we combine suction with a feedback strategy whereby the suction velocity is adjusted in response to either the kinetic energy or the shear stress at the bottom wall. The control procedure is then used in an attempt to stabilize invariant solutions and to carry out direct numerical simulations with a prescribed value of the friction coefficient.

preprint2019arXiv

Nonuniversal transition to condensate formation in two-dimensional turbulence

The occurrence of system-scale coherent structures, so-called condensates, is a well-known phenomenon in two-dimensional turbulence. Here, the transition to condensate formation is investigated as a function of the magnitude of the force and for different types of forcing. Random forces with constant mean energy input lead to a supercritical transition, while forcing through a small-scale linear instability results in a subcritical transition with bistability and hysteresis. That is, the transition to condensate formation in two-dimensional turbulence is nonuniversal. For the supercritical case we quantify the effect of large-scale friction on the value of the critical exponent and the location of the critical point.

preprint2016arXiv

Bypass transition and spot nucleation in boundary layers

The spatio-temporal aspects of the transition to turbulence are considered in the case of a boundary layer flow developing above a flat plate exposed to free-stream turbulence. Combining results on the receptivity to free-stream turbulence with the nonlinear concept of a transition threshold, a physically motivated model suggests a spatial distribution of spot nucleation events. To describe the evolution of turbulent spots a probabilistic cellular automaton is introduced, with all parameters directly fitted from numerical simulations of the boundary layer. The nucleation rates are then combined with the cellular automaton model, yielding excellent quantitative agreement with the statistical characteristics for different free-stream turbulence levels. We thus show how the recent theoretical progress on transitional wall-bounded flows can be extended to the much wider class of spatially developing boundary-layer flows.

preprint2016arXiv

Edge states as mediators of bypass transition in boundary-layer flows

The concept of edge state is investigated in the asymptotic suction boundary layer in relation with the receptivity process to noisy perturbations and the nucleation of turbulent spots. Edge tracking is first performed numerically, without imposing any discrete symmetry, in a large computational domain allowing for full spatial localisation of the perturbation velocity. The edge state is a three-dimensional localised structure recurrently characterised by a single low-speed streak that experiences erratic bursts and planar shifts. This recurrent streaky structure is then compared with predecessors of individual spot nucleation events, triggered by non-localised initial noise. The present results suggest a nonlinear picture, rooted in dynamical systems theory, of the nucleation process of turbulent spots in boundary-layer flows, in which the localised edge state play the role of state-space mediator.

preprint2016arXiv

Streamwise decay of localized states in channel flow

Channel flow, the pressure driven flow between parallel plates, has exact coherent structures that show various degrees of localization. For states which are localized in streamwise direction but extended in spanwise direction, we show that they are exponentially localized, with decay constants that are different on the upstream and downstream side. We extend the analysis of Brand and Gibson [J. Fluid Mech. 750, R1 (2014)] for stationary states to the case of advected structures that is needed here, and derive expressions for the decay in terms of eigenvalues and eigenfunctions of certain second order differential equations. The results are in very good agreement with observations on exact coherent structures of different transversal wave lengths.

preprint2016arXiv

Transition in the asymptotic suction boundary layer over a heated plate

The asymptotic suction boundary layer (ASBL) is a parallel shear flow that becomes turbulent in a bypass transition in parameter regions where the laminar profile is stable. We here add a temperature gradient perpendicular to the plate and explore the interaction between convection and shear in determining the transition. We find that the laminar state becomes unstable in a subcritical bifurcation and that the critical Rayleigh number and wavenumber depend strongly on the Prandtl number. We also track several secondary bifurcations and identify states that are localized in two directions, showing different symmetries. In the subcritical regime, transient turbulent states which are connected to exact coherent states and follow the same transition scenario as found in linearly stable shear flows are identified and analyzed. The study extends the bypass transition scenario from shear flows to thermal boundary layers and highlights the intricate interactions between thermal and shear forces.

preprint2015arXiv

Bypass transition and subcritical turbulence in plane Poiseuille flow

Plane Poiseuille flow shows turbulence at a Reynolds number that is lower than the critical one for the onset of Tollmien-Schlichting waves. The transition to turbulence follows the same route as the by-pass transition in boundary layers, i.e. finite amplitude perturbations are required and the flow is dominated by downstream vortices and streaks in the transitional regime. In order to relate the phenomenology in plane Poiseuille flow to our previous studies of plane Couette flow (Kreilos & Eckhardt, 2012), we study a symmetric subspace of plane Poiseuille flow in which the bifurcation cascade stands out clearly. By tracing the edge state, which in this system is a travelling wave, and its bifurcations, we can trace the formation of a chaotic attractor, the interior crisis that increase the phase space volume affected by the flow, and the ultimate transition into a chaotic saddle in a crisis bifurcation. After the boundary crisis we can observe transient chaos with exponentially distributed lifetimes.

preprint2015arXiv

Crisis bifurcations in plane Poiseuille flow

Many shear flows follow a route to turbulence that has striking similarities to bifurcation scenarios in low-dimensional dynamical systems. Among the bifurcations that appear, crisis bifurcations are important because they cause global transitions between open and closed attractors, or indicate drastic increases in the range of the state space that is covered by the dynamics. We here study exterior and interior crisis bifurcations in direct numerical simulations of transitional plane Poiseuille flow in a mirror-symmetric subspace. We trace the state space dynamics from the appearance of the first three-dimensional exact coherent structures to the transition from an attractor to a chaotic saddle in an exterior crisis. For intermediate Reynolds numbers, the attractor undergoes several interior crises, in which new states appear and intermittent behavior can be observed. The bifurcations contribute to increasing the complexity of the dynamics and to a more dense coverage of state space.

preprint2015arXiv

Direct and noisy transitions in a model shear flow

The transition to turbulence in flows where the laminar profile is linearly stable requires perturbations of finite amplitude. "Optimal" perturbations are distinguished as extrema of certain functionals, and different functionals give different optima. We here discuss the phase space structure of a 2-d simplified model of the transition to turbulence and discuss optimal perturbations with respect to three criteria: energy of the initial condition, energy dissipation of the initial condition and amplitude of noise in a stochastic transition. We find that the states that trigger the transition are different in the three cases, but show the same scaling with Reynolds number.

preprint2015arXiv

Momentum transport in Taylor-Couette flow with vanishing curvature

We numerically study turbulent Taylor-Couette flow (TCF) between two independently rotating cylinders and the transition to rotating plane Couette flow (RPCF) in the limit of infinite radii. By using the shear Reynolds number $Re_S$ and rotation number $R_Ω$ as dimensionless parameters, the transition from TCF to RPCF can be studied continuously without singularities. Already for radius ratios $η\geq0.9$ we find that the simulation results for various radius ratios and for RPCF collapse as a function of $R_Ω$, indicating a turbulent behaviour common to both systems. We observe this agreement in the torque, mean momentum transport, mean profiles, and turbulent fluctuations. Moreover, the central profiles in TCF and RPCF for $R_Ω>0$ are found to conform with inviscid neutral stability. Intermittent bursts that have been observed in the outer boundary layer and have been linked to the formation of a torque maximum for counter-rotation are shown to disappear as $η\rightarrow 1$. The corresponding torque maximum disappears as well. Instead, two new maxima of different origin appear for $η\geq0.9$ and RPCF, a broad and a narrow one, in contrast to the results for smaller $η$. The broad maximum at $R_Ω=0.2$ is connected with a strong vortical flow and can be reproduced by streamwise invariant simulations. The narrow maximum at $R_Ω=0.02$ only emerges with increasing $Re_S$ and is accompanied by an efficient and correlated momentum transport by the mean flow. Since the narrow maximum is of larger amplitude for $Re_S=2\cdot10^4$, our simulations suggest that it will dominate at even higher $Re_S$.

preprint2015arXiv

Order-disorder transitions in a sheared many body system

Motivated by experiments on sheared suspensions that show a transition between ordered and disordered phases, we here study the long-time behavior of a sheared and overdamped 2-d system of particles interacting by repulsive forces. As a function of interaction strength and shear rate we find transitions between phases with vanishing and large single-particle diffusion. In the phases with vanishing single-particle diffusion, the system evolves towards regular lattices, usually on very slow time scales. Different lattices can be approached, depending on interaction strength and forcing amplitude. The disordered state appears in parameter regions where the regular lattices are unstable. Correlation functions between the particles reveal the formation of shear bands. In contrast to single particle densities, the spatially resolved two-particle correlation functions vary with time and allow to determine the phase within a period. As in the case of the suspensions, motion in the state with low diffusivity is essentially reversible, whereas in the state with strong diffusion it is not.

preprint2015arXiv

Turbulence transition in the asymptotic suction boundary layer

We study the transition to turbulence in the asymptotic suction boundary layer (ASBL) by direct numerical simulation. Tracking the motion of trajectories intermediate between laminar and turbulent states we can identify the invariant object inside the laminar-turbulent boundary, the edge state. In small domains, the flow behaves like a travelling wave over short time intervals. On longer times one notes that the energy shows strong bursts at regular time intervals. During the bursts the streak structure is lost, but it reforms, translated in the spanwise direction by half the domain size. Varying the suction velocity allows to embed the flow into a family of flows that interpolate between plane Couette flow and the ASBL. Near the plane Couette limit, the edge state is a travelling wave. Increasing the suction, the travelling wave and a symmetry-related copy of it undergo a saddle-node infinite-period (SNIPER) bifurcation that leads to bursting and discrete-symmetry shifts. In wider domains, the structures localize in the spanwise direction, and the flow in the active region is similar to the one in small domains. There are still periodic bursts at which the flow structures are shifted, but the shift-distance is no longer connected to a discrete symmetry of the flow geometry. Two different states are found by edge tracking techniques, one where structures are shifted to the same side at every burst and one where they are alternatingly shifted to the left and to the right.

preprint2014arXiv

Groebner basis methods for stationary solutions of a low-dimensional model for a shear flow

We use Groebner basis methods to extract all stationary solutions for the 9-mode shear flow model that is described in Moehlis et al, New J. Phys. 6, 54 (2004). Using rational approximations to irrational wave numbers and algebraic manipulation techniques we reduce the problem of determining all stationary states to finding roots of a polynomial of order 30. The coefficients differ by 30 powers of 10 so that algorithms for extended precision are needed to extract the roots reliably. We find that there are eight stationary solutions consisting of two distinct states that each appear in four symmetry-related phases. We discuss extensions of these results for other flows.

preprint2014arXiv

Long wavelength instability of coherent structures in plane Couette flow

We study the stability of coherent structures in plane Couette flow against long-wavelength perturbations in wide domains that cover several pairs of coherent structures. For one and two pairs of vortices, the states retain the stability properties of the small domains, but for three pairs new unstable modes are found. They are shown to be connected to bifurcations that break the translational symmetry and drive the coherent structures from the spanwise extended state to a modulated one that is a precursor to spanwise localized states. Tracking the stability of the orbits as functions of the spanwise wave length reveals a rich variety of additional bifurcations.

preprint2014arXiv

Periodically bursting edge states in plane Poiseuille flow

We investigate the laminar-turbulent boundary in plane Poiseuille flow by the method of edge tracking. In short and narrow computational domains we find for a wide range in Reynolds number that all states in the boundary converge to a period orbit with a period of the order of $10^{3}$ time units. The attracting states in these small domains are periodically extended in the spanwise and streamwise direction, but always localized to one side of the channel in the normal direction. In short and wide domains the edge states are localized in the spanwise direction. The periodic motion found in the small domains then induces a large variety of dynamical activity. The findings are very similar to the ones in the asymptotic suction boundary layer.

preprint2014arXiv

Streamwise and doubly-localised periodic orbits in plane Poiseuille flow

We study localised exact coherent structures in plane Poiseuille flow that are relative periodic orbits. They are obtained from extended states in smaller, periodically continued domains, by increasing the length to obtain streamwise localization and then by increasing the width to achieve spanwise localisation. The states maintain the travelling wave structure of the extended states, which is then modulated by a localised envelope on larger scales. In streamwise direction, the envelope shows exponential localization, with different exponents on the upstream and downstream side. The upstream exponent increases linearly with Reynolds number Re, but the downstream exponent is essentially independent of Re. In the spanwise direction the decay is compatible with a power-law localisation. As the width increases the localised state undergoes further bifurcations which add additional unstable directions, so that the edge state the relative attractor on the boundary between the laminar and turbulent motions, in the system becomes chaotic.

preprint2014arXiv

Symmetry related dynamics in parallel shear flows

Parallel shear flows come with continuous symmetries of translation in the downstream and spanwise direction. As a consequence, flow states that differ in their spanwise or downstream location but are otherwise identical are dynamically equivalent. In the case of travelling waves, this trivial degree of freedom can be removed by going to a frame of reference that moves with the state, thereby turning the travelling wave in the laboratory frame to a fixed point in the comoving frame of reference. We here discuss a general method by which the translational displacements can be removed also for more complicated and dynamically active states and demonstrate its application for several examples. For flows states in the asymptotic suction boundary layer we show that in the case of the long-period oscillatory edge state we can find local phase speeds which remove the fast oscillations and reveal the slow vortex dynamics underlying the burst phenomenon. For spanwise translating states we show that the method removes the drift but not the dynamical events that cause the big spanwise displacement. For a turbulent case we apply the method to the spanwise shifts and find slow components that are correlated over very long times. Calculations for plane Poiseuille flow show that the long correlations in the transverse motions are not special to the asymptotic suction boundary layer.

preprint2014arXiv

Turbulent states in plane Couette flow with rotation

Shearing and rotational forces in fluids can significantly alter the transport of momentum.A numerical investigation was undertaken to study the role of these forces using plane Couette flow subject to rotation about an axis perpendicular to both wall-normal and streamwise directions. Using a set of progressively higher Reynolds numbers up to Re = 5200, we find that the torque for a given Re is a non-monotonic function of rotation number, Ro. For low-to-moderate turbulent Reynolds numbers we find a maximum that is associated with flow fields that are dominated by downstream vortices and calculations of 2-d vortices capture the maximum also quantitatively. For higher shear Reynolds numbers a second stronger maximum emerges at smaller rotation numbers, closer to non-rotating plane Couette flow. It is carried by flows with a markedly 3-d structure and cannot be captured by 2-d vortex studies. As the Reynolds number increases, this maximum becomes stronger and eventually overtakes the one associated with the 2-d flow state.

preprint2013arXiv

Complexity of localised coherent structures in a boundary-layer flow

We study numerically transitional coherent structures in a boundary-layer flow with homogeneous suction at the wall (the so-called asymptotic suction boundary layer ASBL). The dynamics restricted to the laminar-turbulent separatrix is investigated in a spanwise-extended domain that allows for robust localisation of all edge states. We work at fixed Reynolds number and study the edge states as a function of the streamwise period. We demonstrate the complex spatio-temporal dynamics of these localised states, which exhibits multistability and undergoes complex bifurcations leading from periodic to chaotic regimes. It is argued that in all regimes the dynamics restricted to the edge is essentially low-dimensional and non-extensive.

preprint2013arXiv

Edge states for the turbulence transition in the asymptotic suction boundary layer

We demonstrate the existence of an exact invariant solution to the Navier-Stokes equations for the asymptotic suction boundary layer. The identified periodic orbit with a very long period of several thousand advective time units is found as a local dynamical attractor embedded in the stability boundary between laminar and turbulent dynamics. Its dynamics captures both the interplay of downstream oriented vortex pairs and streaks observed in numerous shear flows as well as the energetic bursting that is characteristic for boundary layers. By embedding the flow into a family of flows that interpolates between plane Couette flow and the boundary layer we demonstrate that the periodic orbit emerges in a saddle-node infinite-period (SNIPER) bifurcation of two symmetry-related travelling wave solutions of plane Couette flow. Physically, the long period is due to a slow streak instability which leads to a violent breakup of a streak associated with the bursting and the reformation of the streak at a different spanwise location. We show that the orbit is structurally stable when varying both the Reynolds number and the domain size.

preprint2013arXiv

Increasing lifetimes and the growing saddles of shear flow turbulence

In linearly stable shear flows turbulence spontaneously decays with a characteristic lifetime that varies with Reynolds number. The lifetime sharply increases with Reynolds number so that a possible divergence marking the transition to sustained turbulence at a critical point has been discussed. We present a mechanism by which the lifetimes increase: in the system's state space, turbulent motion is supported by a chaotic saddle. Inside this saddle a locally attracting periodic orbit is created and undergoes a traditional bifurcation sequence generating chaos. The formed new 'turbulent bubble' is initially an attractor supporting persistent chaotic dynamics. Soon after its creation it collides with its own boundary, by which it becomes leaky and dynamically connected with the surrounding structures. The complexity of the chaotic saddle that supports transient turbulence hence increases by incorporating the remnant of a new bubble. As a a result, the time it takes for a trajectory to leave the saddle and decay to the laminar state is increased. We demonstrate this phenomenon in plane Couette flow and show that characteristic lifetimes vary non-smoothly and non-monotonically with Reynolds number.

preprint2013arXiv

Periodic orbits near onset of chaos in plane Couette flow

We track the secondary bifurcations of coherent states in plane Couette flow and show that they undergo an incomplete periodic doubling cascade that ends with a crisis bifurcation. We introduce a symbolic dynamics for the orbits and show that the ones that exist fall into the universal sequence described by Metropolis, Stein and Stein for unimodal maps. The periodic orbits cover much of the turbulent dynamics in that their temporal evolution overlaps with turbulent motions when projected onto a plane spanned by energy production and dissipation.

preprint2012arXiv

Direct Numerical Simulations of Local and Global Torque in Taylor-Couette Flow up to Re=30.000

The torque in turbulent Taylor-Couette flows for shear Reynolds numbers Re_S up to 3x10^4 at various mean rotations is studied by means of direct numerical simulations for a radius ratio of η=0.71. Convergence of simulations is tested using three criteria of which the agreement of dissipation values estimated from the torque and from the volume dissipation rate turns out to be most demanding. We evaluate the influence of Taylor vortex heights on the torque for a stationary outer cylinder and select a value of the aspect ratio of Γ=2, close to the torque maximum. The connection between the torque and the transverse current J^ωof azimuthal motion which can be computed from the velocity field enables us to investigate the local transport resulting in the torque. The typical spatial distribution of the individual convective and viscous contributions to the local current is analysed for a turbulent flow case. To characterise the turbulent statistics of the transport, PDF's of local current fluctuations are compared to experimental wall shear stress measurements. PDF's of instantaneous torques reveal a fluctuation enhancement in the outer region for strong counter-rotation. Moreover, we find for simulations realising the same shear Re_S>=2x10^4 the formation of a torque maximum for moderate counter-rotation with angular velocities ω_o\approx-0.4ω_i. In contrast, for Re_S<=4x10^3 the torque features a maximum for a stationary outer cylinder. In addition, the effective torque scaling exponent is shown to also depend on the mean rotation state. Finally, we evaluate a close connection between boundary-layer thicknesses and the torque.

preprint2012arXiv

Intermittent boundary layers and torque maxima in Taylor-Couette flow

Turbulent Taylor-Couette flow between counter-rotating cylinders develops intermittently fluctuating boundary layers for sufficient counter-rotation. We demonstrate the phenomenon in direct numerical simulations for radius ratios η=0.5 and 0.71 and propose a theoretical model for the critical value in the rotation ratio. Numerical results as well as experiments show that the onset of this intermittency coincides with the maximum in torque. The variations in torque correlate with the variations in mean Taylor vortex flow which is first enhanced for weak counter-rotation, and then reduced as intermittency sets in. To support the model, we compare to numerical results, experiments at higher Reynolds numbers, and to Wendt's data.

preprint2011arXiv

From travelling waves to mild chaos: a supercritical bifurcation cascade in pipe flow

We study numerically a succession of transitions in pipe Poiseuille flow that leads from simple travelling waves to waves with chaotic time-dependence. The waves at the origin of the bifurcation cascade possess a shift-reflect symmetry and are both axially and azimuthally periodic with wave numbers κ = 1.63 and n = 2, respectively. As the Reynolds number is increased, successive transitions result in a wide range of time dependent solutions that includes spiralling, modulated-travelling, modulated-spiralling, doubly-modulated-spiralling and mildly chaotic waves. We show that the latter spring from heteroclinic tangles of the stable and unstable invariant manifolds of two shift-reflect-symmetric modulated-travelling waves. The chaotic set thus produced is confined to a limited range of Reynolds numbers, bounded by the occurrence of manifold tangencies. The states studied here belong to a subspace of discrete symmetry which makes many of the bifurcation and path-following investigations presented technically feasible. However, we expect that most of the phenomenology carries over to the full state-space, thus suggesting a mechanism for the formation and break-up of invariant states that can sustain turbulent dynamics.

preprint2009arXiv

Extreme vorticity growth in Navier-Stokes turbulence

According to statistical turbulence theory, the ensemble averaged squared vorticity rho_E is expected to grow not faster than drho_E/dt ~ rho_E^{3/2}. Solving a variational problem for maximal bulk enstrophy (E) growth, velocity fields were found for which the growth rate is as large as dE/dt ~ E^3. Using numerical simulations with well resolved small scales and a quasi-Lagrangian advection to track fluid subvolumes with rapidly growing vorticity, we study spatially resolved statistics of vorticity growth. We find that the volume ensemble averaged growth bound is satisfied locally to a remarkable degree of accuracy. Elements with dE/dt ~ E^3 can also be identified, but their growth tends to be replaced by the ensemble-averaged law when the intensities become too large.

preprint2009arXiv

Localised edge states nucleate turbulence in extended plane Couette cells

We study the turbulence transition of plane Couette flow in large domains where localised perturbations are observed to generate growing turbulent spots. Extending previous studies on the boundary between laminar and turbulent dynamics we determine invariant structures intermediate between laminar and turbulent flow. In wide but short domains we find states that are localised in spanwise direction, and in wide and long domains the states are also localised in downstream direction. These localised states act as critical nuclei for the transition to turbulence in spatially extended domains.

preprint2004arXiv

Homogeneous shear turbulence

The addition of suitable volume forces to the Navier-Stokes equation allows to simulate flows in the presence of a homogeneous shear. Because of the explicit form of the driving the flows are accessible to rigorous mathematical treatment and to accurate quantitative modelling in their global properties. The statistics of the fluctuations provide insight into generic behaviour of non-equilibrium systems and into the presistence of anisotropies at small scales. Correlation functions can be used to identify dominant large scale dynamical processes that are relevant for most of the momentum transport across the shear. The numerical studies of homogeneous shear flows complement analytical and experimental investigations and contribute to bridging the gap between ideal homogeneous isotropic turbulence and the more realistic heterogeneous turbulence.

preprint2004arXiv

Non-sequential double ionization of molecules

Double ionization of diatomic molecules by short linearly polarized laser pulses is analyzed. We consider the final stage of the ionization process, that is the decay of a highly excited two electron molecule, which is formed after re-scattering. The saddles of the effective adiabatic potential energy close to which simultaneous escape of electrons takes place are identified. Numerical simulations of the ionization of molecules show that the process can be dominated by either sequential or non-sequential events. In order to increase the ratio of non-sequential to sequential ionizations very short laser pulses should be applied.