2.1. Wind tunnel set-up

The cavity is inserted into the rectangular cross-section duct of a 0.075 m high and 0.30 m wide wind tunnel. The cavity, inserted as a depression to the floor, is $D = 0.05$ m deep, $S=0.30$ m wide and $L = 0.075$ m or $L = 0.0875$ m long following the studied regime. The resulting aspect ratios are $R = L/D = 1.5$ for the narrow-bandwidth regime and $R = L/D = 1.75$ for the mode-switching regime. A schematic of the wind tunnel is depicted in figure 1. The walls are made of anti-reflection treated glass. A Blasius-type boundary layer develops from an elliptical edge located 0.30 m upstream. Laser Doppler velocimetry (LDV) measurements of the velocity upstream of the cavity show that the standard deviation of the incoming flow is less than 1 %.

Figure 1.

Diagram of the cavity with the position of the DBD actuator (in red) and the velocity sensor (in green). The magnified region depicts the velocity profiles of the incoming velocity ($U(y)$ in blue) and the ionic wind produced by the DBD actuator ($U_p(y)$ in red).

An anemometer is located at the exit of the open wind tunnel vein. Measurements show that the free-stream velocity $U_{\infty }$ and the velocity measured at the exit of the tunnel vein are linearly related to the rotation speed of the wind tunnel fan motor. Thus, in this study, $U_{\infty }$ is estimated from the anemometer measurements. For the narrow-bandwidth regime, the incoming velocity is set to $U_{\infty } = 2.13$ m s$^{-1}$, resulting in a Reynolds number equal to $Re_L = 1.04 \times 10^4$. The momentum boundary layer thickness is estimated at $\theta _0/L=1.17 \times 10^{-2}$. Great care has been taken to calibrate and regulate the incoming velocity regarding LDV measurements. However, we have observed a 2 % variation of the incoming velocity for the narrow-bandwidth regime over the 24 hours necessary for the longest learning sessions. The velocity variations are caused by the temperature variation, $T\approx 23.14 \pm 2\,^\circ$C and very low cycle frequencies in the wind tunnel at this low-velocity operating point. The incoming velocity variations reach 5 % for the mode-switching regime. Finally, the flow is in the incompressible range with a Mach number less than $10^{-2}$. A more detailed description of the set-up can be found in Lusseyran, Pastur & Letellier (Reference Lusseyran, Pastur and Letellier2008) and Basley et al. (Reference Basley, Pastur, Delprat and Lusseyran2013).

2.2. Hot-wire sensor

For sensing, we use a constant temperature anemometer (DANTEC hot-wire probe 55P16 and miniCTA54T30 converter) with a single one-dimensional (1-D) hot-wire sensor, $5\ \mathrm {\mu }$m in diameter and 1 mm length. The hot-wire is located at 6 mm above the cavity and 6 mm upstream of the trailing edge, as sketched in green in figures 1 and 2(b). The position of the hot-wire sensor has also been chosen to limit the velocity drops in the mode-switching regime, see § 2.5. The hot-wire output signal $E_{w}(t)$ is converted into streamwise velocity information $u$ according to King's law

(2.1)\begin{equation} E_{w}^2 = A + Bu^n, \end{equation}

where $A = 1.28$, $B = 0.70$ and $n = 0.48$ are determined by calibration of the hot-wire using an LDV anemometer. Before conversion, the signal $E_{w}$ is temperature corrected by the multiplicative factor $(T_w-T_0)/(T_w-T)$, where $T$ is the room temperature, $T_0$ is the calibration temperature and $T_w$ is the wire temperature (Jørgensen Reference Jørgensen2005); $T$ and $T_0$ are both measured with a $P_{t100}$ platinum sensor with $0.02\,^\circ$C accuracy. The velocity measured $u$ is then employed in three ways: first, it serves to compute the performance of the tested controllers (§ 3.1); second, it closes the feedback control loop (§ 3.2); third, it is used to analyse the control mechanisms (§ 3.4). All the following spectra and spectrograms are computed from this velocity measurement.

Figure 2.

The actuation is performed with a DBD actuator placed upstream (a) and the current state of the flow is given by a hot-wire probe downstream (b).

2.3. Plasma actuator

The actuation is carried out with a dielectric barrier discharge (DBD) actuator to locally force the boundary layer at the entrance of the cavity, near the separation edge where the receptivity of the shear layer is maximum (Cattafesta et al. Reference Cattafesta, Garg, Choudhari and Li1997) (see figure 1). The DBD consists of two conductive blades placed on either side of an insulating plate and is subjected to a high alternating voltage. The streamwise shift between the two blades (see figure 2a) creates an electric field parallel to the plate and is responsible for an ionic wind in the streamwise direction. The principle and the adjustment of the parameters for an application as a fluid actuator are thoroughly detailed in Moreau (Reference Moreau2007), Forte et al. (Reference Forte, Jolibois, Pons, Moreau, Touchard and Cazalens2007) and Benard et al. (Reference Benard, Moreau, Griffin and Cattafesta2010). In our experimental set-up, the dielectric is made of 2 mm-thick acrylic glass or poly(methyl methacrylate) (PMMA) and the electrodes are made of 9 mm-wide, 26 cm-long and $200\ \mathrm {\mu }$m-thick copper ribbons. The downstream edge of the lower electrode is placed at $x=4$ mm upstream to the leading edge, see figure 2(a).

To produce an ionic wind, a carrying signal $E(t)$ at high frequency $f_p$ $({\approx }3\ {\rm kHz})$ is sent to the active electrode. The signal $E(t)$ is produced by an Agilent Function Generator and amplified ($\times 3000$) by a Trek high-voltage amplifier. The expression of the carrying signal is

(2.2)\begin{equation} E(t) =A(t) \sin(2{\rm \pi} f_pt), \end{equation}

with $A$ being the amplitude of the carrying signal. The control is then achieved by modulation of the amplitude $A$ through the actuation command $b \in [-1,1]$. In practice, $A$ is an affine function of $b$ such as $A|_{b=-1} = A_{min}$ and $A|_{b=1}=A_{max}$; $A_{min}$ is the ionization voltage, the threshold above which an ionic wind is produced. The generated wind acts then as a localized body force whose intensity increases with the voltage and thus with $b$. The increasing level of the body force results in the reduction of the main peak of the power spectrum until the dynamics is completely modified. A steady actuation forcing study of the open cavity flow is reported in § 4.1; $A_{max}$ is defined as the maximum voltage that ensures that the main resonance of the cavity is present in the power spectrum.

In practice, $A_{min}$ and $A_{max}$ are measured before each experiment as they are sensible to the atmospheric pressure, room temperature, moisture and number of hours of use of the electrode. To make the control robust against these variations, the range of the actuation command $b$ is set independent of $A_{min}$ and $A_{max}$.

Forte et al. (Reference Forte, Jolibois, Pons, Moreau, Touchard and Cazalens2007) and Moreau (Reference Moreau2007) describe the typical velocity profile generated by a DBD actuator with LDV measurements. In particular, Forte et al. (Reference Forte, Jolibois, Pons, Moreau, Touchard and Cazalens2007) show that for a voltage of $20$ kV and a carrying frequency of $1\ {\rm kHz}$ applied between $0.1$ mm-thick, $20$ cm-long aluminium electrodes, the velocity profile displays a maximum at $y=0.5$ mm from the wall. As the velocity profile moves downstream, the value of the maximum velocity decreases, and its height increases up to ${\sim }1$ mm. For our experiment, Pitot measurements indicate that, for a tension equal to $6$ kV, the maximum velocity is around $0.8$ m s$^{-1}$ and is reached at $y=1.25$ mm of the wall. Unfortunately, the tension value is not significant as $A_{min}$ and $A_{max}$ have changed between two experiments.

2.4. Control unit

In our experiment, the signal acquisitions and actuation command are carried out by a dSPACE real-time controller, including a DS1600 4 core processor board and a DS2201 input/output (I/O) board with 12 bits on a ${\pm }10$ V range analogue-to-digital converter. Only two inputs of the I/O board are exploited, one for the hot-wire signal and one for the voltage delivered by the $P_{t100}$ platinum sensor. The hot-wire signal $E_w$ is translated and amplified ($\times 40$) before analogue-to-digital conversion. All signals are sampled at $250$ Hz such that the Nyquist–Shannon sampling theorem is respected up to three times the highest frequency of interest $f^+\approx 40$ Hz, also avoiding aliasing of the second harmonics. One output of the I/O board is employed to send the command signal to the Agilent Function Generator.

The control optimization process includes two loops: a fast evaluation loop and a slow learning loop, see figure 9. The fast evaluation loop is managed by the ControlDesk software and Simulink. For our study, the evaluation loop operates at the sampling frequency ($250$ Hz). For each control law tested, the time series of the actuation command and the hot-wire signal are recorded and post-processed with MATLAB. The slow learning loop includes the post-processing of the control and the control law update; it is automated with Python and MATLAB scripts. Finally, the whole control unit is supervised by a PowerShell script that automates all the steps of the control optimization.

2.5. Unforced dynamics

As described in § 1, the cavity allows a wide range of complex intra-cavity dynamics by tuning the two remaining cavity flow parameters, namely the upstream speed $U_\infty$ and the width $L$. We recall that the width $S$ and the depth $D$ of the cavity are fixed throughout this study and that the flow is incompressible (Mach number ${<}10^{-2}$). In this article, we aim to stabilize two different flow regimes of different dynamical complexity. For both regimes, the power spectrum is mainly organized within five frequency bands: the very low frequencies, not considered here, a low frequency $f_b$ and the three peaks directly reflecting the resonance of the mixing layer $f^-$, $f_a$ and $f^+$, see figure 3. These frequencies are nonlinearly coupled and satisfy the relationships $f^-= f_a - f_b$ and $f^+ = f_a+f_b$. The two regimes studied differs by the power ratios of the frequencies $f_a$ and $f^+$.

