A harmonic balance technique for the reduced-order computation of vortex-induced vibration

We present a harmonic balance (HB) method to model frequency lock-in effect during vortex-induced vibration (VIV) of elastically mounted circular cylinder and a flexible riser section in a freestream uniform flow. The fluid flow and structure are coupled by a fixed-point iteration process through a frequency updating algorithm. By minimizing the structural residual in the standard least-square norm, the convergence of HB-based fixed-point algorithm is achieved for a range of reduced velocity. To begin with, the HB solver is first assessed for a periodic unsteady flow around a stationary circular cylinder. A freely vibrating circular cylinder is then adopted for the reduced-order computation of VIV at low Reynolds numbers of Re=100 and 180 with one- and two-degrees-of-freedom. The coupled VIV dynamics and the frequency lock-in phenomenon are accurately captured. The results show that the HB solver is able to predict the amplitude of vibration, frequency and forces comparable to its time domain counterpart, while providing a significant reduction with regard to overall computational cost. The proposed new scheme is then demonstrated for a fully-coupled three dimensional (3D) analysis of a linear-elastic riser section undergoing vortex-induced vibration in the lock-in range. The results reveal the 3D effects through isosurfaces of streamwise vorticity blobs distributed over the span of flexible riser section. In comparison to time domain results, the 3D flow-structure interactions are accurately predicted while providing a similar speed up rate that of 2D simulations. This further corroborates that the HB solver can be extended to 3D flow-structure dynamics without compromising efficiency and accuracy.     

Leading edge slot blowing on an open cavity in supersonic flow

Supersonic flow over an open cavity can create intense pressure loads on the surfaces within the cavity. In order to combat these loads, the development of a control scheme to reduce them is becoming increasingly important for many engineering applications. The present study implements steady leading edge blowing through various configurations of spanwise-aligned rectangular leading edge slots. The effects of this control on the flow field were examined to determine the suppression mechanisms exploited by the leading edge blowing. The cavity studied here had a length-to-depth ratio of 6 and was placed in a freestream flow with a Mach number of 1.4. Actuators with one continuous slot and three and five segmented slots spanning the width of the cavity were installed at the leading edge. Surface pressure reductions of nearly 45% were achieved on the aft wall of the cavity using the 5-slot configuration. Velocity field measurements acquired through 2-component (streamwise-aligned measurement plane) and 3-component stereoscopic (cross-stream-aligned measurement plane) particle image velocimetry revealed the presence of streamwise-aligned vortices created by the segmented slots. These act to significantly alter the shear layer formed at the mouth of the cavity creating highly three-dimensional flow field features.     

Numerical study of adverse pressure gradient generation over a flat plate using a rotating cylinder

Generating an adverse pressure gradient (APG), using a rotating cylinder in the proximity of a plane wall under a laminar freestream flow, is studied numerically in this work. The magnitude of the generated APG is a function of the gap, G, between the cylinder and the wall, and the rotational speed of the cylinder, Ω. The flow in such a configuration is characterized by periodic transient vortex shedding at high Reynolds number. A numerical model for the computation of the transient flow for this configuration is developed using the ANSYS CFD simulation tool. The model is validated against published experimental and numerical data for similar flow configurations and excellent agreement is observed. A parametric study is carried out for different combinations of G and Ω for two different Reynolds numbers of 200 and 1000 to examine the development of the resulting separation bubble due to the generated APG. The mechanism of the boundary layer separation over the plane wall and the corresponding wake dynamics is investigated. Results are presented in terms of the distribution of the pressure coefficient as well as skin friction coefficient along the wall and flow patterns around and downstream of the cylinder in the proximity of the wall. The results of these computations confirm that using a rotating cylinder over a plane wall in a freestream flow is an effective technique to generate a controlled range of adverse pressure gradients.     

Numerical and experimental investigation of the means for reducing the aeroacoustic loads in an extended rectangular cavity at subsonic and transonic freestream velocities

The results of a comprehensive investigation including numerical calculations and experiments with models in a wind tunnel and a vehicle under flight conditions aiming to find the ways of reducing the pressure fluctuation levels in an extended cavity at subsonic and transonic freestream velocities are presented. It is shown that the reduction of these loads can be achieved using the means which have demonstrated their effectiveness for cavities with the open-type flows, for example, a permeable deflector and the bevelling of the rear wall but only in the case of a given combination of their geometric parameters. The mechanisms of the action of these devices on the flow, thanks to which the intensity of the wave disturbances generated by the rear wall is reduced, the instability wave growth in the mixing layer behind the deflector is limited, and the fluctuation level in the cavity decreases, are investigated. The results of numerical investigations of the flow in a cavity with a permeable deflector are apparently among the first.     

