Assessment of low-altitude atmospheric turbulence models for aircraft aeroelasticity

We investigate the dynamic aeroelastic response of large but slow aircraft in low-altitude atmospheric turbulence. To this end, three turbulence models of increasing fidelity, namely, the one-dimensional von Kármán model, the two-dimensional Kaimal model and full three-dimensional wind fields extracted from large-eddy simulations (LES) are used to simulate ambient turbulence near the ground. Load calculations and flight trajectory predictions are conducted for a representative high-aspect-ratio wing aircraft, using a fully coupled nonlinear flight dynamics/aeroelastic model, when it operates in background atmospheric turbulence generated by the aforementioned models. Comparison of load envelopes and spectral content, on vehicles of varying flexibility, shows strong dependency between the selected turbulence model and aircraft aeroelastic response (e.g. 58% difference in the predicted magnitude of the wing root bending moment between LES and von Kármán models). This is mainly due to the presence of large flow structures at low altitudes that have comparable dimensions to the vehicle, and which despite the relatively small wind speeds within the Earth boundary layer, result in overall high load events for slow-moving vehicles. Results show that one-dimensional models that do not capture those effects provide fairly non-conservative load estimates and are unsuitable for very flexible airframe design.     

Dynamic wind loads and wake characteristics of a wind turbine model in an atmospheric boundary layer wind

An experimental study was conducted to characterize the dynamic wind loads and evolution of the unsteady vortex and turbulent flow structures in the near wake of a horizontal axis wind turbine model placed in an atmospheric boundary layer wind tunnel. In addition to measuring dynamic wind loads (i.e., aerodynamic forces and bending moments) acting on the wind turbine model by using a high-sensitive force-moment sensor unit, a high-resolution digital particle image velocimetry (PIV) system was used to achieve flow field measurements to quantify the characteristics of the turbulent vortex flow in the near wake of the wind turbine model. Besides conducting “free-run” PIV measurements to determine the ensemble-averaged statistics of the flow quantities such as mean velocity, Reynolds stress, and turbulence kinetic energy (TKE) distributions in the wake flow, “phase-locked” PIV measurements were also performed to elucidate further details about evolution of the unsteady vortex structures in the wake flow in relation to the position of the rotating turbine blades. The effects of the tip-speed-ratio of the wind turbine model on the dynamic wind loads and wake flow characteristics were quantified in the terms of the variations of the aerodynamic thrust and bending moment coefficients of the wind turbine model, the evolution of the helical tip vortices and the unsteady vortices shedding from the blade roots and turbine nacelle, the deceleration of the incoming airflows after passing the rotation disk of the turbine blades, the TKE and Reynolds stress distributions in the near wake of the wind turbine model. The detailed flow field measurements were correlated with the dynamic wind load measurements to elucidate underlying physics in order to gain further insight into the characteristics of the dynamic wind loads and turbulent vortex flows in the wakes of wind turbines for the optimal design of the wind turbines operating in atmospheric boundary layer winds.     

Numerical investigation of the turbulent energy budget in the wake of freely oscillating elastically mounted cylinder at low reduced velocities

We present a numerical study of the turbulent kinetic energy budget in the wake of cylinders undergoing Vortex-Induced Vibration (VIV). We show three-dimensional Large Eddy Simulations (LES) of an elastically mounted circular cylinder in the synchronization regime at Reynolds number of Re=8000. The Immersed Boundary Method (IBM) is used to account for the presence of the cylinder. The flow field in the wake is decomposed using the triple decomposition splitting the flow variables in mean, coherent and stochastic components. The energy transfer between these scales of motions are then studied and the results of the free oscillation are compared to those of a forced oscillation. The turbulent kinetic energy budget shows that the maximum amplitude of VIV is defined by the ability of the mean flow to feed energy to the coherent structures in the wake. At amplitudes above this maximum amplitude, the energy of the coherent structures needs to be fed additionally by small scale, stochastic energy in form of backscatter to sustain its motion. Furthermore, we demonstrate that the maximum amplitude of the VIV is defined by the integral length scale of the turbulence in the wake.     

Stability of an incompressible two-dimensional wake

