Application of wake oscillators to two-dimensional vortex-induced vibrations of circular cylinders in oscillatory flows

A nonlinear time-domain simulation model for predicting two-dimensional vortex-induced vibration (VIV) of a flexibly mounted circular cylinder in planar and oscillatory flow is presented. This model is based on the utilization of van der Pol wake oscillators, being unconventional since wake oscillators have typically been applied to steady flow VIV predictions. The time-varying relative flow–cylinder velocities and accelerations are accounted for in deriving the coupled hydrodynamic lift, drag and inertia forces leading to the cylinder cross-flow and in-line oscillations. The system fluid–structure interaction equations explicitly contain the time-dependent and hybrid trigonometric terms. Depending on the Keulegan–Carpenter number (KC) incorporating the flow maximum velocity and excitation frequency, the model calibration is performed, entailing a set of empirical coefficients and expressions as a function of KC and mass ratio. Parametric investigations in cases of varying KC, reduced flow velocity, cylinder-to-flow frequency ratio and mass ratio are carried out, capturing some qualitative features of oscillatory flow VIV and exploring the effects of system parameters on response prediction characteristics. The model dependence of hydrodynamic coefficients on the Reynolds number is studied. Discrepancies and limitations versus advantages of the present model with different feasible solution scenarios are illuminated to inform the implementation of wake oscillators as a computationally efficient prediction model for VIV in oscillatory flows.     

A nonlinear model for aerodynamic configuration of wake behind horizontal-axis wind turbine

Determination of the aerodynamic configuration of wake is the key to analysis and evaluation of the rotor aerodynamic characteristics of a horizontal-axis wind turbine.According to the aerodynamic configuration, the real magnitude and direction of the onflow velocity at the rotor blade can be determined, and subsequently, the aerodynamic force on the rotor can be determined. The commonly employed wake aerodynamic models are of the cylindrical form instead of the actual expanding one. This is because the influence of the radial component of the induced velocity on the wake configuration is neglected. Therefore, this model should be called a "linear model". Using this model means that the induced velocities at the rotor blades and aerodynamic loads on them would be inexact. An approximately accurate approach is proposed in this paper to determine the so-called "nonlinear" wake aerodynamic configuration by means of the potential theory,where the influence of all three coordinate components of the induced velocity on wake aerodynamic configuration is taken into account to obtain a kind of expanding wake that approximately looks like an actual one. First, the rotor aerodynamic model composed of axial(central), bound, and trailing vortexes is established with the help of the finite aspect wing theory. Then, the Biot-Savart formula for the potential flow theory is used to derive a set of integral equations to evaluate the three components of the induced velocity at any point within the wake. The numerical solution to the integral equations is found,and the loci of all elementary trailing vortex filaments behind the rotor are determined thereafter. Finally, to formulate an actual wind turbine rotor, using the nonlinear wake model, the induced velocity everywhere in the wake, especially that at the rotor blade,is obtained in the case of various tip speed ratios and compared with the wake boundary in a neutral atmospheric boundary layer. Hereby, some useful and referential conclusions are offered for the aerodynamic computation and design of the rotor of the horizontal-axis wind turbine.     

均匀、振荡及组合流绕平板流动的尾迹研究

本文利用离散涡模型及改进的新生涡产生机制对三种不同来流绕平板的近尾迹进行数值研究。计算讨论了定常流中平板绕流流动的总体特性和近尾迹流场;对于简谐振荡来流,相应于K_c=2.0、4.0 和10.0 分别得到两种不同的尾迹形态。给出了小 K_c 数平板尾迹涡配对、运动的新模式而相应的阻力、惯性力系数计算比以前涡模拟结果更接近于 U 型管实验结果。对于流向组合来流本文模拟了涡锁定及其动力特性并于实验相符,给出了流向扰动对平板绕流流动的影响。     

Frequency lock-in during nonlinear vibration of an airfoil coupled with van der Pol Oscillator

