Rotating-disk flow has been investigated not only as a simple model of cross flow instability to compare with swept-wing flow but also for industrial flow applications with rotating configurations. However the exact nature of laminar-turbulent transi- tion on the rotating-disk flow is still major problem and further research is required for it to be fully understood, in particular, the laminar-turbulent transition process with absolute instability. In addition the studies of the rotating-disk turbulent boundary- layer flow are inadequate to understand the physics of three-dimensional turbulent boundary-layer flow. In present thesis, a rotating-rotating disk boundary-layer flow has been inves- tigated experimentally using hot-wire anemometry. A glass disk with a flat surface has been prepared to archieve low disturbance rotating-disk environment. Azimuthal velocity measurements using a hot-wire probe have been taken for various conditions. To get a better insight into the laminar-turbulent transition region, a new way to describe the process is proposed using the probability density function (PDF) map of azimuthal fluctuation velocity. The effect of the edge of the disk on the laminar-turbulent transition process has been investigated. The disturbance growth of azimuthal fluctuation velocity as a function of Reynolds number has a similar trend irrespective of the various edge conditions. The behaviour of secondary instability and turbulent breakdown has been in- vestigated. It has been found that the kinked azimuthal velocity associated with secondary instability just before turbulent breakdown became less apparent at a cer- tain wall normal heights. Furthermore the turbulent breakdown of the stationary mode seems not to be triggered by its amplitude, however, depend on the appearance of the travelling secondary instability. Finally, the turbulent boundary layer on a rotating disk has been investigated. An azimuthal friction velocity has been directly measured from the azimuthal velocity profile in the viscous sub-layer. The turbulent statistics normalized by the inner and outer sclaes are presented. / QC 20120529

Transition from laminar to turbulent flow has beentraditionally studied in terms of exponentially growingeigensolutions to the linearized disturbance equations.However, experimental findings show that transition may occuralso for parameters combinations such that these eigensolutionsare damped. An alternative non-modal growth mechanism has beenrecently identified, also based on the linear approximation.This consists of the transient growth of streamwise elongateddisturbances, mainly in the streamwise velocity component,called streaks. If the streak amplitude reaches a thresholdvalue, secondary instabilities can take place and provoketransition. This scenario is most likely to occur in boundarylayer flows subject to high levels of free-stream turbulenceand is the object of this thesis. Different stages of theprocess are isolated and studied with different approaches,considering the boundary layer flow over a flat plate. Thereceptivity to free-stream disturbances has been studiedthrough a weakly non-linear model which allows to disentanglethe features involved in the generation of streaks. It is shownthat the non-linear interaction of oblique waves in thefree-stream is able to induce strong streamwise vortices insidethe boundary layer, which, in turn, generate streaks by thelift-up effect. The growth of steady streaks is followed bymeans of Direct Numerical Simulation. After the streaks havereached a finite amplitude, they saturate and a new laminarflow, characterized by a strong spanwise modulation isestablished. Using Floquet theory, the instability of thesestreaks is studied to determine the features of theirbreakdown. The streak critical amplitude, beyond which unstablewaves are excited, is 26% of the free-stream velocity. Theinstability appears as spanwise (sinuous-type) oscillations ofthe streak. The late stages of the transition, originating fromthis type of secondary instability, are also studied. We foundthat the main structures observed during the transition processconsist of elongated quasi-streamwise vortices located on theflanks of the low speed streak. Vortices of alternating signare overlapping in the streamwise direction in a staggeredpattern. Descriptors:Fluid mechanics, laminar-turbulenttransition, boundary layer flow, transient growth, streamwisestreaks, lift-up effect, receptivity, free-stream turbulence,nonlinear mechanism, streak instability, secondary instability,Direct Numerical Simulation. / QC 20100518

The main aim of this thesis was to understand the mechanism behind the classical transition scenario inside the boundary layer over an airfoil and eventually attempting to control this transition utilizing passive devices for transition delay. The initial objective of analyzing the transition phenomenon based on TS wave disturbance growth was conducted at 90 Hz using LDV and CTA measurement techniques at two different angles of attack. This was combined with the studies performed on two other frequencies of 100 and 110 Hz, in order to witness its impact on the neutral stability curve behavior. The challenges faced in the next phase of the thesis while trying to control the transition location, was to understand and encompass the effect of adverse pressure gradient before setting up the passive control devices, which in this case was miniature vortex generators. Consequently, several attempts were made to optimize the parameters of the miniature vortex generators depending upon the streak strength and stability. Finally, for 90 Hz a configuration of miniature vortex generators have been found to successfully stabilize the TS wave disturbances below a certain forcing amplitude, which also led to transition delay.

