Showing papers on "Flow separation published in 2009"

Detached-Eddy Simulation

Abstract: Detached-eddy simulation (DES) was first proposed in 1997 and first used in 1999, so its full history can be surveyed. A DES community has formed, with adepts and critics, as well as new branches. The initial motivation of high–Reynolds number, massively separated flows remains, for which DES is convincingly more capable presently than either unsteady Reynolds-averaged Navier-Stokes (RANS) or large-eddy simulation (LES). This review discusses compelling examples, noting the visual and quantitative success of DES. Its principal weakness is its response to ambiguous grids, in which the wall-parallel grid spacing is of the order of the boundary-layer thickness. In some situations, DES on a given grid is then less accurate than RANS on the same grid or DES on a coarser grid. Partial remedies have been found, yet dealing with thickening boundary layers and shallow separation bubbles is a central challenge. The nonmonotonic response of DES to grid refinement is disturbing to most observers, as is the absence of...

Drag reduction in turbulent flows over superhydrophobic surfaces

Abstract: In this paper, we demonstrate that periodic, micropatterned superhydrophobic surfaces, previously noted for their ability to provide laminar flowdrag reduction, are capable of reducing drag in the turbulent flow regime. Superhydrophobic surfaces contain micro- or nanoscale hydrophobic features which can support a shear-free air-water interface between peaks in the surface topology. Particle image velocimetry and pressure drop measurements were used to observe significant slip velocities, shear stress, and pressure drop reductions corresponding to drag reductions approaching 50%. At a given Reynolds number,drag reduction is found to increase with increasing feature size and spacing, as in laminar flows. No observable drag reduction was noted in the laminar regime, consistent with previous experimental results for the channel geometry considered. The onset of drag reduction occurs at a critical Reynolds number where the viscous sublayer thickness approaches the scale of the superhydrophobic microfeatures and performance is seen to increase with further reduction in viscous sublayer height. These results indicate superhydrophobic surfaces may provide a significant drag reducing mechanism for marine vessels.

Single dielectric barrier discharge plasma enhanced aerodynamics: physics, modeling and applications

Abstract: The term “plasma actuator” has been a part of the fluid dynamics flow control vernacular for more than a decade. A particular type of plasma actuator that has gained wide use is based on a single dielectric barrier discharge (SDBD) mechanism that has desirable features for use in air at atmospheric pressures. For these actuators, the mechanism of flow control is through a generated body force vector that couples with the momentum in the external flow. The body force can be derived from first principles and the plasma actuator effect can be easily incorporated into flow solvers so that their placement and operation can be optimized. They have been used in a wide range of applications that include bluff body wake control; lift augmentation and separation control on a variety of lifting surfaces ranging from fixed wings with various degrees of sweep, wind turbine rotors and pitching airfoils simulating helicopter rotors; flow separation and tip-casing clearance flow control to reduce losses in turbines, to control flow surge and stall in compressors; and in exciting instabilities in boundary layers at subsonic to supersonic Mach numbers for turbulent transition control. New applications continue to appear through programs in a growing number of US universities and government laboratories, as well as in Germany, France, England, Netherland, Russia, Japan and China. This paper provides an overview of the physics, design and modeling of SDBD plasma actuators. It then presents their use in a number of applications that includes both numerical flow simulations and experiments together.

Three-dimensional flows around low-aspect-ratio flat-plate wings at low Reynolds numbers

Abstract: Three-dimensional flows over impulsively translated low-aspect-ratio flat plates are investigated for Reynolds numbers of 300 and 500, with a focus on the unsteady vortex dynamics at post-stall angles of attack. Numerical simulations, validated by an oil tow-tank experiment, are performed to study the influence of aspect ratio, angle of attack and planform geometry on the wake vortices and the resulting forces on the plate. Immediately following the impulsive start, the separated flows create wake vortices that share the same topology for all aspect ratios. At large time, the tip vortices significantly influence the vortex dynamics and the corresponding forces on the wings. Depending on the aspect ratio, angle of attack and Reynolds number, the flow at large time reaches a stable steady state, a periodic cycle or aperiodic shedding. For cases of high angles of attack, an asymmetric wake develops in the spanwise direction at large time. The present results are compared to higher Reynolds number flows. Some non-rectangular planforms are also considered to examine the difference in the wakes and forces. After the impulsive start, the time at which maximum lift occurs is fairly constant for a wide range of flow conditions during the initial transient. Due to the influence of the tip vortices, the three-dimensional dynamics of the wake vortices are found to be quite different from the two-dimensional von Karman vortex street in terms of stability and shedding frequency.

Microramp Control of Supersonic Oblique Shock-Wave/Boundary-Layer Interactions

