scholarly journals Linear iterative method for closed-loop control of quasiperiodic flows

2019 ◽  
Vol 868 ◽  
pp. 26-65 ◽  
Author(s):  
Colin Leclercq ◽  
Fabrice Demourant ◽  
Charles Poussot-Vassal ◽  
Denis Sipp

This work proposes a feedback-loop strategy to suppress intrinsic oscillations of resonating flows in the fully nonlinear regime. The frequency response of the flow is obtained from the resolvent operator about the mean flow, extending the framework initially introduced by McKeon & Sharma (J. Fluid Mech., vol. 658, 2010, pp. 336–382) to study receptivity mechanisms in turbulent flows. Using this linear time-invariant model of the nonlinear flow, modern control methods such as structured ${\mathcal{H}}_{\infty }$-synthesis can be used to design a controller. The approach is successful in damping self-sustained oscillations associated with specific eigenmodes of the mean-flow spectrum. Despite excellent performance, the linear controller is however unable to completely suppress flow oscillations, and the controlled flow is effectively attracted towards a new dynamical equilibrium. This new attractor is characterized by a different mean flow, which can in turn be used to design a second controller. The method can then be iterated on subsequent mean flows, until the coupled system eventually converges to the base flow. An intuitive parallel can be drawn with Newton’s iteration: at each step, a linearized model of the flow response to a perturbation of the input is sought, and a new linear controller is designed, aiming at further reducing the fluctuations. The method is illustrated on the well-known case of two-dimensional incompressible open-cavity flow at Reynolds number $Re=7500$, where the fully developed flow is initially quasiperiodic (2-torus state). The base flow is reached after five iterations. The present work demonstrates that nonlinear control problems may be solved without resorting to nonlinear reduced-order models. It also shows that physically relevant linear models can be systematically derived for nonlinear flows, without resorting to black-box identification from input–output data; the key ingredient being frequency-domain models based on the linearized Navier–Stokes equations about the mean flow. Applicability to amplifier flows and turbulent dynamics has, however, yet to be investigated.

2002 ◽  
Vol 465 ◽  
pp. 213-235 ◽  
Author(s):  
D. R. GRAHAM ◽  
J. J. L. HIGDON

Oscillatory forcing of a porous medium may have a dramatic effect on the mean flow rate produced by a steady applied pressure gradient. The oscillatory forcing may excite nonlinear inertial effects leading to either enhancement or retardation of the mean flow. Here, in Part 1, we consider the effects of non-zero inertial forces on steady flows in porous media, and investigate the changes in the flow character arising from changes in both the strength of the inertial terms and the geometry of the medium. The steady-state Navier–Stokes equations are solved via a Galerkin finite element method to determine the velocity fields for simple two-dimensional models of porous media. Two geometric models are considered based on constricted channels and periodic arrays of circular cylinders. For both geometries, we observe solution multiplicity yielding both symmetric and asymmetric flow patterns. For the cylinder arrays, we demonstrate that inertial effects lead to anisotropy in the effective permeability, with the direction of minimum resistance dependent on the solid volume fraction. We identify nonlinear flow phenomena which might be exploited by oscillatory forcing to yield a net increase in the mean flow rate. In Part 2, we take up the subject of unsteady flows governed by the full time-dependent Navier–Stokes equations.


2010 ◽  
Vol 660 ◽  
pp. 37-54 ◽  
Author(s):  
OLAF MARXEN ◽  
ULRICH RIST

The mutual interaction of laminar–turbulent transition and mean flow evolution is studied in a pressure-induced laminar separation bubble on a flat plate. The flat-plate boundary layer is subjected to a sufficiently strong adverse pressure gradient that a separation bubble develops. Upstream of the bubble a small-amplitude disturbance is introduced which causes transition. Downstream of transition, the mean flow strongly changes and, due to viscous–inviscid interaction, the overall pressure distribution is changed as well. As a consequence, the mean flow also changes upstream of the transition location. The difference in the mean flow between the forced and the unforced flows is denoted the mean flow deformation. Two different effects are caused by the mean flow deformation in the upstream, laminar part: a reduction of the size of the separation region and a stabilization of the flow with respect to small, linear perturbations. By carrying out numerical simulations based on the original base flow and the time-averaged deformed base flow, we are able to distinguish between direct and indirect nonlinear effects. Direct effects are caused by the quadratic nonlinearity of the Navier–Stokes equations, are associated with the generation of higher harmonics and are predominantly local. In contrast, the stabilization of the flow is an indirect effect, because it is independent of the Reynolds stress terms in the laminar region and is solely governed by the non-local alteration of the mean flow via the pressure.


