5,225 research outputs found
Koopman operator-based model reduction for switched-system control of PDEs
We present a new framework for optimal and feedback control of PDEs using
Koopman operator-based reduced order models (K-ROMs). The Koopman operator is a
linear but infinite-dimensional operator which describes the dynamics of
observables. A numerical approximation of the Koopman operator therefore yields
a linear system for the observation of an autonomous dynamical system. In our
approach, by introducing a finite number of constant controls, the dynamic
control system is transformed into a set of autonomous systems and the
corresponding optimal control problem into a switching time optimization
problem. This allows us to replace each of these systems by a K-ROM which can
be solved orders of magnitude faster. By this approach, a nonlinear
infinite-dimensional control problem is transformed into a low-dimensional
linear problem. In situations where the Koopman operator can be computed
exactly using Extended Dynamic Mode Decomposition (EDMD), the proposed approach
yields optimal control inputs. Furthermore, a recent convergence result for
EDMD suggests that the approach can be applied to more complex dynamics as
well. To illustrate the results, we consider the 1D Burgers equation and the 2D
Navier--Stokes equations. The numerical experiments show remarkable performance
concerning both solution times and accuracy.Comment: arXiv admin note: text overlap with arXiv:1801.0641
Reduced-order modeling of transonic flows around an airfoil submitted to small deformations
A reduced-order model (ROM) is developed for the prediction of unsteady transonic flows past an airfoil submitted to small deformations, at moderate Reynolds number. Considering a suitable state formulation as well as a consistent inner product, the Galerkin projection of the compressible flow Navier–Stokes equations, the high-fidelity (HF) model, onto a low-dimensional basis determined by Proper Orthogonal Decomposition (POD), leads to a polynomial quadratic ODE system relevant to the prediction of main flow features. A fictitious domain deformation technique is yielded by the Hadamard formulation of HF model and validated at HF level. This approach captures airfoil profile deformation by a modification of the boundary conditions whereas the spatial domain remains unchanged. A mixed POD gathering information from snapshot series associated with several airfoil profiles can be defined. The temporal coefficients in POD expansion are shape-dependent while spatial POD modes are not. In the ROM, airfoil deformation is introduced by a steady forcing term. ROM reliability towards airfoil deformation is demonstrated for the prediction of HF-resolved as well as unknown intermediate configurations
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
The LifeV library: engineering mathematics beyond the proof of concept
LifeV is a library for the finite element (FE) solution of partial
differential equations in one, two, and three dimensions. It is written in C++
and designed to run on diverse parallel architectures, including cloud and high
performance computing facilities. In spite of its academic research nature,
meaning a library for the development and testing of new methods, one
distinguishing feature of LifeV is its use on real world problems and it is
intended to provide a tool for many engineering applications. It has been
actually used in computational hemodynamics, including cardiac mechanics and
fluid-structure interaction problems, in porous media, ice sheets dynamics for
both forward and inverse problems. In this paper we give a short overview of
the features of LifeV and its coding paradigms on simple problems. The main
focus is on the parallel environment which is mainly driven by domain
decomposition methods and based on external libraries such as MPI, the Trilinos
project, HDF5 and ParMetis.
Dedicated to the memory of Fausto Saleri.Comment: Review of the LifeV Finite Element librar
Feedback control of unstable steady states of flow past a flat plate using reduced-order estimators
We present an estimator-based control design procedure for flow control,
using reduced-order models of the governing equations, linearized about a
possibly unstable steady state. The reduced models are obtained using an
approximate balanced truncation method that retains the most controllable and
observable modes of the system. The original method is valid only for stable
linear systems, and we present an extension to unstable linear systems. The
dynamics on the unstable subspace are represented by projecting the original
equations onto the global unstable eigenmodes, assumed to be small in number. A
snapshot-based algorithm is developed, using approximate balanced truncation,
for obtaining a reduced-order model of the dynamics on the stable subspace. The
proposed algorithm is used to study feedback control of 2-D flow over a flat
plate at a low Reynolds number and at large angles of attack, where the natural
flow is vortex shedding, though there also exists an unstable steady state. For
control design, we derive reduced-order models valid in the neighborhood of
this unstable steady state. The actuation is modeled as a localized body force
near the leading edge of the flat plate, and the sensors are two velocity
measurements in the near-wake of the plate. A reduced-order Kalman filter is
developed based on these models and is shown to accurately reconstruct the flow
field from the sensor measurements, and the resulting estimator-based control
is shown to stabilize the unstable steady state. For small perturbations of the
steady state, the model accurately predicts the response of the full
simulation. Furthermore, the resulting controller is even able to suppress the
stable periodic vortex shedding, where the nonlinear effects are strong, thus
implying a large domain of attraction of the stabilized steady state.Comment: 36 pages, 17 figure
- …