Abstract: The potential of microramp sub-boundary-layer vortex generators for flow control in supersonic engine inlets is investigated. In particular, the study focuses on the ability of these devices to beneficially affect oblique shock-wave/ boundary-layer interactions. Experiments have been conducted at Mach 2.5 to determine the nature of flow controlled by microramps and to investigate their ability to delay separation in a reflected shock interaction. Various ramp heights between 30 and 90% of the boundary-layer thickness were investigated. The details of the vortical flow generated by such devices were identified. The general flow features were found to scale with device height and it is suggested that smaller devices need to be placed closer to the expected adverse pressure gradients. When applied to a separated oblique shock-wave/boundary-layer interaction generated with a 7 degree wedge, microramps were not able to completely eliminate flow separation, although they were shown to break up separated regions. Other performance indicators across the shock-wave/boundary-layer interaction were also improved through the application of the devices.

On vortex shedding from an airfoil in low-Reynolds-number flows

Abstract: Development of coherent structures in the separated shear layer and wake of an airfoil in low-Reynolds-number flows was studied experimentally for a range of airfoil chord Reynolds numbers, 55 × 10 3 ≤ Re c ≤ 210 × 10 3 , and three angles of attack, α = 0°, 5° and 10°. To illustrate the effect of separated shear layer development on the characteristics of coherent structures, experiments were conducted for two flow regimes common to airfoil operation at low Reynolds numbers: (i) boundary layer separation without reattachment and (ii) separation bubble formation. The results demonstrate that roll-up vortices form in the separated shear layer due to the amplification of natural disturbances, and these structures play a key role in flow transition to turbulence. The final stage of transition in the separated shear layer, associated with the growth of a sub-harmonic component of fundamental disturbances, is linked to the merging of the roll-up vortices. Turbulent wake vortex shedding is shown to occur for both flow regimes investigated. Each of the two flow regimes produces distinctly different characteristics of the roll-up and wake vortices. The study focuses on frequency scaling of the investigated coherent structures and the effect of flow regime on the frequency scaling. Analysis of the results and available data from previous experiments shows that the fundamental frequency of the shear layer vortices exhibits a power law dependency on the Reynolds number for both flow regimes. In contrast, the wake vortex shedding frequency is shown to vary linearly with the Reynolds number. An alternative frequency scaling is proposed, which results in a good collapse of experimental data across the investigated range of Reynolds numbers.

Experimental Study of Shock Oscillation over a Transonic Supercritical Profile

Abstract: Shock wave/turbulent boundary-layer interaction and flow separation may induce self-sustained large-scale oscillations on a profile at transonic Mach number. This phenomenon, known as transonic buffet, is at the origin of intense pressure fluctuations which can have detrimental effects, both in external and internal aerodynamics. The present paper describes a new experiment executed in the ONERA S3Ch transonic wind tunnel on shock oscillations over the OAT15A supercritical profile. These experiments have allowed the precise definition of the conditions for buffet onset and the characterization of the properties of the periodic motion from unsteady surface pressure measurements. The flowfield behavior has been described in great detail thanks to high-speed schlieren cinematography and surveys with a two-component laser Doppler velocimetry along with a conditional sampling technique. The first aim of this study was to provide the computational fluid dynamics community with well-documented test cases to validate advanced computing methods. Concerning the physics of the phenomenon, it is suggested that it is mediated by acoustic waves which are produced at the trailing edge and which travel on the two sides of the airfoil. Also, the experimental results strongly suggest that the phenomenon is essentially two-dimensional, even if three-dimensional effects are also detected.

Global stability of a jet in crossflow

Abstract: A linear stability analysis shows that the jet in crossflow is characterized by self-sustained global oscillations for a jet-to-crossflow velocity ratio of 3. A fully three-dimensional unstable steady-state solution and its associated global eigenmodes are computed by direct numerical simulations and iterative eigenvalue routines. The steady flow, obtained by means of selective frequency damping, consists mainly of a (steady) counter-rotating vortex pair (CVP) in the far field and horseshoe-shaped vortices close to the wall. High-frequency unstable global eigenmodes associated with shear-layer instabilities on the CVP and low-frequency modes associated with shedding vortices in the wake of the jet are identified. Furthermore, different spanwise symmetries of the global modes are discussed. This work constitutes the first simulation-based global stability analysis of a fully three-dimensional base flow.

Shallow and deep dynamic stall for flapping low Reynolds number airfoils

