An aluminum NACA 0018 airfoil testbed was constructed with 95 static pressure taps and 25 embedded microphones to enable novel time-resolved measurements of surface pressure. The main objective of this investigation is to utilize time-resolved surface pressure measurements to estimate salient flow characteristics in the separated flow region over the upper surface of an airfoil. The flow development over the airfoil was examined using hot wire anemometry and mean surface pressure for a range of Reynolds numbers from 80x103 to 200x103 and angles of attack from 0° to 18°. For these parameters, laminar boundary layer separation takes place on the upper surface and two flow regimes occur: (i) separation is followed by flow reattachment, so that a separation bubble forms and (ii) separation occurs without subsequent reattachment. Measurements of velocity and mean surface pressure were used to characterize the separated flow region and its effect on airfoil performance using the lift coefficient. In addition, the transition process and the evolution of disturbances were examined. The lift curve characteristics were found to be linked to the rate of change of the separation, transition, and reattachment locations with the angle of attack. For both flow regimes, transition was observed in the separated shear layer. Specifically, the amplification of disturbances within a band of frequencies in the separated shear layer resulted in laminar to turbulent transition. Validation of time-resolved surface pressure measurements was performed for Rec = 100x103 at α = 8° and α = 12°, corresponding to regimes of flow separation with and without reattachment, respectively. A comparative analysis of simultaneous velocity and time-resolved surface pressure measurements showed that the characteristics and development of velocity fluctuations associated with disturbances in the separated shear layer can be extracted from time-resolved surface pressure measurements. Specifically, within the separated flow region, the amplitude of periodic oscillations in the surface pressure signal associated with disturbances in the separated shear layer grew in the streamwise direction. In addition, the frequency at the spectral peak of the amplified disturbances in the separated shear layer was identified. Based on the results of the validation analysis, time-resolved surface pressure measurement analysis techniques were applied for a Reynolds number range from 60x103 to 130x103 and angles of attack from 6° to 16°. Within the separated flow region, the streamwise growth of surface pressure fluctuations is distinctly different depending on the flow regime. Specifically, within the separation bubble, the RMS surface pressure fluctuations increase in the streamwise direction and reach a peak just upstream of the reattachment location. The observed trend is in agreement with that observed for other separating-reattaching flows on geometries such as the forward and backward facing step and splitter plate with fence. In contrast to the separation bubble formation, when the separated shear layer fails to reattach to the airfoil surface, RMS surface pressure fluctuations increase in the streamwise direction with no maximum and the amplitude is significantly lower than those observed in the separation bubble. Surface pressure signals were further examined to identify the frequency, convective velocity, and spanwise uniformity of disturbances in the separated shear layer. Specifically, for both flow regimes, the fundamental frequency and corresponding Strouhal number exhibit a power-law dependency on the Reynolds number. Based on the available data for which velocity measurements were obtained in the separated flow region, the convective velocity matched the mean velocity at the wall-normal distance corresponding to the maximum turbulence intensity. A distinct increase in the convective velocity of disturbances in the separated shear layer was found when the airfoil was stalled in comparison to that found in the separation bubble. From statistical analysis of surface pressure signals in the spanwise direction, it was found that disturbances are strongly two-dimensional in the laminar portion of the separated shear layer and become three-dimensional through the transition process.

