2026-04-13 | | Total: 6
This paper presents an improved Matlab routine, FO_LE, for the numerical computation of Lyapunov exponents of fractional-order systems modeled by Caputo's derivative. It is conceived as an enhanced version of the former FO_Lyapunov and FO_NC_Lyapunov codes for commensurate and non-commensurate orders, respectively. The proposed approach replaces the Gram-Schmidt orthogonalization procedure with QR-based reorthonormalization and uses the new quadratic LIL predictor-corrector scheme for the integration of the extended variational system. Compared with the former implementations, the present routine benefits from the higher order of the fractional integrator LIL and applies to both commensurate and non-commensurate models. Like the previous code, FO_LE retains the full memory structure of the underlying Caputo model. The Matlab code for the LIL solver and for the computation of Lyapunov exponents with FO_LE are provided, while a fast implementation of LIL for commensurate and non-commensurate orders, LIL_nc, is available on MathWorks File Exchange. A benchmark problem with exact solution is used to compare the LIL-based solver with ABM-type methods, whereas the Rabinovich-Fabrikant system illustrates the computation of Lyapunov exponents in different dynamical regimes. The results indicate that the proposed implementation is a compact, robust, and efficient tool for the numerical study of stability and chaos in fractional-order systems.
We study the nature of chaos in finite and infinite dimensional systems through a comparison between the Kuramoto Sivashinsky (KS) equation, the Lorenz system, and a Lorenz type reduction of the KS equation proposed by Wilczak. Numerical simulations of the KS equation reveal intrinsic spatio temporal chaos, with disorder evolving simultaneously in space and time. In contrast, the Lorenz system and the Wilczak reduction exhibit low dimensional temporal chaos lacking spatial complexity. Lyapunov exponent analysis highlights the finite-dimensional convergence properties of the reduced systems and underscores the fundamentally different dynamical nature of chaos in the KS equation. In particular, we demonstrate that low-dimensional reductions may reproduce transient chaotic signatures but do not necessarily retain the structural properties of infinite-dimensional dissipative systems.
Irreversible transport in time-periodic flows is commonly attributed to vorticity, nonlinear forcing, or symmetry breaking. We show that finite-memory reconstruction of the velocity gradient generates a purely geometric mechanism for transport even when the instantaneous flow remains locally irrotational at all times. Memory promotes the velocity gradient to a history-dependent connection along particle trajectories whose noncommutativity produces a finite curvature over one forcing cycle. The associated holonomy generates a measurable loop displacement controlled solely by the dimensionless parameter ωτ_m, which quantifies the phase mismatch between forcing and reconstruction. The predicted scaling is consistent with independently reported measurements across distinct oscillatory flow configurations, supporting the interpretation of memory-induced curvature as a minimal geometric origin of irreversible transport in periodically driven continua.
In this paper, the Ensemble Kalman Filter is compared with a 4DVAR Data Assimilation System in chaotic dynamics. The Lorenz model is chosen for its simplicity in structure and its dynamical similarities with primitive equation models, such as modern numerical weather forecasting. It was examined whether the Ensemble Kalman Filter and 4DVAR are effective in tracking the control for 10%, 20%, and 40% of error in the initial conditions. With 10% of noise, the trajectories of both methods are almost perfect. With 20% of noise, the differences between the simulated trajectories and the observations, as well as the true trajectories, are rather small for the Ensemble Kalman Filter but almost perfect for 4DVAR. However, the differences become increasingly significant at the later part of the integration period for the Ensemble Kalman Filter, due to the chaotic behavior of the system. For the case with 40% error in the initial conditions, neither the Ensemble Kalman Filter nor 4DVAR could track the control with only three observations ingested. To evaluate a more realistic assimilation application, an experiment was created in which the Ensemble Kalman Filter ingested a single observation at the 180th time step in the X, Y, and Z Lorenz variables, and only in the X variable. The results show a perfect fit of 4DVAR and the control during a complete integration period, but the Ensemble Kalman Filter shows disagreement after the 80th time step. On the other hand, a considerable disagreement between the Ensemble Kalman Filter trajectories and the control is observed, as well as a total failure of 4DVAR. Better results were obtained for the case in which observations cover all the components of the model vector.
The inverse energy cascade in turbulent Taylor-Couette flow is studied in line with the results of the large eddy simulation. The simulation results show that the inverse energy cascade first occurs within the core region of the flow channel of the Taylor-Couette flow at higher Reynolds number. It is uncovered that this phenomenon is induced by the pulsed zero shear stress resulting from the singularities of the Navier-Stokes equation. In the core area between the two cylinders, the shear stress is nearly zero at higher Reynolds number. The turbulence generated there has high turbulent energy due to discontinuity of the tangential velocity. Since the energy transfer between the fluid layers is inhibited due to the low shear stress, the turbulent energy cannot be transferred along the radial direction, and small-scale vortices with high turbulent energy are produced. These small-scale vortices are located with the large-scale vortices and cannot be dissipated owing to low shear stress. A peak in the energy spectrum at middle frequency (or wave number) is formed due to the concentration of the small-scale vortices. As the number of the singular points of the Navier-Stokes equation increases with the increasing Reynolds number, the region with zero shear stress expands along the radial direction, intensifying nonlinear instability and energy accumulation. This, in turn, leads to more prominent peaks in the energy spectrum, resulting in a more pronounced inverse energy cascade.
Reservoir observers provide a data-driven approach to the inference of unmeasured variables from observed ones for nonlinear dynamical systems. While previous studies have demonstrated wide applicability, their performance may vary considerably with different input variables, even compromising reliability in the worst cases. To enhance the performance of inference, we integrate residual calibration and attention mechanism into the reservoir observer design. The residual calibration module leverages information from the estimation residuals to refine the observer output, and the attention mechanism exploits the temporal dependencies of the data to enrich the representation of reservoir internal dynamics. Experiments on typical chaotic systems demonstrate that our method substantially improves inference accuracy, especially for the worst cases resulting from the traditional reservoir observers. We also invoke the notion of transfer entropy to explain the reason for the input-dependent observation discrepancy and the effectiveness of the proposed method.