Abstract: We consider a combined experimental (based on flow visualization, direct force measurement and phaseaveraged 2D particle image velocimetry in a water tunnel), computational (2D Reynolds-averaged Navier-Stokes) and theoretical (Theodorsen’s formula) approach to study the fluid physics of rigid-airfoil pitch-plunge in nominally two-dimensional conditions. Shallow-stall (combined pitch-plunge) and deep-stall (pure-plunge) are compared at a reduced frequency commensurate with flapping-flight in cruise in nature. Objectives include assessment of how well attached-flow theory can predict lift coefficient even in the presence of significant separation, and how well 2D velocimetry and 2D computation can mutually validate one another. The shallow-stall case shows promising agreement between computation and experiment, while in the deepstall case, the computation’s prediction of flow separation lags that of the experiment, but eventually evinces qualitatively similar leading edge vortex size. Dye injection was found to give good qualitative match with particle image velocimetry in describing leading edge vortex formation and return to flow reattachment, and also gave evidence of strong spanwise growth of flow separation after leadingedge vortex formation. Reynolds number effects, in the range of 10,000-60,000, were found to influence the size of laminar separation in those phases of motion where instantaneous angle of attack was well below stall, but have limited effect on post-stall flowfield behavior. Discrepancy in lift coefficient time history between experiment, theory and computation was mutually comparable, with no clear failure of Theodorsen’s formula. This is surprising and encouraging, especially for the deep-stall case, because the theory’s assumptions are clearly violated, while its prediction of lift coefficient remains useful for capturing general trends.

High-Fidelity Simulation of Transitional Flows Past a Plunging Airfoil

Abstract: This investigation addresses the simulation of the unsteady separated flows encountered by a plunging airfoil under low-Reynolds-number conditions (Rec 6 ◊ 10 4 ). The flow fields are computed employing a previously developed and extensively validated high-fidelity implicit large-eddy simulation (ILES) approach. In order to permit comparison with available experimental measurements, calculations are performed first for an SD7003 airfoil section at an angle of attack o = 4 plunging with reduced frequency k = 3.93 and nondimensional amplitude ho = 0.05. Under these conditions, it is demonstrated that for Rec = 10 4 , transitional effects are not significant and that the dynamic-stall vortices remain fairly coherent as they propagate along the airfoil. For Rec = 4 ◊ 10 4 , the dynamic-stall vortex system is laminar at is inception, however shortly afterwards, it experiences an abrupt breakdown associated with the onset of spanwise instability effects. A detailed description of this transition process near the leading edge is provided. The computed phased-averaged structures for both values of Reynolds number are found to be in good agreement with the experimental data. As a second example, the suppression of static stall at high angle of attack ( o = 14 ) is investigated using high-frequency small-amplitude vibrations (k = 10,ho = 0.005). At Rec = 6 ◊ 10 4 , separation is completely eliminated in a time-averaged sense, and the mean drag is reduced by approximately 40%. The instantaneous flow is characterized by the periodic generation of dynamic-stall vortices near the leading edge and by their subsequent transition as they convect close to the airfoil. For Rec = 10 4 , significant reduction of the timeaveraged separation region is still possible with transitional effects present in the aft-portion of the airfoil. For larger forcing amplitude (ho = 0.04,Rec = 10 4 ), a very intriguing regime emerges. The dynamic stall vortex moves around and in front of the leading edge and experiences a dramatic breakdown as it impinges against the airfoil. As a result, the phased-averaged flow displays no coherent vortices propagating along the airfoil upper surface. This new flow structure is also characterized in the mean by the existence of a strong jet in the near wake which manifests in a high value of net thrust. The present study demonstrates the importance of transitional effects for low-Reynolds-number maneuvering airfoils, as well as the suitability of the ILES approch for exploring such flow regime.

Low-frequency dynamics of shock-induced separation in a compression ramp interaction

Abstract: The low-frequency dynamics of the shock-induced separation region in a Mach 2 compression ramp interaction is investigated by performing high-speed particle image velocimetry (HSPIV) measurements, at a rate of 6kHz, in a streamwise–spanwise plane. The HSPIV measurements made in the upstream turbulent boundary layer indicate the presence of spanwise strips of elongated regions of uniform streamwise velocity that extend to lengths greater than 30?, validating previous results based on planar laser scattering measurements obtained by Ganapathisubramani, Clemens & Dolling (J. Fluid Mech., vol. 585, 2007, p. 369). At a wall normal-location of y/?=0.2, a surrogate for separation based on a velocity threshold is found to fluctuate over a streamwise range of ±1.2?, consistent with previous studies. The amplitude of unsteadiness has contributions from at least two sources that are related to the incoming boundary layer. First, the velocity threshold based surrogate separation line exhibits large-scale undulations along the spanwise direction that conform to the passage of elongated low- and high-speed regions in the upstream boundary layer. This motion is classified as the local influence of the upstream boundary layer. Second, the spanwise-averaged surrogate separation is found to respond to the overall change in streamwise velocity in the incoming boundary layer and is classified as the global influence of the upstream boundary layer. However, this global influence includes the contributions from the elongated low- and high-speed regions. Preliminary findings based on statistical analysis suggest that the local influence contributes nearly 50% more than the global influence. Regardless, the low-frequency unsteadiness of the separation-region can be attributed to the local and global influences of the incoming boundary layer.