2017 ◽  
Vol 139 (12) ◽  
Author(s):  
Yuanchao Li ◽  
Huang Chen ◽  
Joseph Katz

Modeling of turbulent flows in axial turbomachines is challenging due to the high spatial and temporal variability in the distribution of the strain rate components, especially in the tip region of rotor blades. High-resolution stereo-particle image velocimetry (SPIV) measurements performed in a refractive index-matched facility in a series of closely spaced planes provide a comprehensive database for determining all the terms in the Reynolds stress and strain rate tensors. Results are also used for calculating the turbulent kinetic energy (TKE) production rate and transport terms by mean flow and turbulence. They elucidate some but not all of the observed phenomena, such as the high anisotropy, high turbulence levels in the vicinity of the tip leakage vortex (TLV) center, and in the shear layer connecting it to the blade suction side (SS) tip corner. The applicability of popular Reynolds stress models based on eddy viscosity is also evaluated by calculating it from the ratio between stress and strain rate components. Results vary substantially, depending on which components are involved, ranging from very large positive to negative values. In some areas, e.g., in the tip gap and around the TLV, the local stresses and strain rates do not appear to be correlated at all. In terms of effect on the mean flow, for most of the tip region, the mean advection terms are much higher than the Reynolds stress spatial gradients, i.e., the flow dynamics is dominated by pressure-driven transport. However, they are of similar magnitude in the shear layer, where modeling would be particularly challenging.


2011 ◽  
Vol 68 (9) ◽  
pp. 2042-2060 ◽  
Author(s):  
David A. Ortland ◽  
M. Joan Alexander ◽  
Alison W. Grimsdell

Abstract Convective heating profiles are computed from one month of rainfall rate and cloud-top height measurements using global Tropical Rainfall Measuring Mission and infrared cloud-top products. Estimates of the tropical wave response to this heating and the mean flow forcing by the waves are calculated using linear and nonlinear models. With a spectral resolution up to zonal wavenumber 80 and frequency up to 4 cpd, the model produces 50%–70% of the zonal wind acceleration required to drive a quasi-biennial oscillation (QBO). The sensitivity of the wave spectrum to the assumed shape of the heating profile, to the mean wind and temperature structure of the tropical troposphere, and to the type of model used is also examined. The redness of the heating spectrum implies that the heating strongly projects onto Hough modes with small equivalent depth. Nonlinear models produce wave flux significantly smaller than linear models due to what appear to be dynamical processes that limit the wave amplitude. Both nonlinearity and mean winds in the lower stratosphere are effective in reducing the Rossby wave response to heating relative to the response in a linear model for a mean state at rest.


2018 ◽  
Vol 859 ◽  
pp. 516-542 ◽  
Author(s):  
Calum S. Skene ◽  
Peter J. Schmid

A linear numerical study is conducted to quantify the effect of swirl on the response behaviour of premixed lean flames to general harmonic excitation in the inlet, upstream of combustion. This study considers axisymmetric M-flames and is based on the linearised compressible Navier–Stokes equations augmented by a simple one-step irreversible chemical reaction. Optimal frequency response gains for both axisymmetric and non-axisymmetric perturbations are computed via a direct–adjoint methodology and singular value decompositions. The high-dimensional parameter space, containing perturbation and base-flow parameters, is explored by taking advantage of generic sensitivity information gained from the adjoint solutions. This information is then tailored to specific parametric sensitivities by first-order perturbation expansions of the singular triplets about the respective parameters. Valuable flow information, at a negligible computational cost, is gained by simple weighted scalar products between direct and adjoint solutions. We find that for non-swirling flows, a mode with azimuthal wavenumber $m=2$ is the most efficiently driven structure. The structural mechanism underlying the optimal gains is shown to be the Orr mechanism for $m=0$ and a blend of Orr and other mechanisms, such as lift-up, for other azimuthal wavenumbers. Further to this, velocity and pressure perturbations are shown to make up the optimal input and output showing that the thermoacoustic mechanism is crucial in large energy amplifications. For $m=0$ these velocity perturbations are mainly longitudinal, but for higher wavenumbers azimuthal velocity fluctuations become prominent, especially in the non-swirling case. Sensitivity analyses are carried out with respect to the Mach number, Reynolds number and swirl number, and the accuracy of parametric gradients of the frequency response curve is assessed. The sensitivity analysis reveals that increases in Reynolds and Mach numbers yield higher gains, through a decrease in temperature diffusion. A rise in mean-flow swirl is shown to diminish the gain, with increased damping for higher azimuthal wavenumbers. This leads to a reordering of the most effectively amplified mode, with the axisymmetric ($m=0$) mode becoming the dominant structure at moderate swirl numbers.


