scholarly journals Adjoint-based parametric sensitivity analysis for swirling M-flames

2018 ◽  
Vol 859 ◽  
pp. 516-542 ◽  
Author(s):  
Calum S. Skene ◽  
Peter J. Schmid
Keyword(s):  
Sensitivity Analysis ◽  
Frequency Response ◽  
Stokes Equations ◽  
Numerical Study ◽  
Computational Cost ◽  
Base Flow ◽  
Sensitivity Analyses ◽  
Perturbation Expansions ◽  
Mean Flow ◽  
Adjoint Solutions

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.

scholarly journals Reinforcement-learning-based control of confined cylinder wakes with stability analyses

2021 ◽  
Vol 932 ◽  
Author(s):  
Jichao Li ◽  
Mengqi Zhang
Keyword(s):  
Reinforcement Learning ◽  
Vortex Shedding ◽  
Stokes Equations ◽  
Synthetic Jets ◽  
Base Flow ◽  
Sensitivity Analyses ◽  
Unstable State ◽  
Mean Flow ◽  
Sensitive Region ◽  
Reward Function

This work studies the application of a reinforcement learning (RL)-based flow control strategy to the flow past a cylinder confined between two walls to suppress vortex shedding. The control action is blowing and suction of two synthetic jets on the cylinder. The theme of this study is to investigate how to use and embed physical information of the flow in the RL-based control. First, global linear stability and sensitivity analyses based on the time-mean flow and the steady flow (which is a solution to the Navier–Stokes equations) are conducted in a range of blockage ratios and Reynolds numbers. It is found that the most sensitive region in the wake extends itself when either parameter increases in the parameter range we investigated here. Then, we use these physical results to help design RL-based control policies. We find that the controlled wake converges to the unstable steady base flow, where the vortex shedding can be successfully suppressed. A persistent oscillating control seems necessary to maintain this unstable state. The RL algorithm is able to outperform a gradient-based optimisation method (optimised in a certain period of time) in the long run. Furthermore, when the flow stability information is embedded in the reward function to penalise the instability, the controlled flow may become more stable. Finally, according to the sensitivity analyses, the control is most efficient when the probes are placed in the most sensitive region. The control can be successful even when few probes are properly placed in this manner.

scholarly journals RIP CURRENTS ON A BARRED BEACH

2012 ◽  
Vol 1 (33) ◽  
pp. 38
Author(s):  
Andrea Ruju ◽  
Pablo Higuera ◽  
Javier L. Lara ◽  
Inigo J. Losada ◽  
Giovanni Coco
Keyword(s):  
Boundary Conditions ◽  
Stokes Equations ◽  
Numerical Study ◽  
Three Dimensional ◽  
Wave Generation ◽  
Dimensional Structure ◽  
Navier Stokes ◽  
Rip Currents ◽  
Mean Flow ◽  
Forcing Mechanisms

This work presents the numerical study of rip current circulation on a barred beach. The numerical simulations have been carried out with the IH-FOAM model which is based on the three dimensional Reynolds Averaged Navier-Stokes equations. The new boundary conditions implemented in IH-FOAM have been used, including three dimensional wave generation as well as active wave absorption at the boundary. Applying the specific wave generation boundary conditions, the model is validated to simulate rip circulation on a barred beach. Moreover, this study addresses the identification of the forcing mechanisms and the three dimensional structure of the mean flow.

A Numerical Study of Vortex Breakdown in Turbulent Swirling Flows

1999 ◽  
Vol 122 (1) ◽  
pp. 179-183 ◽  
Author(s):  
Robert E. Spall ◽  
Blake M. Ashby
Keyword(s):  
Reynolds Stress ◽  
Stokes Equations ◽  
Vortex Breakdown ◽  
Numerical Study ◽  
Reynolds Stress Model ◽  
Navier Stokes ◽  
Stress Model ◽  
Mean Flow ◽  
Navier Stokes Equations ◽  
The Mean