Steady separated flow past a circular cylinder at low Reynolds numbers

Abstract: The steady two-dimensional laminar flow around a stationary circular cylinder has been investigated via a stabilized finite-element method. The Reynolds number Re is based on the cylinder diameter and free-stream speed. The results have been presented for 6 ≤ Re ≤ 40 and the blockages between 0.000125 and 0.80. The blockage B is the ratio of the cylinder diameter to the domain width. There is large scatter in the value of Re s , reported in the literature, marking the onset of the flow separation. From the present study the Re s is found to be 6.29, approximately for B = 0.005. The effect of the blockage on the characteristic flow parameters is found to be insignificant for B ≤ 0.01. The bubble length, separation angle and Re s exhibit non-monotonic variation with the blockage. It is for the first time that such a behaviour for the separation angle and Re s is being reported. Two types of boundary conditions at the lateral walls have been studied: the slip wall and towing tank. In general for high blockage, the results from the slip boundary condition are closer to the ones for the unbounded flow. In that sense, the use of the slip boundary condition as opposed to the towing tank boundary condition on the lateral walls is advocated. The bubble length, separation angle, base suction, total drag, pressure drag, viscous drag and maximum vorticity on the cylinder surface for the steady flow are found to vary as Re, Re -0.5 , Re -1 , Re -0.5 , Re -0.64 , Re -0.60 and Re 0.5 , respectively. The extrapolated results for the steady flow, for higher Re, are found to match quite well with the other results from the literature.

Unsteady fluid-structure interactions of membrane airfoils at low Reynolds numbers

Abstract: Membrane wings are used both in nature and small aircraft as lifting surfaces. Separated flows are common at low Reynolds numbers and are the main sources of unsteadiness. Yet, the unsteady aspects of the fluid-structure interactions of membrane airfoils are largely unknown. An experimental study of unsteady aerodynamics of two-dimensional membrane airfoils at low Reynolds numbers has been conducted. Measurements of membrane shape with a high-speed camera were complemented with the simultaneous measurements of unsteady velocity field with a high frame-rate particle image velocimetry system and flow visualization. Vibrations of the membrane and mode shapes were investigated as a function of angle of attack and free stream velocity. While the mean membrane shape is not very sensitive to angle of attack, the amplitude and mode of the vibrations of the membrane depend on the relative location and the magnitude of the unsteadiness of the separated shear layer. The results indicate strong coupling of unsteady flow with the membrane oscillations. There is evidence of coupling of membrane oscillations with the vortex shedding in the wake, in particular, for the post-stall incidences. Comparison of rigid (but cambered) and flexible membrane airfoils shows that the flexibility might delay the stall. Hence this is a potential passive flow control method using flexibility in nature and engineering applications.

On the origin of the inflectional instability of a laminar separation bubble

Abstract: This is an experimental and theoretical Study of a laminar separation bubble and the associated linear stability mechanisms. Experiments were performed over a flat plate kept in a wind tunnel, with an imposed pressure gradient typical of an aerofoil that would involve a laminar separation bubble. The separation bubble was characterized by measurement of surface-pressure distribution and streamwise velocity using hot-wire anemometry. Single component hot-wire anemometry was also used for a detailed study of the transition dynamics. It was foundthat the so-called dead-air region in the front portion of the bubble corresponded to a region of small disturbance amplitudes, with the amplitude reaching a maximum value close to the reattachment point. An exponential growth rate of the disturbance was seen in the region upstream of the mean maximum height of the bubble, and this was indicative of a linear instability mechanism at work. An infinitesimal disturbance was impulsively introduced into the boundary layer upstream of separation location, and the wave packet was tracked (in an ensemble-averaged sense) while it was getting advected downstream. The disturbance was found to be convective in nature. Linear stability analyses (both the Orr-Sommerfeld and Rayleigh calculations) were performed for mean velocity profiles, starting from an attached adverse-pressure-gradient boundary layer all the way up to the front portion of the separation-bubble region (i.e. up to the end of the dead-air region in which linear evolution of the disturbance could be expected). The conclusion from the present work is that the primary instability mechanism in a separation bubble is inflectional in nature, and its origin can be traced back to upstream of the separation location. In other words, the inviscid inflectional instability of the separated shear layer should be logically seen as an extension of the instability of the upstream attached adverse-pressure-gradient boundary layer. This modifies the traditional view that pegs the origin of the instability in a separation bubble to the detached shear layer Outside the bubble, with its associated Kelvin-Helmholtz mechanism. We contendthat only when the separated shear layer has moved considerably away from the wall (and this happens near the maximum-height location of the mean bubble), a description by the Kelvin-Helmholtz instability paradigm, with its associated scaling principles, Could become relevant. We also propose a new scaling for the most amplified frequency for a wall-bounded shear layer in terms of the inflection-point height and the vorticity thickness and show it to be universal.

