Nonlinear evolution of harmonically forced perturbations on a wingtip vortex

Wingtip vortices are created by flying airplanes due to lift generation. The vortex interaction with the trailing aircraft has sparked researchers’ interest to develop an efficient technique to destroy these vortices. Different models have been used to describe the vortex dynamics and they all show that, under real flight conditions, the most unstable modes produce a very weak amplification. Another linear instability mechanism that can produce high energy gains in short times is due to the non-normality of the system. Recently, it has been shown that these non-normal perturbations also produce this energy growth when they are excited with harmonic forcing functions. In this study, we analyze numerically the nonlinear evolution of a spatially, pointwise and temporally forced perturbation, generated by a synthetic jet at a given radial distance from the vortex core. This type of perturbation is able to produce high energy gains in the perturbed base flow (10^3), and is also a suitable candidate for use in engineering applications. The flow field is solved for using fully nonlinear three-dimensional direct numerical simulation with a spectral multidomain penalty method model. Our novel results show that the nonlinear effects are able to produce locally small bursts of instability that reduce the intensity of the primary vortex.

Vortex Shedding Dynamics in Long Aspect-Ratio Aerodynamics Bodies

The focus of the current dissertation is to study qualitatively the underlying physics of vortex-shedding and wake dynamics in long aspect-ratio aerodynamics in incompressible viscous flow through the use of the KLE method. We carried out a long series of numerical experiments in the cases of flow around the cylinder at low Reynolds numbers. The study of flow at low Reynolds numbers provides an insight in the fluid physics and also plays a critical role when applying to stalled turbine rotors. Many of the conclusions about the qualitative nature of the physical mechanisms characterizing vortex formation, shedding and further interaction analyzed here at low Re could be extended to other Re regimes and help to understand the separation of the boundary layers in airfoils and other aerodynamic surfaces. In the long run, it aims to provide a better understanding of the complex multi-physics problems involving fluid-structure-control interaction through improved mathematical computational models of the multi-physics process. Besides the scientific conclusions produced, the research work on streamlined and bluff-body condition will also serve as a valuable guide for the future design of blade aerodynamics and the placement of wind turbines and hydrakinetic turbines, increasing the efficiency in the use of expensive workforce, supplies, and infrastructure. After the introductory section describing the main fields of application of wind power and hydrokinetic turbines, we describe the main features and theoretical background of the numerical method used here. Then, we present the analysis of the numerical experimentation results for the oscillatory regime right before the onset of vortex shedding for circular cylinders. We verified the wake length of the closed near-wake behind the cylinder and analysed the decay of the wake at the wake formation region, and then studied the St-Re relationship at the Reynolds numbers before the wake sheds compared to the experimental data. We found a theoretical model that describes the time evolution of the amplitude of fluctuations in the vorticity field on the twin vortex wake, which accurately matches the numerical results in terms of the frequency of the oscillation and rate of decay. We also proposed a model based on an analog circuit that is able to interpret the concerning flow by reducing the number of degrees of freedom. It follows the idea of the non-linear oscillator and resembles the dynamics mechanism of the closed near-wake with a common configured sine wave oscillator. This low-dimensional circuital model may also help to understand the underlying physical mechanisms, related to vorticity transport, that give origin to those oscillations.

Vortex evolution in the wake of accelerating low-aspect-ratio plates and three-dimensional bodies of revolution

To find examples of effecient locomotion and manoeuvrability, one need only turn to the elegant solutions natural flyers and swimmers have converged upon. This dissertation is specifically motivated by processes of evolutionary convergence, which have led to the propulsors and body shapes in nature that exhibit strong geometric collapse over diverse scales. These body features are abstracted in the studies presented herein using low-aspect-ratio at plates and a three-dimensional body of revolution (a sphere). The highly-separated vortical wakes that develop during accelerations are systematically characterized as a function of planform shape, aspect ratio, Reynolds number, and initial boundary conditions. To this end, force measurements and time-resolved (planar) particle image velocimetry have been used throughout to quantify the instantaneous forces and vortex evolution in the wake of the bluff bodies. During rectilinear motions, the wake development for the flat plates is primarily dependent on plate aspect ratio, with edge discontinuities and curvature playing only a secondary role. Furthermore, the axisymmetric case, i.e. the circular plate, shows strong sensitivity to Reynolds number, while this sensitivity quickly diminishes with increasing aspect ratio. For rotational motions, global insensitivity to plate aspect ratio has been observed. For the sphere, it has been shown that accelerations play an important role in the mitigation of flow separation. These results - expounded upon in this dissertation - have begun to shed light on the specific vortex dynamics that may be coopted by flying and swimming species of all shapes and sizes towards efficient locomotion.