Figure 3.

Time series (top) and spectral content (bottom) of the velocity measured downstream for the two unforced regimes. The main frequencies of the flow are depicted in red: $f_b$, $f^-$, $f_a$, $f^+$ and $f_a$ harmonics. The vertical axis of the spectral plots are in $\log _{10}$ scale. The cost function (see § 3.2) detects the maximum peak in the green-shaded window $St \in [0.5,1.75]$. (a) Narrow-bandwidth regime and (b) mode-switching regime.

The first regime is referred to as the narrow-bandwidth regime and corresponds to a flow dynamics mainly centred on a single frequency $f_a$ and its harmonics. This regime is achieved with $L=7.50$ cm and with an incoming velocity of $U_{\infty }=2.13$ m s$^{-1}$. For this case, the ratio of the powers associated with $f^+$ and $f_a$ is close to $10^{-3}$, see figure 3(a). The coupling between $f_a$ and $f_b$ is then insignificant. In contrast, for the second regime, referred to as the mode-switching regime, the power ratio between $f^+$ and $f_a$ is greater than $0.22$, see figure 3(b). In this case, the nonlinear couplings between frequencies are strong and lead to a chaotic intermittency between $f_a$ and $f^+$ (Lusseyran et al. Reference Lusseyran, Pastur and Letellier2008). In the mode-switching regime, two modes compete in the flow leading to a switch of the dominant frequency. Such intermittency has been mentioned for the first time by Kegerise et al. (Reference Kegerise, Spina, Garg and Cattafesta2004) for a compressible cavity flow. In this study, the mode-switching regime is obtained for incompressible conditions (${Ma}<0.01$) for $L=8.75$ cm and a slightly higher incoming velocity $U_{\infty }=2.23$ m s$^{-1}$, corresponding to a Reynolds number $Re_L=1.28\times 10^{-4}$. The momentum boundary layer thickness is estimated at $\theta _0/L=9.72\times 10^{-3}$. All the cavity flow parameters and experimental conditions are grouped in table 2.

Table 2.

Cavity flow parameters and experimental conditions for the two studied regimes. The temperature $T$ and kinematic viscosity $\nu$ values are averaged over the learning session.

The episodic velocity drops, observed in the time series of the mode-switching regime (figure 3b), are due to slow vertical undulations of the mixing layer which bring the low velocities of the lower part of the mixing layer to the level of the measurement point. The position of the hot-wire sensor has been chosen to minimize these low-velocity incursions while limiting the damping of the oscillations to be controlled. The undulations of the mixing layer are stronger for the mode-switching regime, and the incursions could not be avoided. The presence of these incursions suggests a strong interaction between the mixing layer and the slower intra-cavitary flow for the mode-switching regime.

The temporal evolution of the frequency content for the two regimes is depicted in figure 4. In particular, figure 4(a) shows a clear line for the frequency $f_a$ and less intense lines for its harmonics, whereas figure 4(b) displays a switching between frequencies $f_a$ and $f^+$ and their harmonics over time. The time between two switches is estimated between $15$ and $20$ s; Exceptionally, this time may exceed $40$ s.

Figure 4.

Spectrograms of the downstream velocity for the two studied regimes. The colour code corresponds to the power level in $\log _{10}$scale. The ticks of the colour bar are the exponent values. (a) Narrow-bandwidth regime and (b) mode-switching regime.

To understand the difference in dynamics between the two regimes, we locate them in the Strouhal vs $L/\theta _0$ map (figure 5). Similar maps have been plotted for different impinging shear flows revealing jumps between the modes and linear-like relationships between the Strouhal number and the dimensionless cavity length $L/\theta _0$ or $L/\delta _0$, $\delta _0$ being the boundary layer thickness (Sarohia Reference Sarohia1977; Rockwell & Naudascher Reference Rockwell and Naudascher1978; Knisely & Rockwell Reference Knisely and Rockwell1982). Indeed, Basley et al. (Reference Basley, Pastur, Delprat and Lusseyran2013) shows that, in such incompressible flow, most main frequencies measured in the downstream shear layer align with lines of locked-on modes such that the Strouhal number based on $L$ is given by

(2.3)\begin{equation} St^{L}_{n}(L/\theta_0)=\frac{f_n L}{U_\infty} =\frac{n-\gamma_n(L/\theta_0)}{2},\end{equation}

where the parameter $n = 1, 2, 3$ can be seen as the number of cycles within the cavity length, and the corrective term, $\gamma _n$, can be interpreted as a wave adaptation to the effective resonance length. The authors also propose a model for $\gamma _n$, linear with respect to the dimensionless cavity length $L/\theta _0$

(2.4)\begin{equation} \gamma_n(L/\theta_0)=\frac{41n-L/\theta_0}{10(17-n)}.\end{equation}

On the other hand, (2.3) presents a resemblance with Rossiter's formula for compressible flows in Rossiter (Reference Rossiter1964) where the corrective term is associated with the propagation time of the acoustic waves.

Figure 5.

Strouhal number ($St^{L}$) – dimensionless cavity length ($L/\theta _0$) map, $\theta _0$ being the momentum boundary layer thickness and $L$ the cavity length. Black lines represent the locus of (2.3) for locked-on modes $n=1,2,3$. Black dashed lines represent the locus of (2.3) without the corrective term $\gamma _n$: $St^{L}_{n}=n/2$. The left green band corresponds to the operating point for the narrow-bandwidth regime ($L/\theta _0=85.83$). The right green band corresponds to the operating point for the mode-switching regime ($L/\theta _0=102.90$). The red dots symbolize the three main oscillation modes of the mixing layer. For more details see Appendix B).

The Strouhal distribution is well described by Basley et al. (Reference Basley, Pastur, Delprat and Lusseyran2013), however, it is worth noting that there is still no consensual overview of the origin of the incommensurable frequencies in incompressible open cavity flows. As a first interpretation, the peaks in the spectrum are the result of nonlinear interactions inside the mixing layer dynamics. Indeed, the resonance of the cavity occurs for regimes beyond a critical Reynolds number $Re_c$ in contrast to the Kelvin–Helmholtz instability of the mixing layer that is unstable for all Reynolds numbers. Beyond $Re_c$, self-sustained oscillations appear. This scenario originating from a two-dimensional (2-D) perspective is, however, a little simplistic when considering a real 3-D cavity, especially in an incompressible regime at moderate Reynolds number. In fact, the mixing layer develops in the streamwise direction and, the flow not being strictly parallel, Squire's theorem fails: it is the transverse centrifugal instabilities that transit first, possibly several times, below $Re_c$. With increasing Reynolds number, the centrifugal Görtler–Taylor-type instabilities feature first several bifurcations and then a strong nonlinear development before the cavity resonance. This transition scenario is valid for the two studied regimes whose spectral signature has been described previously and in Appendix A.

In addition, Sipp & Lebedev (Reference Sipp and Lebedev2007) and Meliga (Reference Meliga2017) have shown, using a linearization around the average flow in 2-D simulations, that the occurrence of self-sustaining instabilities in shear-driven cavities is due to a supercritical Hopf bifurcation. Following a similar method, Bengana et al. (Reference Bengana, Loiseau, Robinet and Tuckerman2019) and Tuerke et al. (Reference Tuerke, Sciamarella, Pastur, Lusseyran and Artana2015) manage to predict two incommensurable frequencies in simulations and experiments, respectively. Another approach based on the identification of two characteristic delay times can predict the two frequencies of the flow (Tuerke et al. Reference Tuerke, Lusseyran, Sciamarella, Pastur and Artana2020). The authors show that the nonlinear interactions between these two frequencies can be captured by the resolution of a Stuart–Landau-type amplitude equation, whose quadratic damping term consists of two delayed amplitude terms. In this equation, the first delay time characterizes the upstream travelling hydrodynamic instability wave and hence the feedback of the reflected shear-layer instability (Tuerke et al. Reference Tuerke, Sciamarella, Pastur, Lusseyran and Artana2015). The second delay time is motivated by the hydrodynamic feedback of the recirculating vortices, also referred to as ‘vortex carousel’ and corresponds to an intra-cavity overturning time (Tuerke et al. Reference Tuerke, Pastur, Fraigneau, Sciamarella, Lusseyran and Artana2017). In the following, the description of Basley et al. (Reference Basley, Pastur, Delprat and Lusseyran2013) that leads to figure 5 is sufficient to guide the choice of parameters leading to the two regimes we have chosen to control.

Figure 5 plots the values of the measured frequencies in Strouhal number vs the dimensionless cavity length, for the two regimes and their relation to the resonance points. The computation details of the momentum boundary layer thickness $\theta _0$ are detailed below. For the first regime, $f_a$ is close to a Strouhal number equal to $n/2|_{n=2}$, i.e. at the intersection between the black line and the dashed line, while $f^+$ is clearly below the resonance at $n=3$. As for the second regime, the Strouhal number corresponding to $f_a$ is above the resonant mode $n=2$ ($\gamma _2=-0.14$), and the one corresponding to $f^+$ is below the resonant mode $n=3$ ($\gamma _3=+0.14$). In practice, $U_\infty$ has been chosen such that the average presence rate of the two frequencies $f_a$ and $f^+$ is equalized. The fact that $\lvert \gamma _2 \rvert \approx \lvert \gamma _3 \rvert$ only appears after calculation shows clearly that the parameter guiding the relative intensity of the two main modes is indeed $\lvert \gamma _n \rvert$. The values of Strouhal number and $\gamma _n$ for each frequency are grouped in table 3.

Table 3.

Characteristics of the regimes dynamics. Values of the measured frequencies $f_b$, $f_a$ and $f^+$ and their corresponding Strouhal numbers $St_n^L$. The corrective term $\gamma _n$ is then computed from (2.3).