In experimental investigations of the wake flow behind a plate, a monochromatic (in time) signal is usually observed immediately behind the end of the plate. Downstream, the signal is distorted and then becomes random, i.e., a turbulent flow regime is realized. Theoretically, a branch point is found at the experimentally observed frequency in the spectrum of three-dimensional perturbations of the problem linearized with respect to the steady solution [1]. Mattingly and Criminale [1] attribute all the characteristics of the observed signal to this point. As in other similar investigations, the mechanism of the appearance of the monochromatic signal in the near wake was not elucidated in [1]. In the present paper, the problem of the characteristic oscillations of the flow in the near wake is studied. The appearance of the monochromatic signal is explained by the presence in the near wake of a standing wave of the required frequency, the wave being formed by two scattering points. The first is the end of the plate, and the second the branch point in the spectrum of linear three-dimensional perturbations.     

An Assessment of Detached-Eddy Simulations of Unsteady Crosswind Aerodynamics of Road Vehicles

In this paper, the availability of Detached-Eddy Simulations (DES) for simulating unsteady crosswind aerodynamics in industrial applications is explored and reported. The flow around a simple car geometry under headwind, steady crosswind and time-dependent wind gust respectively, is analysed for two different types of mesh using a commercial software, STAR-CD. The typical Reynolds number of the flow cases studied is 2.0 · 10⁶ based on the vehicle length. The performances of polyhedral meshes together with a comparable hexahedral mesh are evaluated for capturing the time development of the flow structures during a gust. From the results presented, it can be concluded that the aerodynamic coefficients can be reliably calculated using both polyhedral and hexahedral meshes, although the fine turbulent structures in the wake are better captured with the hexahedral mesh. Using turbulence methods like DES in order to resolve the flow scales provides a significant insight for designing a ground vehicle and, due to the reasonable computational times involved, can hence be foreseen to be incorporated in a design process in a near future.     

Flow/acoustic interactions of two cylinders in cross-flow

The noise generated by two tandem cylinders in a cross-flow (i.e., with the second in the wake of the first) has been investigated. Measurements of turbulence and of fluctuating pressure have been obtained between the two cylinders for different flow velocities and incident levels of turbulence. Although, for a number of cases, up to four peaks related to vortex shedding were evident in the spectrum, most measurements exhibited two peaks, a dominant one at the vortex-shedding frequency, with a secondary peak at twice this value. The measurements show that vortex generated noise is strongest at the mid-point between the cylinders and at the rear cylinder with levels of 130 dB. The harmonic component was strongest at the downstream cylinder where peak values of 110 dB were obtained. The nonlinear flow/acoustic interactions are examined using bispectral analysis to identify the quadratic interactions in the parameters. A novel quadratic modelling method is proposed and shown to be capable of both identifying and quantifying the nonlinear interactions which give rise to noise at harmonics of the vortex-shedding frequency.     

Solid/free-surface juncture boundary layer and wake

 The Reynolds-averaged flow for a solid/free-surface juncture boundary layer and wake is documented. The three mean-velocity components and five of the Reynolds stresses are measured for a surface-piercing flat plate in a towing tank using a laser-Doppler velocimeter system for both boundary-layer and wake planes in regions close to the free surface. The experimental method is described, including the foil-plate model, laser-Doppler velocimeter system, conditions, and uncertainty analysis. The underlying flow data is in excellent agreement with benchmark data. Inner (near the plate and wake centerplane and below the free surface) and outer (near the free surface) regions of high streamwise vorticity of opposite sign are observed, which transport, respectively, high mean velocity and low turbulence from the outer to the inner and low mean velocity and high turbulence from the inner to the outer portions of the boundary layer and wake. For the wake, the inner region of vorticity is relatively weak. The physical mechanism for the streamwise vorticity is analyzed with regard to the Reynolds-averaged streamwise vorticity equation. The anisotropy of the crossplane normal Reynolds stresses closely correlates with the vorticity and, additionally, indicates similarity, i.e., its nature is such that it only depends on the proximity to the plate and free surface boundaries or wake centerplane symmetry plane. Free-surface effects on the Reynolds stresses are analyzed with regard to the behavior close to the free surface of the turbulent kinetic energy and the normal components of the anisotropy tensor and the anisotropy invariants. Close to the free surface, the turbulent kinetic energy is nearly constant and increases for the inner and outer portions, respectively, of the boundary layer and wake and the normal components of the anisotropy tensor and the anisotropy invariants roughly correspond to the limiting values for two-component turbulence. The similarities and differences between the present results and analysis with those from related studies are discussed. The data and analysis should have practical application with regard to the development of turbulence models for computational fluid dynamics methods for the Reynolds-averaged Navier–Stokes equations. Received: 27 May 1997/Accepted: 1 August 1997     