An aluminum NACA 0018 airfoil testbed was constructed with 95 static pressure taps and 25 embedded microphones to enable novel time-resolved measurements of surface pressure. The main objective of this investigation is to utilize time-resolved surface pressure measurements to estimate salient flow characteristics in the separated flow region over the upper surface of an airfoil. The flow development over the airfoil was examined using hot wire anemometry and mean surface pressure for a range of Reynolds numbers from 80x103 to 200x103 and angles of attack from 0° to 18°. For these parameters, laminar boundary layer separation takes place on the upper surface and two flow regimes occur: (i) separation is followed by flow reattachment, so that a separation bubble forms and (ii) separation occurs without subsequent reattachment. Measurements of velocity and mean surface pressure were used to characterize the separated flow region and its effect on airfoil performance using the lift coefficient. In addition, the transition process and the evolution of disturbances were examined. The lift curve characteristics were found to be linked to the rate of change of the separation, transition, and reattachment locations with the angle of attack. For both flow regimes, transition was observed in the separated shear layer. Specifically, the amplification of disturbances within a band of frequencies in the separated shear layer resulted in laminar to turbulent transition. Validation of time-resolved surface pressure measurements was performed for Rec = 100x103 at α = 8° and α = 12°, corresponding to regimes of flow separation with and without reattachment, respectively. A comparative analysis of simultaneous velocity and time-resolved surface pressure measurements showed that the characteristics and development of velocity fluctuations associated with disturbances in the separated shear layer can be extracted from time-resolved surface pressure measurements. Specifically, within the separated flow region, the amplitude of periodic oscillations in the surface pressure signal associated with disturbances in the separated shear layer grew in the streamwise direction. In addition, the frequency at the spectral peak of the amplified disturbances in the separated shear layer was identified. Based on the results of the validation analysis, time-resolved surface pressure measurement analysis techniques were applied for a Reynolds number range from 60x103 to 130x103 and angles of attack from 6° to 16°. Within the separated flow region, the streamwise growth of surface pressure fluctuations is distinctly different depending on the flow regime. Specifically, within the separation bubble, the RMS surface pressure fluctuations increase in the streamwise direction and reach a peak just upstream of the reattachment location. The observed trend is in agreement with that observed for other separating-reattaching flows on geometries such as the forward and backward facing step and splitter plate with fence. In contrast to the separation bubble formation, when the separated shear layer fails to reattach to the airfoil surface, RMS surface pressure fluctuations increase in the streamwise direction with no maximum and the amplitude is significantly lower than those observed in the separation bubble. Surface pressure signals were further examined to identify the frequency, convective velocity, and spanwise uniformity of disturbances in the separated shear layer. Specifically, for both flow regimes, the fundamental frequency and corresponding Strouhal number exhibit a power-law dependency on the Reynolds number. Based on the available data for which velocity measurements were obtained in the separated flow region, the convective velocity matched the mean velocity at the wall-normal distance corresponding to the maximum turbulence intensity. A distinct increase in the convective velocity of disturbances in the separated shear layer was found when the airfoil was stalled in comparison to that found in the separation bubble. From statistical analysis of surface pressure signals in the spanwise direction, it was found that disturbances are strongly two-dimensional in the laminar portion of the separated shear layer and become three-dimensional through the transition process.

Small wind turbines that operate in low-wind environments are prone to suffer performance degradation as they often fail to accelerate to a steady, power-producing condition. The behaviour during this process is called “starting behaviour” and it is the subject of this present work. This thesis evaluates potential benefits that can be obtained from the improvement of starting behaviour, investigates, in particular, small wind turbine starting behaviour (both horizontal- and vertical-axis), and presents aerofoil performance characteristics (both steady and unsteady) needed for the analysis. All of the investigations were conducted using a new set of aerodynamic performance data of six aerofoils (NACA0012, SG6043, SD7062, DU06-W-200, S1223, and S1223B). All of the data were obtained at flow conditions that small wind turbine blades have to operate with during the startup - low Reynolds number (from 65000 to 150000), high angle of attack (through 360◦), and high reduced frequency (from 0.05 to 0.20). In order to obtain accurate aerodynamic data at high incidences, a series of CFD simulations were undertaken to illustrate effects of wall proximity and to determine test section sizes that offer minimum proximity effects. A study was carried out on the entire horizontal-axis wind turbine generation system to understand its starting characteristics and to estimate potential benefits of improved starting. Comparisons of three different blade configurations reveal that the use of mixed-aerofoil blades leads to a significant increase in starting capability. The improved starting capability effectively reduces the time that the turbine takes to reach its power-extraction period and, hence, an increase in overall energy yield. The increase can be as high as 40%. Investigations into H-Darriues turbine self-starting capability were made through the analogy between the aerofoil in Darrieus motion and flapping-wing flow mechanisms. The investigations reveal that the unsteadiness associated with the rotor is key to predicting its starting behaviour and the accurate prediction can be made when this transient aerofoil behaviour is correctly modelled. The investigations based upon the analogy also indicate that the unsteadiness can be exploited to promote the turbine ability to self-start. Aerodynamically, this exploitation is related to the rotor geometry itself.

