skip to main content

Explore Research Products in the NSF-PAR

Title: Feedback control of transitional flows: A framework for controller verification using quadratic constraints
The dynamics of incompressible fluid flows are governed by a non-normal linear dynamical system in feedback with a static energy-conserving nonlinearity. These dynamics can be altered using feedback control but verifying performance of a given control law can be challenging. The conventional approach is to perform a campaign of high-fidelity direct numerical simulations to assess performance over a wide range of parameters and disturbance scenarios. In this paper, we propose an alternative simulation-free approach for controller verification. The incompressible Navier-Stokes equations are modeled as a linear system in feedback with a static and quadratic nonlinearity. The energy conserving property of this nonlinearity can be expressed as a set of quadratic constraints on the system, which allows us to perform a nonlinear stability analysis of the fluid dynamics with minimal complexity. In addition, the Reynolds number variations only influence the linear dynamics in the Navier-Stokes equations. Therefore, the fluid flow can be modeled as a parameter-varying linear system (with Reynolds number as the parameter) in feedback with a quadratic nonlinearity. The quadratic constraint framework is used to determine the range of Reynolds numbers over which a given flow will be stable, without resorting to numerical simulations. We demonstrate the framework on a reduced-order model of plane Couette flow. We show that our proposed method allows us to determine the critical Reynolds number, largest initial disturbance, and a range of parameter variations over which a given controller will stabilize the nonlinear dynamics.  more » « less
Award ID(s):
1943988
NSF-PAR ID:
10309420
Author(s) / Creator(s):
 ;  ;  
