Variational optimization and data assimilation in chaotic time-delayed systems with automatic-differentiated shadowing sensitivity
arXiv:2011.08794
Abstract
In this computational paper, we perform sensitivity analysis of long-time (or ensemble) averages in the chaotic regime using the shadowing algorithm. We introduce automatic differentiation to eliminate the tangent/adjoint equation solvers used in the shadowing algorithm. In a gradient-based optimization, we use the computed shadowing sensitivity to minimize different long-time averaged functionals of a chaotic time-delayed system by optimal parameter selection. In combined state and parameter estimation for data assimilation, we use the computed sensitivity to predict the optimal trajectory given information from a model and data from measurements beyond the predictability time. The algorithms are applied to a thermoacoustic model. Because the computational framework is rather general, the techniques presented in this paper may be used for sensitivity analysis of ensemble averages, parameter optimization and data assimilation of other chaotic problems, where shadowing methods are applicable.
50 pages, under review
References in corpus (13)
- Characterizing dynamics with covariant Lyapunov vectors
- Forward and Adjoint Sensitivity Computation of Chaotic Dynamical Systems
- Sensitivity analysis on chaotic dynamical system by Non-Intrusive Least Square Shadowing (NILSS)
- DiffEqFlux.jl - A Julia Library for Neural Differential Equations
- Convergence of the least squares shadowing method for computing derivative of ergodic averages
- Don't Unroll Adjoint: Differentiating SSA-Form Programs
- Stability analysis of thermo-acoustic nonlinear eigenproblems in annular combustors. Part II. Uncertainty quantification
- Sensitivity of the Rayleigh criterion in thermoacoustics
- Adjoint sensitivity analysis on chaotic dynamical systems by Non-Intrusive Least Squares Adjoint Shadowing (NILSAS)
- Feasibility analysis of ensemble sensitivity computation in turbulent flows
- A computable realization of Ruelle's formula for linear response of statistics in chaotic systems
- Adjoint shadowing directions in hyperbolic systems for sensitivity analysis
- Approximating linear response by nonintrusive shadowing algorithms