It is well known that when a flexible or flexibly-mounted structure is placed perpendicular to the flow of a Newtonian fluid, it can oscillate due to the shedding of separated vortices at high Reynolds numbers. If the same flexible object is placed in non-Newtonian flows, however, the structure's response is still unknown. The main objective of this thesis is to introduce a new field of viscoelastic fluid-structure interactions by showing that the elastic instabilities that occur in the flow of viscoelastic fluids can drive the motion of a flexible structure placed in its path. Unlike Newtonian fluids, the flow of viscoelastic fluids can become unstable at infinitesimal Reynolds numbers due to the onset of a purely elastic flow instability. This instability occurs in the absence of nonlinear effects of fluid inertia and the Reynolds number of the flows studied here are in the order of 10-4. When such an elastic flow instability occurs in the vicinity of a flexible structure, the fluctuating fluid forces exerted on the structure grow large enough to cause a structural instability which in turn feeds back into the fluid resulting in a flow instability. Nonlinear periodic oscillations of the flexible structure are observed which have been found to be coupled to the time-dependent growth and decay of viscoelastic stresses in the wake of the structure. Presented in this thesis are the results of an investigation of the interaction occurring in the flow of a viscoelastic wormlike micelle solution past a flexible rectangular sheet. The structural geometries studied include: flexible sheet inclinations at 20°, 45° and 90° and flexible sheet widths of 5mm and 2.5mm. By varying the flow velocity, the response of the flexible sheet has been characterized in terms of amplitude and frequency of oscillations. Steady and dynamic shear rheology and filament stretching extensional rheology measurements are conducted in order to characterize the viscoelastic wormlike micelle solution. Bright field images show the deformation of the flexible sheet during an unstable oscillation while flow-induced birefringence images highlight the viscoleastic fluid stresses produced in the wake of the flexible sheet.

No description available.

No description available.

Une étude expérimentale et numérique traitant de l'influence de structures tourbillonnaires sur la sédimentation de nuage de particules sphériques sous l'effet de la gravité est présentée dans une première partie de ce manuscrit. L'écoulement est créé par électro-convection, ce qui permet de générer un réseau de vortex contrôlés en vitesse et de taille constante qui imite un écoulement tourbillonnaire. Des techniques de PIV (Particle image-velocimetry) et de suivi de particules sont utilisés pour étudier la sédimentation du nuage.Le nuage est modélisé comme un ensemble de particules ponctuelles pour lesquelles les forces d'interaction hydrodynamiques entre particules sont prépondérantes. Le comportement du nuage est comparé aux prédictions obtenues avec des modèles numériques. Dans une seconde partie est présentée une étude expérimentale et numérique concernant la sédimentation à faible nombre de Reynolds de fibres flexibles dans un fluide visqueux au repos. L'état d'équilibre atteint par la fibre flexible est étudié. Nous identifions trois régimes ayant des signatures différentes sur l'état stationnaire de la fibre flexible: un régime de faibles déformations dans lequel la force de traînée est proportionnelle à celle d'une fibre sédimentant horizontalement par rapport à la gravité; un régime de grandes déformations dans lequel la force de traînée est aussi proportionnelle à la vitesse de la fibre, mais avec un coefficient de traînée qui est celui d'une fibre chutant parallèlement à la gravité; et un régime de reconfiguration élastique où le filament se déforme avec une traînée plus faible qui n'est plus proportionnelle à sa vitesse, mais à la racine carrée de celle-ci / In the first part, a jointed experimental and numerical study examining the influence of vortical structures on the settling of a cloud of solid spherical particles under the action of gravity at low Stokes numbers is presented. We use electro-convection to generate a two-dimensional array of controlled vortices which mimics a simplified vortical flow. Particle image-velocimetry and tracking are used to examine the motion of the cloud within this vortical flow. The cloud is modeled as a set of point-particles for which the hydrodynamic interaction is preponderant. The cloud behavior (trajectory, velocity, aspect ratio, break-up time …) is compared to the predictions of a two-way-coupling numerical simulation. In the second part, a jointed experimentally and numerical study on the dynamics of slender flexible filaments settling in a viscous fluid at low Reynolds number is presented. The equilibrium state of a flexible fiber settling in a viscous fluid is examined using a combination of macroscopic experiments, numerical simulations and scaling arguments. We identify three regimes having different signatures on this equilibrium configuration of the elastic filament: a weak deformation regime wherein the drag is proportional to the fiber velocity settling perpendicular to the gravity; a large deformation regime wherein the drag is proportional to the fiber velocity settling parallel to the gravity and an intermediate elastic reconfiguration regime where the filament deforms to adopt a shape with a smaller drag which is no longer linearly proportional to the velocity but to the square root of the velocity