2018 ◽  
Vol 84 (4) ◽  
Author(s):  
A. Alexakis ◽  
S. Fauve ◽  
C. Gissinger ◽  
F. Pétrélis

We discuss the effect of different types of fluctuations on dynamos generated in the limit of scale separation. We first recall that the magnetic field observed in the VKS (von Karman flow of liquid sodium) experiment is not the one that would be generated by the mean flow alone and that smaller scale turbulent fluctuations therefore play an important role. We then consider how velocity fluctuations affect the dynamo threshold in the framework of mean-field magnetohydrodynamics. We show that the detrimental effect of turbulent fluctuations observed with many flows disappears in the case of helical flows with scale separation. We also find that fluctuations of the electrical conductivity of the fluid, for instance related to temperature fluctuations in convective flows, provide an efficient mechanism for dynamo action. Finally, we conclude by describing an experimental configuration that could be used to test the validity of mean-field magnetohydrodynamics in strongly turbulent flows.


Author(s):  
Roland Gårdhagen ◽  
Jonas Lantz ◽  
Fredrik Carlsson ◽  
Matts Karlsson

Low and/or oscillatory Wall Shear Stress (WSS) has been correlated with elevated risk for increased intima media thickness and atherosclerosis in several studies during the last decades [1, 2]. Most of the studies have addressed laminar flows, in which the oscillations mainly are due to the pulsating nature of blood flow. Turbulent flows however show significant spatial and temporal fluctuations although the mean flow is steady.


2019 ◽  
Vol 104 (2-3) ◽  
pp. 533-552
Author(s):  
J. Casacuberta ◽  
K. J. Groot ◽  
Q. Ye ◽  
S. Hickel

AbstractMicro-ramps are popular passive flow control devices which can delay flow separation by re-energising the lower portion of the boundary layer. We compute the laminar base flow, the instantaneous transitional flow, and the mean flow around a micro-ramp immersed in a quasi-incompressible boundary layer at supercritical roughness Reynolds number. Results of our Direct Numerical Simulations (DNS) are compared with results of BiLocal stability analysis on the DNS base flow and independent tomographic Particle Image Velocimetry (tomo-PIV) experiments. We analyse relevant flow structures developing in the micro-ramp wake and assess their role in the micro-ramp functionality, i.e., in increasing the near-wall momentum. The main flow feature of the base flow is a pair of streamwise counter-rotating vortices induced by the micro-ramp, the so-called primary vortex pair. In the instantaneous transitional flow, the primary vortex pair breaks up into large-scale hairpin vortices, which arise due to linear varicose instability of the base flow, and unsteady secondary vortices develop. Instantaneous vortical structures obtained by DNS and experiments are in good agreement. Matching linear disturbance growth rates from DNS and linear stability analysis are obtained until eight micro-ramp heights downstream of the micro-ramp. For the setup considered in this article, we show that the working principle of the micro-ramp is different from that of classical vortex generators; we find that transitional perturbations are more efficient in increasing the near-wall momentum in the mean flow than the laminar primary vortices in the base flow.


1975 ◽  
Vol 42 (1) ◽  
pp. 51-54 ◽  
Author(s):  
N. W. Wilson ◽  
R. S. Azad

A single set of equations is developed to predict the mean flow characteristics in long circular pipes operating at laminar, transitional, and turbulent Reynolds numbers. Generally good agreement is obtained with available data in the Reynolds number range 100 < Re < 500,000.


Author(s):  
S M Fraser ◽  
Y Zhang

Three-dimensional turbulent flow through the impeller passage of a model mixed-flow pump has been simulated by solving the Navier-Stokes equations with an improved κ-ɛ model. The standard κ-ɛ model was found to be unsatisfactory for solving the off-design impeller flow and a converged solution could not be obtained at 49 per cent design flowrate. After careful analysis, it was decided to modify the standard κ-ɛ model by including the extra rates of strain due to the acceleration of impeller rotation and geometrical curvature and removing the mathematical ill-posedness between the mean flow turbulence modelling and the logarithmic wall function.


Sign in / Sign up

Export Citation Format

Share Document