Source author record

Tobias Grafke

Tobias Grafke 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

14works
10topics
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

14 published item(s)

preprint2022arXiv

Dynamical landscape of transitional pipe flow

The transition to turbulence in pipes is characterized by a coexistence of laminar and turbulent states. At the lower end of the transition, localized turbulent pulses, called puffs, can be excited. Puffs can decay when rare fluctuations drive them close to an edge state lying at the phase-space boundary with laminar flow. At higher Reynolds numbers, homogeneous turbulence can be sustained, and dominates over laminar flow. Here we complete this landscape of localized states, placing it within a unified bifurcation picture. We demonstrate our claims within the Barkley model, and motivate them generally. Specifically, we suggest the existence of an antipuff and a gap-edge -- states which mirror the puff and related edge state. Previously observed laminar gaps forming within homogeneous turbulence are then naturally identified as antipuffs nucleating and decaying through the gap edge.

preprint2022arXiv

Mechanism for turbulence proliferation in subcritical flows

The subcritical transition to turbulence, as occurs in pipe flow, is believed to generically be a phase transition in the directed percolation universality class. At its heart is a balance between the decay rate and proliferation rate of localized turbulent structures, called puffs in pipe flow. Here we propose the first-ever dynamical mechanism for puff proliferation -- the process by which a puff splits into two. In the first stage of our mechanism, a puff expands into a slug. In the second stage, a laminar gap is formed within the turbulent core. The notion of a split-edge state, mediating the transition from a single puff to a two puff state, is introduced and its form is predicted. The role of fluctuations in the two stages of the transition, and how splits could be suppressed with increasing Reynolds number, are discussed. Using numerical simulations, the mechanism is validated within the stochastic Barkley model. Concrete predictions to test the proposed mechanism in pipe and other wall bounded flows, and implications for the universality of the directed percolation picture, are discussed.

preprint2022arXiv

Spontaneous Symmetry Breaking for Extreme Vorticity and Strain in the 3D Navier-Stokes Equations

We investigate the spatio-temporal structure of the most likely configurations realising extremely high vorticity or strain in the stochastically forced 3D incompressible Navier-Stokes equations. Most likely configurations are computed by numerically finding the highest probability velocity field realising an extreme constraint as solution of a large optimisation problem. High-vorticity configurations are identified as pinched vortex filaments with swirl, while high-strain configurations correspond to counter-rotating vortex rings. We additionally observe that the most likely configurations for vorticity and strain spontaneously break their rotational symmetry for extremely high observable values. Instanton calculus and large deviation theory allow us to show that these maximum likelihood realisations determine the tail probabilities of the observed quantities. In particular, we are able to demonstrate that artificially enforcing rotational symmetry for large strain configurations leads to a severe underestimate of their probability, as it is dominated in likelihood by an exponentially more likely symmetry broken vortex-sheet configuration.

preprint2022arXiv

Symmetries and zero modes in sample path large deviations

Sharp large deviation estimates for stochastic differential equations with small noise, based on minimizing the Freidlin-Wentzell action functional under appropriate boundary conditions, can be obtained by integrating certain matrix Riccati differential equations along the large deviation minimizers or instantons, either forward or backward in time. Previous works in this direction often rely on the existence of isolated minimizers with positive definite second variation. By adopting techniques from field theory and explicitly evaluating the large deviation prefactors as functional determinant ratios using Forman's theorem, we extend the approach to general systems where degenerate submanifolds of minimizers exist. The key technique for this is a boundary-type regularization of the second variation operator. This extension is particularly relevant if the system possesses continuous symmetries that are broken by the instantons. We find that removing the vanishing eigenvalues associated with the zero modes is possible within the Riccati formulation and amounts to modifying the initial or final conditions and evaluation of the Riccati matrices. We apply our results in multiple examples including a dynamical phase transition for the average surface height in short-time large deviations of the one-dimensional Kardar-Parisi-Zhang equation with flat initial profile.

preprint2021arXiv

Extreme events and instantons in Lagrangian passive scalar turbulence models