Effects of micro-ramps on a shock wave/turbulent boundary layer interaction

Abstract: Stereoscopic particle image velocimetry is used to investigate the effects of micro-ramp sub-boundary layer vortex generators, on an incident shock wave/boundary layer interaction at Mach 1.84. Single- and double-row arrangements of micro-ramps are considered. The micro-ramps have a height of 20% of the unperturbed boundary layer thickness and the measurement planes are located 0.1 and 0.6 boundary layer thicknesses from the wall. The micro-ramps generate packets of individual vortex pairs downstream of their vertices, which produce counter-rotating longitudinal streamwise vortex pairs in a time-averaged view. These structures induce a pronounced spanwise variation of the flow properties, namely the mixing across the boundary layer interface. The probability of reversed-flow occurrence is decreased by 20 and 30% for the single- and double-row configurations, respectively. Both configurations of micro-ramps stabilize the shock motion by reducing the length of its motion by about 20% in the lower measurement plane. The results are summarized by a conceptual model describing the boundary layer’s and interaction’s flow pattern under the effect of the micro-ramps.

Flow separation and resuspension beneath shoaling nonlinear internal waves

Abstract: [1] Laboratory observations are presented showing the structure and dynamics of the turbulent bottom boundary layer beneath nonlinear internal waves (NLIWs) of depression shoaling upon sloping topography. The adverse pressure gradient beneath the shoaling waves causes the rear face to steepen, flow separation to occur, and wave-induced near-bottom vortices to suspend bed material. The resuspension is directly attributed to the near-bed viscous stress and to near-bed patches of elevated positive Reynolds stress generated by the vortical structures. These results are consistent with published field observations of resuspension events beneath shoaling NLIWs. Elevated near-bed viscous stresses are found throughout the domain at locations that are not correlated to the resuspension events. Near-bed viscous stress is thus required for incipient sediment motion but is not necessarily a precursor for resuspension. Resuspension is dependent on the vertical velocity field associated with positive Reynolds stress and is also found to occur where the mean (wave-averaged) vertical velocity is directed away from the bed. The results are interpreted by analogy to the eddy-stress and turbulent bursting resuspension models developed for turbulent channel flows.

Generalized reynolds number for non-newtonian fluids

Abstract: An extended version of the generalized Reynolds number was derived to characterize the duct flow of non-Newtonian gelled fluids of the Herschel-Bulkley-Extended (HBE) Type. This number allows also estimating the transition from laminar to turbulent flow conditions. An experimental investigation was conducted with a capillary rheometer for several non-Newtonian gelled fluids to evaluate the introduced HBE-generalized Reynolds number Regen HBE. A good correlation between the experimental results and the theory could be found for laminar flow conditions. For one of the examined gelled fuels, the necessary high Reynolds numbers could be realized so that the transition from the laminar to the turbulent flow regime could be measured. Because of its general description, the HBE-generalized Reynolds number can also be applied to Newtonian liquids as well as to non-Newtonian fluids of the Herschel-Bulkley (HB), Ostwald-de-Waele (power-law, PL), and Bingham type.

Effects of wavelength and amplitude of a wavy cylinder in cross-flow at low Reynolds numbers

Abstract: Three-dimensional numerical simulations of laminar flow around a circular cylinder with sinusoidal variation of cross-section along the spanwise direction, named ‘wavy cylinder’, are performed. A series of wavy cylinders with different combinations of dimensionless wavelength (λ/Dm) and wave amplitude (a/Dm) are studied in detail at a Reynolds number of Re = U∞Dm/ν = 100, where U∞ is the free-stream velocity and Dm is the mean diameter of a wavy cylinder. The results of variation of mean drag coefficient and root mean square (r.m.s.) lift coefficient with dimensionless wavelength show that significant reduction of mean and fluctuating force coefficients occurs at optimal dimensionless wavelengths λ/Dm of around 2.5 and 6 respectively for the different amplitudes studied. Based on the variation of flow structures and force characteristics, the dimensionless wavelength from λ/Dm = 1 to λ/Dm = 10 is classified into three wavelength regimes corresponding to three types of wake structures. The wake structures at the near wake of different wavy cylinders are captured. For all wavy cylinders, the flow separation line varies along the spanwise direction. This leads to the development of a three-dimensional free shear layer with periodic repetition along the spanwise direction. The three-dimensional free shear layer of the wavy cylinder is larger and more stable than that of the circular cylinder, and in some cases the free shear layer even does not roll up into a mature vortex street behind the cylinder. As a result, the mean drag coefficients of some of the typical wavy cylinders are less than that of a corresponding circular cylinder with a maximum drag coefficient reduction up to 18%. The r.m.s. lift coefficients are greatly reduced to practically zero at optimal wavelengths. In the laminar flow regime (60 ≤ Re ≤ 150), the values of optimal wavelength are Reynolds number dependent.