Solutions to the incompressible Reynolds-averaged Navier–Stokes equations have been obtained for turbulent vortex breakdown within a slightly diverging tube. Inlet boundary conditions were derived from available experimental data for the mean flow and turbulence kinetic energy. The performance of both two-equation and full differential Reynolds stress models was evaluated. Axisymmetric results revealed that the initiation of vortex breakdown was reasonably well predicted by the differential Reynolds stress model. However, the standard K-ε model failed to predict the occurrence of breakdown. The differential Reynolds stress model also predicted satisfactorily the mean azimuthal and axial velocity profiles downstream of the breakdown, whereas results using the K-ε model were unsatisfactory. [S0098-2202(00)01601-1]

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
Keyword(s):  
Turbulent Flows ◽  
Linear Models ◽  
Stokes Equations ◽  
Linear Time ◽  
Base Flow ◽  
Nonlinear Flow ◽  
Mean Flow ◽  
Loop Control ◽  
Linear Controller ◽  
The Mean

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.

A self-consistent formulation for the sensitivity analysis of finite-amplitude vortex shedding in the cylinder wake

2016 ◽  
Vol 800 ◽  
pp. 327-357 ◽  
Author(s):  
P. Meliga ◽  
E. Boujo ◽  
F. Gallaire
Keyword(s):  
Sensitivity Analysis ◽  
Limit Cycle ◽  
Stokes Equations ◽  
Finite Amplitude ◽  
Cylinder Wake ◽  
Mean Flow ◽  
Cycle Frequency ◽  
Flow Perturbation ◽  
Self Consistent ◽  
The Mean

We use the adjoint method to compute sensitivity maps for the limit-cycle frequency and amplitude of the Bénard–von Kármán vortex street in the wake of a circular cylinder. The sensitivity analysis is performed in the frame of the semi-linear self-consistent model recently introduced by Mantič et al. (Phys. Rev. Lett., vol. 113, 2014, 084501), which allows us to describe accurately the effect of the control on the mean flow, but also on the finite-amplitude fluctuation that couples back nonlinearly onto the mean flow via the formation of Reynolds stress. The sensitivity is computed with respect to arbitrary steady and synchronous time-harmonic body forces. For a small amplitude of the control, the theoretical variations of the limit-cycle frequency predict well those of the controlled flow, as obtained from either self-consistent modelling or direct numerical simulation of the Navier–Stokes equations. This is not the case if the variations are computed in the simpler mean flow approach overlooking the coupling between the mean and fluctuating components of the flow perturbation induced by the control. The variations of the limit-cycle amplitude (that falls out the scope of the mean flow approach) are also correctly predicted, meaning that the approach can serve as a relevant and systematic guideline to control strongly unstable flows exhibiting non-small, finite amplitudes of oscillation. As an illustration, we apply the method to control by means of a small secondary control cylinder and discuss the obtained results in the light of the seminal experiments of Strykowski & Sreenivasan (J. Fluid Mech., vol. 218, 1990, pp. 71–107).

The acoustic analogy in an annular duct with swirling mean flow

2013 ◽  
Vol 726 ◽  
pp. 439-475 ◽  
Author(s):  
H. Posson ◽  
N. Peake
Keyword(s):  
Green’S Function ◽  
Green's Function ◽  
Stokes Equations ◽  
Swirling Flow ◽  
Base Flow ◽  
Order Differential Operator ◽  
Source Terms ◽  
Acoustic Analogy ◽  
Mean Flow ◽  
Annular Duct

AbstractThis paper is concerned with modelling the effects of swirling flow on turbomachinery noise. We develop an acoustic analogy to predict sound generation in a swirling and sheared base flow in an annular duct, including the presence of moving solid surfaces to account for blade rows. In so doing we have extended a number of classical earlier results, including Ffowcs Williams & Hawkings’ equation in a medium at rest with moving surfaces, and Lilley’s equation for a sheared but non-swirling jet. By rearranging the Navier–Stokes equations we find a single equation, in the form of a sixth-order differential operator acting on the fluctuating pressure field on the left-hand side and a series of volume and surface source terms on the right-hand side; the form of these source terms depends strongly on the presence of swirl and radial shear. The integral form of this equation is then derived, using the Green’s function tailored to the base flow in the (rigid) duct. As is often the case in duct acoustics, it is then convenient to move into temporal, axial and azimuthal Fourier space, where the Green’s function is computed numerically. This formulation can then be applied to a number of turbomachinery noise sources. For definiteness here we consider the noise produced downstream when a steady distortion flow is incident on the fan from upstream, and compare our results with those obtained using a simplistic but commonly used Doppler correction method. We show that in all but the simplest case the full inclusion of swirl within an acoustic analogy, as described in this paper, is required.