In fact, we observe a slight discrepancy between the natural frequencies measured and the predictions of Basley et al. (Reference Basley, Pastur, Delprat and Lusseyran2013), which we attribute to a change in the free development of the boundary layer and especially a reduction of the boundary layer thickness. This reduction of the boundary layer thickness can be attributed to the planing effect of the $200\ \mathrm {\mu }$m thick upper electrode, glued just before the leading edge. Therefore, in this work, $L/\theta _0$ was not obtained from the Blasius law ($\theta _0=\kappa \sqrt {2 \nu l_x/U_\infty }$ with $\kappa =0.4696$, $l_x=0.3$ m) and $\nu$ the kinematic viscosity, nor by a direct measurement of $\theta _0$, for lack of optical access, but deduced from (2.3) and (2.4), using the observed frequency (figure 3) and the regime parameters (table 2) for the two considered regimes. First, the value of $\gamma _n$ is computed from $f_n$, $n$, $L$ and $U_{\infty }$ and (2.3), then $L/\theta _0$ is deduced from (2.4). The resulting dimensionless cavity lengths are $L/\theta _0=85.83$ for the narrow-bandwidth regime and $L/\theta _0=102.90$ for the mode-switching regime.

Finally, we have investigated the deviation obtained with the Blasius law. From the values of $L/\theta _0$ and assuming the same expression as the Blasius law, we fit the corresponding kappa for our cases: $\kappa =0.4209$ for the narrow-bandwidth regime and $\kappa =0.4205$ for the mode-switching regime. Both values are close to the one of the Blasius law ($\kappa =0.4696$) but slightly lower, which supports the hypothesis of boundary layer thinning by the presence of the DBD electrode.

The low ($\,f_b$) and very low frequencies ($\,f<1$ Hz) constitute a challenge for automatic learning as their rare occurrences require longer time windows for converged statistics and thus slow down the overall learning process. First, we have chosen to alleviate this difficulty by controlling the narrow-bandwidth regime where the very low frequencies ($\,f<1$ Hz) are around two orders of magnitude lower than $f_a$ in terms of power. Then, we fully embrace the effect of the low frequencies with the mode-switching regime where the nonlinear interactions between $f_a$ and $f^+$ give rise to $f_b$ and especially the very low frequencies $f<1$ Hz: $f_b$ is caused by the triadic interaction between $f_a$ and $f^+$ and the low frequencies ($\,f<1$ Hz) are responsible for the frequency switches in the mode-switching regime. Indeed, the power associated with the very low frequencies ($< 1$ Hz) is more than one order of magnitude greater for the mode-switching regime than for the narrow-bandwidth regime, see figures 3(a) and 3(b). The control of the low frequencies is then performed indirectly by controlling the two other frequencies $f_a$ and $f^+$. Moreover, following Basley et al. (Reference Basley, Pastur, Lusseyran, Soria and Delprat2014), the energetic contribution of the very low frequencies is also due to the coupling between the mixing layer instability and the centrifugal instabilities originating in the spanwise direction within the cavity.

To conclude this description of the cavity dynamics, we recall that the goal we set for the control is to reduce the oscillation of the mixing layer by penalizing the peaks of power in the frequency range that includes $f^-$, $f_a$ and $f^+$, as indicated by the green shaded area of the figure 3.

3.1. Cost function and optimization problem

The aim of this study is to stabilize the open cavity flow in two regimes of different complexity, in particular, the mitigation of the self-sustaining oscillations of the mixing layer. For this, a cost function is built based on the velocity data provided by the hot-wire downstream. The oscillations of the mixing layer are reflected in the oscillations of the velocity signal, thus the goal translates into the reduction of the highest peak of the associated power spectrum. Moreover, the power invested and the power saved by the control must be balanced. In that respect, two terms are considered in the cost function

(3.1)\begin{equation} J = J_a+\gamma J_b.\end{equation}

The term $J_a$ accounts for the peak reduction and $J_b$ for the actuation power invested. Here, $J_{a}$ is defined as the value of the power spectral density maximum in a given frequency window. The value is normalized by the value for the unforced case. Hence, the performance of control law $K$ is given by

(3.2) \begin{equation} J_a(K) = \frac{\mathop{\max}\limits_{St \in [0.5,1.75]} {PSD}(u)}{\mathop{\max}\limits_{St \in [0.5,1.75]} {PSD}(u_0)}, \end{equation}

where ${PSD}(u)$ is the power spectral density of the velocity $u$ measured by the hot-wire for the flow forced with the control law $K$ and $u_0$ is the velocity measured for the unforced flow. A steady actuation forcing study (§ 4.1) shows that the actuation affects both $f_a$ and $f^+$ so the detection window for the maximum of the $PSD$ is set such that it comprises both $f_a$ and $f^+$: $St \in [0.5, 1.75]$. Only the frequencies $f_a$ and $f^+$ are take into account as they are the leading modes of the dynamics; the remaining high-power frequencies ($2f_a$, $3f_a$ in figure 3) are harmonics of the fundamental, i.e. slaved to $f_a$. The detection window is set in Strouhal such that it is independent of the studied regime. The normalization of the cost function $J_a$ by the value of the peak for the unforced flow allows us to have a direct measure of the reduction of the peak.

The $PSD$ is computed over $T_{ev}=40$ s. This choice is motivated for three reasons: first, it allows a good convergence of the statistics; second, the time is short enough to evaluate 1000 individuals in a few hours of experiment, limiting potential drifts and staying close to real-life applications with limited testing budget; third, the mode-switching regime may include one or two switches during this period of time, which is enough to have a record of both frequencies $f_a$ and $f^+$ in the spectrum. Hence, the evaluation time balances practicality and good characterization of the flow dynamics. Anticipating on the results, the value chosen for $T_{ev}$ happened to be enough for the control of the two main frequencies in the mode-switching regime. The control of the mode switching is realized indirectly by the control of the two frequencies involved $f_a$ and $f^+$. It is worth noting that a direct control of the intermittency requires a much longer evaluation time due to its very low frequencies.

The actuation penalization term $J_b$ is estimated from the actuation command $b\in [-1,1]$, as the effective power supplied is not directly accessible in the experiment. Here $J_b$ is based on the square of the actuation command averaged over $T_{ev}$ so that it is an analogue to energy. To simplify the interpretation, $J_b$ is normalized by the range of the actuation so that $J_b=0$ when there is no actuation ($A=A_{min}$) and $J_b=1$ when the controller acts steadily at maximum level ($A=A_{max}$). Therefore

(3.3)\begin{equation} J_b = \frac{\langle (A-A_{min})^2 \rangle}{(A_{max}-A_{min})^2} = \frac{\langle (b+1)^2 \rangle}{4}, \end{equation}

with $\langle {\cdot } \rangle$ denoting the mean value over $T_{ev}=40$ s. The choice of the penalization parameter $\gamma$ is based on the open-loop steady forcing study presented in § 4.1. We show, in particular, that a high level of steady actuation is enough to reduce the cost $J_a$ by at least $90\,\%$. The penalization parameter $\gamma$ is chosen such that the cost for the unforced flow ($J_0=1$) is similar to the cost of the high-level steady actuation ($b =1$), thus the optimal solution aimed needs to efficiently reduce $J_a$ with minimal actuation power. As both cost function components $J_a$ and $J_b$ are normalized, we choose then the penalization parameter to be $\gamma =1$. This choice results in setting the cost of the high-level steady actuation ($b=1$) to $J(b=1)=J_a+J_b \approx 1.1\approx J_0$.

Finally, the normalized standard deviation $\tilde {\sigma }$ of the velocity signal is computed for the best control laws, a posteriori, to characterize the controlled flow. Indeed, effective mitigation of the self-sustained oscillations of the mixing layer results in a reduction of the standard deviation defined as

(3.4)\begin{equation} \tilde{\sigma}(K) = \frac{\sigma (u)}{\sigma (u_0)}, \end{equation}

with $\sigma (u)$ being the standard deviation of the velocity $u$ computed over $T_{ev}=40$ s; $T_{ev}$ is also chosen such that the standard deviation is sufficiently converged. The standard deviation is normalized by the standard deviation of the natural unforced flow so as to have a direct measure of the gain.

3.2. Control problem

As stated previously, the control objective is to stabilize the cavity flow by mitigating the oscillations of the mixing layer downstream. To achieve this goal, the flow is forced with a DBD actuator located at the cavity leading edge. The result of the action is an unsteady body force whose intensity is commanded by the signal $b$ sent to the terminals of the DBD actuator; $b$, also referred as the actuation command, is determined by the control law $K$. The control may operate in an open-loop or closed-loop manner. In this study, the considered open-loop actuations are only steady forcing, and closed-loop control includes the unique velocity sensor and time-delayed records. Thus, the control law reads

(3.5)\begin{equation} b =K(\boldsymbol{a}), \end{equation}

with $\boldsymbol {a}$ being the feature vector comprising flow state information. Then, the control problem to solve can be reformulated as an optimization problem where the goal is to derive the optimal control law $K^*$ that minimizes the cost function $J$

(3.6)\begin{equation} K^* = \mathop{\operatorname{arg\,min}}\limits_{K\in \mathcal{K}} J(K),\end{equation}

with $\mathcal {K}:A \mapsto B$ being the space of all possible control laws, $A$ is the control input domain and $B$ is the range of the actuation command. Deriving the optimal control law $K^*$ without any a priori information on the cost function $J$ is a challenging non-convex optimization problem presenting presumably several minima.

3.3. Gradient-enriched machine learning control

