pith. sign in

arxiv: 2605.23812 · v1 · pith:2MLC35UTnew · submitted 2026-05-22 · ⚛️ physics.flu-dyn

An Ensemble Variational approach for High-Dimensional Open-Loop Flow Control

Pith reviewed 2026-05-25 02:40 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn
keywords ensemble variational methodopen-loop flow controlchaotic flowscavity flowgradient approximationnon-intrusive optimizationfluid dynamics
0
0 comments X

The pith

Finite ensemble of perturbed controls approximates gradients for open-loop optimization in chaotic cavity flows.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper establishes that an ensemble-variational framework can estimate cost-function gradients by applying finite differences across a small set of perturbed control inputs, bypassing the adjoint equations that become unstable in highly nonlinear and chaotic flows. This non-intrusive method is demonstrated on two-dimensional cavity flows where a steady forcing is optimized to reduce kinetic energy fluctuations, working across quasi-periodic and fully chaotic regimes. A reader would care because many practical flow-control tasks involve high-dimensional parameter spaces and complex dynamics where traditional gradient methods break down, and a parallelizable alternative could expand the range of solvable problems. In the quasi-periodic regime the controls match those from adjoint optimization and drive the flow toward a limit cycle; in the chaotic regime the same framework still produces usable gradient estimates and reduces fluctuations.

Core claim

The ensemble-variational framework approximates cost-function gradients through a finite ensemble of perturbed control vectors combined with a finite-difference approximation performed directly in ensemble space. When applied to optimize steady forcing in two-dimensional cavity flows across Reynolds numbers that produce quasi-periodic to chaotic dynamics, the method yields control strategies that reduce kinetic energy fluctuations in both regimes, including cases where adjoint-based optimization encounters convergence difficulties.

What carries the argument

The ensemble-variational (EnVar) framework, which estimates gradients from finite differences across a small ensemble of perturbed control vectors.

If this is right

  • In quasi-periodic cavity flows the method recovers control strategies that match adjoint-based results and drive the system toward a periodic limit cycle.
  • In chaotic cavity flows the same ensemble-based gradient estimates remain usable and still reduce flow fluctuations where adjoint methods typically diverge.
  • The approach requires only forward simulations and is therefore parallelizable across ensemble members without modification of the underlying flow solver.
  • The framework scales to high-dimensional control spaces because the finite-difference step operates in the low-dimensional ensemble space rather than the full parameter space.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same ensemble construction could be tested on three-dimensional or experimentally measured flows to check whether the required ensemble size grows with problem dimension.
  • Because the method never forms an adjoint operator, it could be combined with black-box or legacy flow solvers that lack adjoint capability.
  • The finite-difference step in ensemble space suggests a natural link to ensemble Kalman filter techniques already used for state estimation, potentially allowing joint control and assimilation.

Load-bearing premise

A modest number of perturbed control vectors can supply gradient estimates accurate enough to drive optimization even when the flow is chaotic and the underlying dynamics are highly nonlinear.

What would settle it

Apply the EnVar procedure to a documented chaotic cavity flow, obtain an optimized forcing, and observe whether the resulting reduction in kinetic energy fluctuations is smaller than the reduction produced by a converged adjoint run on the same case.

Figures

Figures reproduced from arXiv: 2605.23812 by Jean-Camille Chassaing, Matthieu Queguineur, Riccardo Maranelli, Taraneh Sayadi, Vincent Mons.