The advection and mixing of a scalar quantity by fluid flow is an important problem in engineering and natural sciences. If the fluid is turbulent, the statistics of the passive scalar exhibit complex behavior. This paper is concerned with two Lagrangian scalar turbulence models based on the recent fluid deformation model that can be shown to reproduce the statistics of passive scalar turbulence for a range of Reynolds numbers. For these models, we demonstrate how events of extreme passive scalar gradients can be recovered by computing the instanton, i.e., the saddle-point configuration of the associated stochastic field theory. It allows us to both reproduce the heavy-tailed statistics associated with passive scalar turbulence, and recover the most likely mechanism leading to such extreme events. We further demonstrate that events of large negative strain in these models undergo spontaneous symmetry breaking.

preprint2015arXiv

Efficient Computation of Instantons for Multi-Dimensional Turbulent Flows with Large Scale Forcing

Extreme events play a crucial role in fluid turbulence. Inspired by methods from field theory, these extreme events, their evolution and probability can be computed with help of the instanton formalism as minimizers of a suitable action functional. Due to the high number of degrees of freedom in multi-dimensional fluid flows, traditional global minimization techniques quickly become prohibitive in their memory requirements. We outline a novel method for finding the minimizing trajectory in a wide class of problems that typically occurs in turbulence setups, where the underlying dynamical system is a non-gradient, non-linear partial differential equation, and the forcing is restricted to a limited length scale. We demonstrate the efficiency of the algorithm in terms of performance and memory by computing high resolution instanton field configurations corresponding to viscous shocks for 1D and 2D compressible flows.

preprint2015arXiv

Large Deviations in Fast-Slow Systems

The incidence of rare events in fast-slow systems is investigated via analysis of the large deviation principle (LDP) that characterizes the likelihood and pathway of large fluctuations of the slow variables away from their mean behavior -- such fluctuations are rare on short timescales but become ubiquitous eventually. This LDP involves an Hamilton-Jacobi equation whose Hamiltonian is related to the leading eigenvalue of the generator of the fast process, and is typically non-quadratic in the momenta -- in other words, the LDP for the slow variables in fast-slow systems is different in general from that of any stochastic differential equation (SDE) one would write for the slow variables alone. It is shown here that the eigenvalue problem for the Hamiltonian can be reduced to a simpler algebraic equation for this Hamiltonian for a specific class of systems in which the fast variables satisfy a linear equation whose coefficients depend nonlinearly on the slow variables, and the fast variables enter quadratically the equation for the slow variables. These results are illustrated via examples, inspired by kinetic theories of turbulent flows and plasma, in which the quasipotential characterizing the long time behavior of the system is calculated and shown again to be different from that of an SDE.

preprint2015arXiv

The instanton method and its numerical implementation in fluid mechanics

A precise characterization of structures occurring in turbulent fluid flows at high Reynolds numbers is one of the last open problems of classical physics. In this review we discuss recent developments related to the application of instanton methods to turbulence. Instantons are saddle point configurations of the underlying path integrals. They are equivalent to minimizers of the related Freidlin-Wentzell action and known to be able to characterize rare events in such systems. While there is an impressive body of work concerning their analytical description, this review focuses on the question on how to compute these minimizers numerically. In a short introduction we present the relevant mathematical and physical background before we discuss the stochastic Burgers equation in detail. We present algorithms to compute instantons numerically by an efficient solution of the corresponding Euler-Lagrange equations. A second focus is the discussion of a recently developed numerical filtering technique that allows to extract instantons from direct numerical simulations. In the following we present modifications of the algorithms to make them efficient when applied to two- or three-dimensional fluid dynamical problems. We illustrate these ideas using the two-dimensional Burgers equation and the three-dimensional Navier-Stokes equations.

preprint2015arXiv

Time-irreversibility of the statistics of a single particle in a compressible turbulence