In this section, we present the gMLC algorithm (Cornejo Maceda et al. Reference Cornejo Maceda, Li, Lusseyran, Morzyński and Noack2021) employed to solve the optimization problem (3.6). Gradient-enriched MLC is an iterative function optimizer to derive control laws directly from the plant. The method is based on machine learning control (Duriez et al. Reference Duriez, Brunton and Noack2017, MLC) and is augmented with downhill simplex steps to accelerate the learning; MLC has already been employed to control dozens of experiments outperforming previous control laws with unexpected frequency cross-talk (Noack Reference Noack2019). The choice of downhill simplex algorithm relies on its fast convergence and its easy implementation as it does not require an analytical expression of the cost function but only its evaluation. In the past, downhill simplex has been successful in deriving an adaptive closed-loop control for lift-to-drag ratio optimization over a NACA 0025 airfoil (Tian, Cattafesta & Mittal Reference Tian, Cattafesta and Mittal2006; Cattafesta, Tian & Mittal Reference Cattafesta, Tian and Mittal2009) and reducing the net drag power of the fluidic pinball and a slanted Ahmed body (Li et al. Reference Li, Cui, Jia, Li, Yang, Morzyński and Noack2022). In Cornejo Maceda et al. (Reference Cornejo Maceda, Li, Lusseyran, Morzyński and Noack2021), the gMLC is introduced and employed to stabilize the fluidic pinball. The authors show, in particular, that gMLC outperforms MLC, managing to derive better performing control laws with greater learning speed. It is now applied for the first time in an experiment. Anticipating the results (see Appendix A), the superiority of gMLC over MLC is also verified for the control of the cavity in experimental conditions. The benefits of gMLC compared with MLC comes from the combination of stochastic optimization for exploration of the search space and deterministic optimization for a fast convergence towards the minimum. The method consists of the generation of candidate solutions to (3.6), evaluation of them and systematic recombining of the best ones to improve their performances.

The starting point of gMLC is MLC based on linear genetic programming (Brameier & Banzhaf Reference Brameier and Banzhaf2006, LGP). Following the genetic programming terminology, the candidate solutions are also referred to as individuals. Like the MLC method, gMLC makes no assumptions on the structure of the relationship between the inputs and the outputs of the controller. The optimal solution needs, however, to be computable, meaning it can be expressed by a finite number of mathematical operations with finite memory. Indeed, the candidate solutions are internally represented by matrices inherited from linear genetic programming. Each matrix resembles a computer program that unequivocally codes a control law. Each line of the matrix is an instruction pointing to basic operations ($+$, $-$, $\times$, $\div$, $\cos$, $\sin$, $\tanh$, etc.) and registers containing constant random numbers and variables ($a_1$, $a_2$, $a_3$, etc.). The $N_{\rm inst}$ lines of the matrix are then read linearly yielding the control commands as outputs of the first registers. We refer to Li et al. (Reference Li, Borée, Noack, Cordier and Harambat2019) for more information on the internal representation of the control laws.

The gMLC algorithm starts with a broad exploration of the control law space with a Monte Carlo sampling (MCS) phase. The MCS generates $N_{MCS}$ random matrices that represent the first set of individuals. The individuals are evaluated and added to the database of all individuals. Then, the algorithm alternates between exploration phases carried out by genetic programming and exploitation phases performed by downhill simplex iterations until a stopping criterion is reached. The role of exploration is to locate new and better minima in the space of control laws with stochastic recombination of the best-performing individuals. The stochastic recombination is achieved with the genetic operations crossover and mutation. This exploration is much like the evolution phase in the LGP method, however, in the case of gMLC, the concept of population that evolves through generations is generalized by considering all the individuals evaluated so far and stored in the database. Thus, during the exploration phase new individuals are generated by recombining the best among all the previously evaluated individuals. This ensures that no crucial information is irretrievably lost. The best individuals to recombine are selected following their cost function. The selection is carried out by the tournament method with a tournament size equal to 7 for 100 individuals, following the Duriez et al. (Reference Duriez, Brunton and Noack2017) recommendation. The tournament size is scaled with the number of individuals in the database in order to keep a 7/100 ratio; $N_{p}$ new individuals are built at each exploration phase by recombining the best individuals of the database.

Each exploration phase is followed by an exploitation phase. This step exploits the local gradient information to slide down towards the neighbouring minimum. This is carried out with a variant of downhill simplex for infinite-dimensional spaces introduced by Rowan (Reference Rowan1990) and referred to as downhill subplex. In the following, we do not differentiate between the downhill simplex and subplex as the algorithm steps are similar and only applied to different spaces. The principle of downhill simplex is to linearly combine the $N_{sub}$ best-performing control laws following the gradient of the cost space to derive more performing individuals. In contrast to the exploration phase, the new individuals are built in a deterministic way. First, the $N_{sub}$ best individuals are selected to describe a simplex that lives in the subspace generated by the $N_{sub}$ best individuals. The simplex then crawls in the subspace according to geometric operations (reflection, expansion, contraction and shrink) following the local gradients. Each geometric operation yields one or several new individuals that are linear combinations of the original $N_{sub}$ individuals. After each downhill simplex iteration, the simplex is updated by replacing the least-performing individuals. The downhill simplex steps are iterated until at least $N_{p}$ individuals are generated. The newly generated individuals are then added to the database of all individuals. We emphasize that all the new individuals belong to the subspace defined by the original $N_{sub}$ individuals. The algorithm returns the best-performing control law once the stopping criterion is reached. Otherwise, a new iteration of exploration and exploitation is carried out. The stopping criterion may be a performance threshold or a total number of evaluations when the testing budget is limited.

We note the critical intermediate phase of reconstruction between each exploitation and exploration occurrence. Indeed, the new individuals generated by the downhill simplex are linear combinations of individuals without a matrix representation, which is essential for genetic recombination during the exploration phase. Thus, a matrix reconstruction is performed for each linearly combined individual by solving a secondary optimization problem. The goal is to derive a matrix that translates into a control law that has the same response as the linearly combined one. Such a problem is similar to a surface fitting problem, which we solve with linear genetic programming. The reconstruction phase builds a matrix representation for the linearly combined individuals so that genetic operation can be performed. For more information on the gMLC algorithm, we refer the readers to Cornejo Maceda et al. (Reference Cornejo Maceda, Li, Lusseyran, Morzyński and Noack2021). Figure 6 schematically illustrates the different phases of the gMLC algorithm and the learning principle in the control law space. The MATLAB implementation of gMLC employed for this study is freely available from the website https://github.com/gycm134/gMLC.

Figure 6.

Schematic of the gMLC algorithm. On the left are depicted the three phases: MCS, exploration and exploitation. Each one of the phases requires several experiment evaluations. On the right is depicted a representation of the control law space and the gMLC learning process. The individuals generated during the MCS are depicted in black. The individuals resulting from genetic operations are depicted in blue and genetic operations by dashed blue arrows. The yellow region represents the subspace built from the best individuals. The individuals resulting from downhill simplex are depicted in yellow, and the red arrows symbolize the simplex sub-steps (only one sub-step is depicted). The reconstruction phase is not displayed for clarity.

3.4. Control law investigation

In this section, we propose two methodologies to analyse the actuation mechanisms of optimized control laws. Firstly, an analytical approximation of the control is performed with an affine mapping between the inputs (components of $\boldsymbol a$) and the actuation command ($b$). Such mapping aims to reveal the most relevant component of the feature vector. The affine approximation $\tilde {K}$ of the control law $K$ reads

(3.7)\begin{equation} \tilde{K}(\boldsymbol a) = k_0 + \sum_{i}k_i a_i, \end{equation}

where $k_i$ are gains determined by linear regression between the components of $\boldsymbol a$ and the actuation command $b$. The quality of the fitting is measured by the coefficient of determination $R^2$, measuring the relative reduction of the residual variance. The closer $R^2$ is to 1, the better $\tilde {K}$ fits the original control law $K$.

Secondly, we propose a visualization of the control laws based on the clustering of the feature vector $\boldsymbol a$ to reconstruct the phase portrait. Cluster-based methods have been successful in reproducing key characteristics of fluid flow dynamics such as temporal evolution and fluctuation levels (Fernex, Noack & Semaan Reference Fernex, Noack and Semaan2021; Li et al. Reference Li, Fernex, Semaan, Tan, Morzyński and Noack2021). For this analysis, all the states of the feature vector are grouped in 10 clusters to reconstruct the dynamics. The cluster centroids, $c_k$, are defined as the average state of all the states in a given cluster. Clustering is performed with the $k$-means algorithm (Lloyd Reference Lloyd1982), and the metric employed is the one induced by the $L^2$ norm. The dynamics of the feature vector is then encapsulated in a probability transition matrix where its elements $p_{ij}$ are the transition probabilities from cluster $i$ to cluster $j$. The probability $p_{ij}$ is defined as $p_{ij}=n_{ij}/n_i$ with $n_{ij}$ being the number of states transitioning from cluster $i$ to cluster $j$ and $n_i$ the total number of states in cluster $i$.

Then all feature vector states and centroids are projected on a 2-D space with classical multidimensional scaling (Kaiser, Li & Noack Reference Kaiser, Li and Noack2017; Li et al. Reference Li, Cui, Jia, Li, Yang, Morzyński and Noack2022, MDS). The MDS is a dimensional reduction method that consists of extracting the two main features of the flow ($\gamma _1$ and $\gamma _2$) by applying a proper orthogonal decomposition on the distance matrix of the feature vector $\boldsymbol a$. The vectors $\gamma _1$ and $\gamma _2$ spawn a 2-D space where all the data are projected. It is the optimal projection, in the $L^2$ norm sense, that preserves the distances between the states. Such representation is referred to as a proximity map.

Adding the probability transitions to the proximity map allows us to build a network model reproducing the phase portrait. The centroids constitute representative states of the flow where the system transitions ergodically, meaning that from any centroid, one can reach any other centroid. Finally, the mean forcing level is computed for each cluster and associated with their corresponding centroid. Such representation allows us to partition the states of strong and low forcing level and to reveal actuation mechanisms. Figure 7 summarizes the two approximation methodologies employed.