Planar velocity measurements of a two-dimensional compressible wake

The present study describes the application of particle image velocimetry (PIV) to investigate the compressible flow in the wake of a two-dimensional blunt base at a freestream Mach number M_X=2. The first part of the study addresses specific issues related to the application of PIV to supersonic wind tunnel flows, such as the seeding particle flow-tracing fidelity and the measurement spatial resolution. The seeding particle response is assessed through a planar oblique shock wave experiment. The measurement spatial resolution is enhanced by means of an advanced image-interrogation algorithm. In the second part, the experimental results are presented. The PIV measurements yield the spatial distribution of mean velocity and turbulence. The mean velocity distribution clearly reveals the main flow features such as expansion fans, separated shear layers, flow recirculation, reattachment, recompression and wake development. The turbulence distribution shows the growth of turbulent fluctuations in the separated shear layers up to the reattachment location. Increased velocity fluctuations are also present downstream of reattachment outside of the wake due to unsteady flow reattachment and recompression. The instantaneous velocity field is analyzed seeking coherent flow structures in the redeveloping wake. The instantaneous planar velocity and vorticity measurements return evidence of large-scale turbulent structures detected as spatially coherent vorticity fluctuations. The velocity pattern consistently shows large masses of fluid in vortical motion. The overall instantaneous wake flow is organized as a double row of counter-rotating structures. The single structures show vorticity contours of roughly elliptical shape in agreement with previous studies based on spatial correlation of planar light scattering. Peak vorticity is found to be five times higher than the mean vorticity value, suggesting that wake turbulence is dominated by the activity of large-scale structures. The unsteady behavior of the reattachment phenomenon is studied. Based on the instantaneous flow topology, the reattachment is observed to fluctuate mostly in the streamwise direction suggesting that the unsteady separation is dominated by a pumping-like motion.     

Investigation on 3Dt wake flow structures of swimming bionic fish

A bionic experimental platform was designed for the purpose of investigating time accurate three-dimensional flow field, using digital particle image velocimetry (DSPIV). The wake behind the flapping trail of a robotic fish model was studied at high spatial resolution. The study was performed in a water channel. A robot fish model was designed and built. The model was fixed onto a rigid support framework using a cable-supporting method, with twelve stretched wires. The entire tail of the model can perform prescribed motions in two degrees of freedom, mainly in carangiform mode, by driving its afterbody and lunate caudal fin respectively. The DSPIV system was set up to operate in a translational manner, measuring velocity field in a series of parallel slices. Phase locked measurements were repeated for a number of runs, allowing reconstruction of phase average flow field. Vortex structures with phase history of the wake were obtained. The study reveals some new and complex three-dimensional flow structures in the wake of the fish, including “reverse hairpin vortex” and “reverse Karman S-H vortex rings”, allowing insight into physics of this complex flow.     

Studies on gas turbulence and particle fluctuation in dense gas-particle flows

Particle fluctuation and gas turbulence in dense gas-particle flows are less studied due to complexity of the phenomena. In the present study, simulations of gas turbulent flows passing over a single particle are carried out first by using RANS modeling with a Reynolds stress equation turbulence model and sufficiently fine grids, and then by using LES. The turbulence enhancement by the particle wake effect is studied under various particle sizes and relative gas velocities, and the turbulence enhancement is found proportional to the particle diameter and the square of velocity. Based on the above results, a turbulence enhancement model for the particle-wake effect is proposed and is incorporated as a sub-model into a comprehensive two-phase flow model, which is then used to simulate dilute gas-particle flows in a horizontal channel. The simulation results show that the predicted gas turbulence by using the present model accounting for the particle wake effect is obviously in better agreement with the experimental results than the prediction given by the model not accounting for the wake effect. Finally, the proposed model is incorporated into another two-phase flow model to simulate dense gasparticle flows in a downer. The results show that the particle wake effect not only enhances the gas turbulence, but also amplifies the particle fluctuation.     