The frequency lock-in during the nonlinear vibration of a turbomachinery blade is modeled using a spring-mounted airfoil coupled with a van der Pol Oscillator (VDP) oscillator. The proposed reduced-order model uses the nonlinear VDP oscillator to represent the oscillatory nature of wake dynamics caused by the vortex shedding. The damping term in the VDP oscillator is assumed to be nonlinear. The coupled equations governing the pitch and plunge motion of an airfoil are used to approximate the vibration of a turbomachinery blade. Springs having cubic-order nonlinearity for their stiffnesses are used to mount the airfoil. The unsteady lift acting on the blade is modeled using a self-excited nonlinear wake oscillator. The model for wake dynamics takes into account the influence of blade inertia. The nonlinear coupled three degrees of freedom (dof) aeroelastic system is studied for instability resulting in the frequency lock-in phenomenon. The equations are transformed into non-dimensional form, and then the frequencies of the coupled system are plotted to demonstrate the frequency lock-in. Further, the method of multiple scales is used to derive modulation equations which represent the amplitude and phase of the oscillation. The results obtained using the method of multiple scales are compared with direct numerical solutions to verify the present modeling method. The steady-state amplitudes of the response are plotted against the detuning parameter, which represents the frequency response curve. Further, the sensitivity of non-dimensional parameters such as coupling coefficients, mass ratio, reduced velocity, static unbalance, structural damping coefficient and the ratio of uncoupled pitch and plunge natural frequencies on the frequency response is investigated. The study revealed that parameters such as mass ratio, reduced velocity, structural damping coefficient, and coupling coefficients have a stronger influence in suppressing the amplitude of vibration. Meanwhile, parameters such as the frequency ratio, static unbalance, reduced velocity, and mass ratio significantly affect the range of frequency in which the lock-in phenomenon happens. Further, linear perturbation analysis is done to understand the qualitative effect of the system parameters such as coupling coefficients, mass ratio, frequency ratio, and static unbalance on the range of lock-in.     

Cluster-based control of a separating flow over a smoothly contoured ramp

The ability to manipulate and control fluid flows is of great importance in many scientific and engineering applications. The proposed closed-loop control framework addresses a key issue of model-based control: The actuation effect often results from slow dynamics of strongly nonlinear interactions which the flow reveals at timescales much longer than the prediction horizon of any model. Hence, we employ a probabilistic approach based on a cluster-based discretization of the Liouville equation for the evolution of the probability distribution. The proposed methodology frames high-dimensional, nonlinear dynamics into low-dimensional, probabilistic, linear dynamics which considerably simplifies the optimal control problem while preserving nonlinear actuation mechanisms. The data-driven approach builds upon a state space discretization using a clustering algorithm which groups kinematically similar flow states into a low number of clusters. The temporal evolution of the probability distribution on this set of clusters is then described by a control-dependent Markov model. This Markov model can be used as predictor for the ergodic probability distribution for a particular control law. This probability distribution approximates the long-term behavior of the original system on which basis the optimal control law is determined. We examine how the approach can be used to improve the open-loop actuation in a separating flow dominated by Kelvin–Helmholtz shedding. For this purpose, the feature space, in which the model is learned, and the admissible control inputs are tailored to strongly oscillatory flows.     

Stability of buoyancy-driven convection in an Oldroyd-B fluid-saturated anisotropic porous layer

The nonlinear stability of thermal convection in a layer of an Oldroyd-B fluid-saturated Darcy porous medium with anisotropic permeability and thermal diffusivity is investigated with the perturbation method. A modified Darcy-Oldroyd model is used to describe the flow in a layer of an anisotropic porous medium. The results of the linear instability theory are delineated. The thresholds for the stationary and oscillatory convection boundaries are established, and the crossover boundary between them is demarcated by identifying a codimension-two point in the viscoelastic parameter plane. The stability of the stationary and oscillatory bifurcating solutions is analyzed by deriving the cubic Landau equations. It shows that these solutions always bifurcate supercritically. The heat transfer is estimated in terms of the Nusselt number for the stationary and oscillatory modes. The result shows that, when the ratio of the thermal to mechanical anisotropy parameters increases, the heat transfer decreases.     

A NUMERICAL STUDY OF CYLINDERS IN WAVES AND CURRENTS

The research on combination flow of planar oscillatory flow plus an in-line steady stream is of importance to the situation of structures in waves and current. The combination flow has not been studied extensively. There is still little information about the effect of current and about combined effect of current and waves on hydrodynamic loading of the structures. The present study investigates the combination flow around a circular cylinder using a vortex-based method incorporating vortex moving particles (discrete vortices) with a finite-difference scheme for the vorticity diffusion. The main attention is paid to the effects of a small current on in-line fluid forces and vortex patterns in the wake. Morison's equation and an equation with two drag terms are examined. The results show that the presence of a small current in an oscillatory flow can reduce the drag coefficient significantly. Morison's equation gives reasonably good predictions for the in-line forces for an oscillatory flow plus a small current. The current tends to bring the whole vortex wake downstream and tries to form the stable Karman asymmetrical form in the downstream wake. The present results show certain agreement with some previous experimental results.     

Investigations of wake flows of a flat plate in steady,oscillatory and combined flows