Figure 7.

Control law investigation methodology for the cavity control. A feature vector $\boldsymbol a$ is built from a direct measurement of the flow $u$. A feature vector of dimension three is depicted for simplicity. The elements of the feature vector are grouped in clusters. For clarity, only three clusters and their centroids (1,2,3) are displayed. The proximity map is defined by $\gamma _1$ and $\gamma _2$ the two flow features extracted with classical multidimensional scaling. For the transition matrix, darker squares symbolize higher transition probabilities. In the control network model, the actuation magnitudes are represented by rectangles; yellow denote the actuation range, red (blue) for a positive (negative) actuation with respect to the mean value.

The latter data-driven methodology for control visualization is expected to aid the human interpretability of machine-learned controls. We study a single-input single-output system; however, we believe that the methodology will be beneficial for the analysis of more complex control systems, including a high number of actuators and sensors.

4.1. Open-loop steady forcing

In this section, the response of the flow to steady actuation is described. For this study, the amplitude of the carrying signal is set to constant values. The flow is excited with 23 levels of actuation equally distant from the ionization level ($A=A_{min}=6.9$ kV) to $A=A_{max}=12$ kV. Figure 8 presents the velocity power spectra for the two regimes. For the narrow-bandwidth regime, figure 8(a) shows that the second peak $f^{+}$ rises and the first peak $f_a$ decreases as the actuation level increases. When the actuation is too strong, the noise level increases and the two peaks are at the same level. The maximum peak reduction is achieved for $A=A_{max}$ with a cost reduction of $97\,\%$. The associated standard deviation slightly decreases to $\tilde {\sigma }=96\,\%$.

Figure 8.

Response of the cavity flow to increasing levels of steady actuations. Only 12 actuation levels among the 23 tested are depicted for clarity. The vertical axes are in $\log _{10}$ scale, and the ticks correspond to the exponent values. The spectra have been successively shifted upwards of 4 orders of magnitude to help the comparison. The difference of actuation intensity between two spectra is $\Delta V=0.6$ kV. The spectrum of the unforced flow is plotted in red and in transparency for comparison with the other levels. (a) Narrow-bandwidth regime and (b) mode-switching regime.

For the mode-switching regime, figure 8(b) shows that a strong actuation level is needed to reduce the amplitude of the peaks associated with $f_a$ and $f^+$, although the broadband noise level also increases. The maximum cost reduction is achieved for the nearly maximum actuation level $A=11.4$ kV ($88\,\%$ of $A_{max}$) and the maximum peak power decreases by $90\,\%$. Also, the standard deviation slightly decreases to $\tilde {\sigma }=95\,\%$. We note that for the third spectrum starting from the bottom ($V=1.2$ kV), the accidental absence of mode switching during the measurement led to a lower $f^+$ peak.

This open-loop analysis shows that a strong steady actuation is able to reduce the fluctuations of the shear layer. However, in both cases the background noise level increases. Therefore, in the following, in order to exclude power demanding controllers, we consider the cost function described in § 3.1 that includes two terms: a measure of the maximum amplitude of the spectrum and an actuation penalization term.

4.2. Implementation of gMLC

The parameters chosen for the gMLC algorithm are similar to the ones chosen in Cornejo Maceda et al. (Reference Cornejo Maceda, Li, Lusseyran, Morzyński and Noack2021). The MCS phase generates $N_{MCS}=100$ individuals. The exploration and exploitation phases both produce $N_{p}=50$ new individuals at each iteration. The exploration and exploitation phases alternate until 1000 individuals are evaluated. We recall that each individual is evaluated over $T_{ev}=40$ s. A relaxation time of $2$ s is intercalated between two consecutive control law evaluations. The experiment time needs also to include the time needed to solve the reconstruction problem, but this constraint can be lifted with additional computation power. The limit of 1000 individuals is then chosen such that all the individuals are evaluated in one day. Cornejo Maceda et al. (Reference Cornejo Maceda, Li, Lusseyran, Morzyński and Noack2021) show also that 1000 evaluations is enough to converge for a multiple-input multiple-output problem. The subplex space is generated by $N_{sub}=10$ control laws to balance speed and performance, as in Cornejo Maceda et al. (Reference Cornejo Maceda, Li, Lusseyran, Morzyński and Noack2021). For the evolution during the exploration phase, the crossover and mutation probabilities are both set to $P_c=P_m=0.5$. The control laws are built from nine mathematical operations ($+$, $-$, $\times$, $\div$, $\sin$, $\cos$, $\tanh$, $\exp$ and $\log$), ten flow features $\{a_i\}_{i \ldots 10}$ and $N_{cr}=10$ random constants. As suggested by Duriez et al. (Reference Duriez, Brunton and Noack2017), the $\div$ and $\log$ operations are protected, allowing them to be defined for all the real numbers; $N_{vr}=14$ registers are employed to derive the control laws. Finally, the maximum number of instructions to be coded in the matrix representation is $N_{inst, max}=50$. The flow features employed for feedback control laws are the velocity signal and nine time-delayed velocity signals. Time-delayed sensor signals are introduced as control inputs to enrich the search space and allow, in principle, autoregressive moving average with extra input (ARMAX)-type controllers (Hervé et al. Reference Hervé, Sipp, Schmid and Samuelides2012), linear and nonlinear combinations of them. The resulting feature vector $\boldsymbol a$ reads

(4.1)\begin{equation} \boldsymbol{a}(t) = [u(t),u(t-\tau),\ldots,u(t-9\tau)]^{\intercal}, \end{equation}

where the time delay is $\tau = 0.008$ s. The choice of the delays allows reconstruction of the phase of the main frequencies of the flow between $25$ and $40$ Hz (Cornejo Maceda et al. Reference Cornejo Maceda, Li, Lusseyran, Morzyński and Noack2021). Only half of the delays are necessary but nine have been taken into account to enrich the phase space and get closer to full-state control. The presence of time-delayed information in the control can play, for example, the role of an embedding process in the new dynamical system consisting of the flow and the closed-loop control. Indeed, we have opted for a single measurement point separated from the actuator by a convective time that intrinsically varies over time. Implicitly, the system under closed-loop control is hence reduced to a purely temporal dynamical system and the spatial information can be interpreted as an embedding of this dynamic into a larger phase space.

Table 4 summarizes the parameters for the gMLC optimization process. Figure 9 displays the experiment set-up to control the open cavity with machine learning control and specifically with gMLC.

Table 4.

Gradient-enriched MLC parameters to control the open cavity.

Figure 9.

Diagram of the machine learning control process. The evaluation of one individual constitutes the fast evaluation loop, operating at $250$ Hz. The slow learning loop updates the control laws following the gMLC algorithm; it operates at $\sim O(10^{-2}\ {\rm Hz})$.

4.3. Closed-loop control of the narrow-bandwidth regime

In this section, we describe the best control law derived by gMLC that mitigates the self-sustaining oscillations of the mixing layer for the narrow-bandwidth regime. In the following, the notations for the control law, cost and standard deviation associated with the learning on the narrow-bandwidth regime are marked by the superscript $\textrm {I}$.

Figure 10 depicts the cost for the 1000 evaluated individuals during the optimization process. The individuals from the MCS and exploration phases are sorted following their cost, as there is no direct causal relationship between two successive individuals, while the individuals generated during the exploitation phase are depicted in the order of their evaluation as each individual depends on the previous one. We recall that the cost function is defined such that the cost of the unforced flow is $J_0=1$. We note that a random sampling of 100 control laws already reduces the cost function to $J=0.2410$. Then, the first exploration phase (individuals $i=101,\ldots,150$) slightly reduces the cost function to $J=0.2276$. In the following exploitation phase (individuals $i=151,\ldots,206$), the downhill subplex individuals gradually reduce the cost function to $J=0.1882$. We note that, at first, the individuals are scattered along the vertical axis and then go down, close to the lowest cost so far. This behaviour shows that the individuals progress towards a minimum in the control law space. However, this descent is interrupted by the next exploration phase (individuals $i=207,\ldots,257$), where a better-performing individual is found; its cost is $J=0.1329$. This new individual replaces the worst-performing individual in the simplex, allowing the algorithm to explore beyond the initial subspace. It is worth noting that the dimension of the subspace remains the same as new individuals replace the worst-performing ones. From there on and until individual $i=512$, only the exploitation phases built better individuals. The cost of the best individual after 512 individuals is $J=0.0402$. Interestingly, we note that downhill simplex can build worse-performing individuals. Indeed, we notice that all exploitation phases starting from the $4$th one include individuals whose costs are scattered, up to $J=1$. The next exploration phase (individuals $i=513,\ldots,563$) finds a better control, whose cost is $J=0.0311$. As the simplex includes now poor performing individuals, four better-performing individuals generated by the exploration phase are introduced in the simplex. From there on, progress is made only with exploitation steps: the cost of the best individual reaches a plateau after $707$ evaluations and slightly improves after 913 evaluations. After 1000 evaluations, the cost of the final control law $K^{\textrm {I}}$ is $J^{\textrm {I}} = 0.0192$. The corresponding peak reduction is ${J_a}^{\textrm {I}}=0.0129$, i.e. 99 % of the fluctuation level In other terms, the amplitude of the oscillations is reduced by a factor of 9 or by 19 dB. The structure and components of $K^{\textrm {I}}$ is thoroughly described in Appendix C.

Figure 10.

Distribution of the costs for the 1000 evaluated individuals during the gMLC optimization process for the narrow-bandwidth regime. Each dot represents the cost $J$ of one individual. The colour of the dots symbolizes how the individuals are generated: black dots for randomly generated individuals (MCS phase), blue dots for individuals generated by a genetic operator (exploration phase) and yellow dots for the individuals arising from the subplex method (exploitation phase). The individuals from the MCS and exploration phases are sorted following their costs. The red line follows the evolution of the best cost. The vertical axis is in log scale.