This thesis deals with the instabilities of the incompressible boundary-layer flow that is induced by a disk rotating in otherwise still fluid. The results presented are mostly limited to linear instabilities derived from direct numerical simulations (DNS) but with the objective that further work will focus on the nonlinear regime, providing greater insights into the transition route to turbulence. The numerical code Nek5000 has been chosen for the DNS using a spectral-element method in an effort to reduce spurious effects from low-order discretizations. Large-scale parallel simulations have been used to obtain the present results. The known similarity solution of the Navier–Stokes equation for the rotating-disk flow, also called the von Karman flow, is investigated and can be reproduced with good accuracy by the DNS. With the addition of small roughnesses on the disk surface, convective instabilities appear and data from the DNS are analysed and compared with experimental and theoretical data. A theoretical analysis is also presented using a local linear-stability approach, where two stability solvers have been developedbased on earlier work. A good correspondence between DNS and theory is found and the DNS results are found to explain well the behaviour of the experimental boundary layer within the range of Reynolds numbers for small amplitude (linear) disturbances. The comparison between the DNS and experimental results, presented for the first time here, shows that the DNS allows (for large azimuthal domains) a range of unstable azimuthal wavenumbers β to exist simultaneously with the dominantβ varying, which is not accounted for in local theory, where β is usually fixed for each Reynolds number at which the stability analysis is applied. Furthermore, the linear impulse response of the rotating-disk boundary layer is investigated using DNS. The local response is known to be absolutely unstable. The global response is found to be stable if the edge of the disk is assumed to be at infinity, and unstable if the domain is finite and the edge of the domain is placed such that there is a large enough pocket region for the absolute instability to develop. The global frequency of the flow is found to be determined by the edge Reynolds number. / <p>QC 20140708</p>

Há a necessidade de prevenir formação de gelo nas asas e nos estabilizadores de aeronaves, pois as formas de gelo podem causar a degradação do desempenho aerodinâmico, o aumento de peso, bem como dificuldades de controle e manobra que, em casos críticos, leva a uma diminuição da margem de segurança operacional. Quando as aeronaves atravessam nuvens com gotículas de água sub-resfriadas, ou seja, em equilíbrio metaestável, o crescimento de gelo ocorre nas superfícies não protegidas. Usualmente, os sistemas antigelo térmicos de aerofólios são projetados, desenvolvidos e certificados com o auxílio de programas de simulação numérica. O presente trabalho visa desenvolver e implementar um modelo matemático para prever a transferência de calor e massa no escoamento bidimensional bifásico em torno de aerofólios de uso aeronáuticos, equipados com sistema de antigelo térmico operando em regime permanente. Em condições de formação de gelo, é necessário aquecer o bordo de ataque e controlar a temperatura da região protegida para que não ocorra formação de gelo. O sistema de aquecimento compensa os efeitos do resfriamento imposto principalmente pelos mecanismos acoplados de evaporação e transferência de calor por convecção, que são causados pelo escoamento do ar carregado de gotículas sub-resfriadas e pelo escoamento da água líquida residual. O modelo deverá estimar a distribuição de temperaturas de superfície e o coeficiente de transferência de calor com precisão ao uso em aplicações aeronáuticas. O presente trabalho implementou novos submodelos para: 1) estimar a molhabilidade da superfície do aerofólio por meio de um modelo matemático para caracterizar o escoamento da água líquida residual na padrão de filme e de filetes; 2) avaliar o comportamento dinâmico e térmico da camada-limite laminar e turbulenta por meio de análises integral e diferencial, que considera efeitos do gradiente de pressão, da transição laminar-turbulenta, da transpiração e da não uniformidade de temperatura da superfície e 3) estimar o início e o término da região de transição laminar-turbulenta. O presente trabalho seguiu um processo de desenvolvimento de código numérico que: verificou os resultados de cada submodelo separadamente para depois implementados no modelo do antigelo; validou os resultados da simulação de desempenho do sistema antigelo com os novos submodelos implementados. Os resultados obtidos foram considerados satisfatórios para o modelo do antigelo que utilizou os submodelos de ruptura de filme e formação de filetes pelo critério da Energia Mecânica Total Mínima, de camada-limite diferencial compressível e de previsão da transição laminar-turbulenta por correlações algébricas, que consideraram efeitos do gradiente de pressão e do nível de turbulência ao longe. / It is required to prevent ice accretion on wings and horizontal stabilizers because it may cause aerodynamic performance degradation, weight increase, flight control difficulties and, in critical cases, may lead to operational safety margins reduction. When aircraft flies through clouds containnig supercooled water droplets, which are in metastable equilibrium, ice will form in all non-protected surfaces. Usually, anti-ice protection systems are designed, developed and certified with a support from a numerical tool. The present describes the development and implementation of a mathematical model for prediction of heat and mass transfer in two-phase flow around airfoils, which are equipped with thermal anti-ice system and operating in steady state regime. Under icing conditions, it is necessary to heat and control the temperature of the airfoil surface at leading edge region to prevent ice formation. The heating system balances the evaporative cooling effects, which are caused by the coupled heat and mass convection transfer, imposed by the air flow loaded with supercooled water droplets and the runback water flow around the airfoil. The present work implemented submodels to: 1) estimate airfoil surface wetness factor by adopting a liquid water film flow model as well as a rivulet formation and flow model; 2) evaluate laminar and turbulent boundary layers with pressure gradient and laminar-turbulent transition over nonisothermal and permeable airfoil surface by implementing differential boundary layer analysis and 3) predict the onset position and length of laminar-turbulent transition region. The present paper followed a validation and verification process during the numerical code development. All sub-models results were verified separately against experimental data before their inclusion in anti-ice model.The results of anti-ice model with selected submodels were validated against reference cases. The results were considered suficiently accurate when solving the film breakdown and rivulets formation by total mechanical energy method, compressible boundary layer by differential analysis and laminar-turbulent transition prediction by algebraic correlations, which considered pressure gradient and freestream turbulence level.

