212 research outputs found
Efficient Solution of Large-Scale Algebraic Riccati Equations Associated with Index-2 DAEs via the Inexact Low-Rank Newton-ADI Method
This paper extends the algorithm of Benner, Heinkenschloss, Saak, and
Weichelt: An inexact low-rank Newton-ADI method for large-scale algebraic
Riccati equations, Applied Numerical Mathematics Vol.~108 (2016), pp.~125--142,
doi:10.1016/j.apnum.2016.05.006 to Riccati equations associated with Hessenberg
index-2 Differential Algebratic Equation (DAE) systems. Such DAE systems arise,
e.g., from semi-discretized, linearized (around steady state) Navier-Stokes
equations. The solution of the associated Riccati equation is important, e.g.,
to compute feedback laws that stabilize the Navier-Stokes equations. Challenges
in the numerical solution of the Riccati equation arise from the large-scale of
the underlying systems and the algebraic constraint in the DAE system. These
challenges are met by a careful extension of the inexact low-rank Newton-ADI
method to the case of DAE systems. A main ingredient in the extension to the
DAE case is the projection onto the manifold described by the algebraic
constraints. In the algorithm, the equations are never explicitly projected,
but the projection is only applied as needed. Numerical experience indicates
that the algorithmic choices for the control of inexactness and line-search can
help avoid subproblems with matrices that are only marginally stable. The
performance of the algorithm is illustrated on a large-scale Riccati equation
associated with the stabilization of Navier-Stokes flow around a cylinder.Comment: 21 pages, 2 figures, 4 table
Relaminarisation of Re_{\tau} = 100 channel flow with globally stabilising linear feedback control
The problems of nonlinearity and high dimension have so far prevented a
complete solution of the control of turbulent flow. Addressing the problem of
nonlinearity, we propose a flow control strategy which ensures that the energy
of any perturbation to the target profile decays monotonically. The
controller's estimate of the flow state is similarly guaranteed to converge to
the true value. We present a one-time off-line synthesis procedure, which
generalises to accommodate more restrictive actuation and sensing arrangements,
with conditions for existence for the controller given in this case. The
control is tested in turbulent channel flow () using full-domain
sensing and actuation on the wall-normal velocity. Concentrated at the point of
maximum inflection in the mean profile, the control directly counters the
supply of turbulence energy arising from the interaction of the wall-normal
perturbations with the flow shear. It is found that the control is only
required for the larger-scale motions, specifically those above the scale of
the mean streak spacing. Minimal control effort is required once laminar flow
is achieved. The response of the near-wall flow is examined in detail, with
particular emphasis on the pressure and wall-normal velocity fields, in the
context of Landahl's theory of sheared turbulence
A Snapshot Algorithm for Linear Feedback Flow Control Design
The control of fluid flows has many applications. For micro air vehicles, integrated flow control designs could enhance flight stability by mitigating the effect of destabilizing air flows in their low Reynolds number regimes. However, computing model based feedback control designs can be challenging due to high dimensional discretized flow models. In this work, we investigate the use of a snapshot algorithm proposed in Ref. 1 to approximate the feedback gain operator for a linear incompressible unsteady flow problem on a bounded domain. The main component of the algorithm is obtaining solution snapshots of certain linear flow problems. Numerical results for the example flow problem show convergence of the feedback gains
- …