The learned control law $K^{\textrm {I}}$ is re-evaluated 20 times afterwards to test its efficiency outside the learning loop. We note that the performances slightly dropped as the averaged cost reduction of $J^{\textrm {I}}$ went from $-98\,\%$ to $-94\,\%$ and the standard deviation increased from ${\tilde {\sigma }}^{\textrm {I}}=61\,\%$ to $\tilde {\sigma }=65\,\%$. Such discrepancy is expected as the experimental conditions are always evolving. Indeed, the temperature in the room, the evolution over time of the DBD actuator and the incoming velocity variation are all possible sources of fluctuations on the measured velocity. However, the discrepancies are small, and the overall performance of the control law is retained.

Figure 11 depicts the flow response controlled by the best control law $K^{\textrm {I}}$ derived by gMLC for the narrow-bandwidth regime. We observe that the control law $K^{\textrm {I}}$ effectively answers to the control objective as the oscillation amplitude is reduced compared with the unforced flow, see figure 11(a). This goes along with the decrease of the standard deviation to ${\tilde {\sigma }}^{\textrm {I}} = 61\,\%$. Such a feature was consistently observed in several learning realizations. The power spectra in figure 11(b) show that $K^{\textrm {I}}$ effectively reduces the highest peak of the spectrum at frequency $f_a$ by almost two orders of magnitude. The effect of the control is also observed beyond the observation window between $St \in [0.5,1.75]$ as the harmonics of $f_a$ are also mitigated. Note that the peak $f^+$ associated with the mode $n=3$ increases with the control. This behaviour is not surprising as the frequency is associated with a quasi-stable mode of the flow, which then grows when energy is supplied to the system. Moreover, the frequency $f^+$ is split into two peaks. This peak-splitting phenomenon, referred as spillover, is well known when building closed-loop transfer functions with unstable zeros and poles (Rowley & Williams Reference Rowley and Williams2006; Rowley et al. Reference Rowley, Williams, Colonius, Murray and Macmynowski2006). Despite being beyond a linear framework, one can suspect that a similar mechanism is behind the observed peak splitting. The learned law $K^{\textrm {I}}$ fails to control $f^+$ as its power level is around one order of magnitude lower than $f_a$ under control.

Figure 11.

Characteristics of the best control law $K^{\textrm {I}}$ controlling the narrow-bandwidth regime. (a) Time series of the velocity measured downstream and voltage level at the terminals of the DBD actuator. (b) Power spectra of the velocity measured downstream. Black is for the unforced flow, red for the flow controlled with $K^{\textrm {I}}$ and blue for the flow controlled with the recorded actuation command $A^{\textrm {I}}(t)$ employed as an open-loop control. The spectra for the unforced flow (black) and the closed-loop controlled flow (red) are averaged over 20 realizations. The horizontal dashed lines denote the maximum of the spectra in the observation window (green background). The vertical axis is in $\log _{10}$ scale.

In addition, such control has been achieved with a minimum actuation power. Indeed the associated cost is ${J_b}^{\textrm {I}}=0.0063$, less than $1\,\%$ of the maximum actuation power. The actuation command, plotted in figure 11(a), shows that the control combines low-level steady actuation and low-amplitude feedback control. To demonstrate the effectiveness of the low-amplitude level, we use the closed-loop actuation command as an open-loop control. The closed-loop actuation command recorded during the learning process has been employed as open-loop control signal to force the flow. The spectrum of the resulting flow (blue spectrum in figure 11b) shows that that the main frequency of the flow $f_a$ resurfaces and also that the mode $n=3$ associated with the frequency $f^+$ is also excited. This open-loop test reveals that despite the low-amplitude level, feedback plays a crucial role in stabilizing the flow. Moreover, when we compare the results of the open-loop control with a steady actuation of equivalent level ($\sim 10\,\% A_{max}$) in figure 8(a), we note that both frequencies $f_b$ and $f^+$ are excited, indicating the effect of the unsteady component of the command.

Now we describe the learned law $K^{\textrm {I}}$ with an analytical approximation and a cluster-based visualization. For the analytical approximation, the determination coefficient for the affine reconstruction $R^2=0.87$ indicates an acceptable reconstruction. The gains associated with each feature component are presented in table 5. We note that, aside the mean component, the dominant term is $a_1=u$, as its gain ($k_1=0.51$) is more than 7 times higher than the second highest gain. The strong correlation between $K^{\textrm {I}}$ and $a_1$ reveals that phasor control or direct feedback of the system's state plays a fundamental role for this control. However, figure 12 shows that the relationship between $b$ and $a_1$ is not fully affine as two regions of significant width are displayed. This analysis is in agreement with the cluster-based investigation of the control law.

Figure 12.

Actuation command $b$ vs $a_1$ for the case: $K^{\textrm {I}}$ controlling the narrow-bandwidth regime.

Table 5.

Gains for the affine reconstruction of $K^{\textrm {I}}$ controlling the narrow-bandwidth regime.

The control visualization, following the method described in § 3.4, allows us to interpret the dynamics of the control achieved by $K^{\textrm {I}}$. The proximity map (figure 13a) shows that the centroids are arranged in a circular manner around the origin where centroid 1 is located. The transition probability matrix (figure 13b) gives the probability transition from one cluster to the other for each time step ${\rm d} t=0.004$ s. The transition probability matrix also shows one fixed point and 2 cycles: a small cycle including centroids 2, 3, 4 and 5 and a large one including centroids 6, 7, 8, 9 and 10. Combining the proximity map and the transition probability matrix, we can reconstruct the phase space of the dynamics by deriving a control network model (figure 13c). The identified cycles in the transition probability matrix are represented by limit cycles in the phase space. A frequency analysis based on a Poincaré section and angular first return map of the dynamics similar to Lusseyran et al. (Reference Lusseyran, Pastur and Letellier2008) reveals that the frequencies associated with the large limit and small limit cycles are around $28.48$ Hz and $41.14$ Hz, respectively. We can then assume that the large limit cycle is associated with the dynamics of mode $n=2$ (frequency $f_a=28.81$) and the small limit cycle is associated with mode $n=3$ (frequency $f^+=38.72$). Following Lusseyran et al. (Reference Lusseyran, Pastur and Letellier2008), such organization frequencies in the phase space shows that the dynamics may be structured around a fixed point of the stable spiral type where the oscillation period increases with the distance to the fixed point, here represented by cluster 1. Figure 13(c) also depicts the actuation regions in the phase space, revealing two mains regions of opposite actuation sign separated by a straight line. Interestingly, the sign of the actuation changes when approaching the fixed point. Such actuation map shows that the main stabilization mechanism exploited by $K^{\textrm {I}}$ is phasor control. The change of actuation sign when approaching the fixed point is explained by a phase shift due to the frequency change. The resulting control is similar to the one obtained with linear control by Yan et al. (Reference Yan, Debiasi, Yuan, Little, Ozbay and Samimy2006) where there is a rapid switching between two modes competing for the available energy, thus mitigating any resonance. Moreover, a spectral analysis of the actuation command $b$ and $a_1=u$ shows that they share the same frequency peaks supporting the phase relation between the control and the state of the system, see figure 26(a) in Appendix D.

Figure 13.

Investigation of the law $K^{\textrm {I}}$ controlling the narrow-bandwidth regime. Clusters are painted in different colours. The centroids are marked with numbers from 1 to 10. Panel (c) represents the control network model. The actuation amplitudes are denoted with bars, red (blue) for positive (negative) amplitude with respect to the mean actuation. The yellow boxes underneath indicate 25 % of the maximum amplitude. The black arrows serve as the most probable transition from one cluster to the other ($p>0.40$). The grey arrows are for less probable transitions ($0.40 \geq p \geq 0.05$). Lower probability transitions and self-transitions are omitted for clarity. The red (blue) background denotes the supposed regions of positive (negative) amplitude. The dotted black lines are the supposed limits that separate the regions. (a) Proximity map of the feature vector $\boldsymbol a$, (b) transition probability matrix and (c) control network.

Finally, a comparison between MLC and gMLC has been performed (see Appendix B) and reveals that gMLC outperforms MLC in terms of speed and final solution. In total, the learning has been accelerated by one order of magnitude.

Gradient-enriched MLC has been successfully applied to the stabilization of the open-cavity flow experiment. Exploration and exploitation phases both participated in the fast learning of a feedback control law. The evolution phases managed to discover new minima in the search space, and the simplex steps succeeded in converging towards a new minimum. A feedback control law is built, outperforming the steady actuation and allowing a similar reduction of the level of the maximum peak of the spectrum but with small actuation power. Both the analytical approximation and the cluster-based analysis hint at a control combining phasor control and nonlinear interactions. Hence, the control achieved is close to an ideal stabilization scenario, where some base state or fixed point is stabilized with a vanishing cost. This interpretation is certainly to be considered heuristically, given the complexity of the real dynamics of the 3-D intra-cavity flow and its nonlinear temporal and spatial interactions with the mixing layer. However, it aims at capturing the remarkable properties of the control law learned by the gMLC, whose mode of action is radically opposed to that of a control by steady forcing. In this section, the learning has been realized for a single-frequency regime; in the next section, a more challenging regime is controlled where two modes compete, strengthening nonlinear coupling.

4.4. Closed-loop control of the mode-switching regime

In this section, gMLC is employed to control a flow regime with strong nonlinear coupling, which leads, in particular, to intermittency between the two main instability modes of the mixing layer. The dynamics of intermittency is chaotic (Lusseyran et al. Reference Lusseyran, Pastur and Letellier2008) with the appearance of long time scales that are demanding from the point of view of machine learning. The goal is again to stabilize the flow by reducing the oscillations of the mixing layer, but this time in the case where two modes compete, as described in § 2.5. This constitutes a challenging problem as gMLC needs to learn a control law able to control two modes simultaneously. For the gMLC optimization, the same parameters as for the narrow-bandwidth regime have been employed, see § 4.2. In the following, the notations for the control law, cost and standard deviation associated with the learning on the mode-switching regime are marked by the superscript $\textrm {II}$.