Drag and lift reduction of a 3D bluff-body using active vortex generators

In this study, a passive flow control experiment on a 3D bluff-body using vortex generators (VGs) is presented. The bluff-body is a modified Ahmed body (Ahmed in J Fluids Eng 105:429–434 1983) with a curved rear part, instead of a slanted one, so that the location of the flow separation is no longer forced by the geometry. The influence of a line of non-conventional trapezoïdal VGs on the aerodynamic forces (drag and lift) induced on the bluff-body is investigated. The high sensitivity to many geometric (angle between the trapezoïdal element and the wall, spanwise spacing between the VGs, longitudinal location on the curved surface) and physical (freestream velocity) parameters is clearly demonstrated. The maximum drag reduction is ?12%, while the maximum global lift reduction can reach more than ?60%, with a strong dependency on the freestream velocity. For some configurations, the lift on the rear axle of the model can be inverted (?104%). It is also shown that the VGs are still efficient even downstream of the natural separation line. Finally, a dynamic parameter is chosen and a new set-up with motorized vortex generators is proposed. Thanks to this active device. The optimal configurations depending on two parameters are found more easily, and a significant drag and lift reduction (up to ?14% drag reduction) can be reached for different freestream velocities. These results are then analyzed through wall pressure and velocity measurements in the near-wake of the bluff-body with and without control. It appears that the largest drag and lift reduction is clearly associated to a strong increase of the size of the recirculation bubble over the rear slant. Investigation of the velocity field in a cross-section downstream the model reveals that, in the same time, the intensity of the longitudinal trailing vortices is strongly reduced, suggesting that the drag reduction is due to the breakdown of the balance between the separation bubble and the longitudinal vortices. It demonstrates that for low aspect ratio 3D bluff-bodies, like road vehicles, the flow control strategy is much different from the one used on airfoils: an early separation of the boundary layer can lead to a significant drag reduction if the circulation of the trailing vortices is reduced.     

Bypass transition receptivity modes

The prediction of bypass transition remains an important problem in many engineering applications. This is largely because there is no suitable theoretical model for bypass transition and predictions are made using empirical models. This paper presents numerical results for the receptivity of a zero pressure gradient boundary layer subjected to simple freestream waveforms which are the constituent parts of a turbulent flow field. Significant receptivities are only obtained for a minority of freestream waveforms and these lead to two types of flow structure in the boundary layer. The first type of flow structure is essentially two dimensional in nature and consists of two rows of counter-rotating spanwise vortices and is induced by freestream waves of large normal and spanwise wavelength and streamwise wavelengths approximately equal to the boundary layer thickness. The second type of flow structure are the streamwise streaks frequently observed in flow visualisation experiments. These streaks are induced by freestream waves of long streamwise and normal wavelength and spanwise wavelengths in the range of 14.5-46 θ (1.7-5.4δ). The freestream waves can be formed of velocity components in any direction, however the boundary layer is most receptive to fluctuations that lie in a plane perpendicular to the streamwise direction. The overall receptivity to a full spectrum of waves typical of freestream turbulence is considered and is shown to have similar characteristics to those from experiments.     

Unsteadiness of the shock wave structure in attached and separated compression ramp flows

Wall pressure fluctuations have been measured upstream of the corner-line in several two dimensional, adiabatic, compression ramp flows. The nominal freestream Mach number was 3 and the Reynolds number, based on boundary layer thickness, was 1.4 million. The measurements show that the shockwave structure is unsteady in both separated and attached flows, resulting in a region in which the wall pressure signal is intermittent. Statistical properties of this intermittent region, and of the separated flow, are presented and correlated with results from other studies.     

The Laminar Boundary Layer over a Permeable Wall

