Linear feedback control of transient energy growth
and control performance limitations
in subcritical plane Poiseuille flow
Abstract
Suppression of the transient energy growth in subcritical plane Poiseuille flow via feedback control is addressed. It is assumed that the time derivative of any of the velocity components can be imposed at the walls as control input, and that fullstate information is available. We show that it is impossible to design a linear statefeedback controller that leads to a closedloop flow system without transient energy growth. In a subsequent step, fullstate feedback controllers – directly targeting the transient growth mechanism – are designed, using a procedure based on a Linear Matrix Inequalities approach. The performance of such controllers is analyzed first in the linear case, where comparison to previously proposed linearquadratic optimal controllers is made; further, transition thresholds are evaluated via Direct Numerical Simulations of the controlled threedimensional Poiseuille flow against different initial conditions of physical interest, employing different velocity components as wall actuation. The present controllers are effective in increasing the transition thresholds in closed loop, with varying degree of performance depending on the initial condition and the actuation component employed.
I Introduction
Transient energy growth has been recognized as a possible mechanism explaining subcritical transition in wallbounded flows; in fact, subcritical flows may experience large transient amplifications of the energy of perturbations, that could trigger nonlinear mechanisms and eventually lead to transition to turbulencebutlerfarrel1992 (); reddyhenningson1993 (); schmidhenningson1994 ().
In viscous shear flows, transient energy growth is related to the nonnormality of the linearized NavierStokes operator with respect to the energy inner productschmidhenningson2001 (); schmid2007 (). In the last few years, several investigators attempted to reduce the transient growth phenomenon in Poiseuille and boundary layer flows by employing wall actuation and applying linear control theory to an appropriate discretization of the linearized equations. In their seminal work on feedback control of instabilities in twodimensional Poiseuille flow, Joshi et al.joshispeyerkim1997 () employed a compensator in the form of a constantgain integral feedback, and demonstrated stabilization of the linearly unstable flow as well as attenuation of finite amplitude disturbances. They further pointed out that transient amplifications in the flow energy may not be properly detected by the sensors, and also that the control itself may trigger nonlinear mechanisms by introducing transient disturbances on short times. Leveraging a statespace formulation obtained after discretization of the boundarycontrolled OrrSommerfeldSquire equations, optimal and robust control theory was applied to transitional channel flows by Bewley & LiuBewley1998305 () for a single wavenumber pair and by Högberg et al.Hogberg2003149 () for a large array of wavenumber pairs, leading to a reduction of the maximum transient energy growth as well as an increase in transition thresholds. It has been recently shownkimlim2000 () that the linear coupling term in the OrrSommerfeldSquire equations plays a role not only in the nonnormal behavior of the small perturbation dynamics, but also in the selfsustaining process of nearwall, low Reynolds number turbulence. This evidence led investigators to test in turbulent channel flows the optimal controllers designed on linearized flow models, and encouraging results have been obtained in terms of drag reductionleeetal2001 (); bewleymointemam2001 (); kim2003 (); bewley2001 (). Feedback control of nonmodal disturbances in boundary layer flows has been recently considered by Corbett & Bottarocorbettbottaro2001 () in the framework of optimal control theory, while Zuccher et al.zuccherluchinibottaro2004 () applied steady suction in the attenuation of the growth of given optimal disturbances in a Blasius boundary layer.
Although it has been demonstratedbewley2001 () that optimal and robust control laws are well suited for reducing the nonnormal behavior of fluid flow systems, to date no feedback control law has been devised with the capability of ensuring closedloop monotonic stability, when boundary actuation is employed. It is therefore natural to ask whether such performance can be obtained with a linear feedback, and further which control techniques are available to directly target the transient growth mechanism. In the controls literature, the transient amplification of certain norms of the state bears the name of peaking phenomenon and the monotonic stability requirement is generally referred to as strict dissipativity. Active controllers with the capability of targeting transient norm amplifications have received attention in the analysis of a class of partially linear cascade systemssussmannkokotovic1991 (), and more recently in conjunction with a Linear Matrix Inequality (LMI) approachplischkewirth2004 (); plischke2005 (). In a very recent paper, Whidborne & McKernanWhidborne20071762 () extended these results giving conditions on the existence of a feedback controller ensuring the strict dissipativity of the closedloop system. These results were then exploited by Whidborne et. alWhidborneMckernanPapadakisJSCE2008 (), who considered the feedback control of a single wavenumber pair in plane Poiseuille flow via LMI design and wallnormal blowing and suction, and by Martinelli et al.martinellietal2009 (), where LMIbased feedback controllers have been designed and tested for an array of wavenumber pairs. The present paper aims at expanding and completing these recent results, by showing first that it is impossible to design a linear statefeedback controller ensuring the plane Poiseuille flow – controlled via wall transpiration with any velocity component – to be strictly dissipative. In a second step, feedback control laws are designed using an LMI technique, for an array of wavenumber pairs; their performance is compared against that of optimal controllers in the linear case, and furthermore closedloop transition thresholds are evaluated for optimal initial conditions in the form of a pair of oblique waves and antisymmetric streamwise vortices, at the Reynolds number , using different velocity components as wall actuators.
Ii Model of the system
We consider the dynamics of threedimensional small perturbations to the laminar Poiseuille solution in a plane channel. A Cartesian coordinate system is introduced, where , and denote the streamwise, wallnormal and spanwise directions, and , , denote the corresponding perturbation velocity components. The NavierStokes equations, linearized about the laminar solution , are nondimensionalized with the centerline velocity and the channel halfwidth , and rewritten in the form of a single equation for oneway coupled to an equation for the wallnormal vorticity . Fourier transformation in and direction yields the well known OrrSommerfeldSquire form:
(1) 
at the wavenumber pair . Here, the tilde denotes Fourier coefficients, the dot denotes time derivative, the prime denotes differentiation, , is , and .
We select boundary conditions representing timevarying wall transpiration on any of the velocity components at the two channel walls (“vectorized transpiration”). In turn, this results in inhomogeneous Dirichlet and Neumann conditions on , as well as inhomogeneous Dirichlet conditions on :
(2)  
A standard lifting procedureHogberg2003149 () is then employed, i.e. the unknowns are rewritten as homogeneous components (satisfying homogeneous boundary conditions) plus inhomogeneous components thus:
(3) 
where are polynomials in chosen to satisfy unitary boundary conditions for each lifted component appropriately.
We discretise the homogeneous components and in wallnormal direction using a modified Chebyshev series cardinal function basis
(4)  
where the modified Chebyshev functions and implicitly enforce the required homogeneous Dirichlet and Neumann boundary conditions, lead to good conditioning of the discrete Laplacian operator, and no spurious modes are generated thesis (). The OrrSommerfeldSquire equations are then evaluated on a set of GaussLobatto collocation points in direction and rearranged to have the time rate of change of actuation velocity as an inputHogberg2003149 (). This results in the the linear timeinvariant system modelSkogestad01 ()
where , are constant system and input matrices, and u, are respective input and state vectors. As noted in previous workHogberg2003149 (), the natural outcome of this procedure is an augmented statespace form, where additional integrators associated to the values of the input velocity components at the walls are explicitly introduced. Consequently, the openloop dynamics of the system above (i.e. setting ) is different from the dynamics of the original, unactuated system (where the velocity at the walls is fixed by the noslip condition). This terminology will be used throughout the paper to distinguish between the two cases.
The kinetic energy per unit mass of flow perturbations in the volume
can be expressed as a function of the state vector using the continuity equation, the definition of , and Parseval’s identity:
where the matrix is an Hermitian, positive definite matrix; here, the superscript denotes conjugate transpose. In the following, we have transformed the state vector via the change of variable (C being the Cholesky factor of Q) in order to rewrite the system dynamics as
(5) 
such that the system energy is directly given by the Euclidean norm Skogestad01 (); schmid2007 ().
Iii Closedloop monotonic stability
We consider the linear timeinvariant system model (5) and further assume that , that is B has full column rank (i.e. all the actuators are independent – a condition which is trivially satisfied in the present problem). Contraction analysis for this kind of system has been presented by Whidborne & McKernanWhidborne20071762 () and, in the general case of nonlinear systems, by Lohmiller & Slotinelohmillerslotine1998 (). In particular, referring to the linear case, it has been shownWhidborne20071762 () that there exists a static statefeedback controller , where K is a constant matrix, such that the closedloop system has strict dissipativity (i.e. energy decays monotonically from all initial conditions ), if and only if
(6) 
where is the left null space of B. AdditionallyWhidborne20071762 () if no static controller that achieves strict dissipativity exists, then no dynamic statefeedback controller, where u is given from x by the dynamic system
(7)  
(8) 
where , and are constant matrices and are controller states, exists either.
It is immediate to verify that the second criterion in (6) is never satisfied in the present system, as the Hermitian matrix is never positive definite but it is always positive semidefinite, because the dimension of the input vector is always smaller than that of the state vector. In order to have , a number of independent actuators equal to the number of flow states is required; this is a situation that is unlikely to occur in practical flow control problems, where normally actuators are placed at the walls. Even when volume forcing is available, this condition is unlikely to be satisfied, since practical volume forces are not as flexible as to enforce an arbitrary force distribution in the entire flow domain at any time instant.
The first algebraic criterion in Eq. (6) is equivalent to requiring that the portion of the system dynamics which is not accessible by the controls must be dissipative. Verifying this criterion is not trivial. Here we evaluate it numerically, in order to identify those regions in the parametric space where subcritical Poiseuille flow may be rendered monotonically stable by feedback transpiration. To this aim, the statespace model (5) is first obtained on a fine grid and, as suggested by Reddy & Henningsonreddyhenningson1993 (), a limited number of eigenfunctions is retained, discarding those corresponding to highly damped and poorly resolved eigenvalues. Properly rescaling the variables such that energy is written as an Euclidean norm leads to a reduced order model ,, and the negativedefiniteness of the corresponding matrix in (6) is verified by computing its maximum (real) eigenvalue . Figure 1 shows the present result on the plane for , along with the wellknown result on the transient growth dependence in plane Poiseuille flowreddyhenningson1993 () (i.e. the unactuated case). The white area corresponds to the domain where the unactuated system is monotonically stable, while the shaded area is the region where the unactuated system admits transient energy growth. Solid lines correspond to isocontours of , and it appears that the contour lies on the very boundary between the shaded and white area, implying that the Hermitian matrix is indefinite when the unactuated system is not monotonically stable. An analogous result is reported in Fig. 2, where isocontours of are reported on the plane , at ; again, the contour lies on the boundary between the regions of monotonic and nonmonotonic stability. From the aforestated theoremWhidborne20071762 (), this implies that it is not possible to design a statefeedback controller that ensures the closed loop Poiseuille flow to be monotonically stable, when the corresponding unactuated flow is not.
This result shows an inherent limitation in the feedback control of the transient growth mechanism, when vectorized wall transpiration in terms of time rate of change of zero netmass flux blowing/suction is employed. Note that vectorized transpiration, although being rather idealized, exploits all the degrees of freedom available for boundary control in the present problem; therefore, the present result is representative of a limiting situation. In a practical setting, actuatordependent constraints may introduce additional mechanisms restricting the control authority even more. It is also worth mentioning that strict dissipativity is a rather tough requirement for a controller bound to operate on a largely underactuated flow, and in fact an approach based on strict dissipativity is in general quite conservative, i.e. some energy growth is tolerable in transition control. Finally, it should be emphasized that the present analysis is limited to linear feedback laws, and that performance of nonlinear controllers may be more promising. For example, it has been shownhogbergbewleyhenningson2003 () that introducing a nonlinearity in the form of gainscheduling on fullstate feedback laws led to relaminarization of low turbulence even employing wall actuation only, and that adjointbased optimization on the nonlinear turbulent flow can be successfully employed in feedback relaminarization bewleymointemam2001 ().
Iv Upperbound minimizing feedback controller
In order to design a statefeedback controller with the capability of targeting the transient growth mechanism directly, an estimate of the maximum transient growth is required. Such estimate is obtained as an upper bound on the maximum growth via Lyapunov theory. For the linear, time invariant, asymptotically stable system:
it can be shown that an upper bound on the maximum transient growth is given by Whidborne20071762 (); WhidborneMckernanPapadakisJSCE2008 () :
where satysfies the Lyapunov inequality
A minimal upper bound can be obtained by solving the following minimization problem boydetal1994 () :
(9)  
where the last inequality ensures . The problem stated in Eq. (9) is a LMI generalized eigenvalue problem, and standard solution methods based on interior point algorithms are availablenesterovnemirovskii1994 ().
An analogous problem to that stated in Eq. (9) can be obtained if the feedback minimization of the upper bound is of interest. Indeed, let us consider the system (5) along with a statefeedback control law in the form , so that in closedloop the system dynamics is described by:
Leveraging the additional degrees of freedom due to the controller gains , we move to minimizing the closedloop upper bound on the maximum transient growth. The associated Lyapunov inequality now reads:
this inequality can be rewritten in the LMI form by recalling that a similarity transformation preserves the eigenvalues. Therefore, defining and , the closedloop upperbound minimization problem can be written as boydetal1994 ()
(10)  
where the last, additional inequality ensures a limit in the control effort in the form . The problem (10) has to be solved for , and ; controller gains are obtained from . This problem is again a LMI generalized eigenvalue problem, that can be solved using standard methods nesterovnemirovskii1994 ().
V Results and discussion
LMI controllers are designed wavenumberwise, using the model of the system (5) and the design equations (10). In particular, we consider each actuation component (, or on both walls) independently; a limit on the control effort , kept constant in wavenumber space, is used in the design of all controllers. In the design procedure, the linear equations pertaining to each wavenumber pair are discretized using Chebyshev polynomials, and modal truncation at (the maximum affordable size of the computational problem) is performed prior to the actual solution of the generalized eigenvalue problem (10). In addition to removing poorly resolved dynamics, modal truncation proved to be necessary due to the exacting memory requirements of the existing LMI solvers (scaling as ); in performing modal truncation, it was thoroughly verified that the reduced order model preserves the linear transient energy growth of the unactuated system.
v.1 Linear analysis
The performance of LMI controllers is evaluated first in the linear setting at ; in particular, perturbations at two representative wavenumber pairs are considered, namely, an oblique wave and a streamwise vortex . In the results reported here, the effectiveness of different actuation components is also addressed.
The analysis reported in sec. III shows that the closedloop system will have a nonnormal behavior; it is therefore natural to contrast the maximum closedloop transient growth with the openloop one, in order to verify that a consistent reduction in is obtained via the minimization in (10). Results, obtained for the same truncated system used in the design and for the two wavenumber pairs considered, are reported in table 1. It is shown that, in both the oblique wave case and the streamwise vortex case, solution of (10) leads to a closedloop system experiencing a reduced maximum transient energy growth. In particular, for the oblique wave case, is the most effective actuation component, whereas actuating with is most effective in the streamwise vortex case.
Oblique waves  Streamwise vortices  