Modes of vortex formation and transition to three-dimensionality in the wake of a freely vibrating cylinder

The present work is aimed to give some insight into the relation between vortex shedding modes and transition to three-dimensionality in the wake of a freely vibrating cylinder by establishing a numerical model and analyzing the relevant results of two- and three-dimensional simulations. The compressible flow past an elastically-mounted cylinder is solved by using the immersed boundary method (IB method). The cylinder is free to vibrate in the transverse direction with zero structure damping. The response of displacement amplitude is studied with the variation of reduced velocity. Whether P+S mode exists in three-dimensional flow and the occurrence of 2P mode is caused by flow transition from two-dimensional to three-dimensional are problems of concern. Both 2P and P+S wake modes are observed in two- and three-dimensional simulations. The numerical results indicate that the flow transition from two-dimensional to three-dimensional is coupled with the cylinder vibration in the synchronization/lock-in regime. The wake formation given by three-dimensional simulations suggests that the P+S mode might exist in reality when the flow is reverted to two-dimensional by vortex induced vibration (VIV) at Re=300–350. When Reynolds number increases to 425, the wake formation undergoes transition to three-dimensionality and 2P mode is observed. The effect of mass ratio on the flow transition to three-dimensionality is studied. The relationship between wake modes and aerodynamic forces is discussed.     

Drag reduction of motor vehicles by active flow control using the Coanda effect

A test facility has been constructed to realistically simulate the flow around a two dimensional car shaped body in a wind tunnel. A moving belt simulator has been employed to generate the relative motion between model and ground. In a first step, the aerodynamic coefficients c _L and c _D of the model are determined using static pressure and force measurements. LDA-measurements behind the model show the large vortex and turbulence structures of the near and far wake. In a second step, the ambient flow around the model is modified by way of an active flow control which uses the Coanda effect, whereby the base-pressure increases by nearly 50% and the total drag can be reduced by 10%. The recirculating region is completely eliminated. The current work reveals the fundamental physical phenomena of the new method by observing the pressure forces on the model surface as well as the time averaged velocities and turbulence distributions for the near and far wake. A theory resting on this empirical information is developed and provides information about the effectiveness of the blowing method. For this, momentum and energy equations were applied to the flow around the vehicle to enable a validation of the theoretical results using experimental values. Received: 9 June 1998 / Accepted: 20 July 1999     

DIRECT NUMERICAL SIMULATIONS ON THE TURBULENT FLOW PAST A CONFINED CIRCULAR CYLINDER WITH THE INFLUENCE OF THE STREAMWISE MAGNETIC FIELDS1)

绕流是托卡马克装置中液态包层内常见的流动形态,对流场与热量分布有着重要的影响.本文通过直接数值模拟(DNS),研究了不同磁场强度下$Re=3900$的圆柱绕流,分析了磁场强度对于湍流尾迹的影响.无磁场情况下,直接数值模拟的结果与前人的实验及模拟结果吻合很好.圆柱下游的尾迹中,随着流向距离的增大, 流向速度剖面逐渐从U型进化呈V型, 并慢慢趋于平缓,这表明尾迹中的流动结构受圆柱影响逐渐减小.圆柱后方两侧的剪切层中,由于Kelvin-Helmholtz不稳定性的影响,可以清晰地看到小尺度剪切层涡的脱落.通过对无磁场的计算结果施加流向磁场,本文计算了哈特曼数($Ha$)分别为20, 40和80的工况,以研究磁场效应对于湍流的影响.结果表明磁场较弱时,流动依然呈三维湍流状态.随着磁场增强, 近圆柱尾流区受磁场抑制明显,回流区被拉长,剪切层失稳位置向下游转移.圆柱后方的涡结构由于受到竖直方向洛伦兹力的挤压作用,随着哈特曼数的增加尾迹区域逐渐变窄.相比于无磁场情况的涡结构,由于磁场的耗散作用,相应的涡结构尺度变小.该研究不仅扩展了现有磁场下湍流运动的参数范围,对于液态包层的设计及安全运行同样具有重要的理论指导意义和工程应用价值.     

