Search

Invariant solutions of the Navier-Stokes equations play an important role in the spatiotemporally chaotic dynamics of turbulent shear flows. Despite the significance of these solutions, their identification remains a computational challenge, rendering many solutions inaccessible and thus hindering progress towards a dynamical description of turbulence in terms of invariant solutions. We compute equilibria of three-dimensional wall-bounded shear flows using an adjoint-based matrix-free variational approach. To address the challenge of computing pressure in the presence of solid walls, we develop a formulation that circumvents the explicit construction of pressure and instead employs the influence matrix method. Together with a data-driven convergence acceleration technique based on dynamic mode decomposition, this yields a practically feasible alternative to state-of-the-art Newton methods for converging equilibrium solutions. We compute multiple equilibria of plane Couette flow starting from inaccurate guesses extracted from a turbulent time series. The variational method outperforms Newton(-hookstep) iterations in successfully converging from poor initial guesses, suggesting a larger convergence radius.

An important feature of turbulent boundary layers are persistent large-scale coherent structures in the flow. Here, we use Dynamic Mode Decomposition (DMD), a data-driven technique designed to detect spatio-temporal coherence, to construct optimal low-dimensional representations of such large-scale dynamics in the asymptotic suction boundary layer (ASBL). In the ASBL, fluid is removed by suction through the bottom wall, resulting in a constant boundary layer thickness in streamwise direction. That is, the streamwise advection of coherent structures by the mean flow ceases to be of dynamical importance and can be interpreted as a continuous shift symmetry in streamwise direction. However, this results in technical difficulties, as DMD is known to perform poorly in presence of continuous symmetries. We address this issue using symmetry-reduced DMD (Marensi et al., J. Fluid Mech. 721, A10 (2023)), and find the large-scale dynamics of the ASBL to be low-dimensional indeed and potentially self-sustained, featuring ejection and sweeping events at large scale. Interactions with near-wall structures are captured when including only a few more modes.

In a dynamical systems description of spatiotemporally chaotic PDEs including those describing turbulence, chaos is viewed as a trajectory evolving within a network of non-chaotic, dynamically unstable, time-invariant solutions embedded in the chaotic attractor of the system. While equilibria, periodic orbits and invariant tori can be constructed using existing methods, computations of heteroclinic and homoclinic connections mediating the evolution between the former invariant solutions remain challenging. We propose a robust matrix-free variational method for computing connecting orbits between equilibrium solutions of a dynamical system that can be applied to high-dimensional problems. Instead of a common shooting-based approach, we define a minimization problem in the space of smooth state space curves that connect the two equilibria with a cost function measuring the deviation of a connecting curve from an integral curve of the vector field. Minimization deforms a trial curve until, at a global minimum, a connecting orbit is obtained. The method is robust, has no limitation on the dimension of the unstable manifold at the origin equilibrium, and does not suffer from exponential error amplification associated with time-marching a chaotic system. Owing to adjoint-based minimization techniques, no Jacobian matrices need to be constructed and the memory requirement scales linearly with the size of the problem. The robustness of the method is demonstrated for the one-dimensional Kuramoto-Sivashinsky equation.

It has recently been speculated that statistical properties of chaos may be captured by weighted sums over unstable invariant tori embedded in the chaotic attractor of hyperchaotic dissipative systems; analogous to sums over periodic orbits formalized within periodic orbit theory. Using a novel numerical method for converging unstable invariant 2-tori in a chaotic PDE, we identify many quasiperiodic, unstable, invariant 2-torus solutions of a modified Kuramoto-Sivashinsky equation exhibiting hyperchaotic dynamics with two positive Lyapunov exponents. The set of tori covers significant parts of the chaotic attractor and weighted averages of the properties of the tori -- with weights computed based on their respective stability eigenvalues -- approximate statistics for the chaotic dynamics. These results are a step towards including higher-dimensional invariant sets in a generalized periodic orbit theory for hyperchaotic systems including spatio-temporally chaotic PDEs.

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.