Date Published:
Journal Name:
AIAA AVIATION 2021 FORUM
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ( , Fluids)
    null (Ed.)
    A closed-loop control algorithm for the reduction of turbulent flow separation over NACA 0015 airfoil equipped with leading-edge synthetic jet actuators (SJAs) is presented. A system identification approach based on Nonlinear Auto-Regressive Moving Average with eXogenous inputs (NARMAX) technique was used to predict nonlinear dynamics of the fluid flow and for the design of the controller system. Numerical simulations based on URANS equations are performed at Reynolds number of 106 for various airfoil incidences with and without closed-loop control. The NARMAX model for flow over an airfoil is based on the static pressure data, and the synthetic jet actuator is developed using an incompressible flow model. The corresponding NARMAX identification model developed for the pressure data is nonlinear; therefore, the describing function technique is used to linearize the system within its frequency range. Low-pass filtering is used to obtain quasi-linear state values, which assist in the application of linear control techniques. The reference signal signifies the condition of a fully re-attached flow, and it is determined based on the linearization of the original signal during open-loop control. The controller design follows the standard proportional-integral (PI) technique for the single-input single-output system. The resulting closed-loop response tracks the reference value and leads to significant improvements in the transient response over the open-loop system. The NARMAX controller enhances the lift coefficient from 0.787 for the uncontrolled case to 1.315 for the controlled case with an increase of 67.1%. 
    more » « less
  2. ; ( , AIAA SciTech)
    null (Ed.)
    In this paper, we started by summarizing our recently developed viscous unsteady theory based on coupling potential flow with the triple deck boundary layer theory. This approach provides a viscous extension of potential flow unsteady aerodynamics. As such, a Reynolds- number-dependence could be determined. We then developed a finite-state approximation of such a theory, presenting it in a state space model. This novel nonlinear state space model of the viscous unsteady aerodynamic loads is expected to serve aerodynamicists better than the classical Theodorsen’s model, as it captures viscous effects (i.e., Reynolds number dependence) as well as nonlinearity and additional lag in the lift dynamics; and allows simulation of arbitrary time-varying airfoil motions (not necessarily harmonic). Moreover, being in a state space form makes it quite convenient for simulation and coupling with structural dynamics to perform aeroelasticity, flight dynamics analysis, and control design. We then proceeded to develop a linearization of such a model, which enables analytical results. So, we derived an analytical representation of the viscous lift frequency response function, which is an explicit function of, not only frequency, but also Reynolds number. We also developed a state space model of the linearized response. We finally simulated the nonlinear and linear models to a non- harmonic, small-amplitude pitching maneuver at 100 , 000 Reynolds number and compared the resulting lift and pitching moment with potential flow, in reference to relatively higher fidelity computations of the Unsteady Reynolds-Averaged Navier-Stokes equations. 
    more » « less
  3. ; ; ; ( , International journal for numerical methods in biomedical engineering)
    Computational fluid dynamics (CFD) is increasingly used to study blood flows in patient-specific arteries for understanding certain cardiovascular diseases. The techniques work quite well for relatively simple problems but need improvements when the problems become harder when (a) the geometry becomes complex (eg, a few branches to a full pulmonary artery), (b) the model becomes more complex (eg, fluid-only to coupled fluid-structure interaction), (c) both the fluid and wall models become highly nonlinear, and (d) the computer on which we run the simulation is a supercomputer with tens of thousands of processor cores. To push the limit of CFD in all four fronts, in this paper, we develop and study a highly parallel algorithm for solving a monolithically coupled fluid-structure system for the modeling of the interaction of the blood flow and the arterial wall. As a case study, we consider a patient-specific, full size pulmonary artery obtained from computed tomography (CT) images, with an artificially added layer of wall with a fixed thickness. The fluid is modeled with a system of incompressible Navier-Stokes equations, and the wall is modeled by a geometrically nonlinear elasticity equation. As far as we know, this is the first time the unsteady blood flow in a full pulmonary artery is simulated without assuming a rigid wall. The proposed numerical algorithm and software scale well beyond 10 000 processor cores on a supercomputer for solving the fluid-structure interaction problem discretized with a stabilized finite element method in space and an implicit scheme in time involving hundreds of millions of unknowns. 
    more » « less
  4. ; ; ; ( , International Journal for Numerical Methods in Biomedical Engineering)
    Abstract

    Computational fluid dynamics (CFD) is increasingly used to study blood flows in patient‐specific arteries for understanding certain cardiovascular diseases. The techniques work quite well for relatively simple problems but need improvements when the problems become harder when (a) the geometry becomes complex (eg, a few branches to a full pulmonary artery), (b) the model becomes more complex (eg, fluid‐only to coupled fluid‐structure interaction), (c) both the fluid and wall models become highly nonlinear, and (d) the computer on which we run the simulation is a supercomputer with tens of thousands of processor cores. To push the limit of CFD in all four fronts, in this paper, we develop and study a highly parallel algorithm for solving a monolithically coupled fluid‐structure system for the modeling of the interaction of the blood flow and the arterial wall. As a case study, we consider a patient‐specific, full size pulmonary artery obtained from computed tomography (CT) images, with an artificially added layer of wall with a fixed thickness. The fluid is modeled with a system of incompressible Navier‐Stokes equations, and the wall is modeled by a geometrically nonlinear elasticity equation. As far as we know, this is the first time the unsteady blood flow in a full pulmonary artery is simulated without assuming a rigid wall. The proposed numerical algorithm and software scale well beyond 10 000 processor cores on a supercomputer for solving the fluid‐structure interaction problem discretized with a stabilized finite element method in space and an implicit scheme in time involving hundreds of millions of unknowns.

     
    more » « less
  5. ; ; ; ( , Physics of Fluids)
    Human phonation involves the flow-induced vibrations of the vocal folds (VFs) that result from the interaction with airflow through the larynx. Most voice dysfunctions correspond with the fluid–structure interaction (FSI) features as well as the local changes in perfusion within the VF tissue. This study aims to develop a multiphysics computational framework to simulate the interstitial fluid flow dynamics in vibrating VFs using a biphasic description of the tissue and FSI methodology. The integration of FSI and a permeable VF model presents a novel approach to capture phonation physics' complexity and investigate VF tissue's porous nature. The glottal airflow is modeled by the unsteady, incompressible Navier–Stokes equations, and the Brinkman equation is employed to simulate the flow through the saturated porous medium of the VFs. The computational model provides a prediction of tissue deformation metrics and pulsatile glottal flow, in addition to the interstitial fluid velocity and flow circulation within the porous structure. Furthermore, the model is used to characterize the effects of variation in subglottal lung pressure and VF permeability coefficient by conducting parametric studies. Subsequent investigations to quantify the relationships between these input variables, flow perfusion, pore pressure, and vibration amplitude are presented. A linear relationship is found between the vibration amplitude, pore pressure, and filtration flow with subglottal pressure, whereas a nonlinear dependence between the filtration velocity and VF permeability coefficient is detected. The outcomes highlight the importance of poroelasticity in phonation models. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email