An analysis is given of the laminar boundary layer over a permeable/porous wall. The porous wall is passive in the sense that no suction or blowing velocity is imposed. To describe the flow inside and above the porous wall a continuum approach is employed based on the Volume-Averaging Method (S. Whitaker The method of volume averaging). With help of an order-of-magnitude analysis the boundary-layer equations are derived. The analysis is constrained by: (a) a low wall permeability; (b) a low Reynolds number for the flow inside the porous wall; (c) a sufficiently high Reynolds number for the freestream flow above the porous wall. Two boundary layers lying on top of each other can be distinguished: the Prandtl boundary layer above the porous wall, and the Brinkman boundary layer inside the porous wall. Based on the analytical solution for the Brinkman boundary layer in combination with the momentum transfer model of Ochoa-Tapia and Whitaker (Int. J. Heat Mass Transfer 38 (1995) 2635). for the interface region, a closed set of equations is derived for the Prandtl boundary layer. For the stream function a power series expansion in the perturbation parameter is adopted, where is proportional to ratio of the Brinkman to the Prandtl boundary-layer thickness. A generalization of the Falkner–Skan equation for boundary-layer flow past a wedge is derived, in which wall permeability is incorporated. Numerical solutions of the Falkner–Skan equation for various wedge angles are presented. Up to the first order in wall permeability causes a positive streamwise velocity at the interface and inside the porous wall, but a wall-normal interface velocity is a second-order effect. Furthermore, wall permeability causes a decrease in the wall shear stress when the freestream flow accelerates, but an increase in the wall shear stress when the freestream flow decelerates. From the latter it follows that separation, as indicated by zero wall shear stress, is delayed to a larger positive pressure gradient.     

有限长正方形棱柱绕流的双稳态现象

运用大涡模拟(large eddy simulation, LES) 方法对有限长正方形棱柱绕流进行了研究. 棱柱高宽比为5,一端固定于平板上,另一端为自由端. 平板表面边界层厚度可忽略不计. 基于自由流速度和棱柱宽度的雷诺数为3 900. 由计算结果发现,有限长棱柱的尾流在自由端后下扫流的作用下具有显著的三维性. 尾流中会交替出现两种典型的流动状态:一种是柱体两侧所脱落的涡在尾流中形成类似卡门涡街的交错排列状态;另一种是柱体尾流中涡团对称分布且无明显的周期性. 两种典型流动状态的交替出现,对有限长棱柱气动力特性有直接的影响.     

Controlling the near vortex wake of a swept half-wing by means of mini-flaps

The effect of mini-flaps on the flow pattern in the near vortex wake behind a model swept half-wing is investigated. The distributions of the time-average flow velocity were measured in a subsonic wind tunnel, in a section normal to the freestream velocity vector located at a distance of 3.8 wing half-spans from its trailing edge. When mini-flaps are mounted on both upper and lower wing surfaces, two vortices (tip and auxiliary) of the same sign are observable in the above-mentioned flow section; they are separated by an extended region of vorticity of the opposite sign. The model angle-of-attack effect on the intensities of the tip and auxiliary vortices is estimated.     

侧加热腔体内重力波演化过程的数值模拟

利用二维数值模拟的方法研究了侧加热腔体内的自然对流．基于数值模拟结果,描述了水平热入侵流(intrusion)的整个演化过程,并对该过程的物理机制进行了讨论．结果表明：当热入侵流抵达腔体冷壁后,由于冷壁无法卷入所有的热入侵流,热入侵流在冷上角堆积并产生一个反向流动,在冷壁边界层附近形成一个顺时针涡,该涡在浮力效应驱动下可返回热壁,并在腔体的冷热壁之间形成了腔体尺度的流体振荡,即内重力波．     

Existence of supersonic zones in three-dimensional separation flows

Up till now the region of three-dimensional separation flows which occur with supersonic flow past obstacles has received insufficient study. Supersonic flow with a Mach number of 2.5 past a cylinder mounted on a plate was studied in [1]. A local zone with supersonic velocities was found in the reverse subsonic flow region ahead of the cylinder. Its presence is explained by the three-dimensional nature of the flow. Similar supersonic zones are not observed in the case of supersonic flow over plane and axisymmetric steps.The present paper presents the results of experimental studies whose objective was refinement of the flow pattern ahead of a cylinder on a plate and the study of the local supersonic zones.The experiments were performed in a supersonic wind tunnel with a freestream Mach number M₁=3.11. The 24-mm-diameter cylinder with pressure taps along the generating line was mounted perpendicular to the surface of a sharpened plate. The distance from the plate leading edge to the cylinder axis wasl ₀=140 mm. The plate was pressure tapped along the flow symmetry axis. The Reynolds number was Rl ₀=u₀ l ₀/v ₁, Rl ₀=1.87.10⁷, where u₁ andv ₁ are the freestream velocity and the kinematic viscosity, respectively. The pressures were measured using a Pilot probe with internal and external diameters of 0.15 and 0.9 mm, respectively.The probe was displaced in the flow symmetry plane at a distance of 1.6 mm from the plate surface and at a distance of 1.1 mm along the leading generator of the cylinder. The flow on the surface of the plate and cylinder was studied with the aid of a visualization composition and the flow past the model was photographed with a schlieren instrument. Typical patterns of the visualization composition distribution and the pressure distribution curves over the plate surface, and also photographs of the flow past the model, are shown in [1].     