High-fidelity simulations of moving and flexible airfoils at low Reynolds numbers

Abstract: The present paper highlights results derived from the application of a high-fidelity simulation technique to the analysis of low-Reynolds-number transitional flows over moving and flexible canonical configurations motivated by small natural and man-made flyers. This effort addresses three separate fluid dynamic phenomena relevant to small fliers, including: laminar separation and transition over a stationary airfoil, transition effects on the dynamic stall vortex generated by a plunging airfoil, and the effect of flexibility on the flow structure above a membrane airfoil. The specific cases were also selected to permit comparison with available experimental measurements. First, the process of transition on a stationary SD7003 airfoil section over a range of Reynolds numbers and angles of attack is considered. Prior to stall, the flow exhibits a separated shear layer which rolls up into spanwise vortices. These vortices subsequently undergo spanwise instabilities, and ultimately breakdown into fine-scale turbulent structures as the boundary layer reattaches to the airfoil surface. In a timeaveraged sense, the flow displays a closed laminar separation bubble which moves upstream and contracts in size with increasing angle of attack for a fixed Reynolds number. For a fixed angle of attack, as the Reynolds number decreases, the laminar separation bubble grows in vertical extent producing a significant increase in drag. For the lowest Reynolds number considered \((Re_c = 10^4)\), transition does not occur over the airfoil at moderate angles of attack prior to stall. Next, the impact of a prescribed high-frequency small-amplitude plunging motion on the transitional flow over the SD7003 airfoil is investigated. The motioninduced high angle of attack results in unsteady separation in the leading edge and in the formation of dynamic-stalllike vortices which convect downstream close to the airfoil. At the lowest value of Reynolds number \((Re_c = 10^4)\), transition effects are observed to be minor and the dynamic stall vortex system remains fairly coherent. For \(Re_c = 4 \times 10^4\), the dynamic-stall vortex system is laminar at is inception, however shortly afterwards, it experiences an abrupt breakdown associated with the onset of spanwise instability effects. The computed phased-averaged structures for both values of Reynolds number are found to be in good agreement with the experimental data. Finally, the effect of structural compliance on the unsteady flow past a membrane airfoil is investigated. The membrane deformation results in mean camber and large fluctuations which improve aerodynamic performance. Larger values of lift and a delay in stall are achieved relative to a rigid airfoil configuration. For \(Re_c = 4.85 \times 10^4\), it is shown that correct prediction of the transitional process is critical to capturing the proper membrane structural response.

Passive separation control using a self-adaptive hairy coating

Abstract: A model of hairy medium is developed using a homogenized approach, and the fluid flow around a circular cylinder partially coated with hair is analyzed by means of nu- merical simulations. The capability of this coating to adapt to the surrounding flow is investigated, and its benefits are discussed in the context of separation control. This fluid-structure interaction problem is solved with a partitioned approach, based on the direct resolution of the Navier-Stokes equations together with a non-linear set of equa- tions describing the dynamics of the coating. A volume force, arising from the presence of a cluster of hair, provides the link between the fluid and the structure problems. For the structure part, a subset of reference elements approximates the whole layer. The dy- namics of these elements is governed by a set of equations based on the inertia, elasticity, interaction and losses effects of articulated rods. The configuration chosen is that of the two-dimensional flow past a circular cylinder at Re = 200, a simple and well documented test case. Aerodynamics performances quantified by the Strouhal number, the drag and the maximum lift in the laminar unsteady regime are modified by the presence of the coating. A set of parameters corresponding to a realistic coating (length of elements, porosity, rigidity) is found, yielding an average drag reduction of 15% and a decrease of lift fluctuations by about 40%, associated to a stabilization of the wake.

Large Eddy Simulation of Turbulent Flow Through Submerged Vegetation