Openloop  Closedloop  Openloop  Closedloop  
 actuation  
 actuation  
 actuation 
The performance of the present controllers is further compared against fullstate controllers designed with the LQR approachSkogestad01 (), considered in a similar transition problem by Hogberg et al.Hogberg2003149 (). The aim of the LQR control is the minimization of the time integral of the perturbation energy, while keeping the time integral of the control effort as low as possible; in fact, the control objective is given in terms of the closedloop minimization of a functional in the form
This is substantially different from the control objective of the present LMI formulation (10), which considers bounds on the disturbance energy and control expenditure; therefore – at a fixed control expenditure – LMI controllers can be used to estimate a possible best performance (in terms of peaking suppression) of other control strategies. In order to present a fair comparison between the LQR and the LMI formulation, we iteratively design and test a LQR controller keeping fixed (the same used in LMI design) and as a free parameter, and we evaluate the integral of the control energy in closed loop until it matches the value computed for the LMI controller. The closedloop systems, controlled via both LQR and LMI gains, are tested against the respective optimal perturbations, using actuation for the oblique wave case and actuation for the streamwise vortex case (the best performance cases reported in table 1); further, for the streamwise vortex we consider the antisymmetric (with respect to the plane) optimal perturbation. The time evolution of the perturbation energy is displayed in fig. 3 and 4. Results show that, at a given control expenditure, the worstcase initial condition for the LMIcontrolled system experiences a lower amplification than the corresponding perturbation for the LQRcontrolled system; the peak for the LMIcontrolled system occurs at later times for the streamwise vortex case. Despite the mild reduction in maximum amplification, the results shown here suggest that a control design technique directly targeting the growth mechanism is able to better exploit the degrees of freedom in the controller to achieve a minimal transient peaking of the energy.
The performance of LMI controllers is also evaluated against optimal initial conditions for the unactuated flow, and results are reported in fig. 5 and 6. Since the closedloop system has additional state equations associated with the dynamics of wall velocity components used as actuators, these values are set to zero, assuming that at initial time the openloop flow satisfies the noslip and notranspiration condition at the walls. In particular, for the oblique wave case, it is shown that actuating with reduces the maximum amplification of the optimal disturbance by a factor , whereas a less effective reduction (by a factor ) is obtained using or . Further, the growth curves in these latter cases are very close to each other, as the effect of actuators on the oblique wave is symmetric. In the antisymmetric streamwise vortex case, the most effective components are and (reduction by a factor and , respectively), whereas has a quite poor performance (amplification reduced by a factor ). The differences in performance between and may be interpreted with a geometric argument, as for a streamwiseinvariant perturbation the component acts in a weakly controllable direction. It is also noteworthy that, in all these cases, the control action reduces the time interval after which the initial disturbance gets to its maximum amplification.
It is finally worth emphasizing that, when the closedloop system experiences initial conditions in the form of optimal perturbations for the unactuated flow, the LMI and LQR performance is practically equivalent, for a given global control effort. Considering the best performing LMI controllers ( and actuation, respectively), results are given in fig. 7 and 8, for the oblique wave and streamwise vortex. For the optimal oblique wave (fig.7), the linear evolution of the perturbation energy using LQR control matches almost perfectly that obtained with the LMI controller. In the case of antisymmetric streamwise vortex, we obtain a slightly larger maximum amplification for the LQR, that is however followed by a faster transient to zero if compared to the LMI case. The results reported in fig. 7 and 8 are substantially independent on further decrease of the value of : no significant changes in the time evolution of the perturbation energy are obtained, but at a far larger expense. This indicates that these results are close to the limit of small control weight.
v.2 Closedloop transition thresholds
After these numerical experiments in the linear setting, the performance of LMIbased controllers in terms of transition delay capabilities has been verified using Direct Numerical Simulations (DNS) of transitional Poiseuille flow at , using an existing computer code and computing system luchiniquadrio2006 (). Controllers are tested against initial conditions in the form of:

a pair of oblique waves , in a box of size ;

antisymmetric streamwise vortices , in a box of size .
These initial conditions are obtained by computing the optimal perturbations for the unactuated flow at the corresponding wavenumber pairs, and specifically for the streamwise vortices the optimal initial condition having antisymmetric (with respect to the centerplane ) distribution of wallnormal velocity is considered, as it provides the optimal transition time in the nonlinear casereddyetal1998 (). Random noise, in the form of a random combination of the first Stokes modes – ordered by decreasing real part of the corresponding eigenvalues – is added on the wavenumber array and , and the noise energy is chosen as of the total perturbation energy. The resulting optimal perturbations are identical to those reported in previous worksreddyetal1998 (); Hogberg2003149 ().
In order to reduce the computational problem of control design to an affordable size, LMI controllers are designed on the same array of wavenumber pairs where random noise is introduced. Furthermore, the control effort tuning parameter is set at , a value which is derived from preliminary tests and previous work on the control of the linearized dynamics of streamwise vorticesthesis (). The value of the parameter has been kept constant in wavenumber space; however, it should be emphasized that this parameter could be a function of the wavenumber pair – thus providing room for optimization of the control performance.
The performance of LMI controllers is quantified by evaluating the closedloop transition thresholdreddyetal1998 () of a given initial condition, for all the wall actuation components. The mixed spatial discretization (fourthorder, compact finite differences in direction and Fourier expansion in and directions) employs grid points in and modes in and . Time integration is performed via the usual semiimplicit approach, where nonlinear terms are advanced explicitly using a lowstorage RungeKutta algorithm whereas linear diffusion terms are advanced implicitly via a CrankNicholson scheme. Each simulation was run over a time window of nondimensional time units, that proved sufficiently long to ensure that a laminar or a turbulent state was reached after the initial transient growth of the perturbation energy. In order to obtain the thresholds reported in table 2, a bisection algorithm was employed; this procedure requires a large number of simulations, corresponding to approximately months of CPU time.
The openloop transition thresholds reported in table 2 agree with previous findingsreddyetal1998 (). Results summarized in the table indicate that the LMI controller is able to increase the transitional energy of the initial conditions considered. In particular, a synthetic performance measure is indicated in the table as improvement factor (I.F.), corresponding to the ratio between the threshold energy computed in the controlled case over that corresponding to the unactuated flow. In general, for both the oblique waves and the streamwise vortices, actuation with the wallnormal velocity outperforms actuation with the other components. This behavior is expected, as forcing with wallparallel components affects the flow by means of viscous diffusion only, whereas forcing with introduces an additional nonzero momentum flux at the boundary. The improvement factors associated with the oblique wave case are higher than those pertaining to the streamwise vortex case, when using and actuation. For the component case, this is coherent with previous worksHogberg2003149 (), and can be interpreted physically by the argument that targeting oblique waves mitigates the subsequent development of streamwise vortices, therefore reducing the strength of the associated streak instability. The component provides the overall worst performance in the streamwise vortex case, a result that can be interpreted as a consequence of the particular geometrical configuration. In fact, with respect to a streamwiseinvariant spatial structure, actuation works in an approximately null direction, whereas actuation with is bettersuited, as shown by its improvement factor. A similar geometrical interpretation can be given for the almost equal improvement factors obtained with and , when an pair of oblique waves is given as initial condition.
Oblique waves  Streamwise vortices  