Modélisation pharmacocinétique en imagerie par résonance magnétique et en tomographie d’émission par positrons appliquée à un modèle de glioblastome chez le rat

Résumé : En imagerie médicale, il est courant d’associer plusieurs modalités afin de tirer profit des renseignements complémentaires qu’elles fournissent. Par exemple, la tomographie d’émission par positrons (TEP) peut être combinée à l’imagerie par résonance magnétique (IRM) pour obtenir à la fois des renseignements sur les processus biologiques et sur l’anatomie du sujet. Le but de ce projet est d’explorer les synergies entre l’IRM et la TEP dans le cadre d’analyses pharmacocinétiques. Plus spécifiquement, d’exploiter la haute résolution spatiale et les renseignements sur la perfusion et la perméabilité vasculaire fournis par l’IRM dynamique avec agent de contraste afin de mieux évaluer ces mêmes paramètres pour un radiotraceur TEP injecté peu de temps après. L’évaluation précise des paramètres de perfusion du radiotraceur devrait permettre de mieux quantifier le métabolisme et de distinguer l’accumulation spécifique et non spécifique. Les travaux ont porté sur deux radiotraceurs de TEP (18F-fluorodésoxyglucose [FDG] et 18F-fluoroéthyle-tyrosine [FET]) ainsi que sur un agent de contraste d’IRM (acide gadopentétique [Gd DTPA]) dans un modèle de glioblastome chez le rat. Les images ont été acquises séquentiellement, en IRM, puis en TEP, et des prélèvements sanguins ont été effectués afin d’obtenir une fonction d’entrée artérielle (AIF) pour chaque molécule. Par la suite, les images obtenues avec chaque modalité ont été recalées et l’analyse pharmacocinétique a été effectuée par régions d’intérêt (ROI) et par voxel. Pour le FDG, un modèle irréversible à 3 compartiments (2 tissus) a été utilisé conformément à la littérature. Pour la FET, il a été déterminé qu’un modèle irréversible à 2 tissus pouvait être appliqué au cerveau et à la tumeur, alors qu’un modèle réversible à 2 tissus convenait aux muscles. La possibilité d’effectuer une conversion d’AIF (sanguine ou dérivée de l’image) entre le Gd DTPA et la FET, ou vice versa, a aussi été étudiée et s’est avérée faisable dans le cas des AIF sanguines obtenues à partir de l’artère caudale, comme c’est le cas pour le FDG. Finalement, l’analyse pharmacocinétique combinée IRM et TEP a relevé un lien entre la perfusion du Gd-DTPA et du FDG, ou de la FET, pour les muscles, mais elle a démontré des disparités importantes dans la tumeur. Ces résultats soulignent la complexité du microenvironnement tumoral (p. ex. coexistence de divers modes de transport pour une même molécule) et les nombreux défis rencontrées lors de sa caractérisation chez le petit animal.

Desenvolvimento e validação de método espectrofotométrico para determinação de corante à base de luteína adicionado em iogurte desnatado

A simple analytical method for extraction and quantification of lutein colorant added to yogurt was developed and validated. The method allowed complete extraction of carotenoids using tetrahydrofuran in vortex, followed by centrifugation, partition to diethyl ether/petroleum ether, and drying. The carotenoids dissolved in ethanol were quantified by UV-Vis spectrophotometry. This method showed linearity in the range tested (1.41-13.42 µg g-1), limits of detection and quantification of 0.42 and 1.28 µg g-1, respectively, low relative standard deviation (3.4%) and recovery ranging from 95 to 103%. The method proved reliable for quantification of lutein added to yogurt.