Abstract: Large Eddy Simulations (LES) are performed for an open channel flow through idealized submerged vegetation with a water depth (h) to plant height (hp) ratio of h/hp = 1.5 according to the experimental configuration of Liu et al. (J Geophys Res Earth Sci, 2008). They used a 1D laser Doppler velocimeter (LDV) to measure longitudinal and vertical velocities as well as turbulence intensities along several verticals in the flow and the data are used for the validation of the present simulations. The code MGLET is used to solve the filtered Navier–Stokes equations on a Cartesian non-uniform grid. In order to represent solid objects in the flow, the immersed boundary method is employed. The computational domain is idealized with a box containing 16 submerged circular cylinders and periodic boundary conditions are applied in both longitudinal and transverse directions. The predicted streamwise as well as vertical mean velocities are in good agreement with the LDV measurements. Furthermore, fairly good agreement is found between calculated and measured streamwise and vertical turbulence intensities. Large-scale flow structures of different shapes are present in the form of vortex rolls above the vegetation tops as well as locally generated trailing and von- Karman-type vortices due to flow separation at the free end and the sides of the cylinders. In this paper, the flow field is analyzed statistically and evidence is provided for the existence of these structures based on the LES.

Spatial development of particle-laden turbulent pipe flow

Abstract: The inhomogeneity of turbulence in wall bounded flows induces the phenomenology called turbophoresis whereby inertial particles of suitable mass accumulate at the solid wall. Particles injected near the axis of a fully turbulent pipe flow, after an initial spreading phase, undergo a segregation process which eventually leads to a pseudoequilibrium distribution sufficiently downstream. Wall densities up to thousand times the reference value can be easily achieved. The process is discussed here by analyzing the direct numerical simulation (DNS) data of a spatially developing particle laden pipe flow under the assumption of dilute suspension. Development phase and asymptotic state are addressed in quantitative terms. A Shannon-like entropy is introduced to quantify the level of spreading/segregation achieved by the particle distributions along the pipe. This allows to define on a physically sound basis the length of the developing region and to summarize in a single indicator the accumulation level as a function of the particle response time. By conditional statistics, it is unequivocally shown that particles approach the wall dragged by relatively fast yet comparatively rare events where highly accumulating particles follow the fluid in-rush toward the wall. On the contrary, the outward particle flux takes place in the form of much more frequent and gentle motions away from the wall. The analysis of DNS data and a simple argument highlight the role of the elongated clusters of particles at the wall as essential features responsible for the eventual asymptotic equilibrium.

Structures in turbulent boundary layers subjected to adverse pressure gradients

Abstract: The effects of adverse pressure gradients on turbulent structures were investigated by carrying out direct numerical simulations of turbulent boundary layers subjected to adverse and zero pressure gradients. The equilibrium adverse pressure gradient flows were established by using a power law free-stream distribution U∞ ~ xm. Two-point correlations of velocity fluctuations were used to show that the spanwise spacing between near-wall streaks is affected significantly by a strong adverse pressure gradient. Low-momentum regions are dominant in the middle of the boundary layer as well as in the log layer. Linear stochastic estimation was used to provide evidence for the presence of low-momentum regions and to determine their statistical properties. The mean width of such large-scale structures is closely associated with the size of the hairpin-like vortices in the outer layer. The conditionally averaged flow fields around events producing Reynolds stress show that hairpin-like vortices are the structures associated with the production of outer turbulence. The shapes of the vortices beyond the log layer were found to be similar when their length scales were normalized according to the boundary layer thickness. Estimates of the conditionally averaged velocity fields associated with the spanwise vortical motion were obtained by using linear stochastic estimation. These results confirm that the outer region of the adverse pressure gradient boundary layer is populated with streamwise-aligned vortex organizations, which are similar to those of the vortex packet model proposed by Adrian, Meinhart & Tomkins (J. Fluid Mech., vol. 422, 2000, pp. 1–54). The adverse pressure gradient augments the inclination angles of the packets and the mean streamwise spacing of the vortex heads in the packets.

Experimental study of flow separation in laminar falling liquid films

