Thomas Burton, Sean Symon, Ati S Sharma, Davide Lasagna
{"title":"Resolvent-based optimization for approximating the statistics of a chaotic Lorenz system.","authors":"Thomas Burton, Sean Symon, Ati S Sharma, Davide Lasagna","doi":"10.1103/PhysRevE.111.025104","DOIUrl":null,"url":null,"abstract":"<p><p>We propose a 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 [D. Lasagna, Phys. Rev. E 102, 052220 (2020)2470-004510.1103/PhysRevE.102.052220] 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 optimization to minimize the violation of the projected governing equations, producing trajectories that approximate, rather than exactly solve, the system dynamics. An attempt at utilizing 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 behavior, as the trajectories obtained from partial optimization provide a sufficient \"sketch\" of the attractor in state space.</p>","PeriodicalId":48698,"journal":{"name":"Physical Review E","volume":"111 2-2","pages":"025104"},"PeriodicalIF":2.2000,"publicationDate":"2025-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Physical Review E","FirstCategoryId":"101","ListUrlMain":"https://doi.org/10.1103/PhysRevE.111.025104","RegionNum":3,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q2","JCRName":"PHYSICS, FLUIDS & PLASMAS","Score":null,"Total":0}
引用次数: 0
Abstract
We propose a 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 [D. Lasagna, Phys. Rev. E 102, 052220 (2020)2470-004510.1103/PhysRevE.102.052220] 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 optimization to minimize the violation of the projected governing equations, producing trajectories that approximate, rather than exactly solve, the system dynamics. An attempt at utilizing 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 behavior, as the trajectories obtained from partial optimization provide a sufficient "sketch" of the attractor in state space.
期刊介绍:
Physical Review E (PRE), broad and interdisciplinary in scope, focuses on collective phenomena of many-body systems, with statistical physics and nonlinear dynamics as the central themes of the journal. Physical Review E publishes recent developments in biological and soft matter physics including granular materials, colloids, complex fluids, liquid crystals, and polymers. The journal covers fluid dynamics and plasma physics and includes sections on computational and interdisciplinary physics, for example, complex networks.