Effects of notches on the surface pressure field downstream of a wall-mounted spoiler

Notched spoilers have been observed to be more effective than uniform spoilers to suppress the flow-induced cavity resonance of vehicles with open sunroofs. In this study, a few mechanisms possibly involved in buffeting suppression from notched spoilers were investigated experimentally and numerically. One objective was to investigate the spatial coherence and phase of the wall pressure fluctuations downstream of notched spoilers in comparison with the same quantities for uniform spoilers. Another objective was to gather detailed measured data to allow the verification of computer simulations of the flow over the notched spoiler. Experiments were performed to measure the velocity and wall pressure fields downstream of spoilers mounted on the rigid floor of a closed test section wind tunnel. Efforts were made to reproduce the spoiler and wind tunnel geometry and boundary conditions of the experimental setup in the numerical simulations. The numerical investigation used the Lattice Boltzmann Method (LBM), with the so-called Very Large Eddy Simulation (VLES) viscosity turbulence model. The results of the numerical investigation were in satisfactory agreement with measured data at low frequencies, where buffeting is expected to occur. The results suggested that the notches break down the homogeneity of the leading edge cross-stream vortices predominantly responsible for the cavity excitation. This decreased the cross-stream coherence of the surface pressure field, thereby reducing the magnitude of the net equivalent excitation force acting over the surface downstream.     

Passively controlled supersonic cavity flow using a spanwise cylinder

An experimental investigation of a passively controlled open cavity with a length to depth ratio of six and freestream Mach number of 1.4 was conducted to investigate the mechanisms responsible for the observed surface pressure reductions. The passive control comes from placing a spanwise aligned cylinder in the boundary layer near the leading edge of the cavity. The two control configurations were isolated from previous experiments of the fluctuating surface pressure and correspond to a larger diameter rod near the top of the boundary layer and a smaller diameter rod placed near the wall. These were further analyzed using particle image velocimetry in an attempt to elicit the responsible mechanism for the flow control. The use of two-point statistics revealed the wall normal turbulent velocity correlation’s evolution became elongated in the wall normal direction. This suggests that the shear layer may be less-organized and consists of smaller-scale structures. The disturbance of the feedback receptivity loop is clearly demonstrated for the controlled configurations evidenced by weakened correlation signals between the aft wall sensor and positions on the cavity floor. The presence of the rod is shown to decrease the mean shear gradient, more effectively for the large rod placed at the top of the boundary layer, throughout the shear layer. The efficacy of the control leads to an initially thicker shear layer which spreads more rapidly and is clearly demonstrated by vorticity growth rates, mean, and turbulent flowfield statistics.     

The Effect of Pressure Gradient on Boundary Layer Receptivity

This paper presents numerical results for the receptivity of three laminar boundary layers with zero (ZPG), adverse (APG) and favourable (FPG) pressure gradients. Each boundary layer is subjected to a series of simple freestream waveforms which can be considered as constituent parts of either an isotropic or a non-isotropic turbulent freestream. Each freestream waveform has a single frequency in each spatial direction and is divided into two mutually perpendicular components. The first component has a zero spanwise velocity and hence lies in the streamwise normal plane whereas the second component lies in a plane which is perpendicular both to this plane and the spatial frequency vector. High boundary layer receptivities are only obtained for a minority of these waveforms and so only the resulting flow structures for these waveforms are considered in detail. The dominant flow structures are identified as either Tollmien Schlichting (T-S) waves or streaky structures. The streaky structures can be induced by both freestream components, but the response to the second component, which results in streamwise vortices in the freestream, is considerably stronger and occurs over a much larger streamwise frequency range. The boundary layer is only receptive to a relatively narrow band of spanwise wavelengths ranging from approximately one to four times the local boundary layer thickness. The APG leads to receptivities which are more than double those for the FPG case. The ratio of the freestream fluctuation streamwise wavelength to the distance from the plate leading edge is identified as an important influential parameter for receptivity leading to streaks. Significant T-S activity is only observed for APG, but is also detected for ZPG.     

Numerical study of turbulent flow over two‐dimensional surface‐mounted ribs in a channel

This study accurately predicts the cases of turbulent flow around a surface‐mounted two‐dimensional rib with varying lengths. The numerical method employs a differencing scheme for integrating the elliptic Reynolds‐averaged Navier–Stokes equations and the continuity equation. A two‐equation k–ε turbulence model is employed to simulate the turbulent transport quantities and close the solving problem. The near‐wall regions of the separated sides of the rib are resolved by a near‐wall model of a two‐layer approach instead of the wall function approximation. Computations for flow over a surface‐mounted rectangular rib are conducted for the variations in the rib lengths. Results indicate that upstream of the obstacle, the length of the recirculating region remains unchanged with varying rib lengths; while the downstream length of the recirculating region is a strong function of rib length and changes nearly linearly for the varying lengths of B/H=0.1 to B/H=4.0. Reattachment on top of the rib, owing to its increasing length, affects the downstream boundary layer development. Copyright © 1999 John Wiley & Sons, Ltd.     