Abstract: In a previous publication, Dietze, Leefken & Kneer (J Fluid Mech, vol 595, 2008, p 435) showed that flow separation takes place in the capillary wave region of falling liquid films That investigation focused on the mechanistic explanation of the phenomenon mainly on the basis of numerical data The present publication for the first time provides clear experimental evidence of the phenomenon obtained by way of highly resolving velocity measurements in a specifically designed optical test set-up Characteristically, the refractive index of the working fluid was matched to that of the glass test section to provide optimal access to the cross-section of the film for the employed optical velocimetry techniques, namely, laser doppler velocimetry (LDV) and particle image velocimetry (PIV) Using LDV, time traces of the streamwise velocity component were recorded in high spatial (0025 mm) and temporal resolutions (04 ms) showing negative velocity values in the capillary wave region In addition, simultaneous film thickness measurements were performed using a Confocal Chromatic Imaging (CCI) technique enabling the correlation of velocity data and wave dynamics Further, using PIV the spatio-temporal evolution of the velocity field in the cross-section of the film was measured with high spatial (002 mm) and temporal (05 ms) resolutions yielding insight into the topology of the flow Most importantly these results clearly show the existence of a separation eddy in the capillary wave region Due to the high temporal resolution of the PIV measurements, enabled by the use of a high-speed camera with a repetition rate of up to 4500 Hz, the effect of wave dynamics on the velocity field in all regions of the wavy film was elucidated All experiments were performed using a dimethylsulfoxide (DMSO)―water solution and focused on laminar vertically falling liquid films with externally excited monochromatic surface waves Systematic variations of both the Reynolds number (Re = 86―150) and the excitation frequency (f = 16―24 Hz) were performed Results show that an increase in the wavelength of large wave humps, produced either by an increase in the Reynolds number or a decrease in the excitation frequency, leads to an increase in the size of the capillary separation eddy (CSE) Thereby, the CSE is shown to grow larger than the local film thickness, assuming an open shape with streamlines ending at the free surface

Two-dimensional unsteady laminar flow of a power law fluid across a square cylinder

Abstract: The two-dimensional and unsteady free stream flow of power law fluids past a long square cylinder has been investigated numerically in the range of conditions 60 ≤ R e ≤ 160 and 0.5 ≤ n ≤ 2.0 . Over this range of Reynolds numbers, the flow is periodic in time. A semi-explicit finite volume method has been used on a non-uniform collocated grid arrangement to solve the governing equations. The global quantities such as drag coefficients, Strouhal number and the detailed kinematic variables like stream function, vorticity and so on, have been obtained for the above range of conditions. While, over this range of Reynolds number, the flow is known to be periodic in time for Newtonian fluids, a pseudo-periodic flow regime displaying more than one dominant frequency in the lift is observed for shear-thinning fluids. This seems to occur at Reynolds numbers of 120 and 140 for n = 0.5 and 0.6, respectively. Broadly speaking, the smaller the value of the power law index, lower is the Reynolds number of the onset of the pseudo-periodic regime. This work is concerned only with the fully periodic regime and, therefore, the range of Reynolds numbers studied varies with the value of the power law index. Not withstanding this aspect, in particular here, the effects of Reynolds number and of the power law index have been elucidated in the unsteady laminar flow regime. The leading edge separation in shear-thinning fluids produces an increase in drag values with the increasing Reynolds number, while shear-thickening fluid behaviour delays this separation and shows the lowering of the drag coefficient with the Reynolds number. Also, the preliminary results suggest the transition from the steady to unsteady flow conditions to occur at lower Reynolds numbers in shear-thinning fluids than that in Newtonian fluids.

Direct numerical simulation of the early development of a turbulent mixing layer downstream of a splitter plate

Abstract: A direct numerical simulation is carried out of the initial stages of development of a mixing layer with a velocity ratio of ten, a fast stream Mach number of 0.6 and equal free-stream temperatures. The fast stream is a fully developed turbulent boundary layer with a trailing-edge displacement thickness Reynolds number of 2300, while the slow stream is laminar. The computations include a splitter plate with zero thickness. The initial flow development is dominated by the rapid spreading of an internal shear layer formed as the viscous sublayer of the upstream turbulent boundary layer crosses the trailing edge. A tendency towards spanwise-coherent structures is observed very early in the shear layer development, within five displacement thicknesses of the trailing edge, despite such structures not being present in the upstream boundary layer. A numerical search for a global mode in the vicinity of the splitter plate trailing edge found only convective growth of disturbances. Instead, a convective mechanism...

Role of surface roughness characterized by fractal geometry on laminar flow in microchannels.

Abstract: A three-dimensional model of laminar flow in microchannels is numerically analyzed incorporating surface roughness effects as characterized by fractal geometry. The Weierstrass-Mandelbrot function is proposed to characterize the multiscale self-affine roughness. The effects of Reynolds number, relative roughness, and fractal dimension on laminar flow are all investigated and discussed. The results indicate that unlike flow in smooth microchannels, the Poiseuille number in rough microchannels increases linearly with the Reynolds number, Re, and is larger than what is typically observed in smooth channels. For these situations, the flow over surfaces with high relative roughness induces recirculation and flow separation, which play an important role in single-phase pressure drop. More specifically, surfaces with the larger fractal dimensions yield more frequent variations in the surface profile, which result in a significantly larger incremental pressure loss, even though at the same relative roughness. The accuracy of the predicted Poiseuille number as calculated by the present model is verified using experimental data available in the literature.