O estudo do comportamento de um perfil de uma seção aerolástica típica, com Reynolds na faixa de microaeronaves, constitui o principal foco deste trabalho, tomando como objetivo a estimativa de parâmetros do fenômeno de flutter. A pesquisa analisa o escoamento de um fluido sobre um corpo (cilindro e perfil de aerofólio) em estado estacionário e oscilante, em escoamento de velocidade constante, e, especificamente, o fenômeno de flutter. As equações de Navier-Stokes, com termo de força, são resolvidas pelo método da fronteira virtual para modelagem da interface escoamento/estrutura, representada pela geometria de um corpo de geometria complexa. Na discretização das equações governantes foi utilizado o método de diferenças finitas, sobre malhas deslocadas, com avanço temporal das velocidades do escoamento por meio de um esquema de Runge-Kutta de ordem 4. Os códigos computacionais, para as simulações das diretrizes e a lógica de cálculo, foram criados no contexto deste trabalho. A verificação foi feita através do método da solução manufaturada por meio de um problema fictício, que tem uma solução analítica conhecida, e que preenche as condições de contorno implementadas no código. O modelo da fronteira virtual é testado para os casos de escoamento sobre cilindro de base quadrada, cilindro de base circular e perfil de aerofólio tipo NACA0012, com malhas regular e não regular, e para condições estacionária e sob oscilação forçada. Foi estudado o comportamento de formação de vórtices, provocados por escoamento uniforme sobre o perfil de aerofólio, através dos coeficientes de arrasto, sustentação e pressão com visualização por meio da vorticidade e linhas de corrente, para vários ângulos de ataque e oscilação forçada com elevação e rotação em torno de um pivô posicionado no centro geométrico do perfil (50% da corda). Finalmente, é apresentada uma determinação numérica das características aeroelásticas para o perfil de aerofólio NACA0012, em escoamento de número de Reynolds ultra baixo (Re = 1.000), e parâmetros de flutter para um caso de baixa frequência de oscilação. / The behavior study of a profile of a typical aerolastic section, with Reynolds in range of micro aerial vehicle, is the main focus of this work, taking as objective the estimation of parameters of flutter phenomenon. The research analyzes of the flow of a incompressible fluid on a body (cylinder and airfoil profile) at steady state and oscillating with constant speed and, specifically, the flutter phenomenon. The Navier-Stokes equations, with force term, are solved by virtual boundary method for modeling interface flow/structure, represented by the geometry of a body of complex geometry. In discretization of the governing equations, the method of finite differences on staggered grid, with temporal advancement of discharge velocity through a Runge-Kutta of order 4. The computer codes, for simulations guidelines and logic calculation, were created in the context of this work. The verification was done by method ofmanufactured solution through a fictional problem, which has a known analytical solution, and satisfies the boundary conditions implemented in code. The model of the virtual boundary is tested for cases of flow over a square cylinder, circular cylinder and profile of a NACA0012 airfoil type, with regular and non-regular meshes, over stationary and forced oscillation conditions. We studied the behavior of vortex formation, caused by uniform flow over the airfoil profile, by the drag, lift and pressure coefficients with view through the vorticity and streamlines for various attack angles and forced oscillation with plunge and pich around a pivot witch was positioned at the geometric airfoil profile (half chord). Finally, it is presented a numerical determination of aeroelastic characteristics for the NACA0012 airfoil profile, flow under ultra low Reynolds number, and flutter parameters for a case of low oscillation frequency.