Turbulent diffusion near a free surface

2000 ◽  
Vol 407 ◽  
pp. 145-166 ◽  
Author(s):  
LIAN SHEN ◽  
GEORGE S. TRIANTAFYLLOU ◽  
DICK K. P. YUE
Keyword(s):  
Boundary Layer ◽  
Free Surface ◽  
Turbulent Diffusion ◽  
Similarity Solution ◽  
Stokes Equations ◽  
Numerical Study ◽  
Turbulent Shear ◽  
Surface Boundary ◽  
Mean Flow ◽  
The Mean

We study numerically and analytically the turbulent diffusion characteristics in a low-Froude-number turbulent shear flow beneath a free surface. In the numerical study, the Navier–Stokes equations are solved directly subject to viscous boundary conditions at the free surface. From an ensemble of such simulations, we find that a boundary layer develops at the free surface characterized by a fast reduction in the value of the eddy viscosity. As the free surface is approached, the magnitude of the mean shear initially increases over the boundary (outer) layer, reaches a maximum and then drops to zero inside a much thinner inner layer. To understand and model this behaviour, we derive an analytical similarity solution for the mean flow. This solution predicts well the shape and the time-scaling behaviour of the mean flow obtained in the direct simulations. The theoretical solution is then used to derive scaling relations for the thickness of the inner and outer layers. Based on this similarity solution, we propose a free-surface function model for large-eddy simulations of free-surface turbulence. This new model correctly accounts for the variations of the Smagorinsky coefficient over the free-surface boundary layer and is validated in both a priori and a posteriori tests.

scholarly journals Computation of Second-order Design Sensitivities for Steady State Incompressible Laminar Flows Using the Extended Complex Variables Method

2020 ◽  
Author(s):  
Mahdi Hassanzadeh
Keyword(s):  
Sensitivity Analysis ◽  
Steady State ◽  
Stokes Equations ◽  
General Procedure ◽  
Finite Difference Methods ◽  
Computational Cost ◽  
Complex Variables ◽  
Second Order ◽  
Wide Range ◽  
Extended Complex

In the current paper, the general procedure of the first and second-order sensitivity analysis is presented using the extended complex variables method (ECVM). In the traditional complex variables method, only the imaginary step is used for sensitivity analysis. However, in the ECVM, both of the real and imaginary parts are employed to improve the efficiency of the method. To show this, the ECVM is applied to the steady state incompressible laminar flow around a cylinder. The governing Navier-Stokes equations are solved by the finite element method and then the ECVM is employed. The results are validated through comparing with those of obtained by an analytical as well as the finite difference methods and the convergence rate is investigated. It is illustrated that the first-order sensitivity analysis is not influenced by the change of the step length for both of the traditional and extended complex variables methods. However, it is shown that unlike the traditional complex variables method, the ECVM is less dependent on the step size for calculating the second-order sensitivity. This can be considered as an enhancement in the efficiency of this method. Hence, the ECVM is suggested as an appropriate technique for calculating simultaneously the first and second-order sensitivities with high accuracy as well as low computational cost. The proposed method is applicable to a wide range of problems having simple or complex parameters.

Mean flow deformation in a laminar separation bubble: separation and stability characteristics

2010 ◽  
Vol 660 ◽  
pp. 37-54 ◽  
Author(s):  
OLAF MARXEN ◽  
ULRICH RIST
Keyword(s):  
Flat Plate ◽  
Stokes Equations ◽  
Separation Bubble ◽  
Plate Boundary ◽  
Base Flow ◽  
Laminar Separation Bubble ◽  
Mean Flow ◽  
Laminar Separation ◽  
The Mean ◽  
Flow Deformation

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.

Sign in / Sign up

Export Citation Format

Share Document