The University of Southampton
University of Southampton Institutional Repository

Resolvent-based optimization for approximating the statistics of a chaotic Lorenz system

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.
2470-0045
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
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).

Record type: Article

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
Download (5MB)

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
ORCID for Davide Lasagna: ORCID iD orcid.org/0000-0002-6501-6041
ORCID for Ati Sharma: ORCID iD orcid.org/0000-0002-7170-1627

Catalogue record

Date deposited: 24 Feb 2025 17:38
Last modified: 15 May 2025 01:46

Export record

Altmetrics

Contributors

Author: Thomas Burton
Author: Davide Lasagna ORCID iD
Author: Sean Symon
Author: Ati Sharma ORCID iD

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

Atom RSS 1.0 RSS 2.0

Contact ePrints Soton: eprints@soton.ac.uk

ePrints Soton supports OAI 2.0 with a base URL of http://eprints.soton.ac.uk/cgi/oai2

This repository has been built using EPrints software, developed at the University of Southampton, but available to everyone to use.

We use cookies to ensure that we give you the best experience on our website. If you continue without changing your settings, we will assume that you are happy to receive cookies on the University of Southampton website.

×