ENERGY DISTRIBUTION IN MODES IN THE WAKE OF A FINITE-LENGTH CYLINDER BEFORE AND AFTER TRANSITION

The flow in the wake of a finite-length cylinder has been studied experimentally both before and after the transition to turbulence. This instability occurs at Reynolds numbers around 180–190. One end of the cylinder was fixed to the bottom of the test-section of a wind tunnel, whilst the other terminated in the open flow (free end). For these boundary conditions four main frequency modes within the wake can be identified. These are a centre-cell mode at the Strouhal frequency, end-cell modes with a frequency below the Strouhal frequency, a mode exhibiting the difference frequency between the centre-cell and end-cell modes, and a low-frequency mode (appearing only after transition to turbulence). In this work the energy content of these four modes has been determined throughout the wake, both before and after transition to turbulence. For three of the modes, the energy content is the same before and after transition, whereas the low-frequency mode exhibits energy two to four orders of magnitude greater after transition than before. Hence it is clear that the additional turbulence energy appearing in the wake after transition is located predominantly in this low-frequency mode. The appearance of this low-frequency mode is characterized by the simultaneous appearance of a peak in the power spectra of the velocity fluctuations centred about zero frequency (but with finite width). Consequently, the appearance of this zero frequency peak can be taken as the signature of the onset of turbulence. By considering the downstream growth rates of this low-frequency mode, evidence is presented which suggests that transition to turbulence may occur as a result of wake transition in the downstream central plane of the cylinder.     

Three-component velocity field measurements of propeller wake using a stereoscopic PIV technique

A stereoscopic PIV (Particle Image Velocimetry) technique was used to measure the three-dimensional flow structure of the turbulent wake behind a marine propeller with five blades. The out-of-plane velocity component was determined using two CCD cameras with an angular displacement configuration. Four hundred instantaneous velocity fields were measured for each of four different blade phases, and ensemble averaged in order to find the spatial evolution of the propeller wake in the region from the trailing edge up to one propeller diameter (D) downstream. The influence of propeller loading conditions on the wake structure was also investigated by measuring the velocity fields at three advance ratios (J=0.59, 0.72 and 0.88). The phase-averaged velocity fields revealed that a viscous wake formed by the boundary layers developed along the blade surfaces. Tip vortices were generated periodically and the slipstream contracted in the near-wake region. The out-of-plane velocity component and strain rate had large values at the locations of the tip and trailing vortices. As the flow moved downstream, the turbulence intensity, the strength of the tip vortices, and the magnitude of the out-of-plane velocity component at trailing vortices all decreased due to effects such as viscous dissipation, turbulence diffusion, and blade-to-blade interaction.     

Non-equilibrium turbulent phenomena in the ?ow over a backward-facing ramp

Non-equilibrium turbulence phenomena have raised great interests in recent years. Significant efforts have been devoted to non-equilibrium turbulence properties in canonical flows, e.g., grid turbulence, turbulent wakes, and homogeneous isotropic turbulence(HIT). The non-equilibrium turbulence in non-canonical flows, however, has rarely been studied due to the complexity of the flows. In the present contribution, a directnumerical simulation(DNS) database of a turbulent flow is analyzed over a backwardfacing ramp, the flow near the boundary is demonstrated, and the non-equilibrium turbulent properties of the flow in the wake of the ramp are presented by using the characteristic parameters such as the dissipation coefficient C and the skewness of longitudinal velocity gradient Sk, but with opposite underlying turbulent energy transfer properties. The equation of Lagrangian velocity gradient correlation is examined, and the results show that non-equilibrium turbulence is the result of phase de-coherence phenomena, which is not taken into account in the modeling of non-equilibrium turbulence. These findings are expected to inspire deeper investigation of different non-equilibrium turbulence phenomena in different flow conditions and the improvement of turbulence modeling.     