We investigate time-irreversibility from the point of view of a single particle in Burgers turbulence. Inspired by the recent work for incompressible flows [Xu et al., PNAS 111.21 (2014) 7558], we analyze the evolution of the kinetic energy for fluid markers and use the fluctuations of the instantaneous power as a measure of time irreversibility. For short times, starting from a uniform distribution of markers, we find the scaling $\left<[E(t)-E(0)]^n\right>\propto t$ and $\left<p^n\right> \propto \textrm{Re}^{n-1}$ for the power as a function of the Reynolds number. Both observations can be explained using the "flight-crash" model, suggested by Xu et al. Furthermore, we use a simple model for shocks which reproduces the moments of the energy difference including the pre-factor for $\left<E(t)-E(0)\right>$. To complete the single particle picture for Burgers we compute the moments of the Lagrangian velocity difference and show that they are bi-fractal. This arises in a similar manner to the bi-fractality of Eulerian velocity differences. In the above setting time irreversibility is directly manifest as particles eventually end up in shocks. We additionally investigate time irreversibility in the long-time limit when all particles are located inside shocks and the Lagrangian velocity statistics are stationary. We find the same scalings for the power and energy differences as at short times and argue that this is also a consequence of rare "flight-crash" events related to shock collisions.

preprint2014arXiv

Relevance of instantons in Burgers turbulence

Instanton calculations are performed in the context of stationary Burgers turbulence to estimate the tails of the probability density function (PDF) of velocity gradients. These results are then compared to those obtained from massive direct numerical simulations (DNS) of the randomly forced Burgers equation. The instanton predictions are shown to agree with the DNS in a wide range of regimes, including those that are far from the limiting cases previously considered in the literature. These results settle the controversy of the relevance of the instanton approach for the prediction of the velocity gradient PDF tail exponents. They also demonstrate the usefulness of the instanton formalism in Burgers turbulence, and suggest that this approach may be applicable in other contexts, such as 2D and 3D turbulence in compressible and incompressible flows.

preprint2012arXiv

Finite-Time Euler singularities: A Lagrangian perspective

We address the question whether a singularity in a three-dimensional incompressible inviscid fluid flow can occur in finite time. Analytical considerations and numerical simulations suggest high-symmetry flows being a promising candidate for a finite-time blowup. Utilizing Lagrangian and geometric non-blowup criteria, we present numerical evidence against the formation of a finite-time singularity for the high-symmetry vortex dodecapole initial condition. We use data obtained from high resolution adaptively refined numerical simulations and inject Lagrangian tracer particles to monitor geometric properties of vortex line segments. We then verify the assumptions made by analytical non-blowup criteria introduced by Deng et. al [Commun. PDE 31 (2006)] connecting vortex line geometry (curvature, spreading) to velocity increase to rule out singular behavior.

preprint2012arXiv

Instanton filtering for the stochastic Burgers equation

We address the question whether one can identify instantons in direct numerical simulations of the stochastically driven Burgers equation. For this purpose, we first solve the instanton equations using the Chernykh-Stepanov method [Phys. Rev. E 64, 026306 (2001)]. These results are then compared to direct numerical simulations by introducing a filtering technique to extract prescribed rare events from massive data sets of realizations. Using this approach we can extract the entire time history of the instanton evolution which allows us to identify the different phases predicted by the direct method of Chernykh and Stepanov with remarkable agreement.

preprint2012arXiv

Lagrangian and geometric analysis of finite-time Euler singularities

We present a numerical method of analyzing possibly singular incompressible 3D Euler flows using massively parallel high-resolution adaptively refined numerical simulations up to 8192^3 mesh points. Geometrical properties of Lagrangian vortex line segments are used in combination with analytical non-blowup criteria by Deng et al [Commun. PDE 31 (2006)] to reliably distinguish between singular and near-singular flow evolution. We then apply the presented technique to a class of high-symmetry initial conditions and present numerical evidence against the formation of a finite-time singularity in this case.

preprint2012arXiv

Turbulence properties and global regularity of a modified Navier-Stokes equation

We introduce a modification of the Navier-Stokes equation that has the remarkable property of possessing an infinite number of conserved quantities in the inviscid limit. This new equation is studied numerically and turbulence properties are analyzed concerning energy spectra and scaling of structure functions. The dissipative structures arising in this new equation are curled vortex sheets contrary to vortex tubes arising in Navier-Stokes turbulence. The numerically calculated scaling of structure functions is compared with a phenomenological model based on the She-Lévêque approach. Finally, for this equation we demonstrate global well-posedness for sufficiently smooth initial conditions in the periodic case and in $\mathbb R^3$. The key feature is the availability of an additional estimate which shows that the $L^4$-norm of the velocity field remains finite.