Source author record

Tobias M. Schneider

Tobias M. Schneider 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

18works
6topics
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

18 published item(s)

preprint2026arXiv

Identifying recurrent flows in high-dimensional dissipative chaos from low-dimensional embeddings

Unstable periodic orbits (UPOs) are the non-chaotic, dynamical building blocks of spatio-temporal chaos, motivating a first-principles based theory for turbulence ever since the discovery of deterministic chaos. Despite their key role in the ergodic theory approach to fluid turbulence, identifying UPOs is challenging for two reasons: chaotic dynamics and the high-dimensionality of the spatial discretization. We address both issues at once by proposing a loop convergence algorithm for UPOs directly within a low-dimensional embedding of the chaotic attractor. The convergence algorithm circumvents time-integration, hence avoiding instabilities from exponential error amplification, and operates on a latent dynamics obtained by pulling back the physical equations using automatic differentiation through the learned embedding function. The interpretable latent dynamics is accurate in a statistical sense, and, crucially, the embedding preserves the internal structure of the attractor, which we demonstrate through an equivalence between the latent and physical UPOs of both a model PDE and the 2D Navier-Stokes equations. This allows us to exploit the collapse of high-dimensional dissipative systems onto a lower dimensional manifold, and identify UPOs in the low-dimensional embedding.

preprint2020arXiv

Adjoint-based variational method for constructing periodic orbits of high-dimensional chaotic systems

Chaotic dynamics in systems ranging from low-dimensional nonlinear differential equations to high-dimensional spatio-temporal systems including fluid turbulence is supported by non-chaotic, exactly recurring time-periodic solutions of the governing equations. These unstable periodic orbits capture key features of the turbulent dynamics and sufficiently large sets of orbits promise a framework to predict the statistics of the chaotic flow. Computing periodic orbits for high-dimensional spatio-temporally chaotic systems remains challenging as known methods either show poor convergence properties because they are based on time-marching of a chaotic system causing exponential error amplification; or they require constructing Jacobian matrices which is prohibitively expensive. We propose a new matrix-free method that is unaffected by exponential error amplification, is globally convergent and can be applied to high-dimensional systems. The adjoint-based variational method constructs an initial value problem in the space of closed loops such that periodic orbits are attracting fixed points for the loop-dynamics. We introduce the method for general autonomous systems. An implementation for the one-dimensional Kuramoto-Sivashinsky equation demonstrates the robust convergence of periodic orbits underlying spatio-temporal chaos. Convergence does not require accurate initial guesses and is independent of the period of the respective orbit.

preprint2020arXiv

Interaction and breakup of droplet pairs in a microchannel Y-junction

We combine theory, numerical simulation and experiments to investigate the breakup of two identical droplets entering a Y-junction with controlled spatial offset by which the second droplet trails the first. Based on fully resolved 3D simulations, we describe the flow physics leading to breakup. Scaling arguments, numerical simulation and experiments consistently show that for small initial offset, breakup always occurs with the droplet fragment volume depending linearly on the offset. Above a critical offset, which increases with the capillary number, the droplets enter the constriction sequentially without breakup. Our results are relevant for understanding the breakup behavior in a dense emulsion flowing through a linearly converging channel leading to a constriction. Such geometry is commonly used for the serial interrogation of droplet content in droplet microfluidic applications, where droplet breakup can limit the maximum throughput for such process. For capillary numbers up to $Ca\simeq10^{-2}$, the results from the two-droplet system in a Y-junction are consistent with breakup observations in dense emulsions flowing through a linearly converging channel. The deterministic relation between initial offset and resulting breakup in the two-droplet system suggests that the stochasticity that is observed in the breakup of a dense emulsion arises from multi-droplet interactions. The numerical value of the prefactor in the linear relation between initial offset and droplet fragment volume determined from experiments differs slightly from the one extracted from fully resolved numerical simulations. This discrepancy suggests that even at very high bulk surfactant concentrations, the rate-limiting surfactant adsorption kinetics allows for Marangoni stresses to develop and modify the droplet dynamics.

preprint2020arXiv

Invariant states in inclined layer convection. Part 1. Temporal transitions along dynamical connections between invariant states

