A new class of shell models is proposed, where the shell variables are defined on a recurrent sequence of integer wave-numbers such as the Fibonacci or the Padovan series, or their variations including a sequence made of square roots of Fibonacci numbers rounded to the nearest integer. Considering the simplest model, which involves only local interactions, the interaction coefficients can be generalized in such a way that the inviscid invariants, such as energy and helicity, can be conserved even though there is no exact self-similarity. It is shown that these models basically have identical features with standard shell models, and produce the same power law spectra, similar spectral fluxes and analogous deviation from self-similar scaling of the structure functions implying comparable levels of turbulent intermittency. Such a formulation potentially opens up the possibility of using shell models, or their generalizations along with discretized regular grids, such as those found in direct numerical simulations, either as diagnostic tools, or subgrid models. It also allows to develop models where the wave-number shells can be interpreted as sparsely decimated sets of wave-numbers over an initially regular grid. In addition to conventional shell models with local interactions that result in forward cascade, a particular helical shell model with long range interactions is considered on a similarly recurrent sequence of wave numbers, corresponding to the Fibonacci series, and found to result in the usual inverse cascade.

We revisit a time-dependent, oval-shaped billiard to investigate a phase transition from bounded to unbounded energy growth. In the static case, the phase space exhibits a mixed structure. The chaotic sea in the static scenario leads to average energy growth for a time-dependent boundary. However, inelastic collisions between the particle and the boundary limit this unbounded energy increase. This transition displays properties similar to continuous phase transitions in statistical mechanics, including scale invariance, interrelated critical exponents governed by scaling laws, and an order parameter/susceptibility approaching zero/infinity at the transition. Furthermore, the system exhibits an elementary excitation that promotes particle diffusion and lacks topological defects that provide modifications to the probability distribution function.

This work explores the relationship between state space methods and Koopman operator-based methods for predicting the time-evolution of nonlinear dynamical systems. We demonstrate that extended dynamic mode decomposition with dictionary learning (EDMD-DL), when combined with a state space projection, is equivalent to a neural network representation of the nonlinear discrete-time flow map on the state space. We highlight how this projection step introduces nonlinearity into the evolution equations, enabling significantly improved EDMD-DL predictions. With this projection, EDMD-DL leads to a nonlinear dynamical system on the state space, which can be represented in either discrete or continuous time. This system has a natural structure for neural networks, where the state is first expanded into a high dimensional feature space followed by a linear mapping which represents the discrete-time map or the vector field as a linear combination of these features. Inspired by these observations, we implement several variations of neural ordinary differential equations (ODEs) and EDMD-DL, developed by combining different aspects of their respective model structures and training procedures. We evaluate these methods using numerical experiments on chaotic dynamics in the Lorenz system and a nine-mode model of turbulent shear flow, showing comparable performance across methods in terms of short-time trajectory prediction, reconstruction of long-time statistics, and prediction of rare events. We also show that these methods provide comparable performance to a non-Markovian approach in terms of prediction of extreme events.

This work explores the intersection of time-delay embeddings, periodic orbit theory, and symbolic dynamics. Time-delay embeddings have been effectively applied to chaotic time series data, offering a principled method to reconstruct relevant information of the full attractor from partial time series observations. In this study, we investigate the structure of the unstable periodic orbits of an attractor using time-delay embeddings. First, we embed time-series data from a periodic orbit into a higher-dimensional space through the construction of a Hankel matrix, formed by arranging time-shifted copies of the data. We then examine the influence of the width and height of the Hankel matrix on the geometry of unstable periodic orbits in the delay-embedded space. The right singular vectors of the Hankel matrix provide a basis for embedding the periodic orbits. We observe that increasing the length of the delay (e.g., the height of the Hankel matrix) leads to a clear separation of the periodic orbits into distinct clusters within the embedded space. Our analysis characterizes these separated clusters and provides a mathematical framework to determine the relative position of individual unstable periodic orbits in the embedded space. Additionally, we present a modified formula to derive the symbolic representation of distinct periodic orbits for a specified sequence length, extending the Poly\'a-Redfield enumeration theorem.

We present a novel methodology for modeling the influence of the unresolved scales of turbulence for sub-grid modeling. Our approach employs the differentiable physics paradigm in deep learning, allowing a neural network to interact with the differential equation evolution and performing an a posteriori optimization by incorporating the solver into the training iteration (an approach known as solver-in-the-loop), thus departing from the conventional a priori instantaneous training approach. Our method ensures that the model is exposed to equations-informed input distributions, accounting for prior corrections and often leading to more accurate and stable time evolution. We present results of our methodology applied to a shell model of turbulence, and we discuss further potential applications to Navier-Stokes equations.