Figure 14 depicts the cost of the individuals evaluated along the learning process. For this specific experiment, most of the learning is realized during the MCS phase, reducing the cost function to $J=0.0713$. From there on, the only improvements are carried out with the simplex steps, bringing the cost function to the final value $J^{\textrm {II}}=0.0565$. Nonetheless, the second exploration phase introduced a new control law ($\#11$) in the simplex (see table 8 of Appendix C). As the gMLC algorithm is partially stochastic, it is possible to fall close to the global minimum by pure luck; however, it usually takes several iterations of the learning phases to converge, as in § 4.3. Such a learning process has been observed in other realizations of the same experiment where a combination of exploration and exploitation have been necessary to reach similar levels of performance.

Figure 14.

Distribution of the costs for the 1000 evaluated individuals during the gMLC optimization process for the mode-switching regime. Each dot represents the cost $J$ of one individual. The colour of the dots symbolize how the individuals have been generated: black dots for randomly generated individuals (MCS phase), blue dots for individuals generated from a genetic operator (exploration phase) and yellow dots for the individuals arising from the subplex method (exploitation phase). The individuals from the MCS and exploration phases are sorted following their costs. The red line follows the evolution of the best cost. The vertical axis is in $\log _{10}$ scale.

After 1000 evaluations, the final control $K^{\textrm {II}}$ is a feedback control law, thoroughly described in Appendix C. The spectra of the flow under control (figure 15b) reveal a drastic decrease in the level of the maximum peak at frequency $f_a$. The dominant frequency is then close to the one associated with the second main mode ($\,f^+$). The relative reduction of the maximum peak in the spectrum is ${J_a}^{\textrm {II}}=0.0335$, i.e. 0.97 % of the fluctuation level. This corresponds to a reduction of the oscillation amplitude by a factor of 5 or by 15 dB. The standard deviation associated decreases to ${\tilde {\sigma }}^{\textrm {II}} = 97\,\%$. Like for the narrow-bandwidth regime, the control is achieved with small actuation power, using around $2\,\%$ of the maximum actuation power. Figure 16 reveals that the controlled flow no longer shows a mode-switching behaviour like in figure 4(b). The actuation command plotted in figure 15(a) shows short-time intermittent high-amplitude spikes emerging from the minimum actuation level ($A_{min}$). As for the narrow-bandwidth regime, the control law $K^{\textrm {II}}$ is re-evaluated 20 times and a small performance drop is observed: $J^{\textrm {II}}$ went from $-94\,\%$ to $-89\,\%$ and the standard deviation from ${\tilde {\sigma }}^{\textrm {II}}=97\,\%$ to $\tilde {\sigma }=112\,\%$. Finally, like in § 4.3, the time series of the actuation command has been employed as a signal input for open-loop control. Surprisingly, the equivalent open-loop control performs as good as the close-loop control law, suggesting that feedback was not at play in the reduction of the power amplitude for the mode-switching regime. However, anticipating the next section (§ 5), controlling the narrow-bandwidth regime with $K^{\textrm {II}}$ reveals that feedback is still a feature selected by gMLC.

Figure 15.

Characteristics of the best control law $K^{\textrm {II}}$ controlling the mode-switching regime. (a) Time series of the velocity measured downstream and voltage level at the terminals of the DBD actuator. A short window is enlarged to display the behaviour close to $A_{min}$ (the proportions are not kept). (b) Power spectra of the velocity measured downstream. Black is for the unforced flow, red for the flow controlled with $K^{\textrm {II}}$ and blue for the flow controlled with the recorded actuation command $A^{\textrm {II}}(t)$ employed as an open-loop control. The spectra for the unforced flow (black) and the closed-loop controlled flow (red) are averaged over 20 realizations. The horizontal dashed lines denote the maximum of the spectra in the observation window (green background). The vertical axis is in $\log _{10}$ scale.

Figure 16.

Spectrogram of the velocity for the mode-switching regime controlled with $K^{\textrm {II}}$.

For the polynomial approximation of the $K^{\textrm {II}}$ control of the mode-switching regime, a linear regression was unable to derive an affine reconstruction of the actuation command; the corresponding determination coefficient is $R^2=0.13$. Even with expanding the affine reconstruction with quadratic and cubic terms, $R^2$ is less than 0.25, implying that no linear control can be inferred from the data.

Regarding the cluster-based analysis of the control, figure 17 reveals a complex structure for the dynamics. The transition matrix (not plotted) shows that the cluster self-transitions are dominant. The reconstructed phase space reveals two regions of opposite actuation signs. The complex interactions between the centroids show that strong nonlinearities are at play. Interestingly, the spectrum of the actuation command (figure 26(b) in Appendix D) does not include any significant peak, except for the very low frequencies. The actuation command corresponds then to a random noise without any correlation with the measured velocity.

Figure 17.

Visualization of the control law $K^{\textrm {II}}$ controlling the mode-switching regime. The control network is depicted like in figure 13. The actuation amplitudes are denoted with bars, red (blue) for positive (negative) amplitude with respect to the mean actuation. The yellow boxes indicate 10 % of the maximum amplitude. The black arrows serve as the most probable transition from one cluster to the other ($p>0.15$). The grey arrows are for less probable transitions ($0.15 \geq p \geq 0.10$). Lower probability transitions and self-transitions are omitted for clarity. The red (blue) background denotes the supposed regions of positive (negative) amplitude. The dotted black line is the supposed limit that separates the regions.

In less than 1000 evaluations, gMLC builds feedback control laws that reduced the maximum peak of the power spectrum with small actuation power in two regimes: the narrow-bandwidth regime and the mode-switching regime. We proposed a visual representation of the control laws to aid the interpretability of the actuation mechanisms that enabled such efficient controls. However, a real analysis of the controlled flow is not done yet and constitutes a study in itself. The identification of the involved control mechanism requires a study of the short transient that leads to the stabilized state. Nonetheless, we can affirm that the control impacts the linear amplifier of the shear layer as the DBD actuator modifies its thickness on average. This effect has been demonstrated by studying the difference between the unforced and forced mean flow even for low actuation levels, i.e. near the ionization threshold. Such a mechanism is expected to remain valid for turbulent flows and even for often-studied transonic cavity flows. Moreover, Cornejo Maceda et al. (Reference Cornejo Maceda, Li, Lusseyran, Morzyński and Noack2021) show that gMLC surpasses MLC in terms of final solution performance and learning speed for a 2-D numerical simulation. In Appendix A, we demonstrate that gMLC also surpasses MLC in experiments, establishing gMLC as a keystone for fast learning of feedback control laws directly on the plant. We foresee that gMLC will greatly contribute to the learning of control laws for MIMO control.

5.1. Robustness of the control

In this study, feedback control laws are optimized for steady conditions: fixed Reynolds number and incoming velocity during the learning process. To achieve robustness, the control laws learned should also perform for a large range of parameters. Thus, to test the robustness of the learned laws, they are re-evaluated for operating conditions different from the learning ones: the law $K^{\textrm {I}}$ optimized for the narrow-bandwidth regime is now employed to control the mode-switching regime and vice versa. Such a test is demanding as the dynamics is partially different from one regime to the other. In particular, the unforced dynamics of the mode-switching regime includes the main frequency of the narrow-bandwidth regime but also displays intermittency with another frequency.

Figure 26 displays the spectra for each controlled case. For the control of the narrow-bandwidth regime with $K^{\textrm {I}}$ (figure 18(a), green line), we note that the law $K^{\textrm {II}}$ manages to reduce the main peak $f_a$ of the spectrum to the same level as the law $K^{\textrm {I}}$ (red line). Moreover, for the frequency range around $f^+$, $K^{\textrm {II}}$ performs better than $K^{\textrm {I}}$ as it mitigates the peak-splitting phenomenon described in § 4.3. Such feature is expected as $K^{\textrm {II}}$ has been optimized to control the mode $f^+$. In fact, the whole spectral range is reduced in amplitude. This non-transfer of energy to other frequencies is a remarkable feature of $K^{\textrm {II}}$.

Figure 18.

Spectral response of the flow controlled by $K^{\textrm {I}}$ and $K^{\textrm {II}}$ for the narrow-bandwidth regime (top) and the mode-switching regime (bottom). Black spectra correspond to the unforced dynamics of each regime; red spectra correspond to the flow controlled by law learned in the given regime, i.e. $K^{\textrm {I}}$ ($K^{\textrm {II}}$) for the narrow-bandwidth (mode-switching) regime. Green corresponds to the flow controlled by the law learned in the other regime and blue to the open-loop equivalent of the latter one. The horizontal dashed lines denote the maximum of each spectrum in the observation window (shaded green section). The vertical axes are in $\log _{10}$ scale. (a) Power spectra for the narrow-bandwidth regime and (b) power spectra for the mode-switching regime.

As for the control of the mode-switching regime with $K^{\textrm {I}}$ (figure 18(b), green line), the power of the main frequency $f_a$ is drastically reduced, reaching the same power level as the control with $K^{\textrm {II}}$ (red line). However, the controller is less efficient than $K^{\textrm {II}}$, as $K^{\textrm {I}}$ fails to reduce the power associated with the frequency range surrounding the frequency $f^+$ of mode $n=3$. Again, like for the control of the narrow-bandwidth regime with $K^{\textrm {I}}$, we observe the spillover effect with a splitting of the main frequency into two frequencies on either side and with a lower power level.

In summary, both learned control laws $K^{\textrm {I}}$ and $K^{\textrm {II}}$ are able to retain their efficiency when controlling regimes that are out of the learning conditions. Expectedly, $K^{\textrm {I}}$ was unable to reduce the peak of the frequency $f^+$ in the mode-switching regime. Nonetheless, it manages to reduce the $f_a$, despite only appearing intermittently. On the other hand, for the narrow-bandwidth regime, the control with $K^{\textrm {II}}$ was more significant than $K^{\textrm {I}}$ as it prevents the peak splitting (or spillover) of the third mode. Thus this test also reveals that $K^{\textrm {II}}$ is not only able to control the frequency $f_a$ and $f^+$ but also to prevent the rise of both modes simultaneously. Of course, from the point of view of the cost function, both control laws $K^{\textrm {I}}$ and $K^{\textrm {II}}$ are similar when controlling the narrow-bandwidth regime, and gMLC could have converged towards any of the two control laws. Yet, $K^{\textrm {II}}$ is the control law that answers the best to the control objective, which is the stabilization of the flow. Therefore, it is the learning conditions that make the difference. This analysis reveals that learning a control law in complex and rich conditions is beneficial for the robustness and the overall efficiency of the control as the richness of the dynamics will be reflected in the control law.

5.2. Interpretation of the resulting controlled flow

Now, we propose an interpretation of the controls performed in the previous section using both analytical approximation and cluster-based visualization of the control laws. Firstly, we analyse the case where $K^{\textrm {II}}$ controls the narrow-bandwidth regime. Like for § 4.4, a linear regression is unable to reconstruct the actuation command despite being in a less complex regime. The determination coefficient is $R^2=0.13$. The addition of quadratic terms brings the coefficient no higher than 0.78. A complex nonlinear control is expected as $K^{\textrm {II}}$ manages to control both frequencies $f_a$ and $f^+$.

As for the control law visualization, figure 19 depicts a similar control network as in § 4.3. However, there is only one large cycle composed of centroids 1, 2, 3, 4, 5, 6 and 7. The limit cycle presents four phases regularly alternating between positive and negative actuation, suggesting that the control operates at twice the frequency of the flow. A spectral analysis of $a_1=u$ and $b$ (see figure 26(c) in Appendix D) shows that the main peaks are respectively $f_a=29.03$ Hz and $f=58.23\ {\rm Hz} \approx 2f_a$, confirming the phase relation between the flow dynamics and the actuation command. Interestingly, the second harmonic $2f_a$ is slightly excited but does not resonate for the controlled flow (figure 18(a), green line), while it resonates for the open-loop equivalent of $K^{\textrm {II}}$ (blue line). In summary, $K^{\textrm {II}}$ controls the flow at twice the main frequency but avoids the resonance of the second harmonic. The efficiency of a controller at twice the main frequency has been previously reported by Schumm, Berger & Monkewitz (Reference Schumm, Berger and Monkewitz1994). The authors note the stabilization of a cylinder wake by transverse vibration of the cylinder at 1.8 times the natural shedding frequency. In particular, they declare that the control effect is due to a nonlinear interaction between the instability and the forcing input. Thus, both the analytical and cluster-based analyses point towards a nonlinear actuation mechanism for the control of the narrow-bandwidth regime with $K^{\textrm {II}}$.

Figure 19.

Visualization of the control law $K^{\textrm {II}}$ controlling the narrow-bandwidth regime. The control network is depicted as in figure 13. The actuation amplitudes are denoted with bars, red (blue) for positive (negative) amplitude with respect to the mean actuation. The yellow boxes indicate 25 % of the maximum amplitude. The black arrows serve as the most probable transition from one cluster to the other ($p>0.50$). The grey arrows are for less probable transitions ($0.50 \geq p \geq 0.5$). Lower probability transitions and self-transitions are omitted for clarity. The red (blue) background denotes the supposed regions of positive (negative) amplitude. The dotted black lines are the supposed limits that separate the regions.

Secondly, we interpret the case where $K^{\textrm {I}}$ controls the mode-switching regime. This time, the linear regression builds an affine approximation of the control as the determination coefficient is $R^2=0.92$. The gains associated with each feature component are displayed in table 6. As for § 4.3, the most relevant feature is $a_1=u(t)$, but again the control cannot be reduced to an affine relationship as figure 20 displays a nonlinear curve. Such observation is in agreement with the spectral analysis of $b$ and $a_1=u$ (see figure 26(d) in Appendix D) showing the peaks at the same frequency. The complexity of the flow translates into a complex control network as in § 4.4. Figure 21 depicts a reconstructed phase space divided into two main regions: one on the left (centroids 1, 2 and 3) with positive actuation amplitude; and one on the right (centroids 4, 5, 6, 7, 9 and 10) with negative actuation amplitude. The role of centroid 8 may be small as its associated actuation is close to the mean value. Interestingly, the overall structure of the control network is similar to the one in figure 21, suggesting that the control mechanism is also a complex nonlinear one.

Table 6.

Gains for the affine reconstruction of $K^{\textrm {I}}$ controlling the mode-switching regime.

Figure 20.

Actuation command $b$ vs $a_1$ for the case: $K^{\textrm {I}}$ controlling the mode-switching regime.

Figure 21.

Visualization of the control law $K^{\textrm {I}}$ controlling the mode-switching regime. The control network is depicted as in figure 13. The actuation amplitudes are denoted with bars, red (blue) for positive (negative) amplitude with respect to the mean actuation. The yellow boxes indicate 25 % of the maximum amplitude. The black arrows serve as the most probable transition from one cluster to the other ($p>0.15$). The grey arrows are for less probable transitions ($0.15 \geq p \geq 0.9$). Lower probability transitions and self-transitions are omitted for clarity. The red (blue) background denotes the supposed regions of positive (negative) amplitude. The dotted black lines are the supposed limits that separate the regions.

5.3. The need of feedback

We have shown with the control of the narrow-bandwidth regime (§ 4.3) that without feedback the control law $K^{\textrm {I}}$ is unable to stabilize the flow and even excites the frequencies $f_a$ and $f^+$ (figure 11b). Moreover, $K^{\textrm {I}}$ is able to partially control the mode-switching regime (§ 5.1) and, as expected, the same controller applied in an open-loop manner is no longer able to control the main mode $f_a$ (blue spectrum in figure 18b). We note, nonetheless, a small shift of the spectrum towards the higher frequencies.

On the other hand, surprisingly, it has been shown in (§ 4.4) that $K^{\textrm {II}}$ performs in open loop as well as in closed loop. This result is consistent with the absence of correlation between the actuation command $b$ and the velocity measure $a_1=u$, see figure 26(b) in Appendix D. So, it seems that the learning process has selected the flow information ($a_i$, see Appendix C) without any improvement of the cost. However, $K^{\textrm {II}}$ applied in closed-loop manner to the narrow-bandwidth regime performs even better than $K^{\textrm {I}}$, while in an open-loop manner, $K^{\textrm {II}}$ fails to achieve any control. Indeed, the corresponding spectrum (blue spectrum in figure 18b) is almost similar to the unforced flow. Therefore, it should be noticed that the flow state information $a_i$ in $K^{\textrm {II}}$ is truly functional and, in fact, gives the controller the ability to remain effective well away from the learning conditions.

This analysis shows the extent to which feedback is a key feature for control. We believe that the ability of gMLC to learn effective and efficient feedback control laws in experiments will greatly benefit future MIMO control experiments.

5.4. Stabilization of the open cavity flow

In this study, the flow is monitored by a single hot-wire sensor downstream of the actuator. The sensor signal is employed for sensor feedback and to characterize the controlled flow. The achieved stabilization near the sensor extends in the spanwise and streamwise directions, assuming a significant portion of the finite aspect ratio cavity.

For large aspect ratios $S/D$, e.g. $O(100)$ or more, the effect of 2-D actuation along the whole span can be expected to depend on the sensor location. The feedback stabilization in the sensor plane will become a non-stabilizing open-loop actuation far beyond the spanwise coherence length. An interesting example has been reported for the stabilization of a large aspect ratio cylinder wake. Roussopoulos (Reference Roussopoulos1993) forced the cylinder wake with a pair of loudspeakers driven in opposite phase and significantly reduced the fluctuations at the downstream sensor location. Far away from the sensor in spanwise direction, no stabilization was observed. For our open cavity flow, we expect a loss of control authority at a distance from the hot-wire position greater than the transverse coherence length. This transverse coherence of the mixing layer instabilities is at least of the order of the cavity depth $D$. Our cavity has an aspect ratio of $S/D=6$, i.e. we control at least one third of the cavity spanwise ($1D$ on either side of the hot-wire plane). The remaining lateral thirds are in the Ekman layers and, therefore, probably less oscillatory. The global stabilization could be augmented with multiple actuators and multiple sensors. As for the spanwise homogeneity of our one-piece DBD actuator, measurements performed for actuation levels close to the ionization threshold, when the ionization is still quite inhomogeneous along the electrode, show that the flow response is independent of the spanwise location of the measurement point. It should be noted that the feedback is only provided by a sensor in the symmetry plane. Yet, the global stabilization of the cavity flow is validated with additional measurements at several spanwise locations as elaborated in Appendix E.

Concerning the streamwise direction, the source of the oscillation is related to the mixing layer through the Kelvin–Helmholtz instability. At the level of the mixing layer, the hypothesis of an oscillation along the length of the cavity which would be cancelled at the downstream edge corresponds to the idea of a standing wave with a node at the hot-wire location. However, with a Kelvin–Helmholtz instability, we are in the case of convective instability which can only be killed by cancelling the disturbances at the source. Otherwise, any oscillations close to the most dangerous frequency would inevitably be amplified and be present at all $x$ and especially at the hot-wire location. This is well observed in the visualization of the natural flow where the oscillations reach maximum nonlinear amplitudes and break on the downstream corner. Hence, the downstream stabilization of the flow results necessarily in the control of the mixing layer in the streamwise direction.