Análise de um ciclone semi-estacionário na costa sul do Brasil associado a bloqueio atmosférico

Na primeira semana de maio de 2008, durante quatro dias, um ciclone em superfície permaneceu semi-estacionário na costa da região sul do Brasil. Este sistema foi responsável por chuvas e ventos fortes no Rio Grande do Sul e Santa Catarina, os quais causaram muitos danos (queda de árvores, enchentes e desabamentos). O objetivo deste trabalho é avaliar o processo de formação e entender os mecanismos responsáveis pelo lento deslocamento do ciclone, já que a maioria dos ciclones nesta região possui deslocamento mais rápido. A equação de desenvolvimento de Sutcliffe mostrou que a advecção de vorticidade absoluta ciclônica na média troposfera e a advecção de ar quente na camada entre 1000-500 hPa foram mecanismos importantes para a ciclogênese. Neste período, o intenso aquecimento diabático também contribuiu para a ciclogênese, à medida que se contrapôs ao intenso resfriamento adiabático devido aos movimentos verticais ascendentes. A advecção de vorticidade absoluta ciclônica que favoreceu a ciclogênese esteve associada a um Vórtice Ciclônico em Altos Níveis (VCAN), que se formou numa região de anomalia de vorticidade potencial. O VCAN se manteve semi-estacionário e compôs o setor norte de um bloqueio do tipo dipolo. Tal bloqueio intensificou um anticiclone em superfície, situado a sul/leste do ciclone, o que contribuiu para o ciclone se manter semi-estacionário. O movimento atípico e lento do ciclone para sul, e em alguns períodos para sudoeste, esteve associado com advecções de vorticidade absoluta ciclônica na média troposfera e de ar quente no seu setor sul. Somente quando o bloqueio em níveis médios e a anomalia de vorticidade potencial em níveis médios/altos se enfraqueceram, o ciclone em superfície se afastou da costa sul do Brasil.

Convecção linearmente organizada na área de Petrolina, semi-árido do Nordeste do Brasil: aspectos em meso e grande escala

É apresentado um estudo sobre sistemas convectivos linearmente organizados e observados por um radar meteorológico banda-C na região semi-árida do Nordeste do Brasil. São analisados três dias (27 a 29) de março de 1985, com ênfase na investigação do papel desempenhado por fatores locais e de grande escala no desenvolvimento dos sistemas. No cenário de grande escala, a área de cobertura do radar foi influenciada por um cavado de ar superior austral no dia 27 e por um vórtice ciclônico de altos níveis no dia 29. A convergência de umidade próxima à superfície favoreceu a atividade convectiva nos dias 27 e 29, enquanto que divergência de umidade próxima à superfície inibiu a atividade convectiva no dia 28. No cenário de mesoescala, foi observado que o aquecimento diurno é um fator importante para a formação de células convectivas, somando-se a ele o papel determinante da orografia na localização dos ecos. De maneira geral, as imagens de radar mostram os sistemas convectivos linearmente organizados em áreas elevadas e núcleos convectivos intensos envolvidos por uma área de precipitação estratiforme. Os resultados indicam que convergência do fluxo de umidade em grande escala e aquecimento radiativo, são fatores determinantes na evolução e desenvolvimento dos ecos na área de estudo.

Numerical investigation of the three-dimensional secondary instabilities in the time-developing compressible mixing layer