Figure 1
Figure 1. Figure 1: Sketch of the open-cavity flow configuration. The filled yellow circles denote the sensors. The slip portions of the lower [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Schematic EnVar loop: ensemble generation, forward propagation, optimisation, optimum update and iteration [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Selection of an actuation region within the domain where the perturbation vector (spatial component of the forcing) is [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Spectrum and leading eigenvectors of the covariance matrix [PITH_FULL_IMAGE:figures/full_fig_p012_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Re = 6250. Streamwise velocity component of mean flow 5a; vorticity field 5b. 0 5 10 15 20 25 30 35 40 Frequency 10 5 10 4 10 3 10 2 10 1 10 0 10 1 Amplitude 2.50 2.75 3.00 3.25 3.50 3.75 4.00 Time 1e2 0.5 1.0 1.5 2.0 Energy 1e 3 (a) Fourier spectrum of the velocity fluctuations with respect to the mean flow. The inset reports the corresponding time history. 1.0 0.5 0.0 0.5 1.0 1.5 uy(3/5, 0) 1e 1 2 1 0 1 … view at source ↗
Figure 6
Figure 6. Figure 6: Characterization of the quasi-periodic dynamic. [PITH_FULL_IMAGE:figures/full_fig_p014_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Evolution and decomposition of the cost function for the EnVar Gauss–Newton ( [PITH_FULL_IMAGE:figures/full_fig_p015_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Optimisation of the forcing shape, adjoint(a), EnVar Gauss–Newton(b), EnVar L-BFGS-B(c). [PITH_FULL_IMAGE:figures/full_fig_p016_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Comparative characterization of quasi-periodic flow dynamics under controlled and uncontrolled conditions. The un [PITH_FULL_IMAGE:figures/full_fig_p017_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Reynolds stresses, ⟨u ′u ′ ⟩ component. Uncontrolled flow (a) and controlled flow (b). 1.0 0.5 0.0 0.5 1.0 1.5 2.0 2.5 x 1.0 0.5 0.0 0.5 y a) 1.0 0.5 0.0 0.5 1.0 1.5 2.0 2.5 x b) 4 2 0 2 4 1e 2 [PITH_FULL_IMAGE:figures/full_fig_p017_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Reynolds stresses, ⟨v ′ v ′ ⟩ component. Uncontrolled flow (a) and controlled flow (b). 17 [PITH_FULL_IMAGE:figures/full_fig_p017_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Evolution of the cost function for the EnVar-based optimisation algorithm for di [PITH_FULL_IMAGE:figures/full_fig_p018_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Optimisation of the forcing shape, EnVar using 5 ensemble members (a), EnVar using 10 ensemble members (b), EnVar [PITH_FULL_IMAGE:figures/full_fig_p019_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Re = 14000. Streamwise velocity component of mean flow; vorticity field. high rotational motion and shear, corresponding to regions susceptible to Kelvin–Helmholtz instabilities. The averaged flow organization and the spatial locations where unsteady structures develop highlights the increased complexity of the flow compared to lower Reynolds numbers. With respect to the previous case, the same actuation … view at source ↗
Figure 15
Figure 15. Figure 15: Evolution of the cost function. − EnVar Gauss–Newton; − · − EnVar L-BFGS-B. 1.0 0.5 0.0 0.5 1.0 1.5 2.0 2.5 x 1.0 0.5 0.0 0.5 y a) 1.0 0.5 0.0 0.5 1.0 1.5 2.0 2.5 x b) 4 2 0 2 4 1e 1 [PITH_FULL_IMAGE:figures/full_fig_p021_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: Optimisation of the forcing shape, obtained with EnVar Gauss–Newton (a) and the EnVar L-BFGS-B (b). [PITH_FULL_IMAGE:figures/full_fig_p021_16.png] view at source ↗
Figure 17
Figure 17. Figure 17: Comparative characterization of quasi-periodic flow dynamics under controlled and uncontrolled conditions. The [PITH_FULL_IMAGE:figures/full_fig_p022_17.png] view at source ↗
Figure 18
Figure 18. Figure 18: Reynolds stresses, ⟨u ′u ′ ⟩ component. Uncontrolled flow (a) and controlled flow (b). Figure17c shows trajectories constructed from the wall-normal velocity component measured at stream￾wise positions x = 3 5 , x = 4 5 in the shear layer. Uncontrolled flow exhibits predominantly chaotic dynamics with complex, irregular trajectories. Under steady forcing, however, the trajectories become more compact. Thi… view at source ↗
Figure 19
Figure 19. Figure 19: Reynolds stresses, ⟨v ′ v ′ ⟩ component. Uncontrolled flow (a) and controlled flow (b). simulation and the corresponding backward adjoint solve. The cost of the latter depends on the implemen￾tation, in particular on the strategy used to access the forward trajectory during the backward integration. In the present adjoint solver, which relies on checkpointing, this leads to an approximate iteration cost o… view at source ↗
read the original abstract

Designing effective optimisation strategies for unsteady flows in the presence of complex dynamics is challenging. Gradient-based optimisation algorithms that rely on gradient information obtained from adjoint equations are efficient for high-dimensional control problems such as those considered here. However, they can be prone to numerical sensitivities when the underlying physics is complex, i.e. when it is highly nonlinear, non-differentiable and chaotic. This work proposes an ensemble-variational (EnVar) framework, which provides a non-intrusive alternative to classical, adjoint-based approaches for flow control applications. This framework approximates cost-function gradients through a finite ensemble of perturbed control vectors. A formulation based on a finite-difference approximation in the ensemble space is employed to address high-dimensional parameter spaces. The methodology is evaluated on two-dimensional cavity flows across Reynolds regimes spanning quasi-periodic to chaotic dynamics, where a steady forcing is optimised. In the quasi-periodic regime, the method identifies control strategies consistent with adjoint-based optimization and achieves a significant reduction of kinetic energy fluctuations, driving the flow toward a periodic limit cycle. In the chaotic regime, the framework remains effective in estimating gradients and mitigating flow fluctuations in situations where adjoint-based approaches typically exhibit convergence issues. This work demonstrates that the EnVar method serves as a computationally efficient, parallelizable, and non-intrusive alternative for high-dimensional optimization problems in complex fluid dynamic regimes.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The paper proposes an ensemble-variational (EnVar) framework as a non-intrusive, adjoint-free alternative for gradient-based optimization of steady forcing in high-dimensional unsteady flows. Gradients are approximated via finite-difference operations on a finite ensemble of perturbed control vectors. The approach is evaluated on 2D cavity flows spanning quasi-periodic to chaotic Reynolds regimes, with claims that it recovers adjoint-consistent controls in the former and remains effective for fluctuation reduction in the latter where adjoints typically fail.

Significance. If the ensemble gradient estimates are shown to be sufficiently accurate, the method would supply a parallelizable, non-intrusive route to high-dimensional control in chaotic fluid systems, addressing a recognized limitation of adjoint methods.

major comments (2)
  1. [Abstract] Abstract: the central claim that the framework 'remains effective in estimating gradients' in chaotic regimes rests on empirical success without any reported ensemble-size convergence study, gradient-error metric, or comparison against adjoint gradients (even in the quasi-periodic regime).
  2. [Methodology (finite-difference formulation)] The finite-difference approximation in ensemble space is load-bearing for all reported performance; the manuscript supplies no analysis of how sampling error or the free parameters (ensemble size, perturbation amplitude) affect gradient quality in exponentially sensitive chaotic cost landscapes.
minor comments (1)
  1. [Abstract] Abstract: the specific Reynolds numbers, ensemble sizes, and perturbation amplitudes used in each regime are not stated, making it difficult to assess reproducibility.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments, which highlight opportunities to strengthen the quantitative support for our claims. We address each major comment below and will revise the manuscript accordingly.

read point-by-point responses
  1. Referee: [Abstract] Abstract: the central claim that the framework 'remains effective in estimating gradients' in chaotic regimes rests on empirical success without any reported ensemble-size convergence study, gradient-error metric, or comparison against adjoint gradients (even in the quasi-periodic regime).

    Authors: We agree that explicit quantitative validation would strengthen the central claim. In the revised manuscript we will add an ensemble-size convergence study together with gradient-error metrics obtained by direct comparison against adjoint gradients in the quasi-periodic regime. For the chaotic regime, where adjoint methods are known to diverge, we will clarify that effectiveness is demonstrated through the achieved fluctuation reduction and physical consistency of the resulting controls rather than through direct gradient comparison. revision: yes

  2. Referee: [Methodology (finite-difference formulation)] The finite-difference approximation in ensemble space is load-bearing for all reported performance; the manuscript supplies no analysis of how sampling error or the free parameters (ensemble size, perturbation amplitude) affect gradient quality in exponentially sensitive chaotic cost landscapes.

    Authors: We acknowledge that a dedicated sensitivity analysis is warranted. The revised manuscript will include additional results (new subsection and/or appendix) that quantify the influence of ensemble size and perturbation amplitude on gradient estimates, including sampling-error estimates obtained from repeated ensemble realizations in the chaotic regime. revision: yes

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The paper defines the EnVar framework directly via finite ensemble perturbations and finite-difference approximation in ensemble space as a non-intrusive alternative to adjoint methods. No equations or steps reduce the claimed gradient estimates or optimization performance to fitted parameters, self-citations, or prior results by the same authors. The abstract and description present an independent methodological construction evaluated on cavity flows, with no evidence of self-definitional loops, renamed known results, or load-bearing self-citation chains. This is the common case of a self-contained empirical method proposal.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

The central claim rests on the ability of a modest ensemble to approximate gradients in high-dimensional spaces; this introduces choices for ensemble size and perturbation scale that are not derived from first principles.

free parameters (2)
  • ensemble size
    Number of perturbed control vectors used to span the approximation space; must be chosen to balance accuracy and cost.
  • perturbation amplitude
    Magnitude of the control perturbations applied to generate the finite-difference estimates.
axioms (1)
  • domain assumption Numerical simulation of the cavity flow produces repeatable cost-function evaluations that can be differenced across ensemble members.
    The method requires multiple forward simulations whose outputs are treated as reliable inputs to the gradient estimator.

pith-pipeline@v0.9.0 · 5788 in / 1336 out tokens · 28471 ms · 2026-05-25T02:40:07.066060+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

49 extracted references · 49 canonical work pages

  1. [1]

    M. J. Kochenderfer, T. A. Wheeler, Algorithms for Optimization, MIT Press, 2019. 26

  2. [2]

    J. C. Spall, Introduction to stochastic search and optimization: estimation, simulation, and control, Wiley-Interscience, 2003

  3. [3]

    M. B. Giles, N. A. Pierce, An introduction to the adjoint approach to design, Flow, Turbulence and Combustion 65 (3) (2000) 393–415.doi:10.1023/A:1011430410075

  4. [4]

    Lions, Optimal Control of Systems Governed by Partial Differential Equations, no

    J.-L. Lions, Optimal Control of Systems Governed by Partial Differential Equations, no. 170 in Die Grundlehren der mathematischen Wissenschaften in Einzeldarstellungen, Springer, 1971

  5. [5]

    Pironneau, On optimum design in fluid mechanics, Journal of Fluid Mechanics 64 (1) (1974) 97– 110.doi:10.1017/S0022112074002023

    O. Pironneau, On optimum design in fluid mechanics, Journal of Fluid Mechanics 64 (1) (1974) 97– 110.doi:10.1017/S0022112074002023

  6. [6]

    Jameson, Aerodynamic design via control theory, Journal of Scientific Computing 3 (3) (1988) 233–260.doi:10.1007/BF01061285

    A. Jameson, Aerodynamic design via control theory, Journal of Scientific Computing 3 (3) (1988) 233–260.doi:10.1007/BF01061285

  7. [7]

    S. B. Pope, Turbulent Flows, Cambridge University Press, Cambridge, UK, 2000

  8. [8]

    Vishnampet, D

    R. Vishnampet, D. J. Bodony, J. B. Freund, A practical discrete-adjoint method for high-fidelity compressible turbulence simulations, Journal of Computational Physics 285 (2015) 173–192.doi: 10.1016/j.jcp.2015.01.009. URLhttps://www.sciencedirect.com/science/article/pii/S0021999115000133

  9. [9]

    C. Liu, Q. Xiao, B. Wang, An ensemble-based four-dimensional variational data assimilation scheme. part i: Technical formulation and preliminary test, Monthly Weather Review 136 (9) (2008) 3363– 3373.doi:10.1175/2008MWR2312.1

  10. [10]

    Raissi, P

    V . Mons, J.-C. Chassaing, T. Gomez, P. Sagaut, Reconstruction of unsteady viscous flows using data assimilation schemes, Journal of Computational Physics 316 (2016) 255–280.doi:10.1016/j.jcp. 2016.04.022

  11. [11]

    Evensen, Data Assimilation: The Ensemble Kalman Filter, Springer, Berlin, Heidelberg, 2009

    G. Evensen, Data Assimilation: The Ensemble Kalman Filter, Springer, Berlin, Heidelberg, 2009. doi:10.1007/978-3-642-03711-5

  12. [12]

    N. K. Chada, Y . Chen, D. Sanz-Alonso, Iterative ensemble kalman methods: A unified perspective with some new variants, FoDS 3 (3) (2021) 331–369.doi:10.3934/fods.2021011

  13. [13]

    Nakano, Behavior of the iterative ensemble-based variational method in nonlinear problems, Non- linear Processes in Geophysics 28 (1) (2021) 93–109.doi:10.5194/npg-28-93-2021

    S. Nakano, Behavior of the iterative ensemble-based variational method in nonlinear problems, Non- linear Processes in Geophysics 28 (1) (2021) 93–109.doi:10.5194/npg-28-93-2021

  14. [14]

    R. N. Bannister, A review of operational methods of variational and ensemble-variational data as- similation, Quarterly Journal of the Royal Meteorological Society 143 (703) (2017) 607–633.doi: 10.1002/qj.2982

  15. [15]

    Z. Ding, Q. Li, Ensemble kalman inversion: Mean-field limit and convergence analysis, Stat Comput 31 (1) (2021) 9.doi:10.1007/s11222-020-09976-0

  16. [16]

    Jahanbakhshi, T

    R. Jahanbakhshi, T. A. Zaki, Nonlinearly most dangerous disturbance for high-speed boundary-layer transition, Journal of Fluid Mechanics 876 (2019) 87–121.doi:10.1017/jfm.2019.527

  17. [17]

    Jahanbakhshi, T

    R. Jahanbakhshi, T. A. Zaki, Optimal heat flux for delaying transition to turbulence in a high-speed boundary layer, Journal of Fluid Mechanics 916 (2021) A46.doi:10.1017/jfm.2021.210. 27

  18. [18]

    D. A. Buchta, T. A. Zaki, Observation-infused simulations of high-speed boundary-layer transition, Journal of Fluid Mechanics 916 (2021) A44.doi:10.1017/jfm.2021.172

  19. [19]

    Journal of Fluid Mechanics 960, 30 (2023) https://doi.org/10.1017/jfm

    J. Lorente-Macias, Y . Bengana, Y . Hwang, Shape optimisation for a stochastic two dimensional cylin- der wake using ensemble variation, Journal of Fluid Mechanics 959 (2023).doi:10.1017/jfm. 2023.122

  20. [20]

    Zhang, F

    X.-L. Zhang, F. Zhang, Z. Li, X. Yang, G. He, Large-eddy simulation-based shape optimization for mitigating turbulent wakes of a bluffbody using the regularized ensemble kalman method, Journal of Fluid Mechanics 1001 (2024) A31.doi:10.1017/jfm.2024.1090

  21. [21]

    Y . Liu, S. Wang, X. Zhang, G. He, Optimizing flow control with ensemble kalman method for miti- gating flow-induced vibration, AIAA Journal (2025). URLhttps://arc.aiaa.org/doi/10.2514/1.J064850

  22. [22]

    A. Q. Rodríguez, Numerical methods for modeling and optimization of interfacial flows, Ph.D. thesis, Sorbonne Université (2024). URLhttps://theses.hal.science/tel-04842433

  23. [23]

    V . Mons, Y . Du, T. Zaki, Ensemble-variational assimilation of statistical data in large-eddy simulation, Physical Review Fluids 6 (10) (2021) 104607.doi:10.1103/PhysRevFluids.6.104607

  24. [24]

    V . Mons, Q. Wang, T. A. Zaki, Kriging-enhanced ensemble variational data assimilation for scalar- source identification in turbulent environments, Journal of Computational Physics 398 (2019) 108856. doi:10.1016/j.jcp.2019.07.054

  25. [25]

    Bocquet, P

    M. Bocquet, P. Sakov, An iterative ensemble kalman smoother, Quarterly Journal of the Royal Mete- orological Society 140 (682) (2013) 1521–1535.doi:10.1002/qj.2236

  26. [26]

    D. Sipp, A. Lebedev, Global stability of base and mean flows: A general approach and its applications to cylinder and open cavity flows, Journal of Fluid Mechanics 593 (2007) 333–358.doi:10.1017/ S0022112007008907

  27. [27]

    Schlichting, K

    H. Schlichting, K. Gersten, Boundary-Layer Theory, 9th Edition, Springer, Berlin, 2017

  28. [28]

    F. M. White, Fluid Mechanics, 8th Edition, McGraw-Hill, New York, 2016

  29. [29]

    Tani, Separation and reattachment in laminar and turbulent boundary layers, Progress in Aerospace Sciences 5 (1964) 201–228

    I. Tani, Separation and reattachment in laminar and turbulent boundary layers, Progress in Aerospace Sciences 5 (1964) 201–228

  30. [30]

    H. von Helmholtz, Über discontinuierliche flüssigkeits-bewegungen [on the discontinuous movements of fluids], Monatsberichte der Königlichen Preussischen Akademie der Wissenschaften zu Berlin 23 (1868) 215–228

  31. [31]

    L. K. W. Thomson), Hydrokinetic solutions and observations, Philosophical Magazine 42 (1871) 362– 377

  32. [32]

    P. G. Drazin, Kelvin–helmholtz instability of a slowly varying flow, Journal of Fluid Mechanics 65 (4) (1974) 781–797.doi:10.1017/S0022112074001650. 28

  33. [33]

    Hwang, M

    K.-J. Hwang, M. L. Goldstein, M. M. Kuznetsova, Y . Wang, A. F. Viñas, D. G. Sibeck, The first in situ observation of kelvin–helmholtz waves at high-latitude magnetopause during strongly dawnward interplanetary magnetic field conditions, Journal of Geophysical Research: Space Physics 117 (2012) A08233.doi:10.1029/2011JA017256

  34. [34]

    Citro, F

    V . Citro, F. Giannetti, L. Brandt, P. Luchini, Linear three-dimensional global and asymptotic stability analysis of incompressible open cavity flow, Journal of Fluid Mechanics 768 (2015) 113–140.doi: 10.1017/jfm.2015.72

  35. [35]

    Picella, J.-C

    F. Picella, J.-C. Loiseau, F. Lusseyran, J.-C. Robinet, S. Cherubini, L. Pastur, Successive bifurcations in a fully three-dimensional open cavity flow, Journal of Fluid Mechanics 844 (2018) 855–877.doi: 10.1017/jfm.2018.169

  36. [36]

    Barbagallo, D

    A. Barbagallo, D. Sipp, P. J. Schmid, Closed-loop control of an open cavity flow using reduced-order models, Journal of Fluid Mechanics 641 (2009) 1–50.doi:10.1017/S0022112009991418

  37. [37]

    D. Sipp, O. Marquet, P. Meliga, A. Barbagallo, Dynamics and control of global instabilities in open- flows: A linearized approach, Applied Mechanics Reviews 63 (030801) (2010).doi:10.1115/1. 4001478

  38. [38]

    Barbagallo, D

    A. Barbagallo, D. Sipp, P. J. Schmid, Input–output measures for model reduction and closed-loop control: Application to global modes, Journal of Fluid Mechanics 685 (2011) 23–53.doi:10.1017/ jfm.2011.271

  39. [39]

    Sipp, Open-loop control of cavity oscillations with harmonic forcings, Journal of Fluid Mechanics 708 (2012) 439–468.doi:10.1017/jfm.2012.329

    D. Sipp, Open-loop control of cavity oscillations with harmonic forcings, Journal of Fluid Mechanics 708 (2012) 439–468.doi:10.1017/jfm.2012.329

  40. [40]

    Leclercq, F

    C. Leclercq, F. Demourant, C. Poussot-Vassal, D. Sipp, Linear iterative method for closed-loop control of quasiperiodic flows, Journal of Fluid Mechanics 868 (2019) 26–65.doi:10.1017/jfm.2019.112

  41. [41]

    V . Mons, O. Marquet, B. Leclaire, P. Cornic, F. Champagnat, Dense velocity, pressure and eulerian ac- celeration fields from single-instant scattered velocities through navier–stokes-based data assimilation, Measurement Science and Technology 33 (2022) 124004

  42. [42]

    D. C. Liu, J. Nocedal, On the limited memory bfgs method for large scale optimization, Mathematical Programming 45 (1-3) (1989) 503–528.doi:10.1007/BF01589116

  43. [43]

    Hecht, New development in freefem++, J

    F. Hecht, New development in freefem++, J. Numer. Math. 20 (3-4) (2012) 251–265. URLhttps://freefem.org/

  44. [44]

    Codina, Pressure stability in fractional step finite element methods for incompressible flows, Journal of Computational Physics 170 (1) (2001) 112–140.doi:10.1006/jcph.2001.6725

    R. Codina, Pressure stability in fractional step finite element methods for incompressible flows, Journal of Computational Physics 170 (1) (2001) 112–140.doi:10.1006/jcph.2001.6725

  45. [45]

    Codina, S

    R. Codina, S. Badia, On some pressure segregation methods of fractional-step type for the finite ele- ment approximation of incompressible flow problems, Computer Methods in Applied Mechanics and Engineering 195 (23-24) (2006) 2900–2918.doi:10.1016/j.cma.2004.06.048

  46. [46]

    Katzfuss, J

    M. Katzfuss, J. R. Stroud, C. K. Wikle, Understanding the ensemble kalman filter, The American Statistician 70 (4) (2016) 350–357. 29

  47. [47]

    X. Wang, S. Ma, D. Goldfarb, W. Liu, Stochastic quasi-newton methods for nonconvex stochastic optimization, SIAM Journal on Optimization 27 (2) (2017) 927–956

  48. [48]

    Cartis, N

    C. Cartis, N. I. M. Gould, P. L. Toint, On the complexity of steepest descent, newton’s and regularized newton’s methods for nonconvex unconstrained optimization problems, SIAM Journal on Optimiza- tion 20 (6) (2010) 2833–2852

  49. [49]

    Costanzo, T

    S. Costanzo, T. Sayadi, M. Fosas de Pando, P. J. Schmid, P. Frey, Parallel-in-time adjoint-based opti- mization – application to unsteady incompressible flows, Journal of Computational Physics 471 (2022) 111664.doi:10.1016/j.jcp.2022.111664. 30