Thermal convection in an inclined layer between two parallel walls kept at different fixed temperatures is studied for fixed Prandtl number Pr=1.07. Depending on the angle of inclination and the imposed temperature difference, the flow exhibits a large variety of self-organized spatio-temporal convection patterns. Close to onset, these patterns have been explained in terms of linear stability analysis of primary and secondary flow states. At larger temperature difference, far beyond onset, experiments and simulations show complex, dynamically evolving patterns that are not described by stability analysis and remain to be explained. Here we employ a dynamical systems approach. We construct stable and unstable exact invariant states, including equilibria and periodic orbits of the fully nonlinear three-dimensional Oberbeck-Boussinesq equations. These invariant states underlie the observed convection patterns beyond their onset. We identify state-space trajectories that, starting from the unstable laminar flow, follow a sequence of dynamical connections between unstable invariant states until the dynamics approaches a stable attractor. Together, the network of dynamically connected invariant states mediates temporal transitions between coexisting invariant states and thereby supports the observed complex time-dependent dynamics in inclined layer convection.

preprint2020arXiv

Invariant states in inclined layer convection. Part 2. Bifurcations and connections between branches of invariant states

Convection in a layer inclined against gravity is a thermally driven non-equilibrium system, in which both buoyancy and shear forces drive spatio-temporally complex flow. As a function of the strength of thermal driving and the angle of inclination, a multitude of convection patterns is observed in experiments and numerical simulations. Several observed patterns have been linked to exact invariant states of the fully nonlinear 3D Oberbeck-Boussinesq equations. These exact equilibria, traveling waves and periodic orbits reside in state space and, depending on their stability properties, are transiently visited by the dynamics or act as attractors. To explain the dependence of observed convection patterns on control parameters, we study the parameter dependence of the state space structure. Specifically, we identify the bifurcations that modify the existence, stability and connectivity of invariant states. We numerically continue exact invariant states underlying spatially periodic convection patterns at Pr=1.07 under changing control parameters for temperature difference between the walls and inclination angle. The resulting state branches cover various inclinations from horizontal layer convection to vertical layer convection and beyond. The collection of all computed branches represents an extensive bifurcation network connecting 16 different invariant states across control parameter values. Individual bifurcation structures are discussed in detail and related to the observed complex dynamics of individual convection patterns. Together, the bifurcations and associated state branches indicate at what control parameter values which invariant states coexist. This provides a nonlinear framework to explain the multitude of complex flow dynamics arising in inclined layer convection.

preprint2020arXiv

Modified snaking in plane Couette flow with wall-normal suction

A specific family of spanwise-localised invariant solutions of plane Couette flow exhibits homoclinic snaking, a process by which spatially localised invariant solutions of a nonlinear partial differential equation smoothly grow additional structure at their fronts while undergoing a sequence of saddle-node bifurcations. Homoclinic snaking is well understood in the context of simpler pattern forming systems such as the one-dimensional Swift-Hohenberg equation with cubic-quintic nonlinearity, whose solutions remarkably well resemble the snaking solutions of plane Couette flow. We study the structural stability of the characteristic snakes-and-ladders structure associated with homoclinic snaking for flow modifications that break symmetries of plane Couette flow. We demonstrate that wall-normal suction modifies the bifurcation structure of three-dimensional plane Couette solutions in the same way, a symmetry-breaking quadratic term modifies solutions of the one-dimensional Swift-Hohenberg equation. These modifications are related to the breaking of the discrete rotational symmetry. At large amplitudes of the symmetry-breaking wall-normal suction the connected snakes-and-ladders structure is destroyed. Previously unknown solution branches are created and can be parametrically continued to vanishing suction. This yields new localised solutions of plane Couette flow that exist in a wide range of Reynolds number.

preprint2020arXiv

Self-sustained large-scale motions in the asymptotic suction boundary layer