Mixing layers are present in very different types of physical situations such as atmospheric flows, aerodynamics and combustion. It is, therefore, a well researched subject, but there are aspects that require further studies. Here the instability of two-and three-dimensional perturbations in the compressible mixing layer was investigated by numerical simulations. In the numerical code, the derivatives were discretized using high-order compact finite-difference schemes. A stretching in the normal direction was implemented with both the objective of reducing the sound waves generated by the shear region and improving the resolution near the center. The compact schemes were modified to work with non-uniform grids. Numerical tests started with an analysis of the growth rate in the linear regime to verify the code implementation. Tests were also performed in the non-linear regime and it was possible to reproduce the vortex roll-up and pairing, both in two-and three-dimensional situations. Amplification rate analysis was also performed for the secondary instability of this flow. It was found that, for essentially incompressible flow, maximum growth rates occurred for a spanwise wavelength of approximately 2/3 of the streamwise spacing of the vortices. The result demonstrated the applicability of the theory developed by Pierrehumbet and Widnall. Compressibility effects were then considered and the maximum growth rates obtained for relatively high Mach numbers (typically under 0.8) were also presented.

Hidden vorticity in binary Bose-Einstein condensates

We consider a binary Bose-Einstein condensate (BEC) described by a system of two-dimensional (2D) Gross-Pitaevskii equations with the harmonic-oscillator trapping potential. The intraspecies interactions are attractive, while the interaction between the species may have either sign. The same model applies to the copropagation of bimodal beams in photonic-crystal fibers. We consider a family of trapped hidden-vorticity (HV) modes in the form of bound states of two components with opposite vorticities S(1,2) = +/- 1, the total angular momentum being zero. A challenging problem is the stability of the HV modes. By means of a linear-stability analysis and direct simulations, stability domains are identified in a relevant parameter plane. In direct simulations, stable HV modes feature robustness against large perturbations, while unstable ones split into fragments whose number is identical to the azimuthal index of the fastest growing perturbation eigenmode. Conditions allowing for the creation of the HV modes in the experiment are discussed too. For comparison, a similar but simpler problem is studied in an analytical form, viz., the modulational instability of an HV state in a one-dimensional (1D) system with periodic boundary conditions (this system models a counterflow in a binary BEC mixture loaded into a toroidal trap or a bimodal optical beam coupled into a cylindrical shell). We demonstrate that the stabilization of the 1D HV modes is impossible, which stresses the significance of the stabilization of the HV modes in the 2D setting.

Formation of vortices in a dense Bose-Einstein condensate

A relaxation method is employed to study a rotating dense Bose-Einstein condensate beyond the Thomas-Fermi approximation. We use a slave-boson model to describe the strongly interacting condensate and derive a generalized nonlinear Schrodinger equation with a kinetic term for the rotating condensate. In comparison with previous calculations, based on the Thomas-Fermi approximation, significant improvements are found in regions where the condensate in a trap potential is not smooth. The critical angular velocity of the vortex formation is higher than in the Thomas-Fermi prediction.

Absorption of planar waves in a draining bathtub

We present an analysis of the absorption of acoustic waves by a black hole analogue in (2 + 1) dimensions generated by a fluid flow in a draining bathtub. We show that the low-frequency absorption length is equal to the acoustic hole circumference and that the high-frequency absorption length is 4 times the ergoregion radius. For intermediate values of the wave frequency, we compute the absorption length numerically and show that our results are in excellent agreement with the low-and high-frequency limits. We analyze the occurrence of superradiance, manifested as negative partial absorption lengths for corotating modes at low frequencies.

Emergence of Turbulence in an Oscillating Bose-Einstein Condensate

We report on the experimental observation of vortex tangles in an atomic Bose-Einstein condensate (BEC) of (87)Rb atoms when an external oscillatory perturbation is introduced in the trap. The vortex tangle configuration is a signature of the presence of a turbulent regime in the cloud. We also show that this turbulent cloud suppresses the aspect ratio inversion typically observed in quantum degenerate bosonic gases during free expansion. Instead, the cloud expands keeping the ratio between their axis constant. Turbulence in atomic superfluids may constitute an alternative system to investigate decay mechanisms as well as to test fundamental theoretical aspects in this field.

Volumetric reconstruction of the mean flow around circular cylinders fitted with strakes