Experimental investigation of the turbulent wake past real seabed elements for velocity variations characterization in the water column.

In high flow velocity areas like those suitable for marine energy application, bathymetry variations create strong velocity fluctuations in the water column. It is therefore essential to characterize the turbulence evolution in the wake of seabed elements which may impact the loads on tidal turbines. For that purpose, experiments are carried out in a flume tank with R_e as high as achievable in Froude similitude, with bathymetry variations experimentally represented with various wall-mounted square elements of height H: a cylinder or a cube as unitary obstacles and combinations of these elements followed by an inclined floor to resemble smooth bathymetry changes. The onset flow is a simple boundary layer profile with height 1.3 H and a low turbulence intensity. PIV and LDV measurements are used to investigate the wake past all test cases in order to distinguish high floor elevation cases (unitary obstacles) from mean roughness effect (obstacle combinations). Results show that the obstacle combinations produce a wake less extended than for a single wide cylinder that produces an extended wake and very energetic turbulent events. With a single cube, no downstream development of large turbulent events exist and the wake reduces by a factor of 3 compared to the wake cylinder case. An inclined floor downstream of a single wall-mounted obstacle reduces its wake length but does not alter the turbulent structures shed. Turbulent velocity profiles extracted from every wake topology investigated are also compared. The general conclusion is that: for small aspect ratio cases, the obstacle will not affect the water column. On the contrary, strong energetic turbulent events are emitted from large aspect ratio obstacles. Combinations cases stand in-between.     

Modeling the Turbulent Vortex Wake Behind a High-Lift Wing

A mathematical model and the results of a numerical simulation of the multiple vortex structure behind the high-lift wing (with deflected flaps) of a civil aircraft model are presented. The calculations are performed within the framework of the simplified (with allowance for the specific features of the flow under consideration) three-dimensional Reynolds equations using three turbulence models. The results obtained are compared with the experimental data; a considerable influence of the streamline curvature, typical of vortex flows, not only on the turbulence parameters but also on the parameters of the average flow is demonstrated. As a result, only one of the turbulence models considered, namely, that taking the streamline curvature and rotation effects on the turbulence into account, ensures acceptable accuracy of the main wake characteristics.     

Experimental investigation of a blunt trailing edge flow field with application to sound generation

The unsteady lift generated by turbulence at the trailing edge of an airfoil is a source of radiated sound. The objective of the present research was to measure the velocity field in the near wake region of an asymmetric beveled trailing edge in order to determine the flow mechanisms responsible for the generation of trailing edge noise. Two component velocity measurements were acquired using particle image velocimetry. The chord Reynolds number was 1.9 × 10⁶. The data show velocity field realizations that were typical of a wake flow containing an asymmetric periodic vortex shedding. A phase average decomposition of the velocity field with respect to this shedding process was utilized to separate the large scale turbulent motions that occurred at the vortex shedding frequency (i.e., those responsible for the production of tonal noise) from the smaller scale turbulent motions, which were interpreted to be responsible for the production of broadband sound. The small scale turbulence was found to be dependent on the phase of the vortex shedding process implying a dependence of the broadband sound generated by the trailing edge on the phase of the vortex shedding process.     

Towards multi-fidelity simulation of flows around an underwater vehicle with appendages and propeller

Flow around a real-life underwater vehicle often happens at a high Reynolds number with flow structures at different scales from the boundary layer around a blade to that around the hull. This poses a great challenge for large-eddy simulation of an underwater vehicle aiming at resolving all relevant flow scales. In this work, we propose to model the hull with appendages using the immersed boundary method, and model the propeller using the actuator disk model without resolving the geometry of the blade. The proposed method is then applied to simulate the flow around Defense Advanced Research Projects Agency(DARPA) suboff. An overall acceptable agreement is obtained for the pressure and friction coefficients. Complex flow features are observed in the near wake of suboff. In the far wake, the core region is featured by a jet because of the actuator disk, surrounded by an annular region with velocity deficit due to the body of suboff.