Large-scale motions, also known as superstructures, are dynamically relevant coherent structures in a wall-bounded turbulent flow, that span the entire domain in wall-normal direction and significantly contribute to the global energy and momentum transport. Recent investigations in channel and Couette flow, suggest that these large-scale motions are self-sustained, implying they are not driven by small-scale motions at the wall. Whether large-scale motions are self-sustained has however not yet been answered for open boundary layers, which are relevant for many applications. Here, using the asymptotic suction boundary layer flow at the friction Reynolds number $Re_τ= 1\,168$ as a testbed, we show that large-scale motions are self-sustained also in boundary layers. Together with the previous investigations in confined flows, this observation provides strong evidence of the robust and general nature of coherent self-sustained processes in turbulent wall-bounded flows. The dynamics of the large-scale self-sustained process involving the growth, breakdown and regeneration of quasi-streamwise coherent streaks and vortices within the boundary layer shows temporal phase relations reminiscent of bursting as observed in buffer layer structures. The dynamics however differs from the quasi-periodic large-scale streak-vortex regeneration cycle observed in confined flows. Based on the similarity of the dynamics of large-scale motions in boundary layers and of small-scale buffer layer structures, we conjecture that the bursting behaviour is associated with the dynamical relevance of only one wall, while two confining walls lead to the quasi-periodic cycle.

preprint2019arXiv

Self-similar invariant solution in the near-wall region of a turbulent boundary layer at asymptotically high Reynolds numbers

At sufficiently high Reynolds numbers, shear-flow turbulence close to a wall acquires universal properties. When length and velocity are rescaled by appropriate characteristic scales of the turbulent flow and thereby measured in \emph{inner units}, the statistical properties of the flow become independent of the Reynolds number. We demonstrate the existence of a wall-attached non-chaotic exact invariant solution of the fully nonlinear 3D Navier-Stokes equations for a parallel boundary layer that captures the characteristic self-similar scaling of near-wall turbulent structures. The branch of travelling wave solutions can be followed up to $Re=1,000,000$. Combined theoretical and numerical evidence suggests that the solution is asymptotically self-similar and exactly scales in inner units for Reynolds numbers tending to infinity. Demonstrating the existence of invariant solutions that capture the self-similar scaling properties of turbulence in the near-wall region is a step towards extending the dynamical systems approach to turbulence from the transitional regime to fully developed boundary layers.

preprint2016arXiv

Localized travelling waves in the asymptotic suction boundary layer

We present two spanwise-localized travelling wave solutions in the asymptotic suction boundary layer, obtained by continuation of solutions of plane Couette flow. One of the solutions has the vortical structures located close to the wall, similar to spanwise-localized edge states previously found for this system. The vortical structures of the second solution are located in the free stream far above the laminar boundary layer and are supported by a secondary shear gradient that is created by a large-scale low-speed streak. The dynamically relevant eigenmodes of this solution are concentrated in the free stream, and the departure into turbulence from this solution evolves in the free stream towards the walls. For invariant solutions in free-stream turbulence, this solution thus shows that that the source of energy of the vortical structures can be a dynamical structure of the solution itself, instead of the laminar boundary layer.

preprint2016arXiv

Spatio-temporal Patterns in Inclined Layer Convection

This paper reports on a theoretical analysis of the rich variety of spatio-temporal patterns observed recently in inclined layer convection at medium Prandtl number when varying the inclination angle $γ$ and the Rayleigh number $R$. The present numerical investigation of the inclined layer convection system is based on the standard Oberbeck-Boussinesq equations. The patterns are shown to originate from a complicated competition of buoyancy-driven and shear-flow driven pattern forming mechanisms. The former are expressed as \rm{longitudinal} convection rolls with their axes oriented parallel to the incline, the latter as perpendicular \rm{transverse} rolls. Along with conventional methods to study roll patterns and their stability, we employ direct numerical simulations in large spatial domains, comparable with the experimental ones. As a result, we determine the phase diagram of the characteristic complex 3D convection patterns above onset of convection in the $γ-R$ plane, and find that it compares very well with the experiments. In particular we demonstrate that interactions of specific Fourier modes, characterized by a resonant interaction of their wavevectors in the layer plane, are key to understanding the pattern morphologies.

preprint2015arXiv

Homoclinic snaking in plane Couette flow: bending, skewing, and finite-size effects

Invariant solutions of shear flows have recently been extended from spatially periodic solutions in minimal flow units to spatially localized solutions on extended domains. One set of spanwise-localized solutions of plane Couette flow exhibits homoclinic snaking, a process by which steady-state solutions grow additional structure smoothly at their fronts when continued parametrically. Homoclinic snaking is well understood mathematically in the context of the one-dimensional Swift-Hohenberg equation. Consequently, the snaking solutions of plane Couette flow form a promising connection between the largely phenomenological study of laminar-turbulent patterns in viscous shear flows and the mathematically well-developed field of pattern-formation theory. In this paper we present a numerical study of the snaking solutions, generalizing beyond the fixed streamwise wavelength of previous studies. We find a number of new solution features, including bending, skewing, and finite-size effects. We show that the finite-size effects result from the shift-reflect symmetry of the traveling wave and establish the parameter regions over which snaking occurs. A new winding solution of plane Couette flow is derived from a strongly skewed localized equilibrium.