Thres.  I.F.  Thres.  I.F.  
Openloop  
 actuation  
 actuation  
 actuation 
Results reported in table 2 for the thresholds obtained using actuation are in qualitative agreement with previous work using the LQR controllersHogberg2003149 (); however, in quantitative terms the LQR approach outperforms the present LMI approach. In particular, the improvement factors reported with LQR for oblique waves and streamwise vortices at are and , respectivelyHogberg2003149 (); therefore, LQR controllers seem to perform substantially better than LMI controllers in presence of oblique waves as initial conditions. It should be emphasized, however, that such comparison is not entirely appropriate. In fact, even if the same energy norm is used to quantify the magnitude of velocity disturbances, control laws are designed with different parameters constraining the control effort. As a matter of fact, the linear results reported in fig. 5 and 6 show that the LMI controller performs similarly to the LQR controller for a value of the LQR control weight and , respectively. These values are different from the value of , used uniformly in wavenumber space in Högberg et al.Hogberg2003149 (). A fair comparison between the two approaches is not possible in this case; therefore, it is impossible to draw a conclusive statement about the effectiveness of feedback minimization of transient growth versus feedback minimization of the disturbance energy in transition delay.
As the ultimate goal of LMI controllers is that of preventing transition to turbulence, it can be important to quantify the energy efficiency of these controllers in the nonlinear case. For instance, given a transitional initial condition, it is possible to compare the energy expenditure of the controller to prevent transition with the additional energy to be introduced into the unactuated flow to compensate for the increase in friction, over the same time window (i.e. the time necessary for the transient in the controlled flow to die out). Referring to the best performing cases in table 2, we consider actuation with wallnormal velocity and, in the two cases, initial conditions having energy below the corresponding closedloop threshold. Using the present nondimensionalization, a conservative estimate of the energy required for the control action in the time interval can be given bybewleymointemam2001 ():
where is the fluctuating wall pressure and the upper and lower blowing/suction surfaces, whereas the additional energy required to drive the unactuated flow against the increased viscous drag on the same time interval is given by:
The ratio reads about and , for the oblique waves and streamwise vortex case, respectively. Furthermore, linear tests using actuation have shown that can be of the same order of magnitude of the actual reduction in maximum transient energy growth. Hence, these results indicate that, even if the energy expenditure due to the control action is comparable to that experienced in the linear amplification of the optimal disturbance, it is nevertheless negligible if compared to the potential energy saving due to transition prevention.
Vi Conclusions
The present work has considered the design of fullstate feedback controllers specifically targeting the transient energy growth mechanism in laminar channel flow. It has been shown that full transpiration at both walls and fullstate knowledge are not sufficient to ensure a monotonically stable closedloop system via a linear feedback law. Further, an advanced control design technique – based on a LMI formulation – has been employed to design feedback controllers that have been tested in both the linearized setting and in nonlinear, transitional flows.
Linear tests indicated that the LMI strategy allows to obtain a consistent reduction of the maximum openloop transient growth. At a given global control expenditure over a time window sufficiently long for the perturbations to decay to zero, the LMIcontrolled closedloop system experiences a lower transient energy amplification than a LQRcontrolled closedloop flow. However, in presence of an optimal perturbation for the unactuated flow, the performance of the two control strategies is practically equivalent. Results obtained in the linear setting further indicate that, in the case of perturbations in the form of oblique waves and streamwise vortices, the most effective actuation components are and , respectively.
In the nonlinear case, it has been found that these controllers are capable of increasing the threshold energy for transition when initial conditions are given to the flow in the form of oblique waves or streamwise vortices; the effectiveness of different actuation components has been addressed, indicating that wall blowing/suction is most effective in providing a higher closedloop threshold energy. Additionally, in transitional conditions, LMI controllers prove to be energyeffective, as the energy required by the control action is negligible when compared to the energy saving due to avoiding transition.
References
 [1] K.M. Butler and B.F. Farrel. Threedimensional optimal perturbations in viscous shear flow. Phys. Fluids, 4(8):1637–1650, 1992.
 [2] S.C. Reddy and D.S. Henningson. Energy growth in viscous channel flows. J. Fluid Mech., 252:209–238, 1993.
 [3] P.J. Schmid and D.S. Henningson. Optimal energy density growth in HagenPoiseuille flow. J. Fluid Mech., 277:197–255, 1994.
 [4] P. J. Schmid and D. S. Henningson. Stability and transition in shear flows. Springer, 2000.
 [5] P. J. Schmid. Nonmodal stability theory. Ann. Rev. Fluid Mech, 39:129–162, 2007.
 [6] S. S. Joshi, J. L. Speyer, and J. Kim. A systems theory approach to the feedback stabilization of infinitesimal and finiteamplitude disturbances in plane Poiseuille flow. J. Fluid Mech., 332(157), 1997.
 [7] T.R. Bewley and S. Liu. Optimal and robust control and estimation of linear paths to transition. J. Fluid Mech., 365:305–349, Jun 1998.
 [8] M. Högberg, T.R. Bewley, and D.S. Henningson. Linear feedback control and estimation of transition in plane channel flow. J. Fluid Mech., 481:149–175, Apr 2003.
 [9] J. Kim and J. Lim. A linear process in wallbounded turbulent shear flows. Phys. Fluids, 12(8):1885–1888, 2000.
 [10] K. Lee, L. Cortelezzi, J. Kim, and J. L. Speyer. Application of reducedorder controller to turbulent flows for drag reduction. Phys. Fluids, 13(5):1321–1330, 2001.
 [11] T. R. Bewley, P. Moin, and R. Temam. DNSbased predictive control of turbulence: an optimal benchmark for feedback algorithms. J. Fluid Mech., 447:179–225, 2001.
 [12] J. Kim. Control of turbulent boundary layers. Phys. Fluids, 15(5):1093–1105, 2003.
 [13] T. R. Bewley. Flow control: new challenges for a new Renaissance. Prog. Aero. Sci., 37:21–58, 2001.
 [14] P. Corbett and A. Bottaro. Optimal control of nonmodal disturbances in boundary layers. Theor. Comp. Fluid Dyn., 15:65–81, 2001.
 [15] S. Zuccher, P. Luchini, and A. Bottaro. Algebraic growth in a Blasius boundary layer: optimal and robust control by mean suction in the nonlinear regime. J. Fluid Mech., 513:135–160, 2004.
 [16] H.J. Sussmann and P.V. Kokotovic. The peaking phenomenon and the global stabilization of nonlinear systems. IEEE Trans. Autom. Control, 36(4):424–440, 1991.
 [17] E. Plischke and F. Wirth. Stabilization of linear systems with prescribed transient bounds. In Proc. 16th Int. Symp. Math. Theory Net. Syst., Leuven, Belgium, 2004.
 [18] E. Plischke. Transient effects of linear dynamical systems. PhD thesis, Universitat Bremen, Germany, 2005.
 [19] J.F. Whidborne and J. McKernan. On the minimization of maximum transient energy growth. IEEE Trans. Autom. Control, 52(9):1762–1767, 2007.
 [20] J.F. Whidborne, J. McKernan, and G. Papadakis. Minimising transient energy growth in plane Poiseuille flow. Proceedings Inst. Mech. Eng., Part I, J. Syst. Control Eng., 222(15):323–331, 2008.
 [21] F. Martinelli, M. Quadrio, J. McKernan, and J. F. Whidborne. Feedback control of transient energy growth in subcritical plane Poiseuille flow. In D. S. Henningson and P. Schlatter, editors, Seventh IUTAM Symposium on LaminarTurbulent Transition, volume 18, 2009.
 [22] J. McKernan. Control of Plane Poiseuille Flow: A Theoretical and Computational Investigation. PhD thesis, Department of Aerospace Sciences, School of Engineering, Cranfield University, 2006.
 [23] S. Skogestad and I. Postlethwaite. Multivariable Feedback Control. Wiley, Chichester, England, 1996.
 [24] W. Lohmiller and J. J. Slotine. On contraction analysis for nonlinear systems. Automatica, 34, 1998.
 [25] M. Högberg, T. R. Bewley, and D. S. Henningson. Relaminarization of =100 turbulence using gain scheduling and linear statefeedback control. Phys. Fluids, 15(11):3572–3575, 2003.
 [26] S. Boyd, L. El Ghaoui, E. Feron, and V. Balakrishnan. Linear matrix inequalities in system and control theory. SIAM, Philadelphia, 1994.
 [27] Y. Nesterov and A. Nemirovskii. Interiorpoint polynomial algorithms in convex programming. SIAM, 1994.
 [28] P. Luchini and M. Quadrio. A lowcost parallel implementation of direct numerical simulation of wall turbulence. J. Comp. Phys., 211(2):551–571, 2006.
 [29] S. C. Reddy, P. J. Schmid, J. S. Baggett, and D. S. Henningson. On stability of streamwise streaks and transition thresholds in plane channel flows. J. Fluid Mech., 365:269–303, 1998.