We present an effective reduced-order model (ROM) technique to couple an incompressible flow with a transversely vibrating bluff body in a state-space format. The ROM of the unsteady wake flow is based on the Navier–Stokes equations and is constructed by means of an eigensystem realization algorithm (ERA). We investigate the underlying mechanism of vortex-induced vibration (VIV) of a circular cylinder at low Reynolds number via linear stability analysis. To understand the frequency lock-in mechanism and self-sustained VIV phenomenon, a systematic analysis is performed by examining the eigenvalue trajectories of the ERA-based ROM for a range of reduced oscillation frequency $(F_{s})$ , while maintaining fixed values of the Reynolds number ( $Re$ ) and mass ratio ( $m^{\ast }$ ). The effects of the Reynolds number $Re$ , the mass ratio $m^{\ast }$ and the rounding of a square cylinder are examined to generalize the proposed ERA-based ROM for the VIV lock-in analysis. The considered cylinder configurations are a basic square with sharp corners, a circle and three intermediate rounded squares, which are created by varying a single rounding parameter. The results show that the two frequency lock-in regimes, the so-called resonance and flutter, only exist when certain conditions are satisfied, and the regimes have a strong dependence on the shape of the bluff body, the Reynolds number and the mass ratio. In addition, the frequency lock-in during VIV of a square cylinder is found to be dominated by the resonance regime, without any coupled-mode flutter at low Reynolds number. To further discern the influence of geometry on the VIV lock-in mechanism, we consider the smooth curve geometry of an ellipse and two sharp corner geometries of forward triangle and diamond-shaped bluff bodies. While the ellipse and diamond geometries exhibit the flutter and mixed resonance–flutter regimes, the forward triangle undergoes only the flutter-induced lock-in for $30\leqslant Re\leqslant 100$ at $m^{\ast }=10$ . In the case of the forward triangle configuration, the ERA-based ROM accurately predicts the low-frequency galloping instability. We observe a kink in the amplitude response associated with 1:3 synchronization, whereby the forward triangular body oscillates at a single dominant frequency but the lift force has a frequency component at three times the body oscillation frequency. Finally, we present a stability phase diagram to summarize the VIV lock-in regimes of the five smooth-curve- and sharp-corner-based bluff bodies. These findings attempt to generalize our understanding of the VIV lock-in mechanism for bluff bodies at low Reynolds number. The proposed ERA-based ROM is found to be accurate, efficient and easy to use for the linear stability analysis of VIV, and it can have a profound impact on the development of control strategies for nonlinear vortex shedding and VIV.

We present an active feedback blowing and suction (AFBS) procedure via model reduction for unsteady wake flow and the vortex-induced vibration (VIV) of circular cylinders. The reduced-order model (ROM) for the AFBS procedure is developed by the eigensystem realization algorithm (ERA), which provides a low-order representation of the unsteady flow dynamics in the neighbourhood of the equilibrium steady state. The actuation is considered via vertical suction and a blowing jet at the porous surface of a circular cylinder with a body-mounted force sensor. While the optimal gain is obtained using a linear quadratic regulator (LQR), Kalman filtering is employed to estimate the approximate state vector. The feedback control system shifts the unstable eigenvalues of the wake flow and the VIV system to the left half-complex-plane, and subsequently results in suppression of the vortex street and the VIV in elastically mounted structures. The resulting controller designed by a linear low-order approximation is able to suppress the nonlinear saturated state of wake vortex shedding from the circular cylinder. A systematic linear ROM-based stability analysis is performed to understand the eigenvalue distribution for the flow past stationary and elastically mounted circular cylinders. The results from the ROM analysis are consistent with those obtained from full nonlinear fluid–structure interaction simulations, thereby confirming the validity of the proposed ROM-based AFBS procedure. A sensitivity study on the number of suction/blowing actuators, the angular arrangement of actuators and the combined versus independent control architectures has been performed for the flow past a stationary circular cylinder. Overall, the proposed control concept based on the ERA-based ROM and the LQR algorithm is found to be effective in suppressing the vortex street and the VIV for a range of reduced velocities and mass ratios.

This work proposes a novel approach to compute transonic limit-cycle oscillations using high-fidelity analysis. Computational-Fluid-Dynamics based harmonic balance methods have proven to be efficient tools to predict periodic phenomena. This paper’s contribution is to present a new methodology to determine the unknown frequency of oscillations, enabling harmonic balance methods to accurately capture limit-cycle oscillations; this is achieved by defining a frequency-updating procedure based on a coupled computational-fluid-dynamics/computational-structural-dynamics harmonic balance formulation to find the limit-cycle oscillation condition. A pitch/plunge airfoil and delta wing aerodynamic and respective linear structural models are used to validate the new method against conventional time-domain simulations. Results show consistent agreement between the proposed and time-marching methods for both limit-cycle oscillation amplitude and frequency while producing at least a one-order-of-magnitude reduction in computational time.

The co-flow jet is a zero-net-mass-flux active flow control strategy and presents great potential to improve the aerodynamic efficiency of future fuel-efficient aircrafts. The present work is to integrate the co-flow jet technology into aerodynamic shape optimization to further realize the potential of co-flow-jet technology and improve co-flow jet airfoil performance. The optimization results show that the maximum energy efficiency ratio of lift augmentation and drag reduction increased by 203.53% (α = 0°) and 10.25% (α = 10°) at the Power-1 condition (power coefficient of 0.3), respectively. A larger curvature is observed near the leading edge of the optimized aerodynamic shape, which leads to the early onset of flow separation and improves energy transfer efficiency from the jet to the free stream. In addition, the higher mid-chord of the optimized airfoil is characterized by accelerating the flow in the middle of the airfoil, increasing the strength of the negative pressure zone, thus improving the stall margin and enhancing the co-flow jet circulation.

A novel surrogate model is proposed in lieu of computational-fluid-dynamics solvers, for fast nonlinear aerodynamic and aeroelastic modeling. A nonlinear function is identified on selected interpolation points by a discrete empirical interpolation method. The flowfield is then reconstructed using a least-square approximation of the flow modes extracted by proper orthogonal decomposition. The aeroelastic reduced-order model is completed by introducing a nonlinear mapping function between displacements and the discrete empirical interpolation method points. The proposed model is investigated to predict the aerodynamic forces due to forced motions using a NACA 0012 airfoil undergoing a prescribed pitching oscillation. To investigate aeroelastic problems at transonic conditions, a pitch/plunge airfoil and a cropped delta wing aeroelastic models are built using linear structural models. The presence of shock waves triggers the appearance of limit-cycle oscillations, which the model is able to predict. For all cases tested, the new reduced-order model shows the ability to replicate the nonlinear aerodynamic forces and structural displacements and reconstruct the complete flowfield with sufficient accuracy at a fraction of the cost of full-order computational-fluid-dynamics model.