Numerical study on near wake flows of a flat plate in three kinds of oncoming flows is made by using the discrete vortex model and improved vorticity creation method. For steady oncoming flow, both gross and detailed features of the wake flow are calculated and discussed. Then, in harmonic oscillatory oncoming flow two different wake flow patterns withK _c=2,4 and 10 are obtained respectively. Our results present a new wake flow pattern for lowKc numbers (Kc<5) describing vortex shedding, pairing and moving in a period of the oscillatory flow starting from rest. The calculated drag and inertia force coefficients are closer to experimental data from the U-tube than the previous results of vortex simulation. For in-line combined oncoming flow the vortex lock-in and dynamic characteristics are simulated. The results are shown to be in good agreement with experiments. The project supported by National Natural Science Fundation of China and LNM of Institute of Mechanics. CAS     

Modelling the non-linear interaction of wind and tide: Its influence on current profiles

A single-point model in the vertical is used to examine the coupling between tidal currents and wind-driven flows in shallow near-coastal regions. Calculations using both a linear slip and a no-slip condition at the sea bed clearly show that coupling between tidal and wind-driven currents cannot occur in a linear model with a time-independent eddy viscosity. However with a physically more realistic time-varying viscosity related to the flow field, coupling does occur, the magnitude of this non-linear interaction depending upon the change in eddy viscosity over a tidal cycle and the intensity of shear in the vertical. A point model in the vertical with flow induced by an oscillatory pressure gradient and an additional constant wind stress is used to examine the influence of viscosity parametrization and water depth upon this coupling. The solution in the vertical is accomplished using both a functional approach and a finite difference method. Some conclusions as to the relative merits of these approaches, particularly the use of a transformed grid in the case of high-shear surface and bed boundary layers, are made in the paper.     

A microstructural investigation of the nonlinear response of electrorheological suspensions

In the preceding paper, Part 1, the transition from linear to nonlinear behavior for electrorheological (ER) suspensions under start-up of steady shear flow was found to first arise from the slight rearrangement of unstable structures. In this paper, we investigate the transition to nonlinear behavior for ER suspensions under oscillatory shear flow, focusing on the role of the rearrangement of unstable structures, and employing experimental and simulation results. Again, we find that nonlinear deformation first arises from these rearrangements, as opposed to the gross rearrangement or rupture of particulate chains. The Fourier transform of the simulated time-dependent shear stress is employed to quantify the dependence of the critical strain on the deformation frequency and electric field strength. The predicted behavior is consistent with experimental trends. Methods for verifying the predictions are discussed, as well as possible avenues for exploiting this information in improved operating strategies and improved ER fluids.     

The onset of convection in a nanofluid saturated porous layer using Darcy model with cross diffusion

Linear and nonlinear stability analysis for the onset of convection in a horizontal layer of a porous medium saturated by a nanofluid is studied. The model used for the nanofluid incorporates the effects of Brownian motion and thermophoresis. The modified Darcy equation that includes the time derivative term is used to model the momentum equation. In conjunction with the Brownian motion, the nanoparticle fraction becomes stratified, hence the viscosity and the conductivity are stratified. The nanofluid is assumed to be diluted and this enables the porous medium to be treated as a weakly heterogeneous medium with variation, in the vertical direction, of conductivity and viscosity. The critical Rayleigh number, wave number for stationary and oscillatory mode and frequency of oscillations are obtained analytically using linear theory and the non-linear analysis is made with minimal representation of the truncated Fourier series analysis involving only two terms. The effect of various parameters on the stationary and oscillatory convection is shown pictorially. We also study the effect of time on transient Nusselt number and Sherwood number which is found to be oscillatory when time is small. However, when time becomes very large both the transient Nusselt value and Sherwood value approaches to their steady state values.     

Frequency lock-in is caused by coupled-mode flutter

The mechanism underlying the lock-in of frequencies in flow-induced vibrations is analysed using elementary linear dynamics. Considering the case of lock-in in vortex-induced vibrations (VIV), we use a standard wake oscillator model, as in previous studies, but in its simplest form where all nonlinear terms and all dissipative terms are neglected. The stability of the resulting linear system is analysed, and a range of coupled-mode flutter is found. In this range, the frequency of the most unstable mode is found to deviate from the Strouhal law when the frequency of the wake oscillator approaches that of the free cylinder motion. Simultaneously the growth rate resulting from coupled-mode flutter increases, which would lead to higher vibration amplitudes. The extent of the range of lock-in is then compared with experimental data, showing a good agreement. It is therefore stated that the lock-in phenomenon, such as in VIV, is a particular case of linear coupled-mode flutter.     

