Resolvent-based optimization for approximating the statistics of a chaotic Lorenz system
Resolvent-based optimization for approximating the statistics of a chaotic Lorenz system
We propose a novel framework for approximating the statistical properties of turbulent flows by combining variational methods for the search of unstable periodic orbits with resolvent analysis for dimensionality reduction. Traditional approaches relying on identifying all short, fundamental unstable periodic orbits to compute ergodic averages via cycle expansion are computationally prohibitive for high-dimensional fluid systems. Our framework stems from the observation in Lasagna, Phys. Rev. E (2020), that a single unstable periodic orbit with a period sufficiently long to span a large fraction of the attractor captures the statistical properties of chaotic trajectories. Given the difficulty of identifying unstable periodic orbits for high-dimensional fluid systems, approximate trajectories residing in a low-dimensional subspace are instead constructed using resolvent modes, which inherently capture the temporal periodicity of unstable periodic orbits. The amplitude coefficients of these modes are adjusted iteratively with gradient-based optimisation to minimise the violation of the projected governing equations, producing trajectories that approximate, rather than exactly solve, the system dynamics. A first attempt at utilising this framework on a chaotic system is made here on the Lorenz 1963 equations, where resolvent analysis enables an exact dimensionality reduction from three to two dimensions. Key observables averaged over these trajectories produced by the approach as well as probability distributions and spectra rapidly converge to values obtained from long chaotic simulations, even with a limited number of iterations. This indicates that exact solutions may not be necessary to approximate the system's statistical behaviour, as the trajectories obtained from partial optimisation provide a sufficient ``sketch'' of the attractor in state space.
Burton, Thomas
68513172-dfe6-43f5-894b-fa83cf962237
Lasagna, Davide
0340a87f-f323-40fb-be9f-6de101486b24
Symon, Sean
2e1580c3-ba27-46e8-9736-531099f3d850
Sharma, Ati
cdd9deae-6f3a-40d9-864c-76baf85d8718
18 February 2025
Burton, Thomas
68513172-dfe6-43f5-894b-fa83cf962237
Lasagna, Davide
0340a87f-f323-40fb-be9f-6de101486b24
Symon, Sean
2e1580c3-ba27-46e8-9736-531099f3d850
Sharma, Ati
cdd9deae-6f3a-40d9-864c-76baf85d8718
Burton, Thomas, Lasagna, Davide, Symon, Sean and Sharma, Ati
(2025)
Resolvent-based optimization for approximating the statistics of a chaotic Lorenz system.
Physical Review E, 111 (2), [025104].
(doi:10.1103/PhysRevE.111.025104).
Abstract
We propose a novel framework for approximating the statistical properties of turbulent flows by combining variational methods for the search of unstable periodic orbits with resolvent analysis for dimensionality reduction. Traditional approaches relying on identifying all short, fundamental unstable periodic orbits to compute ergodic averages via cycle expansion are computationally prohibitive for high-dimensional fluid systems. Our framework stems from the observation in Lasagna, Phys. Rev. E (2020), that a single unstable periodic orbit with a period sufficiently long to span a large fraction of the attractor captures the statistical properties of chaotic trajectories. Given the difficulty of identifying unstable periodic orbits for high-dimensional fluid systems, approximate trajectories residing in a low-dimensional subspace are instead constructed using resolvent modes, which inherently capture the temporal periodicity of unstable periodic orbits. The amplitude coefficients of these modes are adjusted iteratively with gradient-based optimisation to minimise the violation of the projected governing equations, producing trajectories that approximate, rather than exactly solve, the system dynamics. A first attempt at utilising this framework on a chaotic system is made here on the Lorenz 1963 equations, where resolvent analysis enables an exact dimensionality reduction from three to two dimensions. Key observables averaged over these trajectories produced by the approach as well as probability distributions and spectra rapidly converge to values obtained from long chaotic simulations, even with a limited number of iterations. This indicates that exact solutions may not be necessary to approximate the system's statistical behaviour, as the trajectories obtained from partial optimisation provide a sufficient ``sketch'' of the attractor in state space.
Text
Lorenz_Paper
- Accepted Manuscript
More information
Accepted/In Press date: 16 January 2025
Published date: 18 February 2025
Identifiers
Local EPrints ID: 498609
URI: http://eprints.soton.ac.uk/id/eprint/498609
ISSN: 2470-0045
PURE UUID: 0526b039-1487-4748-b487-80b3eae80b53
Catalogue record
Date deposited: 24 Feb 2025 17:38
Last modified: 15 May 2025 01:46
Export record
Altmetrics
Contributors
Author:
Thomas Burton
Author:
Ati Sharma
Download statistics
Downloads from ePrints over the past year. Other digital versions may also be available to download e.g. from the publisher's website.
View more statistics