Abstract:In this paper we develop the Optimal Control Approach to the rotary control of the cylinder wake. We minimize the functional which represents the sum of the work needed to resist the drag force and the work needed to control the flow, where the rotation rate φ̇(t) is the control variable. Sensitivity of the functional to control is determined using the adjoint equations. We solve them in the “vorticity” form, which is a novel approach and leads to computational advantages. Simulations performed at Re=75 and Re… Show more
“…To supplement this main result, we will also give an estimate of the computational savings that can be obtained by a POD ROM based optimal control approach compared with the more "classical" approach where the Navier-Stokes equations are used for constraints [38][39][40]. Consequently in this study our main concern is not to determine the control law with the maximum energetic efficiency as it can be characterized for example by the Power Saving Ratio (PSR) [40, for a definition or hereafter in § 5.2.3].…”
Section: A Generic Configuration Of Separated Flows: the Cylinder Wakmentioning
confidence: 99%
“…Recently, due to the maturity of control theory, optimization methods and computational fluid dynamics, optimal and suboptimal approaches attracted increased attention in flow control setting [35][36][37]. For example, in [38][39][40] the optimal control theory was used with the twodimensional Navier-Stokes equations as the state equation to control by rotary oscillation the unsteady wake of the cylinder (see table 1 for the characteristics of these approaches). An attractive element of the optimal control approach is that the control design is explicitly based on the cost functional.…”
Section: A Generic Configuration Of Separated Flows: the Cylinder Wakmentioning
In this paper, optimal control theory is used to minimize the total mean drag for a circular cylinder wake flow in the laminar regime (Re = 200). The control parameters are the amplitude and the frequency of the time-harmonic cylinder rotation. In order to reduce the size of the discretized optimality system, a Proper Orthogonal Decomposition (POD) Reduced-Order Model (ROM) is derived to be used as state equation. We then propose to employ the Trust-Region Proper Orthogonal Decomposition (TRPOD) approach, originally introduced by Fahl (2000), to update the reduced-order models during the optimization process. A lot of computational work is saved because the optimization process is now based only on low-fidelity models. A particular care was taken to derive a POD ROM for the pressure and velocity fields with an appropriate balance between model accuracy and robustness. The key enablers are the extension of the POD basis functions to the pressure data, the use of calibration methods for the POD ROM and the addition in the POD expansion of several non-equilibrium modes to describe various operating conditions. When the TRPOD algorithm is applied to the wake flow configuration, this approach converges to the minimum predicted by an open-loop control approach and leads to a relative mean drag reduction of 30% at reduced cost.
“…To supplement this main result, we will also give an estimate of the computational savings that can be obtained by a POD ROM based optimal control approach compared with the more "classical" approach where the Navier-Stokes equations are used for constraints [38][39][40]. Consequently in this study our main concern is not to determine the control law with the maximum energetic efficiency as it can be characterized for example by the Power Saving Ratio (PSR) [40, for a definition or hereafter in § 5.2.3].…”
Section: A Generic Configuration Of Separated Flows: the Cylinder Wakmentioning
confidence: 99%
“…Recently, due to the maturity of control theory, optimization methods and computational fluid dynamics, optimal and suboptimal approaches attracted increased attention in flow control setting [35][36][37]. For example, in [38][39][40] the optimal control theory was used with the twodimensional Navier-Stokes equations as the state equation to control by rotary oscillation the unsteady wake of the cylinder (see table 1 for the characteristics of these approaches). An attractive element of the optimal control approach is that the control design is explicitly based on the cost functional.…”
Section: A Generic Configuration Of Separated Flows: the Cylinder Wakmentioning
In this paper, optimal control theory is used to minimize the total mean drag for a circular cylinder wake flow in the laminar regime (Re = 200). The control parameters are the amplitude and the frequency of the time-harmonic cylinder rotation. In order to reduce the size of the discretized optimality system, a Proper Orthogonal Decomposition (POD) Reduced-Order Model (ROM) is derived to be used as state equation. We then propose to employ the Trust-Region Proper Orthogonal Decomposition (TRPOD) approach, originally introduced by Fahl (2000), to update the reduced-order models during the optimization process. A lot of computational work is saved because the optimization process is now based only on low-fidelity models. A particular care was taken to derive a POD ROM for the pressure and velocity fields with an appropriate balance between model accuracy and robustness. The key enablers are the extension of the POD basis functions to the pressure data, the use of calibration methods for the POD ROM and the addition in the POD expansion of several non-equilibrium modes to describe various operating conditions. When the TRPOD algorithm is applied to the wake flow configuration, this approach converges to the minimum predicted by an open-loop control approach and leads to a relative mean drag reduction of 30% at reduced cost.
“…Since Riesz representation (5) does not apply in nonHilbert spaces, we employ a more general procedure to extract the steepest descent directions in Banach spaces which follows the proposal first made by Neuberger in [15]. This procedure will involve a nonlinear transformation of the adjoint field equivalent to a nonlinear change of variables in iterative procedure (4). Further-more, by extracting the steepest descent directions in a continuous family of nested spaces we will allow this change of the metric to vary in the course of iterations (4).…”
Section: Introductionmentioning
confidence: 99%
“…• shape optimization in aerodynamics (see, e.g., [1,2]), • flow control for drag reduction, (see, e.g., [3,4]), • variational data assimilation in dynamic meteorology known as 4DVAR (see, e.g., [5]), • mixing enhancement (see, e.g., [6]). …”
In this work we investigate a technique for accelerating convergence of adjoint-based optimization of PDE systems based on a nonlinear change of variables in the control space. This change of variables is accomplished in the differentiate-then-discretize approach by constructing the descent directions in a control space not equipped with the Hilbert structure. We show how such descent directions can be computed in general Lebesgue and Besov spaces, and argue that in the Besov space case determination of descent directions can be interpreted as nonlinear wavelet filtering of the adjoint field. The freedom involved in choosing parameters characterizing the spaces in which the steepest descent directions are constructed can be leveraged to accelerate convergence of iterations. Our computational examples involving state estimation problems for the 1D Kuramoto-Sivashinsky and 3D Navier-Stokes equations indeed show significantly improved performance of the proposed method as compared to the standard approaches.
“…• shape optimization with application to aircraft design, e.g., Mohammadi and Pironneau (2001); Martins et al (2004), • flow control for drag reduction, e.g., Bewley et al (2001); Protas and Styczek (2002), • variational data assimilation in dynamic meteorology, e.g., Kalnay (2003),…”
This note discusses certain aspects of computational solution of optimal control problems for fluid systems. We focus on approaches in which the steepest descent direction of the cost functional is determined using the adjoint equations. In the first part we review the classical formulation by presenting it in the context of Nonlinear Programming. In the second part we show some new results concerning determination of descent directions in general Banach spaces without Hilbert structure. The proposed approach is illustrated with computational examples concerning a state estimation problem for the 1D Kuramoto-Sivashinsky equation.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.