preprint2015arXiv

Turbulent-laminar patterns in plane Poiseuille flow

Turbulent-laminar banded patterns in plane Poiseuille flow are studied via direct numerical simulations in a tilted and translating computational domain using a parallel version of the pseudospectral code Channelflow. 3D visualizations via the streamwise vorticity of an instantaneous and a time-averaged pattern are presented, as well as 2D visualizations of the average velocity field and the turbulent kinetic energy. Simulations for Reynolds numbers descending from 2300 to 700 show the gradual development from uniform turbulence to a pattern with wavelength 20 half-gaps near Re=1900, to a pattern with wavelength 40 near Re=1300 and finally to laminar flow near Re=800. These transitions are tracked quantitatively via diagnostics using the amplitude and phase of the Fourier transform and its probability distribution. The propagation velocity of the pattern is approximately that of the mean flux and is a decreasing function of Reynolds number. Examination of the time-averaged flow shows that a turbulent band is associated with two counter-rotating cells stacked in the cross-channel direction and that the turbulence is highly concentrated near the walls. Near the walls, the Reynolds stress force accelerates the fluid through a turbulent band while viscosity decelerates it; advection by the laminar profile acts in both directions. In the center, the Reynolds stress force decelerates the fluid through a turbulent band while advection by the laminar profile accelerates it. These characteristics are compared with those of turbulent-laminar banded patterns in plane Couette flow.

preprint2014arXiv

Infinite Geometric Frustration in a Cubic Dipole Cluster

The geometric arrangement of interacting (magnetic) dipoles is a question of fundamental importance in physics, chemistry and engineering. Motivated by recent progress concerning the self-assembly of magnetic structures, the equilibrium orientation of 8 interacting dipoles in a cubic cluster is investigated in detail. Instead of discrete equilibria we find a new type of ground state consisting of infinitely many orientations. This continuum of energetically degenerate states represents a yet unknown form of magnetic frustration. The corresponding dipole rotations in the flat potential valley of this Goldstone mode enable the construction of frictionless magnetic couplings. Using novel computer-assisted algebraic geometry methods, we moreover completely enumerate all equilibrium configurations. The seemingly simple cubic system allows for exactly 9536 unstable discrete equilibria falling into 183 distinct energy families.

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

Studying edge geometry in transiently turbulent shear flows

In linearly stable shear flows at moderate Re, turbulence spontaneously decays despite the existence of a codimension-one manifold, termed the edge of chaos, which separates decaying perturbations from those triggering turbulence. We statistically analyse the decay in plane Couette flow, quantify the breaking of self-sustaining feedback loops and demonstrate the existence of a whole continuum of possible decay paths. Drawing parallels with low-dimensional models and monitoring the location of the edge relative to decaying trajectories we provide evidence, that the edge of chaos separates state space not globally. It is instead wrapped around the turbulence generating structures and not an independent dynamical structure but part of the chaotic saddle. Thereby, decaying trajectories need not cross the edge, but circumnavigate it while unwrapping from the turbulent saddle.

preprint2011arXiv

Algorithm for a microfluidic assembly line

Microfluidic technology has revolutionized the control of flows at small scales giving rise to new possibilities for assembling complex structures on the microscale. We analyze different possible algorithms for assembling arbitrary structures, and demonstrate that a sequential assembly algorithm can manufacture arbitrary 3D structures from identical constituents. We illustrate the algorithm by showing that a modified Hele-Shaw cell with 7 controlled flowrates can be designed to construct the entire English alphabet from particles that irreversibly stick to each other.

preprint2009arXiv

Snakes and ladders: localized solutions of plane Couette flow

We demonstrate the existence of a large number of exact solutions of plane Couette flow, which share the topology of known periodic solutions but are localized in space. Solutions of different size are organized in a snakes-and-ladders structure strikingly similar to that observed for simpler pattern-forming PDE systems. These new solutions are a step towards extending the dynamical systems view of transitional turbulence to spatially extended flows.