No presente trabalho foi implementado um modelo matemático para simular o sistema antigelo eletrotérmico de um aerofólio. Por meio do programa ONERA2D simulou-se o escoamento potencial completo com velocidade 44,7 m/s (100 mph) e 89,4 m/s (200 mph) em torno de um aerofólio perfil NACA0012 de corda 0,914 m (3 pés) com ângulo de ataque de 0°, e calculou-se a eficiência de coleta local de gotículas de água com diâmetro mediano volumétrico de 20 &#956;m. Foram simuladas quatro condições de teste com diferentes distribuições de fluxo de calor nos aquecedores elétricos do sistema antigelo. O modelo previu a distribuição de temperaturas na superfície sólida do aerofólio e no filme de água líquida, e as distribuições de fluxo de água líquida sobre a superfície do aerofólio (\"runback water\") e de coeficiente de transferência de calor por convecção de calor entre a superfície do aerofólio e o escoamento gasoso. Os resultados da simulação obtidos com o modelo foram comparados com resultados experimentais da NASA e os resultados numéricos dos programas LEWICE/ANTICE (EUA) e CANICE (Canada). Para as regiões molhadas pelo filme de água líquida, obteve-se um desvio máximo de temperatura de 2,6°C entre os resultados do presente modelo e o resultados experimentais. Para as regiões secas, onde não existe o filme de água líquida sobre a superfície do aerofólio, obteve-se um desvio de máximo de temperatura de 8°C. As previsões para distribuição de vazão de \"runback\", posição do término do filme de água líquida foram comparadas com os resultados do programa LEWICE/ANTICE. O modelo desenvolvido simula com adequada aproximação os efeitos da transferência de calor e de massa por convecção entre a superfície não-isotérmica do aerofólio ou do filme de água líquida e o escoamento gasoso, bem como os efeitos da transição entre o escoamento laminar e o turbulento na camada limite dinâmica e térmica e ainda a influência do escoamento do filme de água líquida sobre o desempenho do sistema de antigelo do aerofólio. / An electro-thermal anti-ice system was simulated with a mathematical model developed in the present work. A 44.7 m/s (100 mph) and 89.4 m/s (200 mph) full potential flow around a 0.914 m (3 ft) chord NACA0012 airfoil with 0° angle of attack and the local water catch efficiency of 20 &#956;m median volumetric diameter droplets impingement were calculated by the numerical code ONERA2D. Four test conditions were simulated with four different heat flux distributions of the anti-ice system according to the experimental work developed at NASA. The model predicted distributions of solid surface and liquid water film temperatures, runback water flow and convection heat transfer coefficient between airfoil or water surface and gaseous flow. The simulated results obtained by the mathematical model developed were compared to NASA experimental results and the ones predicted by the numerical codes LEWICE/ANTICE (US) and CANICE (Canada). For the regions wetted by the water film, the present model provided 2.6°C maximum temperature deviations between the predicted results and experimental data. For the dry regions, where there is no liquid water on the airfoil surface, an 8°C maximum temperature deviation was obtained. The runback flow and water film ending point position were compared to LEWICE/ANTICE numerical results. The developed model predicts adequately the convection heat and mass transfer effects between the non-isothermal airfoil or liquid water film surface and the gaseous flow, as well the effects of laminar to turbulent flow transition within dynamic and thermal boundary layer and the influence of the liquid water film flow on the anti-ice system performance.