O escoamento em torno de corpos tem sido objeto de estudo de muitos pesquisadores e é muito explorado experimental e computacionalmente, devido a sua grande aplicabilidade na engenharia. No entanto, simular computacionalmente este tipo de escoamento requer uma atenção especial ao escolher o tipo malha a ser utilizado. Em muitos casos faz-se necessário o uso de uma malha que se adapte ao contorno do obstáculo, o que pode ocasionar um aumento no esforço computacional. Um maneira de contornar este problema é a utilização do Método das Fronteiras Imersas, que possibilita o uso de malha cartesiana na simulação computacional do escoamento em torno de obstáculos. Isso é possível através da adição de um termo forçante nas equações que modelam o escoamento, e assim as forças que agem sobre o contorno do corpo são transferidas diretamente para a malha. O objetivo deste trabalho de mestrado foi implementar o método das Fronteiras Imersas e simular o escoamento em torno de um cilindro circular em repouso, movimentando-se na mesma direção do escoamento, na direção perpendicular ao escoamento, ou rotacionando em torno do próprio eixo. As simulações computacionais possibilitaram a captura do fenômeno de Atrelagem Síncrona, caracterizado pela sincronia entre a frequência de desprendimento natural de vórtices e a frequência de oscilação do mesmo. O Método das Fronteiras Imersas mostrou um ótimo desempenho quando comparado a resultados experimentais e numéricos encontrados na literatura / The flow around bodies have been studied by many researchers. Both experimental and computational approaches have been extensively explored in researches on flow around bodies and have been applied in many engeneering problems. However, to choose an appropriate type of mesh to perform computational simulations of this type of problem requires special attention. In many cases, it is necessary to use a mesh that is able to conform to the boundary if a given obstacle. The need to perform this adaptation may increase the computational effort. The Immersed Boundary Method enables the use of cartesian meshes to perform computational simulations of flows around obstacles. The idea of this method is to add a forcing term in the equations that model the flow. Thus, the forces applied on the body boundaries are directly transfered to the mesh. The aim of this work was to perform a computational implementation of the Immersed Boundary Method to simulate the flow over a oscilating circular cylinder. This oscilation may be inline with the flow, cross-flow, or rotating. The computational simulations enabled the capture of the lock-in phenomena, which consists of the syncronization between the vortex shedding frequency and the cylinder oscilation frequency. The results obtained from the computational simulations using the Immersed Boundary Method were in good agreement with the numerical and experimental results found in the literature

The aerodynamic performance of low aspect ratio wings at low Reynolds numbers applicable to micro air vehicle design was studied in this thesis. There is an overall lack of data for this low Reynolds number range, particularly concerning details of local flow behavior along the span. Experiments were conducted to measure the local pressure distributions on a wing at various spanwise locations in a Reynolds number range 30000 < Re < 90000. The model wing consisted of numerous wing sections and had a rectangular planform with NACA0012 airfoil shape with aspect ratio of one. One wing section, with pressure ports at various chordwise locations, was placed at different spanwise locations on a wing to effectively obtain the local pressure information. Integration of the pressure distributions yielded the local lift coefficients. Comparison of the local lift distributions to optimal elliptic lift distribution was conducted. This comparison showed a sharply peaked lift distribution near the wing tip resulting in a drastic deviation from the equivalent elliptic lift distributions predicted by the finite wing theory. The local lift distributions were further analyzed to determine the total lift coefficients vs angle of attack curves, span efficiency factors and the induced drag coefficients. Measured span efficiency factors, which were lower than predictions of the elliptic wing theory, can be understood by studying deviations of measured lift from the elliptic lift distribution. We conclude that elliptic wing theory is not sufficient to predict these aerodynamic performance parameters. Overall, these local measurements provided a better understanding of the low Reynolds number aerodynamics of the low aspect ratio wings.