Adaptive Control of Nonlinear Dynamical Systems Using a Model Reference Approach

In this paper we consider using a model reference adaptive control approach to control nonlinear systems. We consider the controller design and stability analysis associated with these type of adaptive systems. Then we discuss the use of model reference adaptive control algorithms to control systems which exhibit nonlinear dynamical behaviour using the example of a Duffing oscillator being controlled to follow a linear reference model. For this system we show that if the nonlinearity is small then standard linear model reference control can be applied. A second example, which is often found in synchronization applications, is when the nonlinearities in the plant and reference model are identical. Again we show that linear model reference adaptive control is sufficient to control the system. Finally we consider controlling more general nonlinear systems using adaptive feedback linearization to control scalar nonlinear systems. As an example we use the Lorenz and Chua systems with parameter values such that they both have chaotic dynamics. The Lorenz system is used as a reference model and a single coordinate from the Chua system is controlled to follow one of the Lorenz system coordinates.     

Onset of Darcy–Brinkman convection in a binary viscoelastic fluid saturated porous layer

The onset of Darcy–Brinkman double-diffusive convection in a binary viscoelastic fluid-saturated porous layer is studied using both linear and weakly nonlinear stability analyses. The Oldroyd-B model is employed to describe the rheological behavior of the fluid. An extended form of Darcy–Oldroyd law incorporating the Brinkman’s correction and time derivative is used to describe the fluid flow and the Oberbeck–Boussinesq approximation is invoked. The onset criterion for stationary and oscillatory convection is derived analytically. The effects of rheological parameters, Darcy number, normalized porosity, Lewis number, solute Rayleigh number, and Darcy–Prandtl number on the stability of the system is investigated. The results indicated that there is a competition among the processes of thermal, solute diffusions and viscoelasticity that causes the convection to set in through the oscillatory modes rather than the stationary. The Darcy–Prandtl number has a dual effect on the threshold of oscillatory convection. The nonlinear theory based on the method of truncated representation of Fourier series is used to find the transient heat and mass transfer. Some existing results are reproduced as the particular cases of present study.     

Nonlinear dynamics of an aeroelastic airfoil with free-play in transonic flow

Nonlinear dynamic behaviors of an aeroelastic airfoil with free-play in transonic air flow are studied. The aeroelastic response is obtained by using time-marching approach with computational fluid dynamics (CFD) and reduced order model (ROM) techniques. Several standardized tests of transonic flutter are presented to validate numerical approaches. It is found that in time-marching approach with CFD technique, the time-step size has a significant effect on the calculated aeroelastic response, especially for cases considering both structural and aerodynamic nonlinearities. The nonlinear dynamic behavior for the present model in transonic air flow is greatly different from that in subsonic regime where only simple harmonic oscillations are observed. Major features of the responses in transonic air flow at different flow speeds can be summarized as follows. The aeroelastic responses with the amplitude near the free-play are dominated by single degree of freedom flutter mechanism, and snap-though phenomenon can be observed when the air speed is low. The bifurcation diagram can be captured by using ROM technique, and it is observed that the route to chaos for the present model is via period-doubling, which is essentially caused by the free-play nonlinearity. When the flow speed approaches the linear flutter speed, the aeroelastic system vibrates with large amplitude, which is dominated by the aerodynamic nonlinearity. Effects of boundary layer and airfoil profile on the nonlinear responses of the aeroelastic system are also discussed.     

Reduced‐order controllers for control of flow past an airfoil

Reduced‐order controller design by means of reduced‐order model for control of a wake flow is presented. Reduced‐order model is derived by combining the Galerkin projection with proper orthogonal decomposition (POD) or with other related reduced‐order approaches such as singular value decomposition or reduced‐basis method. In the present investigation, we discuss the applicability of the reduced‐order approaches for fast computation of the optimal control for control of vortex shedding behind a thin airfoil through unsteady blowing on the airfoil surface. Accuracy of the reduced‐order model is quantified by comparing flow fields obtained from the reduced‐order models with those from the full‐order simulations under the same free‐stream conditions. A control of vortex shedding is demonstrated for Reynolds number 100. It is found that downstream directed blowing on the upper surface of the airfoil near the leading edge is more efficient in mitigating flow separation and suppressing the vortex shedding. Copyright © 2005 John Wiley & Sons, Ltd.     

螺旋桨梢涡不稳定性机理与演化模型研究