La thèse vise à intégrer des critères de transition dans le solveur Navier-Stokes non structuré Aether utilisé chez Dassault Aviation. Une méthodologie de prévision de la transition laminaire/turbulent a été élaborée et implémentée dans le solveur RANS Aether. Deux stratégies de calcul de transition ont été testées. Soit Aether est couplé avec le code de couche limite de l’ONERA 3C3D. Soit la position de transition est calculée en utilisant directement les profils de vitesse RANS. Les deux méthodes ont été testées pour des écoulements subsoniques et transsoniques. L’influence des solveurs numériques, des critères de transition et du processus de couplage sont étudiés. L’utilisation de schémas numériques d’ordre élevé est également considérée. / This thesis present the state-of-the-art of the transition prediction numerical chain which has been developed at Dassault Aviation in the RANS solver AETHER. Two strategies for transition location estimations exist. First, AETHER is coupled with the ONERA boundary layer code 3C3D. Second, the transition location is computed by using directly the RANS velocity profiles. Both methods were preliminarily tested in subsonic and transonic. The issues of the influence of the numerical solvers, transition onset criteria and coupling process are studied. The influence of higher order numerical method, are also considered.

Em diversos escoamentos sobre superfícies há a presença de protuberâncias, como por exemplo rebites, parafusos e juntas. Estas protuberâncias podem influenciar a camada limite, acelerando a transição do escoamento do estado laminar para o estado turbulento. Em alguns casos isto pode ser indesejável, já que o escoamento turbulento implica necessariamente em uma força de atrito maior do que aquela referente ao escoamento laminar. Existem alguns aspectos neste tipo de escoamento que ainda não estão bem compreendidos. O objetivo deste trabalho é estudar a influência de uma rugosidade isolada no escoamento sobre uma superfície. Este estudo contribui para se entender o que ocorre em casos de maior complexidade. O estudo é de natureza computacional, em que se utiliza simulação numérica direta das equações de Navier-Stokes. A técnica de fronteiras imersas é utilizada para representar a rugosidade no escoamento sobre a superfície. O código numérico é verificado por meio do método de soluções manufaturadas. Comparações entre resultados experimentais, da teoria de estabilidade linear e numéricos também são utilizados para a validação do código. Resultados obtidos com diferentes alturas de rugosidade e variações no gradiente de pressão permitiram analisar a influência de elemento rugoso tridimensional em escoamentos de camada limite. / The presence of protuberances on surfaces, for example, rivets, screws and gaskets, can influence the boundary layer by accelerating the transition from laminar flow to turbulent flow. In some cases this may be undesirable, since the turbulent flow involves frictional forces greater than the ones at the laminar regime. There are some aspects of the flow in the boundary layer perturbed by a single roughness element that are not well understood. The aim of this work is to study the influence of an isolated roughness on the boundary layer. This study is a step towards to the understanding of what can happen in more complex cases. The nature of this study is computational, therefore a Direct Numerical Simulation code is used. The immersed boundary method is used to represent the roughness in the flow on the surface. The numerical code is verified via theMethod ofManufactured Solutions. Comparisons between experimental data, Linear Stability Theory and numerical results are also used for the validation of the code. Results obtained with different roughness heights and variations in the pressure gradient allowed the analysis of the influence of a three-dimensional roughness element in boundary layer flows.