Simulation of three‐dimensional flows over moving objects by an improved immersed boundary–lattice Boltzmann method

An improved immersed boundary–lattice Boltzmann method (IB–LBM) developed recently [28] was applied in this work to simulate three‐dimensional (3D) flows over moving objects. By enforcing the non‐slip boundary condition, the method could avoid any flow penetration to the wall. In the developed IB–LBM solver, the flow field is obtained on the non‐uniform mesh by the efficient LBM that is based on the second‐order one‐dimensional interpolation. As a consequence, its coefficients could be computed simply. By simulating flows over a stationary sphere and torus [28] accurately and efficiently, the proposed IB–LBM showed its ability to handle 3D flow problems with curved boundaries. In this paper, we further applied this method to simulate 3D flows around moving boundaries. As a first example, the flow over a rotating sphere was simulated. The obtained results agreed very well with the previous data in the literature. Then, simulation of flow over a rotating torus was conducted. The capability of the improved IB–LBM for solving 3D flows over moving objects with complex geometries was demonstrated via the simulations of fish swimming and dragonfly flight. The numerical results displayed quantitative and qualitative agreement with the date in the literature. Copyright © 2011 John Wiley & Sons, Ltd.     

Modeling the Effect of Freestream Turbulence on Unsteady Boundary Layer Flow

We present a new version of the two-equation turbulence model, which makes it possible to calculate continuously the entire flow range from laminar to turbulent, including transition, in the case of a time-periodic, high-turbulence-level freestream. The influence of the parameters characterizing the harmonic fluctuations of the external velocity and the freestream turbulence intensity and scale on the parameters of the flat-plate flow are analyzed. A comparison of the numerical solutions with the experimental and theoretical data indicates the possibility of describing the wall flow properties on the basis of a quasi-stationary turbulence model, as the Reynolds number varies from low to high values.     

A low‐dimensional spectral approach for transient axisymmetric free‐surface flow inside thin cavities of arbitrary shape

A low‐dimensional spectral method is used to solve the transient axisymmetric free surface flow inside thin cavities of arbitrary shape. The flow field is obtained on the basis of the lubrication equations, which are expanded in terms of orthonormal functions over the cavity gap. The formulation accounts for nonlinearities stemming from inertia and front location. The work is of close relevance to the filling stage during die casting, and injection molding, or the flow inside annular (extrusion) dies. Both flows under an imposed flow rate, and an imposed pressure at the cavity entrance are examined. The influence of inertia, aspect ratio, gravity, and wall geometry on the evolution of the front, flow rate, and pressure is assessed particularly in the early stage of flow, when a temporal behavior of the ‘boundary‐layer’ type develops. The multiple‐scale method is applied to obtain an approximate solution at small Reynolds number, Re. Comparison with the exact (numerical) solution indicates a wide range of validity for the multiple‐scale approach, including the moderately small Re range. Copyright © 2002 John Wiley & Sons, Ltd.     

Effects of dual jets distance on mixing characteristics and flow path within a cavity in supersonic crossflow

A rectangular open cavity with upstream dual injectors at a freestream Mach number of 1.9 was investigated experimentally. To evaluate the effect of the distance between the jets, the flow characteristics were investigated using the high-speed schlieren photography, particle image velocimetry, and surface oil flow techniques. The dual jet distances of 18 and 54 mm were used. Unstable flow occurs over the cavity in all cases and is not improved by changing the distance between the dual jets. Although the distance between the dual jets does not influence the flow stability, the flow field varies decidedly depending on the dual jets distance. The enhancement of air mixing depends on the distance between the jets. A long dual jets distance was found to yield better mixing characteristics within the cavity than a short one. When the jets are further apart, the mainstream between two counter-rotating vortex pairs behind the jets flows strongly into the cavity because of the increased blow-down occurring between the vortex pairs. Additionally, a counterflow with a low velocity magnitude occurs behind the jets. Hence, mixing is enhanced within the cavity by effects of the opposed flow. When the jet pairs are closer to each other, the counter-rotating vortex pairs are in contact; as a result, the blow-down effect does not occur between them. The flow drawn into the cavity from the mainstream is supplied from the sides of the test section into the cavity.