The volumetric reconstruction technique presented in this paper employs a two-camera stereoscopic particle image velocimetry (SPIV) system in order to reconstruct the mean flow behind a fixed cylinder fitted with helical strakes, which are commonly used to suppress vortex-induced vibrations (VIV). The technique is based on the measurement of velocity fields at equivalent adjacent planes that results in pseudo volumetric fields. The main advantage over proper volumetric techniques is the avoidance of additional equipment and complexity. The averaged velocity fields behind the straked cylinders and the geometrical periodicity of the three-start configuration are used to further simplify the reconstruction process. Two straked cylindrical models with the same pitch (p = 10d) and two different heights (h = 0.1 and 0.2d) are tested. The reconstructed flow shows that the strakes introduce in the wake flow a well-defined wavelength of one-third of the pitch. Measurements of hydrodynamic forces, fluctuating velocity, vortex formation length, and vortex shedding frequency show the interdependence of the wake parameters. The vortex formation length is increased by the strakes, which is an important effect for the suppression of vortex-induced vibrations. The results presented complement previous investigations concerning the effectiveness of strakes as VIV suppressors and provide a basis of comparison to numerical simulations.

Experimental investigation of flow-induced vibration on isolated and tandem circular cylinders fitted with strakes

The effect of varying the geometric parameters of helical strakes on vortex-induced vibration (VIV) is investigated in this paper. The degree of oscillation attenuation or even suppression is analysed for isolated circular cylinder cases. How a cylinder fitted with strakes behaves when immersed in the wake of another cylinder in tandem arrangement is also investigated and these results are compared to those with a single straked cylinder. The experimental tests are conducted at a circulating water channel facility and the cylindrical models are mounted on a low-damping air bearing elastic base with one degree-of-freedom, restricted to oscillate in the transverse direction to the channel flow. Three strake pitches (p) and heights (h) are tested: p = 5, 10, 15d, and h = 0.1, 0.2, 0.25d. The mass ratio is 1.8 for all models. The Reynolds number range is from 1000 to 10000, and the reduced velocity varies up to 21. The cases with h = 0.1d strakes reduce the amplitude response when compared to the isolated plain cylinder, however the oscillation still persists. On the other hand, the cases with h = 0.2, 0.25d strakes almost completely suppress VIV. Spanwise vorticity fields, obtained through stereoscopic digital particle image velocimetry (SDPIV), show an alternating vortex wake for the p = 10d and h = 0.1d straked cylinder. The p = 10d and h = 0.2d cylinder wake has separated shear layers with constant width and no roll-up close to the body. The strakes do not increase the magnitude of the out-of-plane velocity compared to the isolated plain cylinder. However, they deflect the flow in the out-of-plane direction in a controlled way, which can prevent the vortex shedding correlation along the span. In order to investigate the wake interference effect on the strake efficiency, an experimental arrangement with two cylinders in tandem is employed. The centre-to-centre distance for the tandem arrangement varies from 2 to 6. When the downstream p = 10d and h = 0.2d cylinder is immersed in the wake of an upstream fixed plain cylinder, it loses its effectiveness compared with the isolated case. Although the oscillations have significant amplitude, they are limited, which is a different behaviour from that of a tandem configuration with two plain cylinders. For this particular case, the amplitude response monotonically increases for all gaps, except one, a trait usually found in galloping-like oscillations. SDPIV results for the tandem arrangements show alternating vortex shedding and oscillatory wake. (C) 2010 Elsevier Ltd. All rights reserved.

Secondary instabilities in the flow around two circular cylinders in tandem

Direct stability analysis and numerical simulations have been employed to identify and characterize secondary instabilities in the wake of the flow around two identical circular cylinders in tandem arrangements. The centre-to-centre separation was varied from 1.2 to 10 cylinder diameters. Four distinct regimes were identified and salient cases chosen to represent the different scenarios observed, and for each configuration detailed results are presented and compared to those obtained for a flow around an isolated cylinder. It was observed that the early stages of the wake transition changes significantly if the separation is smaller than the drag inversion spacing. The onset of the three-dimensional instabilities were calculated and the unstable modes are fully described. In addition, we assessed the nonlinear character of the bifurcations and physical mechanisms are proposed to explain the instabilities. The dependence of the critical Reynolds number on the centre-to-centre separation is also discussed.