No presente trabalho foi implementado um modelo matemático para simular o sistema antigelo eletrotérmico de um aerofólio. Por meio do programa ONERA2D simulou-se o escoamento potencial completo com velocidade 44,7 m/s (100 mph) e 89,4 m/s (200 mph) em torno de um aerofólio perfil NACA0012 de corda 0,914 m (3 pés) com ângulo de ataque de 0°, e calculou-se a eficiência de coleta local de gotículas de água com diâmetro mediano volumétrico de 20 &#956;m. Foram simuladas quatro condições de teste com diferentes distribuições de fluxo de calor nos aquecedores elétricos do sistema antigelo. O modelo previu a distribuição de temperaturas na superfície sólida do aerofólio e no filme de água líquida, e as distribuições de fluxo de água líquida sobre a superfície do aerofólio (\"runback water\") e de coeficiente de transferência de calor por convecção de calor entre a superfície do aerofólio e o escoamento gasoso. Os resultados da simulação obtidos com o modelo foram comparados com resultados experimentais da NASA e os resultados numéricos dos programas LEWICE/ANTICE (EUA) e CANICE (Canada). Para as regiões molhadas pelo filme de água líquida, obteve-se um desvio máximo de temperatura de 2,6°C entre os resultados do presente modelo e o resultados experimentais. Para as regiões secas, onde não existe o filme de água líquida sobre a superfície do aerofólio, obteve-se um desvio de máximo de temperatura de 8°C. As previsões para distribuição de vazão de \"runback\", posição do término do filme de água líquida foram comparadas com os resultados do programa LEWICE/ANTICE. O modelo desenvolvido simula com adequada aproximação os efeitos da transferência de calor e de massa por convecção entre a superfície não-isotérmica do aerofólio ou do filme de água líquida e o escoamento gasoso, bem como os efeitos da transição entre o escoamento laminar e o turbulento na camada limite dinâmica e térmica e ainda a influência do escoamento do filme de água líquida sobre o desempenho do sistema de antigelo do aerofólio. / An electro-thermal anti-ice system was simulated with a mathematical model developed in the present work. A 44.7 m/s (100 mph) and 89.4 m/s (200 mph) full potential flow around a 0.914 m (3 ft) chord NACA0012 airfoil with 0° angle of attack and the local water catch efficiency of 20 &#956;m median volumetric diameter droplets impingement were calculated by the numerical code ONERA2D. Four test conditions were simulated with four different heat flux distributions of the anti-ice system according to the experimental work developed at NASA. The model predicted distributions of solid surface and liquid water film temperatures, runback water flow and convection heat transfer coefficient between airfoil or water surface and gaseous flow. The simulated results obtained by the mathematical model developed were compared to NASA experimental results and the ones predicted by the numerical codes LEWICE/ANTICE (US) and CANICE (Canada). For the regions wetted by the water film, the present model provided 2.6°C maximum temperature deviations between the predicted results and experimental data. For the dry regions, where there is no liquid water on the airfoil surface, an 8°C maximum temperature deviation was obtained. The runback flow and water film ending point position were compared to LEWICE/ANTICE numerical results. The developed model predicts adequately the convection heat and mass transfer effects between the non-isothermal airfoil or liquid water film surface and the gaseous flow, as well the effects of laminar to turbulent flow transition within dynamic and thermal boundary layer and the influence of the liquid water film flow on the anti-ice system performance.

In this thesis direct numerical simulations are used to investigate two phenomenain shear flows: laminar-turbulent transition over a flat plate and periodicvortex shedding induced by a jet in cross flow. The emphasis is on understanding and controlling the flow dynamics using tools from dynamical systems and control theory. In particular, the global behavior of complex flows is describedand low-dimensional models suitable for control design are developed; this isdone by decomposing the flow into global modes determined from spectral analysisof various linear operators associated with the Navier–Stokes equations.Two distinct self-sustained global oscillations, associated with the sheddingof vortices, are identified from direct numerical simulations of the jet incrossflow. The investigation is split into a linear stability analysis of the steadyflow and a nonlinear analysis of the unsteady flow. The eigenmodes of theNavier–Stokes equations, linearized about an unstable steady solution revealthe presence of elliptic, Kelvin-Helmholtz and von K´arm´an type instabilities.The unsteady nonlinear dynamics is decomposed into a sequence of Koopmanmodes, determined from the spectral analysis of the Koopman operator. Thesemodes represent spatial structures with periodic behavior in time. A shearlayermode and a wall mode are identified, corresponding to high-frequency andlow-frequency self-sustained oscillations in the jet in crossflow, respectively.The knowledge of global modes is also useful for transition control, wherethe objective is to reduce the growth of small-amplitude disturbances to delaythe transition to turbulence. Using a particular basis of global modes, knownas balanced modes, low-dimensional models that capture the behavior betweenactuator and sensor signals in a flat-plate boundary layer are constructed andused to design optimal feedback controllers. It is shown that by using controltheory in combination with sensing/actuation in small, localized, regionsnear the rigid wall, the energy of disturbances may be reduced by an order of magnitude.