螺旋桨尾流场的涡流特性是一个基础但又十分复杂的流体力学问题, 它的复杂性源于其蕴含复杂的漩涡系统, 且该漩涡系统会在高速的剪切层流动中不断演化, 其流体动力学行为, 如由稳定态演变为不稳定态的机理以及复杂工况环境中的流动现象, 一直是流体力学领域的难点和备受关注的热点问题. 从工程应用的角度看, 桨后梢涡的演化特性与船舶结构物的宏观特性直接相关, 更好地理解多工况下螺旋桨尾流的动力学特性, 将有助于改善与振动、噪声以及结构问题等相关的推进器性能, 对综合性能优良的下一代螺旋桨的设计和优化有着重要的现实意义. 本文基于延迟分离涡模拟、大涡模拟和无湍流模型模拟方法以及粒子图像测速流场测试分别开展了螺旋桨尾流动力学特性的数值与试验研究, 对螺旋桨尾流不稳定性的触发机理进行了揭示. 基于均匀来流中螺旋桨梢涡的演化机理, 提出了螺旋桨梢涡演化模型. 该模型能够较为准确地模拟螺旋桨梢涡的演化过程, 预测螺旋桨梢涡融合的时间和位置, 对螺旋桨流噪声预报和控制以及性能优良的螺旋桨设计具有重要意义.      

Experimental estimation of a D-shaped cylinder wake using body-mounted sensors

The effectiveness of a small array of body-mounted sensors, for estimation and eventually feedback flow control of a D-shaped cylinder wake is investigated experimentally. The research is aimed at suppressing unsteady loads resulting from the von-Kármán vortex shedding in the wake of bluff-bodies at a Reynolds number range of 100–1,000. A low-dimensional proper orthogonal decomposition (POD) procedure was applied to the stream-wise and cross-stream velocities in the near wake flow field, with steady-state vortex shedding, obtained using particle image velocimetry (PIV). Data were collected in the unforced condition, which served as a baseline, as well as during influence of forcing within the “lock-in” region. The design of sensor number and placement was based on data from a laminar direct numerical simulation of the Navier-Stokes equations. A linear stochastic estimator (LSE) was employed to map the surface-mounted hot-film sensor signals to the temporal coefficients of the reduced order model of the wake flow field in order to provide accurate yet compact estimates of the low-dimensional states. For a three-sensor configuration, results show that the estimation error of the first two cross-stream modes is within 20–40% of the PIV-generated POD time coefficients. Based on previous investigations, this level of error is acceptable for a moderately robust controller required for feedback flow control.     

Natural convection in a rotating anisotropic porous layer with internal heat generation

In this article, linear and nonlinear thermal instability in a rotating anisotropic porous layer with heat source has been investigated. The extended Darcy model, which includes the time derivative and Coriolis term has been employed in the momentum equation. The linear theory has been performed by using normal mode technique, while nonlinear analysis is based on minimal representation of the truncated Fourier series having only two terms. The criteria for both stationary and oscillatory convection is derived analytically. The rotation inhibits the onset of convection in both stationary and oscillatory modes. Effects of parameters on critical Rayleigh number has also been investigated. A weak nonlinear analysis based on the truncated representation of Fourier series method has been used to find the Nusselt number. The transient behavior of the Nusselt number has also been investigated by solving the finite amplitude equations using a numerical method. Steady and unsteady streamlines, and isotherms have been drawn to determine the nature of flow pattern. The results obtained during the analysis have been presented graphically.     

Superposition rheometry of a wormlike micellar fluid

Oscillatory measurements are often used to explore the nonlinear response of materials, with recently a strong focus on using large amplitude oscillatory experiments. However, the superposition of an oscillatory motion onto a steady-state shear flow is a method where the kinematic history experienced by the sample is simpler. Such a superposed oscillation can be applied either orthogonal or parallel to the main flow direction. Both superposed deformation modes can now be achieved on rotational rheometers equipped with a force-rebalanced transducer, the orthogonal mode requiring a minor modification to the control loop of the normal force. In the present work, the nonlinear properties of a wormlike micellar (WLM) solution are studied. The results are compared with the predictions of the Giesekus model, which is chosen both for its capability to describe the WLM response and for being one of the simplest continuum models that incorporate an anisotropic microstructure. From the fluid response in the homogeneous flow regime, a rate-dependent relaxation time and a rate-dependent plateau modulus can be derived. The latter provides insight into the structural anisotropy during flow at short length scales, which in this case is isotropic. Further analysis of the superposition moduli can be used to separate and quantify the effects of flow on the reptation and breaking of the chains. In the shear-banding regime, the orthogonal moduli show a weaker dependence on shear rate compared to the predictions of the Giesekus model, yet they remain sensitive to changes in the shear-banded state.